Domestication and Genetics of Papaya: A Review

A wealth of plant species used by humans for different purposes, but mainly as food, originated and domesticated in the Mesoamerican region. Papaya (Carica papaya) is the third most cultivated tropical crop worldwide, and it has been hypothesized that Mesoamerica is the most likely center of its origin and domestication. In support of it, many wild populations of papaya occur throughout Mesoamerica and hence represent the gene pool of genetic variability for further evolution and future crop management. Despite its importance, a dearth of information exists regarding the status of wild populations of papaya, as compared to the extent of knowledge, and interest, on domesticated varieties. We review the evidence on the extant wild populations of papaya, as well as its origin and distribution. Also, we synthetize what is known on the domestication history of the species, including the domestication syndrome that distinguishes wild and domesticated papayas. Moreover, we make an account of the use of genetic markers to assess genetic diversity of wild and domesticated papaya, and discuss the importance of papaya as the first species with a transgenic cultivar to be released for human consumption, and one that has its complete genome sequenced. Evidence from different disciplines strongly suggest that papaya originated and was domesticated in Mesoamerica, and that wild populations in the region possess, still, high genetic diversity compared to the domesticated papaya. Finally, we outline papaya as an excellent model species for genomic studies that will help gain insight into the domestication process and improvement of papaya and other tropical crops.

The Caricaceae members arrived to South America from Central America between 27 and 19 MYA, when the Central America land bridge already had begun to form, facilitating the range expansion from Mexico to South America (Carvalho and Renner, 2012). It is estimated that C. papaya diverged from its sister clade about 25 MYA, and belongs to a small clade restricted to Mexico, Guatemala and El Salvador that includes four species: three perennial herbs, Jarilla chocola, J. heterophylla, and J. nana that occur in seasonal tropical forests, and the treelet Horovitzia cnidoscoloides, endemic to Oaxaca in Mexico (Carvalho and Renner, 2014).
Papaya possess a morphological structure and development according to Corner's model of architecture (Hallé et al., 1978): a monopodial, single, orthotropic and nonbranching trunk constructed by one vegetative meristem, with axillary inflorescences, hence with indeterminate growth. Carica papaya produces huge palmate-shaped leaves. Female flowers are produced at the axils of the leaf petiole, and later the fruits will occupy that position along the stem. Male individuals produce inflorescences. Three sex types of papaya are known: female, male, and hermaphrodite trees. In wild populations, most individuals are diclinous (dioecy), whereas cultivated papayas are dioecious or hermaphrodites (Carvalho and Renner, 2012;Chávez-Pesqueira et al., 2014). Cultivars can inbred, resulting in stable characteristics across generations (Manshardt, 1992). For gynodioecious cultivated papaya, it has been reported that two-thirds of the plants correspond to hermaphrodites and onethird to female plants, though dioecious cultivars exist (VanBuren et al., 2015). Sex determination for the three types of papaya, female, male, and hermaphrodite plants, is genetically regulated by the pairing of sex chromosomes (Carvalho and Renner, 2012), through a sex-linked region that behaves like an XY sex chromosome. For male and hermaphrodite individuals, sex is controlled by slightly different Y chromosome regions: Y h in hermaphrodites and Y in males (VanBuren et al., 2015). Female plants produce flowers and fruits all year round in tropical regions; however, in subtropical areas, despite continuous flowering, fruit set is decreased during drier seasons (Gonsalves, 1998).
Papaya is the third most cultivated tropical crop world-wide, Brazil and India are the largest producers of papaya although Mexico is the main exporter (Evans and Ballen, 2012). Under cultivation papaya trees grow fast, producing mature fruits within 9-12 months after planting. Commercially, a density of 1,500-2,500 trees per hectare, can produce from 125,000 to 300,000 lbs per hectare, per year (Gonsalves, 1998). Among common fruits, papaya is ranked first on nutritional scores for the percentage of vitamin A, vitamin C, potassium, folate, niacin, thiamine, riboflavin, iron and calcium, and fiber (Huerta-Ocampo et al., 2012). Moreover, fruits, stems, leaves and roots of papaya are used in a wide range of medical applications and papain production (Ming et al., 2008). Commercial production of papain is directed for protein digestion, mainly as a red meat tenderizer, for the brewing of beer, and the skin treatment of warts and scars (Ming et al., 2012). Because of its rapid growth, continuous harvest and multiple uses, papaya is widely common in home gardens of tropical regions (Manshardt, 1992).
Wild, undomesticated, populations of papaya still grow in many regions of Mesoamerica, but it is common in naturally disturbed tropical forests (Chávez-Pesqueira and Núñez-Farfán, 2016). On occasion, wild papaya is also found in home gardens of ethnic groups in southern Mexico (Paz and Vázquez-Yanes, 1998) where it is used to make jams and dried crystallized fruit candy (Chávez-Pesqueira and Núñez-Farfán, pers. obs). In the rain forests of southern Mexico, wild papaya trees behave like a typical fast-growing, short-lived, nomadic tree (Van Steenis, 1958); they establish rapidly and reach maturity and reproduction only in recent (1-5 years, Paz and Vázquez-Yanes, 1998) and moderately large canopy gaps in mature rain forest, as well as in early secondary forests or "acahuals, " and man-made clearings. Because of this, wild papaya is rare in the rain forest since it depends on the creation of new canopy gaps to persist and disperse (Chávez-Pesqueira and Núñez-Farfán, 2016). In mature forests, papaya plants mainly die due to shading by other trees that "close" the canopy gaps by lateral growth. Because of its nomadic nature, wild papaya represents a key element in the regeneration dynamic of tropical and sub-tropical forests along its natural distribution (Chávez-Pesqueira et al., 2014).
In southern Mexico, isolated individuals of papaya, and some in home gardens, show intermediate fruit phenotypes between wild and cultivated plants (Paz and Vázquez-Yanes, 1998). These phenotypes suggest introgression due to occasional mating between population types (Chávez-Pesqueira et al., 2014). Since wild C. papaya represents the genetic reservoir of the species (Chávez-Pesqueira and Núñez-Farfán, 2016), exhaustive efforts should be made to know and preserve this important species in its wild form. Moreover, because papaya is one of the most economically important tropical crops in the world, and wild populations still occur naturally, the species represents an ideal system to study, in depth, the process of domestication. Here, we present a review of the existent knowledge on the origin, distribution and domestication of papaya, and discuss the usage of genetic and genomic methods to study the domestication of this valuable species.

