Search

Transgenic papaya is widely publicized for controlling papaya ringspot virus. However, the impact of particle bombardment on the genome remains unknown. The transgenic SunUp and its progenitor Sunset genomes were assembled into 351.5 and 350.3 Mb in nine chromosomes, respectively. We identified a 1.64 Mb insertion containing three transgenic insertions in SunUp chromosome 5, consisting of 52 nuclear-plastid, 21 nuclear-mitochondrial and 1 nuclear genomic fragments. A 591.9 kb fragment in chromosome 5 was translocated into the 1.64 Mb insertion. We assembled a gapless 9.8 Mb hermaphrodite-specific region of the Y h chromosome and its 6.0 Mb X counterpart. Resequencing 86 genomes revealed three distinct groups, validating their geographic origin and breeding history. We identified 147 selective sweeps and defined the essential role of zeta-carotene desaturase in carotenoid accumulation during domestication. Our findings elucidated the impact of particle bombardment and improved our understanding of sex chromosomes and domestication to expedite papaya improvement., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Papaya is a model system for the study of sex chromosome evolution in plants. However, the cytological structures of the papaya chromosomes remain largely unknown and chromosomal features have not been linked with any genetic or genomic data. We constructed a cytogenetic map of the papaya sex chromosome (chromosome 1) by hybridizing 16 microsatellite markers and 2 cytological feature-associated markers on pachytene chromosomes using fluorescence in situ hybridization (FISH). Except for three markers, the order of the markers was concordant to that of marker loci along the linkage map. This discrepancy was likely caused by skewed segregation in the highly heterochromatic or centromeric regions. The papaya sex chromosome is largely euchromatic, its heterochromatin spans about 15 % of the Y chromosome and is mostly restricted to the centromeric and pericentromeric regions. Analysis of the recombination frequency along the papaya sex chromosome revealed a complete suppression of recombination in the centromere and pericentromere region and 60 % higher recombination rate in the long arm than in the short arm. The uneven distribution of recombination events might be caused by differences in sequence composition. Sequence analysis of 18 scaffolds in total length of 15 Mb revealed higher gene density towards the telomeres and lower gene density towards the centromere, and a relatively higher gene density in the long arm than in the short arm. In an opposite trend, the centromeric and pericentromeric region contained the highest repetitive sequences and the long arm showed the lowest repetitive sequences. This cytogenetic map provides essential information for evolutionary study of sex chromosomes in Caricaceae and will facilitate the analysis of papaya sex chromosomes.

Carotenoid pigments in fruits are indicative of the ripening process and potential nutritional value. Papaya (Carica papaya) fruit flesh color is caused by the accumulation of lycopene or beta-carotenoids in chromoplasts. It is a distinct feature affecting nutritional composition, fruit quality, shelf life, and consumer preference. To uncover the molecular basis of papaya flesh color, we took map-based cloning and candidate gene approaches using integrated genetic and physical maps. A DNA marker tightly linked to flesh color colocalized on a contig of the physical map with a cDNA probe of the tomato (Solanum lycopersicum) chromoplast-specific lycopene beta-cyclase, CYC-b. Candidate gene sequences were obtained from amplified fragments and verified by sequencing two bacterial artificial chromosomes containing the two alleles. Sequence comparison revealed a 2-bp insertion in the coding region of the recessive red flesh allele resulting in a frame-shift mutation and a premature stop codon. A color complementation test in bacteria confirmed that the papaya CpCYC-b is the gene controlling fruit flesh color. Sequence analysis of wild and cultivated papaya accessions showed the presence of this frame-shift mutation in all red flesh accessions examined. Evaluation of DNA markers near CpCYC-b revealed a recombination hot spot, showing that CpCYC-b is located in a gene-rich region with a recombination rate at 3.7 kb per centimorgan, more than 100-fold higher than the genome average at 400 kb per centimorgan. Conserved microsynteny of the CpCYC-b region is indicated by colinearity of two to four genes between papaya, Arabidopsis (Arabidopsis thaliana), grape (Vitis vinifera), and tomato. Our results enhanced our understanding of papaya flesh color inheritance and generated new tools for papaya improvement.