ORIGIN AND DISTRIBUTION
Different authors suggest a Mesoamerican origin of C. papaya, comprising southern Mexico to Central America (Vavilov, 1926;Storey, 1976), although no direct archaeological evidence regarding the center of origin of papaya has been reported as yet (Fuentes and Santamaría, 2014). One reason for the lack of archaeological data is the difficulty to identify papaya from phytoliths, and their pollen grains have been hardly found (Carvalho and Renner, 2012). Together, the presence of undomesticated populations in Mexico and Central America and its cultivation in the region predating the Spanish conquest of Mexico, support the Mesoamerican hypothesis (Carvalho and Renner, 2014;Fuentes and Santamaría, 2014).
Phylogenetic evidence of the Caricaceae family, also supports a Mesoamerican origin of papaya. Carvalho and Renner (2012) obtained a molecular phylogeny of the Caricaceae family using chloroplast and nuclear data (4,711 bp) of the 34 species in the family. This is the only phylogeny using members from all the Caricaceae family. Their resulting phylogeny points that Carica papaya is more closely related to the genus Horovitzia, endemic to Mexico, and Jarilla, endemic to Mexico and Guatemala, than to the South American genus Vasconcellea, as previously thought (Carvalho and Renner, 2014). Moreover, Carica, Jarilla and Horovitzia show an unilocular ovary, whereas the remaining South American Caricaceae possess 5-locular ovaries (Carvalho and Renner, 2012). This morphological synapomorphy supports the family phylogeny.
Further evidence is offered by the origin of the Y h chromosome in cultivated hermaphrodite individuals of papaya (VanBuren et al., 2015). By sequencing the entire male-specific region of the Y chromosome and comparing it with the previous sequences of the hermaphrodite-specific region of the Y h chromosome, (Wang et al., 2012),VanBuren et al. (2015 found that the Y h chromosome possess lower nucleotide diversity as compared to the Y. This reduced variability is consistent with a genetic bottleneck scenario possibly brought about by domestication, and suggests that dioecy is ancestral in C. papaya. Given the ubiquity of dioecy in wild populations of papaya in Mesoamerica, an origin for the species in this region is a more parsimonious hypothesis. Yet, genomic studies can aid to determine the origin of some innovations related with domestication. High values of genetic variation of wild populations, and wild relatives, are expected to be found in the centers of origin of crop species (Gepts and Papa, 2003). Additionally, domestication is expected to reduce genetic diversity and provoke selective sweeps in genes associated to characters target of domestication (Purugganan and Fuller, 2009). A recent phylogeographic analysis of 19 wild populations of papaya in Mexico using nuclear and chloroplast markers, revealed high values of genetic diversity (Chávez-Pesqueira and Núñez-Farfán, 2016). The higher genetic diversity was found in locations of southern Mexico, suggesting this region as a genetic reservoir for the species (see below). Summing up, evidence strongly suggests that C. papaya originated in Mesoamerica, likely in southern Mexico.
Regarding the natural distribution of papaya, this has been suggested to range from the northern tropical limit of Mexico to Costa Rica in Central America (Aradhya et al., 1999;Carvalho and Renner, 2012). However, the precise assessment of its natural distribution is still lacking. One reason for this is the scarcity of studies in wild populations of papaya, and the paucity of herbarium specimens indicating whether a specimen belongs to a cultivated or wild individual. Fuentes and Santamaría (2014) performed searches in herbaria around the world to explore the distribution of papaya. Although they did not distinguish between wild and cultivated plants, their results show that most specimens belong to Mexico and to a lesser extent to Central America Unfortunately, the present rates of deforestation and habitat fragmentation within the proposed distribution range of wild papaya is high enough as to endanger the persistence of plant species (Barlow et al., 2016). This, coupled with the lack of information on the state of many wild varieties of important crop species, warn us about the relevance of studying and conserving wild populations and wild relatives of papaya and other crop species.

DOMESTICATION
The limited occurrence of wild populations of C. papaya and its four closest relative species in Mexico and Central America, are in line with the domestication of papaya in the Mesoamerican region (Carvalho and Renner, 2012). Mesoamerica is considered as one of the World's centers of plant domestication (Harlan, 1971;Pohl et al., 1996). The Maya was the most important culture present in that region before the conquest of Mexico by Spain in the sixteenth century, and probably one of the first to cultivate and trade the fruits of C. papaya (Colunga-GarcíaMarín and Zizumbo-Villarreal, 2004). Moreover, it has been suggested that the enzime papain was used by the Mesoamerican people to tenderize meat by wrapping it in papaya leaves and that this knowledge was then taken to Europe after the Spanish colonization (Larqué-Saavedra, 2016). By the time of the conquest of Mexico, it is believed that papaya was cultivated by native people all the way from southern Mexico to the Isthmus of Panama, where it was locally known as olocoton (Storey et al., 1986).
In the sixteenth century the Spaniards were probably the first responsible for the spread of papaya beyond Mesoamerica (Carvalho and Renner, 2014). It was introduced into the Hispaniola island (nowadays Haiti and Santo Domingo) in 1,521. There, it acquired the Carib Indian name ababai, that would then be changed to papaia, papia, papeya, and finally papaya. From this island, it spread to other West Indian islands, such as Jamaica, where it was known as pawpaw; to Cuba where it was named fruta bomba; to Venezuela known as lechosa; and to Brazil and Argentina where it is called mamao and mamón, respectively. In 1526 papaya was taken outside of America to Indonesia, and then spread commercially throughout the East Indies and tropical Asia (Storey et al., 1986). After that, it rapidly spread into other Asian countries and finally to Africa, brought by European colonial powers, such as Portugal, Denmark, Great Britain and France (Manshardt, 2014).
Moreover, because the papaya seeds have a moderate period of longevity, it is likely that this trait aided to its rapid spread throughout the tropics, where it has existed practically since man has recorded modern history (Schroeder, 1958). Nowadays, there are many varieties of papaya cultivated in most tropical and subtropical regions of the world, differing in traits such as fruit size, color, flavour, and tree size (Moore, 2014).
Domesticated plants commonly display a suite of traits, defined as the "domestication syndrome, " that distinguish them from their wild populations or relatives. In general, they differ in traits related to yield, food usage, and cultivation (Martínez-Ainsworth and Tenaillon, 2016). For papaya, domesticated plants have become morphologically and physiologically different in relation to their wild counterparts (Paz and Vázquez-Yanes, 1998). The principal characters of papaya that have been studied as target of selection under domestication are: tree size, fruit size, sex types, and morphology, and germinability of seeds. Whereas tree size has been selected to become smaller to facilitate fruit harvesting (Niklas and Marler, 2007), selection has been focused on enlarging the fruit size and increasing the ovary wall (pulp) for human consumption Manshardt, 2014). As a result, several varieties of papaya are cultivated around the world.
Regarding sex types, in dioecious fruit crops, mutations inducing hermaphroditism have been associated with domestication; this is the case of species like strawberry, grape and papaya (Janick, 2006). In the Caricaceae family, dioecy is the ancestral stateand it has been suggested that hermaphrodite individuals in papaya, resulted from a natural mutation in male plants and were likely selected by humans for its favorable fruit phenotype (Ueno et al., 2015). Moreover, the region of the chromosome that produces hermaphrodite papaya plants, which occurs only in cultivated varieties, arose only ∼4,000 years ago (higher posterior density of 95%: 1,400-6,700 years ago; VanBuren et al., 2015). Although crop plant domestication in Mesoamerica occurred around 6,200 years ago, the estimation coincides with the rise of the Maya civilization. Given that hermaphrodite individuals are rarely found in wild populations of Mesoamerica, this strongly supports that papaya was domesticated by the Mayans or other Mesoamerican cultures thousands of years ago. It is believed that the rare allele for hermaphroditism in wild populations was probably selected by early domesticators favoring mutants with red flesh color over the typical yellow/orange color of wild papayas (Manshardt, 2014).
Finally, seeds from wild and domesticated papayas differ in size, germination rate, dormancy, and light sensitivity. Paz and Vázquez-Yanes (1998) found that seeds of wild papayas are smaller than those of cultivated plants in southern Mexico, and that they can remain dormant if buried for long periods of time, reflecting their strong light requirement for germination. In cultivated papaya, seeds germinate with a higher probability, lessening the importance of specific environmental conditions to germinate (Paz and Vázquez-Yanes, 1998). This domestication syndrome has been reported for other crops (Doebley et al., 2006;Gepts, 2010).

GENETICS AND GENOMICS
Understanding crop domestication is crucial to fulfill the demand for improving yield and quality of crops (Tang et al., 2010). Knowledge of the genetic diversity within crop species is essential to understand their origin, domestication and evolutionary relationships, and to efficiently develop strategies for the conservation of their genetic resources, and effective crop improvement (Moore, 2014). Furthermore, crops represent excellent systems for the study of rapid evolution. For the case of papaya, the species now represents an important model in genetic and genomic studies; papaya is one of the first plant species to have its genome sequenced (Ming et al., 2008) and the first transgenic cultivar released for human consumption (Manshardt, 2014). In this section, we briefly review some advances of genetic and genomic methods in papaya; from what is known about the state of wild populations, to recent advances in domesticated varieties. It must be noted, however, that recent specialized reviews about genomics in papaya have been published elsewhere (Ming et al., 2012;Ming and Moore, 2014;Tripathi et al., 2014).
A recurrent consequence of domestication is the reduction of genetic diversity due to the 2-fold effect of genetic drift and selection that operate during the domestication process (Doebley et al., 2006;Gepts, 2014). Assessing the functional variation of wild genetic resources to expand the usable genetic diversity is of utmost importance for breeding and improvement programs (Martínez-Ainsworth and Tenaillon, 2016). However, little information is available about the genetic diversity of wild varieties of important crop species, with relevance for conservation and management (Chávez-Pesqueira and Núñez-Farfán, 2016).
Early on, the use of genetics in tropical fruit crops was confined to the development of isozyme and dominant PCRbased markers and their use for germplasm diversity analysis and clonal fingerprinting; however, information from dominant markers is of limited use in genomic applications (Litz and Padilla, 2012). Since then, codominant markers, known as simple sequence repeat (SSR) markers or microsatellites, have been developed for several tropical fruit crops and used for parentage analysis, clonal fingerprinting, genetic diversity analysis, and development of genetic linkage maps. Nowadays, the increased capability of DNA sequencing by next generation methods, is leading to increased interest in tropical tree fruit crops and to new opportunities to increase the rate of genetic gain in breeding programs (Litz and Padilla, 2012). Recent approaches in massive sequencing tools, such as de novo genome sequencing, whole genome resequencing, reduction of genome complexity using restriction enzymes, transcriptome analysis, and epigenetic studies are becoming of great importance in domestication studies (Guerra-García and Piñero, 2017).
Genetic knowledge of papaya has been accelerated with the advances in molecular markers, linkage and physical maps, comparative genomics studies, and the sequencing of its genome (Tripathi et al., 2014). For genetic diversity studies, some molecular markers have been developed for papaya, mainly SSRs (Ocampo et al., 2006;Eustice et al., 2008;De Oliveira et al., 2010a;Fang et al., 2016). However, most studies have focused on cultivated papaya, with a dearth of studies addressing the genetic diversity and structure of wild populations.
Our recent surveys about the genetic diversity and structure of wild papaya in its natural distribution in northern Mesoamerica (Mexico) (Chávez-Pesqueira et al., 2014; Chávez-Pesqueira and Núñez-Farfán, 2016), have revealed high levels of genetic diversity for the species in the wild. Using both nuclear SSRs and chloroplast DNA markers in 355 individuals of 19 natural populations, we found a mean observed heterozygosity of 0.681 (range of 0.409-0.783) and a mean haplotidic diversity of 0.701 (range of 0.307-0.934) (Chávez-Pesqueira and Núñez-Farfán, 2016). The area with the higher genetic diversity for both markers occur in southeast Mexico, near the region of the Tehuantepec Isthmus. With the chloroplast DNA, we found a lack of phylogeographic structure but a recent structure with the microsatellite data (R ST = 0.149). At the local scale, we also found a negative effect of habitat fragmentation for wild populations of papaya in a fragmented rainforest in Los Tuxtlas rainforest in southern Mexico (Chávez-Pesqueira et al., 2014). Populations that inhabit forest fragments showed a reduced genetic diversity, higher population differentiation, and less migrants. Together, these studies suggest that wild papaya has maintained genetic connectivity among populations throughout time; however, populations are becoming recently structured probably due to human disturbances of its natural habitat, like habitat fragmentation, rendering important conservation concerns for the species in its wild form. Efforts to conserve the natural genetic resources of the species should be addressed as well to assure the conservation of the domesticated varieties through genetic improvement.
Genetic diversity of cultivated and feral papayas has shown to be, in general, lower than wild populations: (1) Using AFLPs markers, a genetic similarity of 0.880 has been reported for 63 cultivated papaya accessions (Kim et al., 2002); (2) an expected heterozygosity (He) ranging from 0.54 to 0.69 in cultivars from the Caribbean region (Ocampo Pérez et al., 2005); (3)  During domestication, gene flow between wild and domesticated conspecific plants plays an important role in addition to other evolutionary forces such as genetic drift and selection. In many cases, crop plants and their wild progenitors belong to the same biological species (Gepts, 2014;Rendón-Anaya et al., 2017), as in papaya, or cohabit with close wild relatives. There is morphological evidence of gene flow between wild and domesticated forms in papaya in southern Mexico (Paz and Vázquez-Yanes, 1998;Chávez-Pesqueira and Núñez-Farfán, 2016; Figure 1C). However, there are no studies evaluating the ecological and evolutionary consequences of this phenomenon. Considering the introduction of transgenic papayas amid the natural distribution of wild populations of C. papaya (see below), would represent a potential ecological risk that demands consideration.
Papaya is known as the world's first transgenic cultivar to have been released for commercial production (in Hawaii, USA), in 1998 by public institutions in the USA (Tecson Mendoza et al., 2008;Manshardt, 2014). This transgenic variety was developed to control the papaya ringspot virus (PRSV). The PRSV acts destroying the photosynthetic capacity of the leaves, leading first to a decrease in fruit quality and yield, loss of vegetative vigor, and eventual death of plants (Lius et al., 1997). After a spread of the virus in 1994 that nearly killed all cultivars in Puna Island in Hawaii, efforts to develop a resistant variety to the PRSV began. Two varieties successfully controlled the virus infection and saved the papaya production in Hawaii: Rainbow and SunUp. According to Gonsalves (2014), transgenic papaya acreage in the year 2012 was about 85% in Hawaii. However, some concern raised regarding the potential effects of gene flow between transgenic and nontransgenic cultivars. Experiments were carried out to assess this; from 2004 through 2009, Gonsalves et al. (2012) studied transgene flow of "Rainbow" (transgenic) to "Kapoho" (nontransgenic) cultivars in Hawaii. Their results showed that observed transgene pollen dispersal to hermaphrodite "Kapoho" from neighboring plants of the "Rainbow" variety was very low, but from "Rainbow" to female "Kapoho" individuals, it was higher (Gonsalves, 2014), confirming the existence of gene flow in some amount among cultivars. Although cultivars of transgenic papaya are confined to Hawaii, little is known about the ecological, evolutionary and commercial consequences of possible transgene flow in the center of origin of papaya where natural populations occur, and where attempts to release transgenic cultivars have prevailed (Silva-Rosales et al., 2010). Nowadays, transgenic papaya from Hawaii is consumed in USA, Canada, and Japan, but no studies have been made to assess the possible movement of seeds to papaya's origin center. Because transgenic papaya saved Hawaii's papaya industry, it has been recognized as a model species of translational biotechnology and the best-characterized commercial transgenic crop. Given that PRSV is widespread in papaya cultivars around the world, papaya SunUp could serve as a transgenic germplasm source to breed suitable virusresistant cultivars (Ming et al., 2008). However, despite these achievements, many countries have rejected transgenic crops; in 2010, only transgenic papaya and plum had been approved for human consumption (Litz and Padilla, 2012).
With the recent advances of DNA next generation sequencing a great interest has been advocated to crops. Papaya is the fifth flowering plant (after Arabidopsis, rice, poplar, and grape) to have its genome sequenced (Kanchana-Udomkan et al., 2014). A draft genome (approximately 3x coverage) was obtained from a transgenic female cultivated papaya from the SunUp variety from Hawaii, using a whole-genome shotgun approach (Ming et al., 2008). The assembled contigs summed up 73% of the genome (271 Mb) with 370 Mb of spanning scaffolds (Wei and Wing, 2008). Further efforts were made to complete the draft genome of papaya; Yu et al. (2009) constructed a BAC (Bacterial Artificial Chromosome) based physical map to support the sequence assembly with an estimated genome coverage of 95.8%. The papaya genome has shown to be three times the size of the Arabidopsis genome, containing fewer genes, with a minimal angiosperm gene set of 13,311. Amplifications in gene number with roles in the evolution of tree-like habit, deposition and remobilization of starch reserves, attraction of seed dispersal agents, and adaptation to tropical day lengths have been identified (Ming et al., 2008).; Moreover, the small genome and diploid nature (2n = 18) of the species, render papaya as an excellent model for genomics (Ming and Moore, 2014;Tripathi et al., 2014).
Recent genetic and genomic methods have revealed interesting facts about domestication in papaya. For instance, Wu et al. (2017) found evidence of a recent emergence of the y allele and its selection in red fruited cultivation. Analysing the levels and patterns of genetic diversity at the CYC-b locus (y allele) and six loci in a 100-kb region flanking this locus, they found evidence of a strong selective sweep, as genetic diversity showed a reduction at the recessive y allele in comparison with the dominant Y allele present in the yellow-fleshed papayas (wild and cultivated). By means of a haplotype network, the authors suggest that the y allele likely originated in the wild and was introduced into the domesticated varieties during the domestication of papaya. However, they found shared haplotype structure among some wild, feral, and cultivated haplotypes near the y allele, suggesting a successive escape of the y allele from red cultivars back into wild populations, probably through feral intermediates in Costa Rica (Wu et al., 2017). In other study, Porter et al. (2009), using a genome-wide analysis, revealed that papaya has relatively few nucleotide-bindings sites (NBS) encoding genes, but structurally diverse showing a novel subgroup. This group of genes are the most common in plant disease resistance proteins. The authors argue that the absence of recent genome duplication and relatively low gene number in papaya, may explain the apparent scarcity NBS gene number, and that this could function as pathogen surveillance, making papaya suitable for functional studies and a better understanding of plant resistance at the genetic level (Porter et al., 2009).
Furthermore, papaya has been recognized as an excellent model for studies about sex determination in plants (Ming et al., 2012;Aryal and Ming, 2014), sex chromosome evolution (Weingartner and Moore, 2012), origin and evolution of dioecy in the Caricaceae family, and the identification of candidate genes and genome-wide DNA markers for papaya improvement (Ming et al., 2012). These studies have helped to expand the knowledge for tropical fruit tree genomics. Moreover, phylogenetically, papaya belongs to the Brassicales, and can serve as an excellent out-group to study genome evolution in the Brassicaceae family, given that it shares a common ancestor with Arabidopsis, about 72 MYA (Ming et al., 2012). Undoubtedly, the sequencing of the papaya genome has opened avenues to investigate the function of many genes for which little is known. The sequenced genome will facilitate knowledge about the genomic regions associated to diverse aspects of plant defense, plant growth, development of leaves, roots, and fruits, flowering, as well as fruit ripening, and circadian clock (Tripathi et al., 2014).
Regarding its commercial use, a variety of markers have been developed to distinguish between the different sexes in juvenile plants (Somsri et al., 1997;Deputy et al., 2002;Lemos et al., 2002;Urasaki et al., 2002; for a complete review of molecular markers in papaya see Kanchana-Udomkan et al., 2014). Most markers achieved to distinguish between male/hermaphroditic and female individuals. This aided in eliminating females, given that hermaphrodites are preferred for their bigger and elongated fruits, rather than the rounded female fruits that require greater container space for shipping (Tripathi et al., 2014). However, in commercial production, males are useless and with these markers, males couldn't be distinguished. Recently, Liao et al. (2017), used a large male-specific retrotransposon insertion of 8396 bp to develop two papaya male-specific markers. Given that sex determination is only possible around 6 months after germination, these markers largely facilitate the elimination of undesired plants. Finally, sex markers are also useful in ecological and evolutionary studies that require the estimation of effective population sizes which can be obtained with the sex ratio of populations (Chávez-Pesqueira et al., 2014).
In spite of these advances, many questions remain unsolved for papaya, mainly for the wild counterpart. For example, the lack of knowledge about the consequences of gene flow from domesticated and/or transgenic populations. Moreover, few studies have evaluated genetic regions of wild individuals (standing variation) for crop improvement purposes (but see Vázquez Calderón et al., 2014).
The use of genomics in tropical crop breeding will greatly assist in solving many difficult challenges in the face of global warming and climate change (Guerra-García and Piñero, 2017). Compared to annual crops and temperate trees, genetic improvement of tropical fruit trees has been limited (Litz and Padilla, 2012). The use of genomics for the development of new tropical cultivars adapted to higher temperatures and increased droughts will become necessary to alleviate food security in the upcoming years (Abberton et al., 2016). The organization of genetic diversity in wild relatives represent important genetic resources that can facilitate the development of climate-tolerant cultivars (Gepts, 2014).
Finally, the continuous advances of genomic tools in crop species will also serve to gain insights into the ecological and evolutionary knowledge of domestication in tropical crops. For instance, the functional genomics aspects of domestication will help assess how many and which genes show differences in expression between wild and domesticated types (Sarah et al., 2017), which factors affect gene expression, and a better understanding of biotic interactions and how these have affected crop evolution, among many others (Gepts, 2014).
Plant domestication represents one of the most relevant events in human history. Today, papaya is one of the most important economically crop species worldwide. Understanding the ecology, evolutionary history and domestication process of such species is necessary to maintain food security in the future and counteract upcoming threats of overpopulation and climate change.

FINAL REMARKS
The understanding of the origin and domestication of crop species has become a topic of great interest and has advanced importantly through combining approaches from diverse disciplines (Martínez-Ainsworth and Tenaillon, 2016). Papaya represents the third most produced crop in the tropics worldwide and an important source of commercial uses for humans. It was most probably originated and domesticated in Mesoamerica where wild populations still occur and distribute in Mexico and Central America. In nature, wild papayas play an important role in the regeneration of their natural habitat and possess high levels of genetic variation, representing the genetic wealth and evolutionary potential of the species. Finally, papaya is an excellent model for genomic studies as it is one of the first plant species to have its complete genome sequenced. However, there is still much more knowledge needed from this important species, mainly about its wild populations and the evolutionary process of domestication.

AUTHOR CONTRIBUTIONS
MC-P and JN-F conceived the review and the outline, and searched for literature. MC-P wrote the manuscript. JN-F reviewed drafts of the manuscript and contributed to writing the final version. Both authors were involved in the final editing and review of the paper.