A high-density genetic linkage map of papaya, previously developed using an F2 mapping population derived from the intraspecific cross AU9 x SunUp, was enriched with AFLP markers. The comprehensive genetic map presented here spans 945.2 cM and covers 9 major and 5 minor linkage groups containing 712 SSR, 277 AFLP, and 1 morphological markers. The average marker density for the 9 major linkage groups is 0.9 cM between adjacent markers, and the total number of gaps >5 cM was reduced from 48 to 27 in the current map. AFLPs generated by EcoRI/MseI primer combinations were distributed throughout the 14 linkage groups and resulted in several large locus order rearrangements within the 9 major linkage groups. Integration of AFLP markers provided tighter linkage association between loci, leading to a reduction in map distance on LGs 1, 2, and 4, which were inflated in the previous map, and correction of the marker order on LG8. Suppression of recombination in the male-specific Y region (MSY) of LG1 is further validated by the addition of 27 sex co-segregating AFLP markers. A large region of distorted segregation surrounding the MSY spans 54.4 cM and represents approximately 71% of the linkage group. This comprehensive high-density genetic map provides a framework for mapping quantitative trait loci and for fine mapping as well as for comparative genomic studies of crop plant development and evolution.

Papaya, a fruit crop cultivated in tropical and subtropical regions, is known for its nutritional benefits and medicinal applications. Here we report a 3x draft genome sequence of 'SunUp' papaya, the first commercial virus-resistant transgenic fruit tree to be sequenced. The papaya genome is three times the size of the Arabidopsis genome, but contains fewer genes, including significantly fewer disease-resistance gene analogues. Comparison of the five sequenced genomes suggests a minimal angiosperm gene set of 13,311. A lack of recent genome duplication, atypical of other angiosperm genomes sequenced so far, may account for the smaller papaya gene number in most functional groups. Nonetheless, striking amplifications in gene number within particular functional groups suggest roles in the evolution of tree-like habit, deposition and remobilization of starch reserves, attraction of seed dispersal agents, and adaptation to tropical daylengths. Transgenesis at three locations is closely associated with chloroplast insertions into the nuclear genome, and with topoisomerase I recognition sites. Papaya offers numerous advantages as a system for fruit-tree functional genomics, and this draft genome sequence provides the foundation for revealing the basis of Carica's distinguishing morpho-physiological, medicinal and nutritional properties.

In the ABC model of flower development, B function organ-identity genes act in the second and third whorls of the flower to control petal and stamen identity. The trioecious papaya has male, female, and hermaphrodite flowers and is an ideal system for testing the B-class gene expression patterns in trioecious plants. We cloned papaya B-class genes, CpTM6-1, CpTM6-2, and CpPI, using MADS box gene specific degenerate primers followed by cDNA library screening and sequencing of positive clones. While phylogenetic analyses show that CpPI is the ortholog of the Arabidopsis gene PI, the CpTM6-1 and CpTM6-2 loci are representatives of the paralogous TM6 lineage that contain paleoAP3 motifs unlike the euAP3 gene observed in Arabidopsis. These two paralogs appeared to have originated from a tandem duplication occurred approximately 13.4 million year ago (mya) (bootstrap range 13.36 +/- 2.42). In-situ hybridization and RT-PCR showed that the papaya B-class genes were highly expressed in young flowers across all floral organ primordia. As the flower organs developed, all three B-class genes were highly expressed in petals of all three-sex types and in stamens of hermaphrodite and male flowers. CpTM6-1 expressed at low levels in sepals and carpels, whereas CpTM6-2 expressed at a low level in sepals and at a high level in leaves. Our results showed that B-class gene homologs could function as predicted by the ABC model in trioecous flowers but differential expressions of CpTM6-1, and CpTM6-2, and CpPI suggested the diversification of their functions after the duplication events.

A high-density genetic map of papaya (Carica papaya L.) was constructed using microsatellite markers derived from BAC end sequences and whole-genome shot gun sequences. Fifty-four F(2) plants derived from varieties AU9 and SunUp were used for linkage mapping. A total of 707 markers, including 706 microsatellite loci and the morphological marker fruit flesh color, were mapped into nine major and three minor linkage groups. The resulting map spanned 1069.9 cM with an average distance of 1.5 cM between adjacent markers. This sequence-based microsatellite map resolved the very large linkage group 2 (LG 2) of the previous high-density map using amplified fragment length polymorphism markers. The nine major LGs of our map represent papaya's haploid nine chromosomes with LG 1 of the sex chromosome being the largest. This map validates the suppression of recombination at the male-specific region of the Y chromosome (MSY) mapped on LG 1 and at potential centromeric regions of other LGs. Segregation distortion was detected in a large region on LG 1 surrounding the MSY region due to the abortion of the YY genotype and in a region of LG6 due to an unknown cause. This high-density sequence-tagged genetic map is being used to integrate genetic and physical maps and to assign genome sequence scaffolds to papaya chromosomes. It provides a framework for comparative structural and evolutional genomic research in the order Brassicales.

Papaya (Carica papaya L.) cell wall matrix polysaccharides are modified as the fruit starts to soften during ripening and an endoxylanase is expressed that may play a role in the softening process. Endoxylanase gene expression, protein amount and activity were determined in papaya cultivars that differ in softening pattern and in one cultivar where softening was modified by the ethylene receptor inhibitor 1-methylcyclopropene (1-MCP). Antibodies to the endoxylanase catalytic domain were used to determine protein accumulation. The three papaya varieties used in the study, 'Line 8', 'Sunset', and 'Line 4-16', differed in softening pattern, respiration rate, ethylene production and showed similar parallel relationships during ripening and softening in endoxylanase expression, protein level and activity. When fruit of the three papaya varieties showed the respiratory climacteric and started to soften, the level of endoxylanase gene expression increased and this increase was related to the amount of endoxylanase protein at 32 kDa and its activity. Fruit when treated at less than 10% skin yellow stage with 1-MCP showed a significant delay in the respiratory climacteric and softening, and reduced ethylene production, and when ripe was firmer and had a 'rubbery' texture. The 1-MCP-treated fruit that had the 'rubbery' texture showed suppressed endoxylanase gene expression, protein and enzymatic activity. Little or no delay occurred between endoxylanase gene expression and the appearance of activity during posttranslational processing from 65 to 32 kDa. The close relationship between endoxylanase gene expression, protein accumulation and activity in different varieties and the failure of the 1-MCP-treated fruit to fully soften, supported de novo synthesis of endoxylanase, rapid posttranslation processing and a role in papaya fruit softening.

Genetic relationships among Carica papaya cultivars, breeding lines, unimproved germplasm, and related species were established using amplified fragment length polymorphism (AFLP) markers. Seventy-one papaya accessions and related species were analyzed with nine EcoRI-MseI primer combinations. A total of 186 informative AFLP markers was generated and analyzed. Cluster analysis suggested limited genetic variation in papaya, with an average genetic similarity among 63 papaya accessions of 0.880. Genetic diversity among cultivars derived from the same or similar gene pools was smaller, such as Hawaiian Solo hermaphrodite cultivars and Australian dioecious cultivars with genetic similarity at 0.921 and 0.912, respectively. The results indicated that self-pollinated hermaphrodite cultivars were as variable as open-pollinated dioecious cultivars. Genetic diversity between C. papaya and six other Carica species was also evaluated. Carica papaya shared the least genetic similarity with these species, with an average genetic similarity of 0.432; the average genetic similarity among the six other species was 0.729. The results from AFLP markers provided detailed estimates of the genetic variation within and among papaya cultivars, and supported the notion that C. papaya diverged from the rest of Carica species early in the evolution of this genus.