Next Article in Journal
Pelletization of Compost from Different Mixtures with the Addition of Exhausted Extinguishing Powders
Next Article in Special Issue
The GASA Gene Family in Cacao (Theobroma cacao, Malvaceae): Genome Wide Identification and Expression Analysis
Previous Article in Journal
Exploring Substrate Water Capture in Common Greenhouse Substrates through Preconditioning and Irrigation Pulsing Techniques
Previous Article in Special Issue
Marker-Assisted Backcross Breeding for Improvement of Submergence Tolerance and Grain Yield in the Popular Rice Variety ‘Maudamani’
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

A Review of Omics Technologies and Bioinformatics to Accelerate Improvement of Papaya Traits

by
Rabiatul-Adawiah Zainal-Abidin
1,
Insyirah-Hannah Ruhaizat-Ooi
2 and
Sarahani Harun
2,*
1
Biotechnology and Nanotechnology Research Centre, Malaysian Agricultural Research and Development Institute (MARDI), Serdang 43400, Malaysia
2
Centre for Bioinformatics Research, Institute of Systems Biology (INBIOSIS), Universiti Kebangsaan Malaysia, Bangi 43600, Malaysia
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(7), 1356; https://doi.org/10.3390/agronomy11071356
Submission received: 30 May 2021 / Revised: 16 June 2021 / Accepted: 17 June 2021 / Published: 1 July 2021
(This article belongs to the Special Issue Omics Approaches for Crop Improvement)

Abstract

:
Papaya (Carica papaya) is an economically important fruit crop that is mostly planted in tropical and subtropical regions. Major diseases of papaya, such as the papaya dieback disease (PDD), papaya ringspot virus (PRSV) disease, and papaya sticky disease (PSD), have caused large yield and economic losses in papaya-producing countries worldwide. Postharvest losses have also contributed to the decline in papaya production. Hence, there is an urgent need to secure the production of papaya for a growing world population. Integration of omics resources in crop breeding is anticipated in order to facilitate better-designed crops in the breeding programme. In papaya research, the application of omics and bioinformatics approaches are gradually increased and are underway. Hence, this review focuses on addressing omics technologies and bioinformatics that are used in papaya research. To date, four traits of the papaya have been studied using omics and bioinformatics approaches, which include its ripening process, abiotic stress, disease resistance, and fruit quality (i.e., sweetness, fruit shape, and fruit size). This review also highlights the potential of genetics and genomics data, as well as the systems biology approach that can be applied in a papaya-breeding programme in the near future.

1. Introduction

Papaya (Carica papaya) is one of the tropical fruits that is widely grown in tropical and subtropical countries such as Australia, Brazil, Malaysia, Thailand, and South America. In worldwide fruit production, the papaya is ranked in third place after mango and pineapple fruit production [1]. The papaya has a high nutritional content and medicinal value, which make it widely planted and globally popular. For instance, the fruit is rich in vitamin A, vitamin C, folate, potassium, and niacin, and its leaves, stems, and roots are suitable for alternative medicine [2].
Papaya has nine pairs of chromosomes and is a diploid. The papaya, which has a genome size of 372 Mb, belongs to the Brassicales order and Caricaceae family. Among the well-known papaya cultivars in Hawaii are Solo, Sunrise, SunUp, and Rainbow [3], while Eksotika and Sekaki are well-known cultivars in Malaysia [3]. To date, six papaya cultivars, namely Eksotika, Eksotika II, Sekaki, Three Pillars, Frangi, and Viorica have been registered with the Department of Agriculture, Malaysia [4]. Correspondingly, Japan also produces papaya cultivars, namely Ishigaki Sango and Ishigaki Wondrous [5]; in India, the well-known papaya cultivars are Co1 and Co2 [3].
As a perennial herb, the breeding cycle of the papaya takes about seven to nine months until a ripe fruit is produced. The breeding cycle of the papaya is divided into three stages that are represented by seed germination, flowering, and fruit setting to fruit harvesting (Figure 1) [6,7]. The papaya seed germinates within two to four weeks after sowing [8]. Then, after being transplanted, its vegetative state starts to grow into single-stemmed trees, each bearing a rosette of large, deeply lobed leaves at the apex about two to three months after transplantation [6,7,8]. The papaya plant produces its flowers at the age of three to six months, and produces ripe fruits at seven to nine months [6,7].
Papaya is cultivated in fertile and well-drained soils with a pH of 5.5 to 6.5 [8,9]. Nutrients for the papaya are essential in order to increase the productivity and quality of papaya production. Nitrogen (N) and phosphorus (P) are essential nutrients during the early growth of the papaya in order to ensure an optimum growth of foliage, trunks, and roots, as well as induce higher productivity [10]. Potassium (K) is needed at the papaya fruiting stage because it is essential for improving the fruit quality (i.e., sweetness and pulp thickness) [11,12]. Hence, these agronomic aspects are contributed to ensure that breeding of papaya is successful. The breeding of papaya has produced more papaya cultivars with desired agronomic traits such as high yield, fruit shape, fruit size, and sweetness. However, the papaya is easily exposed to pathogens and being infected by diseases, (i.e., the papaya dieback disease (PDD), papaya ringspot virus (PRSV) disease, and papaya sticky disease (PSD)), which have resulted in a decline in papaya production. In addition, high temperature and water stress in papaya plantation are being affected by climate change, which is a challenge to the papaya industry, especially for its production yield [13,14,15]. The high temperature (i.e., above 35 °C) causes flowers to drop and sex changes in female and hermaphrodite flowers [13]. Although high temperature (i.e., 28 °C) promotes the fast growth of papaya, low pollen viability and early maturation result in imperfect quality fruits and a low yield [13]. Interestingly, high temperature coupled with a higher moisture content produces higher total soluble solids (TSS) in papaya [13]. A sufficient amount of water is essential to papaya because it determines the fruit size and fruit quality [15]. For instance, in dry conditions, the fruits are smaller with a hard texture when ripe.
One major way to overcome the constraints due to climate change and papaya diseases is by breeding for new and improved papaya that has been enhanced with desired traits, such as resistance to disease, resistance to abiotic stress, delayed ripening, and sweetness. Recent trends in crop improvements have shown the integration of omics approaches (i.e., genomics, transcriptomics, proteomics, and metabolomics) and bioinformatics in breeding programmes [16,17,18,19]. The use of omics and bioinformatic approaches in crop breeding helps to obtain a holistic understanding of the genetic and genomic bases of the crop, as well as to understand the molecular interaction among genes, proteins, and metabolites, especially regarding complex traits. Thus, the integration of omics and bioinformatic approaches in crop-breeding programmes is anticipated to facilitate the development of climate-resilient crop varieties and efficient germplasm screening, as well as to accelerate the rate of the genetic gain in a crop [18,19].
To date, papaya-breeding programmes have focused on improving yield and quality, resistance to abiotic stress condition, resistance to disease, as well as delayed ripening. Remarkable success has been achieved in the genetic engineering of papaya: the first transgenic papaya cultivar SunUp developed resistance to papaya ring spot virus (PRSV) disease [20]. In 2008, the first papaya genome, from the cultivar SunUp, was sequenced using whole-genome shotgun Sanger sequencing [21]. The papaya cultivar SunUp is a transgenic papaya and was the first papaya genome sequenced [21]. The SunUp genome sequence was annotated, and yielded 27,793 protein-coding transcripts. To date, SunUp has been used as a reference genome for various comparative genomics analyses of the papaya genome. Previous efforts have exploited the genome sequences, molecular markers, and physical and genetic maps for improvement of papaya traits [22,23,24,25]. The domestication and genetics of the papaya have also been discussed [26]; while Dhekney et al. [8] and Palei et al. [27] have summarised the use of biotechnology tools and the progress of research that has been conducted on papaya, they have not covered all the omics and bioinformatic approaches that have been used in the improvement of papaya traits.
Hence, this review paper aims to examine the recent progress of the improvement of papaya traits using the omics and bioinformatic approaches, and their application in breeding for new and improved papaya cultivars (Figure 2). Each omics approach is shown in Figure 2, outlining how most of the generated molecular or omics data are analysed and visualised using bioinformatics approaches. The integration of systems biology with the analysed omics data analysis will cater to the identification of molecular markers, as well as candidate genes, proteins, and metabolites, for application in papaya-breeding programmes. In this review paper, the future perspective of using omics approaches in improving the desired papaya traits are also discussed.

2. The Role of Bioinformatics in Analysing Omics Data

Bioinformatics is an interdisciplinary field of various research backgrounds that comprises biology, computer science, mathematics, and statistics in order to extract, analyse, integrate, and visualise the biological data that are generated from omics platform technologies [28]. Several bioinformatics tools and biological databases have been developed and are accessible to other researchers in similar fields [19,28]. The development of useful bioinformatics tools is also important for crop breeders in order to inform them of the target selection of traits effortlessly [19,28]. The continuous effort in building and managing the integrated databases would contribute to the translational research in which multiomics data analysis has become an integral part of systems biology.

Comparative Genomics Analysis of Carica papaya

The completion of the SunUp papaya reference genome sequence provides the opportunity for researchers to perform a comparative genomics analysis of the papaya genome. Comparative genomics, which is also known as comparative genome-wide analysis, is an area of study in which the structure and function of genomes from different species or varieties are compared [29]. Using computational approaches, the comparison of genomes and their contents among the cultivars or species enable researchers to identify sequence features (i.e., genes and proteins) that are conserved among the species [30].
Such examples include the comparative genome-wide analysis of the papaya genome in order to identify, annotate, and classify the genes into several gene and transcription-factor families that are associated with abiotic stress [31], disease resistance [24,32], ripening [33,34], and flower development [35]. The identification and annotation of the selected gene families in the papaya genome have also provided information on their gene structure, and phylogenetic tree relationships in relation to orthologs or paralogs from other species. Consequently, the potential genes and transcription factors were validated to further investigate their expression patterns under different conditions [31,32,33,34,35].
The papaya is susceptible to multiple pathogens, such as Erwinia mallotivora [36], ringspot virus [37], and papaya meleira virus (PMeV) [38]. Hence, in order to understand the basis of papaya resistance and susceptibility, the method is to perform a comparative genome-wide analysis of the papaya genome. The genome-wide analysis of papaya sequences has revealed that the papaya has fewer (0.2%) disease-resistance genes (i.e., nucleotide binding site (NBS)-containing R genes) than Arabidopsis thaliana (0.68%) and O. sativa (1.38%), which makes conventional breeding for resistance difficult [24]. Similarly, Praza-Echeverria et al. [32] investigated the potential of NPR1, a pathogenesis-related gene, in papaya resistance against a pathogen by comparing the NPR1 in Arabidopsis and tomato genome sequences. This could lead to the application of the identified gene in genetic engineering for crop improvement.
The understanding of the papaya ripening process is important in order to reduce the occurrence of postharvest losses in the papaya industry. To improve the ripening trait in papaya, it is crucial to identify the genes that are responsible for the ripening process. Coupled with bioinformatics analysis, Liu et al. [33] identified 14 potential SQUAMOSA protein-binding protein-like (SPL) genes, whereas Xu et al. [34] investigated 18 potential auxin/indole-3-acetic acid (Aux/IAA) genes in the papaya ripening process. The involvement of SPL and Aux/IAA gene families in the papaya ripening process is not well known. Hence, these efforts provide an opportunity to understand the roles of auxin-responsive genes and SPL in the papaya ripening process.
Using the genome sequence of the papaya, the basic helix–loop–helix (bHLH) transcription factors that were associated with abiotic stress were identified [31]. A total of 73 candidate genes from the bHLH family were detected using comparative genome-wide analysis. The quantitative real time PCT (qRT-PCR) experiment also revealed the role of candidate bHLHs that might be responsible for abiotic stress responses (i.e., salt, drought, and cold stresses) in papaya. Using a similar approach, Liu et al. [35] identified 11 potential genes in the auxin response factor (ARF) transcription factor family, and their role in papaya flower development.
In summary, a comparative genomics analysis of papaya can be initiated by retrieving the papaya genome sequence from public databases. To date, five databases can be used to retrieve the papaya genome sequence; namely, the NCBI (https://ftp.ncbi.nlm.nih.gov/genomes/all/annotation_releases/3649/100/GCF_000150535.2_Papaya1.0/ accessed on 22 June 2021) [39], PlantGDB [40], Phytozome [41], EnsemblPlants [42], and PLAZA [43]. HMMER [44] and BLAST [45] were developed to perform similarity and homolog sequence analyses. InterProScan [46], Pfam [47], and SMART [48] tools were designed to identify the gene family, motif sequence, and domain association, respectively. PLACE [49] and PlantCARE [50] are tools that are used for the analysis of cis-regulatory elements and promoter identification of the genes of interest. The Gene Structure Display Server (GSDS) [51] annotates gene structures; i.e., exon, intron, 3′UTR, and 5′UTR. To construct a phylogenetic tree, multiple sequence alignment must be performed using the gene sequences containing the domain or motif of interest as input data. ClustalX [52] and MAFFT [53] are multiple sequence alignment programmes that are used for aligning the gene sequences, and MEGA is a tool that is used for constructing a phylogenetic tree [54].
Figure 3 summarises the workflow of the comparative genomics analysis of the papaya gene families using bioinformatic approaches that are related to disease resistance, ripening, flower development, and abiotic stress.

3. Application of Omics Technologies in Carica papaya

The suffix ‘omics’ can be described as the screening of the whole molecular data from a living organism for as genes (genomics), mRNA (transcriptomics), proteins (proteomics), and metabolites (metabolomics). As seen by the expanding number of publications over the years, these omics methods, as well as systems and synthetic biology, are becoming increasingly popular [55]. As compared to other model fruit crops such as the tomato and grape, the multiomics study of papaya is still in its infancy. However, various omics data of papaya have been generated from each omics platform, which can lead to the understanding of the complex phenotypic variations that can facilitate papaya breeding strategies in the improvement of papaya traits (Table 1) [56].

3.1. Genomics and Molecular Markers

Advancements in next-generation sequencing technology (NGS) have resulted in the different types of sequencing platforms (i.e., Illumina, PacBio, and Oxford Nanopore Technologies), which produce high-quality sequences, including longer sequence reads and fewer sequence error rates [56]. In addition, NGS has offered an affordable cost, facilitating the sequencing of the genomes of various plant species and cultivars [57]. For instance, the whole-genome resequencing of the papaya cultivars Eksotika and Sekaki from Malaysia [58], as well as the Sunset cultivar from Hawaii [59], have been carried out using the Illumina platform. The whole genome resequencing of papaya cultivars has led to the identification of a large number of variants that have been annotated into coding genes, where the identified SNPs could be developed as molecular markers in the application of a marker-assisted selection of papaya breeding.
The large volumes of data that have been generated using next-generation sequencing technology have increased the efficiency of the development of molecular markers, such as simple sequence repeat (SSR) and single nucleotide polymorphism (SNP). SSR and SNP are types of molecular markers that are widely used to improve the agronomic traits of fruit trees [5]. Molecular markers, such as restriction fragment length polymorphism (RFLP), amplified fragment length polymorphism (AFLP), and random amplified polymorphic DNA (RAPD), are often limited by low reproducibility, laborious techniques, and time-consuming processes [60]; using SSR and SNP will greatly expedite the detection by using a high-throughput genotyping platform.
Previous efforts have identified 116,43 SSRs in the SunUp genome sequence [61]. In addition, 73 SSRs were found in 25 genes that were related to fruit ripening. A similar study focused on the development of high polymorphic SSR from the whole-genome resequencing of Sunrise Solo (Hawaiian cultivar) and RB2 (Australian cultivar) [62]. The developed SSR markers can be used to differentiate the elite cultivars, and can be used in papaya breeding selection.
Genomics variants such as SNP and insertion–deletion (InDel) are highly abundant in plant genomes [63]. Previous efforts have discovered putative SNPs in the genome-wide analysis of papaya genome resequencing [20,21,48]. Zainal-Abidin et al. [58] sequenced the genome of the papaya cultivars Eksotika and Sekaki in order to identify and annotate SNPs. The identified SNPs will be useful for genotyping to develop a papaya SNP panel. Similarly, the whole-genome resequencing of papaya cultivars such as SunUp (transgenic) and SunSet (nontransgenic) has identified SNPs, InDels, and structural variations [59]. These variants’ data have been used to perform a comparative genomics analysis between transgenic and nontransgenic papaya. Recently, Bohry et al. [64] carried out genome-wide identification of SNP and InDels between the two parental lines of UC10 hybrids, the Formosa elite lines Sekati and JS-12, using the Illumina MiSeq (Illumina Inc., San Diego, CA, USA) platform. Interestingly, the putative variants that were located in the ripening-related genes (RRGs) were suggested to be validated through functional analysis and the genotyping platform. Subsequently, these candidate SNPs that are associated with RRGs can be applied in diversity and genetic-mapping studies, as well as for application in marker-assisted selection (MAS) of papaya.
Recent efforts have been conducted to identify the structural variations (SVs) in the genomes of 25 wild and 42 cultivated papaya [65]. The SVs (i.e., insertions, deletions, inversions, transposable elements, and copy number variations) contain relatively long DNA changes as compared to SNPs and InDels [63,65]. These SV data have been used to unravel the effect of SVs on the papaya phenotype and its adaptation during the domestication process [65]. Detailed GO enrichment between the wild and cultivated papaya has identified genes that are artificially selected during papaya domestication. Environmental adaptability, sexual reproduction, and essential traits such as pistil development, embryonic development, flowering duration, crop yield, pedicel elongation, defence response, and disease response are all influenced by these genes. This study would facilitate the understanding of the genes that are involved during the process of papaya domestication, and provide potential SV data to develop molecular markers in papaya breeding programmes.
The establishment of a bioinformatics pipeline in discovering large genomic variants (i.e., SNP, InDel) has made it possible to unravel the genomic variants in the papaya genome from the various cultivars. The identified and annotated SNPs in disease-resistance and ripening-related genes could be applied in the MAS of the papaya as a tool to aid the selection process in papaya germplasm, as well as in the study of its genetic diversity.
The genetic map is useful for dissecting the genetic components of complex traits [66]. The first genetic map of papaya was developed between Sunrise Solo x Line UH536, which comprised 61 RAPD markers and was distributed in 11 linkage groups (LG) [22]. Then, Deputy et al. [67] developed the genetic map of Kapoho x SunUp that comprised 1498 AFLP in 12 LG; while Chen et al. [23] developed the genetic map of AU9 x SunUp, which comprised 706 SSR in 12 LGs. Blas et al. identified 14 QTL controlling fruit size and shape of papaya [68]. However, progress in the development of papaya QTL has been limited by a lack of genetic and genomic information. Genome-wide identification of SNPs in the papaya genome sequence has allowed the development of a high-density genetic map. This approach has been applied for unravelling the fruit quality of papaya, which has been developed in the F2 population from the RB2 x Sunrise Solo cross using the genotyping-by-sequencing approach (GBS) [69]. The highlights include the candidate genes that are related to the fruit quality traits, such as fruit size, fruit shape, sweetness, length, and firmness. These QTL data that are associated with the genes, and which are related to the fruit quality traits, can be used as candidates for gene exploration in the selection of SNPs and InDels using bioinformatics analysis [64]. Notably, candidate genes with associated SNP markers represent a valuable resource for the future of strategic selective breeding of elite papaya cultivars.
Although the genome sequences of papaya cultivars have been determined, no work has been carried out on the pangenome analysis of C. papaya. Pangenome analysis enables us to capture the entire set of genes from papaya cultivars, as well as to overcome the limitations of relying on a single reference genome [70,71]. The identification of SVs can represent pangenome analysis that would enable us to capture the entire set of genes from papaya cultivars, as well as to overcome the limitations of relying on a single reference genome [70,71]. In addition, more new candidate genes (i.e., genes for disease resistance and ripening) or gene pools can be identified from the wild germplasm and molecular markers, and these can be developed to screen for resistant varieties in the field [70,71].

3.2. Transcriptomics

The transcriptome is defined as the whole set of transcripts in a cell and the quantification of its specific developmental stage or physiological condition [72]. Hence, transcriptome plays a role in estimating the expression of genes, as well as in deciphering the regulation of genes in tissues and organs. Unravelling the transcriptome of tissues and organs of the species that are of interest can be carried out using RNA sequencing technology (RNA-seq), which has been shown to be highly reproducible and enables the simultaneous study of expressed gene samples [73,74,75,76,77,78,79]. These RNA-seq features have made RNA-seq experiments widely used in most transcriptome studies, including the transcriptome analysis of papaya plants.
Previous studies on papaya transcriptomes have focused on the analysis of several papaya traits, such as drought effects [73], fruit quality [74,75], sex determination [76,77], and disease mechanisms [78,79] (Table 1). The molecular response of papaya plants can be observed by analysing the tissue-specific differentially expressed genes (DEGs) using a gene-enrichment analysis. Gene-enrichment analysis using AgriGO [80] or ShinyGO [81] provides information on the biological processes that are regulated in response to the desired traits in papaya.
A transcriptome study of the delayed sticky disease symptoms in papaya has revealed the involvement of stress-responses genes in tolerance mechanisms at the pre-flowering stage [79]. In addition, the authors found that the salicylic-activated genes (i.e., PR1, PR2, PR5, WRKY) contributed to the delayed symptoms, while the activation of candidate genes such as NPR1, UDP-glucuronosyltransferase (UGT), and ethylene limit salicylic acid allowed the PSD symptoms to develop. The nutrient transporter gene family (i.e., nitrate, ammonium, potassium, sodium, phosphate, and sulfate) was the upregulated gene in the host during the infection, and had been shown to act as sensors for plant immunity [78,79].
Shen et al. [75] used RNA-seq technology to elucidate the fruit-colouration process in the ripening condition of papaya. A total of 13 candidate genes, including beta-carotene hydroxylase (CHYB), carotene ε-monooxygenase (LUT1), violaxanthin de-epoxidase (VDE), phytoene synthases (PSY1, PSY2), phytoene desaturases (PDS1, PDS2), zeta-carotene desaturae (ZDS), lycopene cyclases (CYCB, LCYB1, LCYB2, LCYE), and zeaxanthin epoxidase (ZEP), were detected in the papaya fruit transcriptome, which showed that these genes were involved in the carotenoid biosynthetic pathway.
A transcriptome study of the papaya male flower and male-to-hermaphrodite sex-reversal flower demonstrated the involvement of 1756 differentially expressed genes in sex determination [76]. Of these, four papaya homologous genes, including three PIN1 and one PIN3, were found to be upregulated in the male-to-hermaphrodite sexual-reversal flowers. In addition, the authors found the phytohormone-biosynthesis and signal-transduction pathways in the male-to-hermaphrodite sex-reversal flowers. Similarly, Zhou et al. [77] used RNA-seq technology to investigate the underlying mechanism in DNA methylation contribution to the male and female flowers. Dissecting the molecular mechanism that regulates sex expression in the papaya provides valuable resources that facilitate an understanding of sexual differentiation in the papaya.
In the papaya disease-resistance trait, the transcriptomes of SunUp and Sunset leaves were sequenced to investigate the expression changes in both of the papaya cultivars, and were compared to the regulated genes in the resistant and susceptible varieties [78]. This comparative transcriptomic analysis found that there were few disease-resistance and hormone-related genes in SunUp, indicating the PRSV resistance in SunUp transgenic papaya. The finding of this study also provided evidence that genetically modified papaya is not harmful.
Using transcriptome data, gene clusters can be observed by performing a gene coexpression network analysis (Figure 4). This analysis identifies gene–gene interaction, regulation of biological processes, and molecular functions that occur within a cluster of genes. In a transcriptome study of papaya leaves, sap, and roots that were under mild and severe drought, several transcription factors (i.e., WRKY, MYB, bHLH) that were commonly linked to abiotic stress conditions were identified [73]. These transcription factors can be potential regulators in the leaves and roots of papaya plants under drought-stress conditions.

3.3. Proteomics

The two-dimensional differential gel electrophoresis (2D-DIGE), isobaric-tags for relative and absolute quantification (iTRAQ), and matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry (MALDI-TOF/TOF MS) platforms have been widely used to estimate the expressed and abundance of proteins in plants, animals, and humans [82]. In papaya, the use of a proteomics platform has focused on the ripening process [83,84] and disease mechanism [85,86] (Table 1).
The use of a proteomics platform to study the underlying mechanism in papaya fruit ripening was first reported in 2011. Nogueira et al. [83] performed a comparative proteomic analysis of climacteric and preclimacteric papaya cultivars using the 2D-DIGE platform. Several proteins that were closely related to metabolic changes (i.e., cell wall, ethylene biosynthesis, climacteric respiratory burst, stress response, and chromoplast differentiation) in the ripening papaya were found, suggesting that these candidate proteins might be involved in the fruit-ripening process. The 2D-DIGE protein analysis, which was performed in order to identify the differentially expressed proteins during papaya ripening, also suggested the role of 1-methylcyclopropene (1-MCP) in affecting the fruit-ripening process [83]. In addition, several expressed proteins were associated with sugar and cell-wall metabolism, signalling, defence and stress responses, folding and protein stability, and ion channels. Proteins related to pectinesterases, SODs, and the dienelactone hydrolase family deserve further attention, as these proteins might be involved in the ripening processes [84]. In a proteomic study that was applied to different ripening stages in papaya, differential accumulated proteins (DAPs) were quantified using HPLC fractionation and LC-MS/MS [85]. Interestingly, unsaturated fatty acids (i.e., methyl palmitoleate and methyl alpha-octadecatrienoic) were increased during the ripening, indicating that they might be associated with volatile formation in papaya fruit.
Rodrigues et al. [86] used 2-DE and DIGE to identify and distinguish proteins that had accumulated during sticky disease infection, which was caused by the papaya meleira virus (PeMV). A proteomics study of the compatible reaction of C. papaya cv. Eksotika in response to E. mallotivora attack was carried out using an iTRAQ mass spectrometry analysis, and showed differentially expressed proteins that were related to metabolic processes, defence response, and response to stress [87]. Similarly, iTRAQ mass spectrometry was used to quantify the effector proteins in E. mallotivora, and suggested the type III secretion system (T3SS) as an important protein that contributed to the bacterial pathogenicity and virulance [88]. A high-throughput proteomic study using an LC-MS/MS-based label-free proteomics approach was used to assess the protein expression between PMeV-infected preflowering C. papaya and control plants [89]. It was suggested that the increased modulation of photosynthesis, the 26S proteasome, and cell-wall remodelling-associated proteins were involved in the initiation of papaya plant immunity.
In general, identifying differentially expressed proteins provides valuable resources for selecting essential proteins in improving C. papaya agronomic traits (i.e., disease resistance and ripening). To encourage reproducibility of proteomics data, the MS-based proteomics data can be uploaded to a public database such as the ProteomeXchange Consortium [90]. This effort adds additional value of the data in efforts towards the improvement of papaya traits.

3.4. Metabolomics

Metabolomics is the study of small molecules such as metabolites, substrates, and metabolic production in an organ, tissue, or cell [5]. Metabolomics has been used to identify and measure differentially expressed metabolites and to gain an insight into biochemical composition under different environmental conditions [91]. Papaya is rich in secondary metabolites that serve as a source of nutrients for human health. High-performance liquid chromatography or gas chromatography tandem mass spectrometry (LC/GC-MS) and liquid chromatography-quadrupole time-of-flight (LC-QTOF) have been widely used in studies of plant metabolomics [92]. Determining the standard chemical compounds of metabolites and their MS spectrum data is essential in a metabolomic-profiling analysis.
Previous studies performed a comparative analysis of metabolites among papaya cultivars [93,94] (Table 1). The identification of metabolites among these papaya cultivars is important in order to identify the cultivar with the essential metabolites that are associated with quality traits (i.e., sweetness and ripening) that subsequently are to be applied in papaya breeding programmes. For example, carotenoids, tocopherols, and glucosinolates are highly abundant in papaya, and all of them are sources of antioxidants [95,96]. The identification of unique metabolites in a papaya cultivar has highlighted the importance of the chemical marker in authenticating papaya-based food products [94].
Chilling injury is a part of the constraint in maintaining the freshness of papaya, especially for exported papaya. A study by Wu et al. [97] identified different metabolite profiles in the papaya peel at a temperature of 4 °C. The metabolites were associated with aroma traits, such as organic acid, amino acids, hexanal, carbonic acid, pentadecyl propyl ester, and methyl geranate, in papaya peel [97]. The elucidation of the metabolite profile that involves chilling stress at 4 °C can be applied to regulate the storage temperature of the papaya to prevent chilling injury, and to extend its storage period.
Apart from the fruit, studies of metabolite profiling in other parts of the papaya, including its leaves and seeds, have been conducted to detect the metabolites’ phytomedicinal properties [98,99]. Papaya leaves are rich in metabolites (i.e., phenolics, flavonoid, saponins, and tocopherol) that have potential antimicrobial, anticancer, antioxidant, and pesticide properties [100]. Hence, papaya leaves have been widely used in the pharmacological industry for drug development.
Metabolomics that has been combined with other high-throughput omics technologies, such as transcriptomics and proteomics, is referred to as integrated metabolomics, which is sometimes used in studies that are aimed at understanding metabolism as the phenotype of genome function [92]. However, little is known on the integration between metabolomics and transcriptomics to understand the mechanism during the fruit-ripening process. It has been reported that metabolites (i.e., flavonoids, terpenoids, organic acids, phenolic acids, and alkaloids) are closely related to the ripening disorder of fruits [101,102,103]. The elucidation of these potential metabolites during papaya-ripening processes will provide valuable information for developing a strategy for postharvest storage and improving the fruit quality.
Table 1. A summary of the recent omics and bioinformatics approaches that are used in the improvement of papaya traits.
Table 1. A summary of the recent omics and bioinformatics approaches that are used in the improvement of papaya traits.
Type of Omics PlatformTraits/ConditionsDescriptionsApproachReference
Genomics-Whole-genome sequences of papaya cultivar SunUp. Development of first papaya reference genome sequences.Whole-genome shotgun Sanger sequencing[2]
-Whole-genome resequencing of papaya cultivars Eksotika and Sekaki to identify putative SNPs. The identified SNPs between Eksotika and Sekaki located in genes of interest could be suggested for validation using a genotyping platform.Whole-genome resequencing using Illumina HiSeq2000 (Illumina Inc., CA, USA)
and bioinformatic analysis
[58]
-Whole-genome resequencing of papaya cultivar SunUp (transgenic) and Sunset (nontransgenic) to identify SNPs and InDels, and used in comparing transgenic and nontransgenic papaya. The identified SNPs and InDels that were located in high-impact genes could be applied in marker-assisted PRSV disease-resistance breeding in papaya.Whole-genome resequencing using Illumina HiSeq2000 (Illumina Inc., CA, USA) and bioinformatic analysis[59]
-Whole-genome resequencing of wild-type and cultivated papaya to detect structural variations in papaya, and used in understanding the process of papaya domestication.Whole-genome resequencing using Illumina HiSeq2500 (Illumina Inc., CA, USA)
and bioinformatic analysis
[65]
RipeningGene-based SSR marker development focusing on genes related to fruit ripening.Bioinformatics and
genotyping
[61]
Polymorphic SSR marker development for marker-assisted breeding in papaya.Whole-genome resequencing using Illumina HiSeq4000 (Illumina Inc., Foster City, CA, USA),
bioinformatics, and
genotyping
[62]
Genome-wide identification of SNPs and InDels using whole-genome resequencing of two papaya cultivars, namely Sekati and JS-12. The SNPs that were located in RRGs are potential SNPs to be converted in PCR markers, and could be applied in papaya genetic mapping and diversity studies, as well as marker-assisted selection.Whole-genome resequencing using Illumina Miseq (Illumina Inc., Foster City, CA, USA)[64]
Abiotic
stress
Genome-wide analysis of basic helix–loop–helix (bHLH) transcription factors. Candidate bHLH genes that might be responsible for abiotic stress.Comparative genomics
and quantitative real-time PCR (qRT-PCR)
[31]
Disease resistanceGenome-wide analysis of NBS resistance gene family. Candidate resistance (R) genes potentially responsible for disease-resistance mechanism.Comparative genomics
and quantitative real-time PCR (qRT-PCR)
[24]
Disease resistanceGenome-wide analysis of NPR1 family. Candidate pathogenesis-related genes that might be responsible for a disease-resistance mechanism.Comparative genomics
and quantitative real-time PCR (qRT-PCR)
[32]
RipeningGenome-wide analysis of SQUAMOSA promoter binding protein-like gene family in papaya. Candidate ripening- and development-related genes.Comparative genomics
and quantitative real-time PCR (qRT-PCR)
[33]
RipeningGenome-wide analysis of Aux/IAA gene family. Candidate ripening-related genes in papaya.Comparative genomics
and quantitative real-time PCR (qRT-PCR)
[34]
Flower developmentGenome-wide analysis of auxin response factor (ARF) family genes related to flower and fruit development in papaya.
Candidate genes related to flower and fruit development.
Comparative genomics
and Quantitative real-time PCR (qRT-PCR)
[35]
TranscriptomicsDrought toleranceCoexpression network analysis to identify genes and transcription factors related to abiotic stress.Transcriptome sequencing using Illumina NextSeq500 (Illumina Inc., Foster City, CA, USA)
and coexpression network analysis
[73]
Ripening mechanismIdentification of potential regulatory genes
during papaya ripening underlying 1-MCP treatment.
Transcriptome sequencing using Hiseq Xten (Illumina Inc., Foster City, CA, USA)[74]
Fruit colourationIdentification of potential TF regulating the carotenoid biosynthetic pathway.Transcriptome sequencing using
Illumina HiSeq2500
[75]
Sex determinationDifferential expressed genes in sex determination of papaya, in male-to-hermaphrodite and male flowers.(Illumina Inc., Foster City, CA, USA)
Transcriptome sequencing using
Illumina HiSeq2500 (Illumina Inc., Foster City, CA, USA)
[76]
Disease resistanceIdentification of disease-resistance genes in PRSV-resistant and susceptible cultivars.Transcriptome sequencing using
Illumina HiSeq2500 (Illumina Inc., Foster City, CA, USA)
[78]
Disease resistanceIdentification of stress-response genes and nutrient upregulated genes in tolerance mechanism of papaya sticky disease.Transcriptome sequencing using Illumina HiSeq2000 (Illumina Inc., Foster City, CA, USA)[79]
ProteomicsRipening mechanismComparative proteomic analysis of climacteric and preclimacteric papaya cultivars.2-DGE and LC-MS/MS[83]
Ripening mechanismDifferentially expressed proteins during papaya ripening.2-DGE and QTRAP hybrid tandem mass spectrometer[84]
Ripening mechanismDifferentially accumulated proteins (DAPs) during papaya ripening.HPLC and LC-MS/MS[85]
Disease mechanismIdentification of differentially expressed proteins in healthy and PMev disease leaf samples in the Golden cultivar. Metabolism-related proteins were downregulated, and stress-responsive proteins were upregulated.MALDI-TOF-MS/MS and DIGE/LC-IonTrap-MS/MS[86]
Disease mechanismDifferentially expressed proteins of compatible reaction between Eksotika papaya and E. mallotivoraiTRAQ mass spectrometry[87]
Disease mechanismProtein expression between PMeV-infected preflowering C. papaya and control plantsLC-MS/MS-based label-free proteomics[88]
MetabolomicsFruit
ripening
Comparative analysis of metabolite profiling between Eksotika and Sekaki cultivars.GC-MS[93]
Fruit
ripening
Profiling analysis of bioactive and volatile compounds in two papaya cultivars, namely Sel-42 and Tainung.HPLC-ESI-MS/MS[94]
Fruit
ripening
Comparative profiling of carotenoids and volatile in yellow and red flashed between Sui huang and Sui hong cultivars.HPLC-ApCI-MS[95]
Fruit
ripening
Identification of genes and metabolites regulating fruit ripening and softening in papaya cultivar Suiyou-2.Transcriptome sequencing using Illumina Hiseq Xten (Illumina Inc., Foster City, CA, USA) and metabolomics profiling using HPLC-ESI-MS/MS[103]
Chilling
injury
Elucidating of primary metabolites and volatile changes in papaya peel in response to chilling stress.GC-MS /MS[97]
Bioactive propertiesMetabolite profiling in papaya leaves.UPLC-ESI-MS
and GC-MS/MS
[98,99,100]

4. Future Perspective

The progress of the improvement of papaya traits has been limited due to the lack of genetic and genomic information on papaya. The outcome of recent omics studies of papaya plants suggests that there is potential for using these valuable genetic and genomic resources as a breeding tool to improve the desired traits in papaya. Incorporating omics data in papaya breeding programmes with a focus on abiotic stress, disease resistance, delayed ripening, and sweetness offers a promising strategy for developing high-quality traits in papaya cultivars without compromising their yields or agronomic traits. A previous study in Mexico conducted a network analysis of the interaction between viruses in papaya orchards. This viral metagenomics study, which was coupled with a network analysis, could contribute to the understanding of the host–pathogen interactions, which would cater to the management strategies against PRSV and non-PRSV symptoms in papaya [104]. The potential of computational approaches in understanding these biological systems has been employed in crop improvement. The computational models that were constructed integrate genome and phenome information, which led to new experimental strategies in improving crop production [105].
In papaya, pangenome analysis has not received much attention. Using whole-genome resequencing of papaya cultivars from diverse germplasms enables researchers to perform pangenome analyses, which would facilitate the identification of core and variable genes in various papaya genomes. This effort has been performed in soy bean [106] and B. rapa [107]. A further area of interest is to screen for favourable alleles of diverse resistance genes sourced from the wild germplasm of papaya.
Another direction in papaya-trait improvement is employing the genome-editing approach. Functional analysis of candidate genes in the papaya–pathogen system can be performed using the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9(Cas9) system. The CRISPR-Cas9 system was successfully applied in banana [108], apple [109], and kiwifruit [110].

5. Conclusions

The use of omics data plays a role as an advanced breeding tool that will enable the faster and more accurate selection of key consumer-driven traits. Integration of various high-throughput omics platforms may accelerate the research on papaya crop improvement. In addition, the application of computational approaches is key in revealing and filling data gaps, which will be valuable in the designing of new experimentation and measurement strategies that would result in enhanced papaya quality, as well as the ability of it to be sustained under various environmental conditions.

Author Contributions

R.-A.Z.-A. and S.H. conceived the idea; R.-A.Z.-A. wrote the original draft; S.H. and I.-H.R.-O. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Government of Malaysia under the 12th Malaysia Plan Development through Malaysian Agricultural Research and Development Institute (MARDI), grant number (PRB502).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Food and Agriculture Organization. Major Tropical Fruits, Market Review; Food and Agriculture Organization: Rome, Italy, 2020; pp. 1–5. [Google Scholar]
  2. Ming, R.; Yu, Q.; Moore, P.H.; Paull, R.E.; Chen, N.J.; Wang, M.L.; Zhu, Y.J.; Schuler, M.A.; Jiang, J.; Paterson, A.H. Genome of papaya, a fast growing tropical fruit tree. Tree Genet. Genomes 2012, 8, 445–462. [Google Scholar] [CrossRef]
  3. Chan, Y. Breeding papaya (Carica papaya L.). In Breeding Plantation Tree Crops. Tropical Species, 1st ed.; Jain, S.M., Priyadarshan, P.M., Eds.; Springer: New York, NY, USA, 2009; Volume 1, pp. 121–159. [Google Scholar]
  4. Plant Variety Protection Malaysia. Available online: http://pvpbkkt.doa.gov.my/ (accessed on 10 June 2021).
  5. Ogata, T.; Yamanaka, S.; Shoda, M.; Urasaki, N.; Yamamoto, T. Current status of tropical fruit breeding and genetics for three tropical fruit species cultivated in Japan: Pineapple, mango, and papaya. Breed. Sci. 2016, 66, 69–81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Storey, W.B. Genetics of papaya. J. Hered. 1953, 44, 70–78. [Google Scholar] [CrossRef]
  7. Tamaki, M.; Urasaki, N.; Motomura, K.; Nakamura, I.; Adaniya, S. Shortening the breeding cycle of papaya (Carica papaya L.) by culturing embryos treated with ethrel. Plant Cell Tissues Organ Cult. 2011, 106, 225–233. [Google Scholar] [CrossRef]
  8. Dhekney, S.A.; Kandel, R.; Bergey, D.R.; Sitther, V.; Soorianathasundaram, K.; Litz, R.E. Advances in papaya biotechnology. Biocatal. Agric. Biotechnol. 2016, 5, 133–142. [Google Scholar] [CrossRef]
  9. Vos, C.; Arancon, N. Soil and plant nutrient management and fruit production of papaya (Carica papaya) in Keaau, Hawaii. J. Plant Nutr. 2019, 43, 384–395. [Google Scholar] [CrossRef]
  10. Cruz, A.F.; de Oliveira, B.F.; de Carvalho Pires, M. Optimum Level of Nitrogen and Phosphorus to Achieve Better Papaya (Carica papaya var. Solo) Seedlings Growth and Mycorrhizal Colonization. Int. J. Fruit Sci. 2019, 17, 1–11. [Google Scholar] [CrossRef]
  11. Fallas-Corralesa, R.G.; van der Zee, S.E.A.T.M. Diagnosis and management of nutrient constraints in papaya. In Fruit Crops: Diagnosis and Management of Nutrient Constraints, 1st ed.; Srivastava, A.K., Hu, C., Eds.; Elsevier: Amsterdam, The Netherlands, 2020; Volume 1, pp. 607–628. [Google Scholar]
  12. Santos, E.M.C.; Cavalcante, Í.H.L.; da Silva Junior, G.B.; Albano, F.G. Impact of nitrogen and potassium nutrition on papaya (pawpaw) fruit quality. Biosci. J. Uberlândia 2015, 31, 1341–1348. [Google Scholar] [CrossRef] [Green Version]
  13. Dinesh, M.R.; Reddy, B.M.C. Physiological Basis of Growth and Fruit Yield Characteristics of Tropical and Sub-tropical Fruits to Temperature. In Tropical Fruit Tree Species and Climate Change; Sthapit, S.R., Scherr, S.J., Eds.; Bioversity International: Rome, Italy, 2012; Volume 1, pp. 45–70. [Google Scholar]
  14. Salinas, I.; Hueso, J.J.; Cuevas, J. Active Control of Greenhouse Climate Enhances Papaya Growth and Yield at an Affordable Cost. Agronomy 2021, 11, 378. [Google Scholar] [CrossRef]
  15. Ramírez, A.G.; Rodríguez, L.M.P.; Erosa, F.E.; Fuentes, G.; Santamaría, J.M. Identification of the SHINE clade of AP2/ERF domain transcription factors genes in Carica papaya; Their gene expression and their possible role in wax accumulation and water deficit stress tolerance in a wild and a commercial papaya genotypes. Environ. Exp. Bot. 2021, 183, 104341. [Google Scholar] [CrossRef]
  16. Shiratake, K.; Suzuki, M. Omics studies of citrus, grape and rosaceae fruit trees. Breed. Sci. 2016, 66, 122–138. [Google Scholar] [CrossRef] [Green Version]
  17. Ali, M.W.; Borrill, P. Applying genomic resources to accelerate wheat biofortification. Heredity 2020, 125, 386–395. [Google Scholar] [CrossRef]
  18. Roorkiwal, M.; Jain, A.; Kale, S.M.; Doddamani, D.; Chitikineni, A.; Thudi, M.; Varshney, R.K. Development and evaluation of high-density Axiom®CicerSNP Array for high-resolution genetic mapping and breeding applications in chickpea. Plant Biotechnol. J. 2018, 16, 890–901. [Google Scholar] [CrossRef] [Green Version]
  19. Raza, A.; Tabassum, J.; Kudapa, H.; Varshney, R.K. Can omics deliver temperature resilient ready-to- grow crops? Crit. Rev. Biotechnol. 2021, 1–24. [Google Scholar] [CrossRef]
  20. Fitch, M.M.M.; Manshardt, R.M.; Gonsalves, D.; Slightom, J.L.; Sanford, J.C. Virus Resistant Papaya Plants Derived from Tissues Bombarded with the Coat Protein Gene of Papaya Ringspot Virus. Bio Technol. 1992, 10, 1466–1472. [Google Scholar] [CrossRef]
  21. Ming, R.; Hou, S.; Feng, Y.; Yu, Q.; Dionne-Laporte, A.; Saw, J.H.; Senin, P.; Wang, W.; Ly, B.V.; Lewis, K.L.T.; et al. The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature 2008, 452, 991–996. [Google Scholar] [CrossRef] [Green Version]
  22. Sondur, S.N.; Manshardt, R.M.; Stiles, J.I. A genetic map of papaya based on random amplified polymorphic DNA markers. Theor. Appl. Genet. 1996, 93, 547–553. [Google Scholar] [CrossRef]
  23. Chen, C.; Yu, Q.; Hou, S.; Li, Y.; Eustice, M.; Skelton, R.L.; Veatch, O.; Herdes, R.E.; Diebold, L.; Saw, J.; et al. Construction of a sequence-tagged high-density genetic map of papaya for comparative structural and evolutionary genomics in Brassicales. Genetics 2007, 177, 2481–2491. [Google Scholar] [CrossRef] [Green Version]
  24. Porter, B.W.; Paidi, M.; Ming, R.; Alam, M.; Nishijima, W.T.; Zhu, Y.J. Genome-wide analysis of Carica papaya reveals a small NBS resistance gene family. Mol. Genet. Genom. 2009, 281, 609–626. [Google Scholar] [CrossRef]
  25. Yu, Q.; Tong, E.; Skelton, R.L.; Bowers, J.E.; Jones, M.R.; Murray, J.E.; Hou, S.; Guan, P.; Acob, R.A.; Luo, M.C.; et al. A physical map of the papaya genome with integrated genetic map and genome sequence. BMC Genom. 2009, 10, 371. [Google Scholar] [CrossRef] [Green Version]
  26. Chávez-Pesqueira, M.; Núñez-Farfán, J. Domestication and genetics of papaya: A review. Front. Ecol. Evol. 2017, 5. [Google Scholar] [CrossRef] [Green Version]
  27. Palei, S.; Dash, D.K.; Rout, G.R. Biology and biotechnology of papaya, an important fruit crop of tropics: A review. Vegetos 2018, 31, 1–15. [Google Scholar] [CrossRef]
  28. Marsh, J.I.; Hu, H.; Gill, M.; Batley, J.; Edwards, D. Crop breeding for a changing climate: Integrating phenomics and genomics with bioinformatics. Theor. Appl. Genet. 2021, 2021, 1–14. [Google Scholar] [CrossRef]
  29. Ellegren, H. Comparative genomics and the study of evolution by natural selection. Mol. Ecol. 2008, 17, 4586–4596. [Google Scholar] [CrossRef]
  30. Morrell, P.L.; Buckler, E.S.; Ross-Ibarra, J. Crop genomics: Advances and applications. Nat. Rev. Genet. 2012, 13, 85–96. [Google Scholar] [CrossRef]
  31. Yang, M.; Zhou, C.; Yang, H.; Kuang, R.; Huang, B.; Wei, Y. Genome-wide analysis of basic helix-loop-helix transcription factors in papaya (Carica papaya L.). PeerJ 2020, 8, e9319. [Google Scholar] [CrossRef]
  32. Praza-Echeverria, S.; Santamaría, J.M.; Fuentes, G.; De los Angeles Menéndez-Cerón, M.; Vallejo-Reyna, M.A.; Herrera-Valencia, V.A. The NPR1 family of transcription cofactors in papaya: Insights into its structure, phylogeny and expression. Genes Genom. 2012, 34, 379–390. [Google Scholar] [CrossRef]
  33. Liu, K.; Yuan, C.; Feng, S.; Zhong, S.; Li, H.; Zhong, J.; Shen, C.; Liu, J. Genome-wide analysis and characterization of Aux/IAA family genes related to fruit ripening in papaya (Carica papaya L.). BMC Genom. 2017, 18, 1–12. [Google Scholar] [CrossRef]
  34. Xu, Y.; Xu, H.; Wall, M.M.; Yang, J. Roles of transcription factor SQUAMOSA promoter binding protein-like gene family in papaya (Carica papaya) development and ripening. Genomics 2020, 112, 2734–2747. [Google Scholar] [CrossRef]
  35. Liu, K.; Yuan, C.; Li, H.; Lin, W.; Yang, Y.; Shen, C.; Zheng, X. Genome-wide identification and characterization of auxin response factor (ARF) family genes related to flower and fruit development in papaya (Carica papaya L.). BMC Genom. 2015, 16, 901. [Google Scholar] [CrossRef] [Green Version]
  36. Amin, M.N.; Bunawan, H.; Redzuan, R.A.; Jaganath, I.B.S. Erwinia mallotivora sp., a new pathogen of papaya (Carica papaya) in peninsular Malaysia. Int. J. Mol. Sci. 2011, 12, 39–45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Purcifull, D.E.; Edwardson, J.R.; Hiebert, E.; Gonsalves, D. Papaya Ringspot Virus. In CMI/AAB Description of Plant Viruses; Coronel, R.E., Ed.; Wageningen University: Wageningen, The Netherlands, 1984; Volume 2. [Google Scholar]
  38. Maciel-Zambolim, E.; Kunieda-Alonso, S.; Matsuoka, K.; De Carvalho, M.G.; Zerbini, F.M. Purification and some properties of Papaya meleira virus, a novel virus infecting papayas in Brazil. Plant Pathol. 2003, 52, 389–394. [Google Scholar] [CrossRef]
  39. National Center for Botechnology Information (NCBI). Available online: https://ftp.ncbi.nlm.nih.gov/genomes/all/annotation_releases/3649/100/GCF_000150535.2_Papaya1.0/ (accessed on 29 May 2021).
  40. Duvick, J.; Fu, A.; Muppirala, U.; Sabharwal, M.; Wilkerson, M.D.; Lawrence, C.J.; Lushbough, C.; Brendel, V. PlantGDB: A resource for comparative plant genomics. Nucleic Acids Res. 2008, 36, 959–965. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Goodstein, D.M.; Shu, S.; Howson, R.; Neupane, R.; Hayes, R.D.; Fazo, J.; Mitros, T.; Dirks, W.; Hellsten, U.; Putnam, N.; et al. Phytozome: A comparative platform for green plant genomics. Nucleic Acids Res. 2012, 40, D1178–D1186. [Google Scholar] [CrossRef]
  42. Yates, A.; Akanni, W.; Amode, M.R.; Barrell, D.; Billis, K.; Carvalho-Silva, D.; Cummins, C.; Clapham, P.; Fitzgerald, S.; Gil, L.; et al. Ensembl 2016. Nucleic Acids Res. 2016, 44, D710–D716. [Google Scholar] [CrossRef]
  43. Van Bel, M.; Diels, T.; Vancaester, E.; Kreft, L.; Botzki, A.; Van De Peer, Y.; Coppens, F.; Vandepoele, K. PLAZA 4.0: An integrative resource for functional, evolutionary and comparative plant genomics. Nucleic Acids Res. 2018, 46, 1190–1196. [Google Scholar] [CrossRef]
  44. Prakash, A.; Jeffryes, M.; Bateman, A.; Finn, R.D. The HMMER Web Server for Protein Sequence Similarity Search. Curr. Protoc. Bioinform. 2017, 60, 3.15.1–3.15.23. [Google Scholar] [CrossRef]
  45. Altschul, S.F.; Gertz, E.M.; Agarwala, R.; Schäffer, A.A.; Yu, Y.K. PSI-BLAST pseudocounts and the minimum description length principle. Nucleic Acids Res. 2009, 37, 815–824. [Google Scholar] [CrossRef] [Green Version]
  46. Jones, P.; Binns, D.; Chang, H.Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [Green Version]
  47. El-Gebali, S.; Mistry, J.; Bateman, A.; Eddy, S.R.; Luciani, A.; Potter, S.C.; Qureshi, M.; Richardson, L.J.; Salazar, G.A.; Smart, A.; et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2019, 47, D427–D432. [Google Scholar] [CrossRef]
  48. Letunic, I.; Bork, P. 20 years of the SMART protein domain annotation resource. Nucleic Acids Res. 2018, 46, D493–D496. [Google Scholar] [CrossRef]
  49. Higo, K.; Ugawa, Y.; Iwamoto, M.; Korenaga, T. Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999, 27, 297–300. [Google Scholar] [CrossRef] [Green Version]
  50. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van De Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef]
  51. Hu, B.; Jin, J.; Guo, A.Y.; Zhang, H.; Luo, J.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef] [Green Version]
  52. Thompson, J.D.; Gibson, T.J.; Higgins, D.G. Multiple Sequence Alignment Using ClustalW and ClustalX. Curr. Protoc. Bioinform. 2003, 2.3.1–2.3.22. [Google Scholar] [CrossRef]
  53. Katoh, K.; Misawa, K.; Kuma, K.I.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [Green Version]
  54. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  55. Aizat, W.M.; Ismail, I.; Noor, N.M. Recent Development in Omics Studies. In Omics Applications for Systems Biology, 1st ed.; Aizat, W., Goh, H.H., Baharum, S., Eds.; Springer International Publishing: Cham, Switzerland, 2018; Volume 1102, pp. 1–9. [Google Scholar] [CrossRef]
  56. Scossa, F.; Alseekh, S.; Fernie, A.R. Integrating multi-omics data for crop improvement. J. Plant Physiol. 2021, 257, 1–16. [Google Scholar] [CrossRef]
  57. Choi, H.K. Translational genomics and multi-omics integrated approaches as a useful strategy for crop breeding. Genes Genom. 2019, 41, 133–146. [Google Scholar] [CrossRef] [Green Version]
  58. Zainal-Abidin, R.; Abu-Bakar, N.; Yusof, M.M.F.; Abdullah, N. Sequence information on single nucleotide polymorphism (SNP) through genome sequencing analysis of Carica papaya variety Eksotika and Sekaki. J. Trop. Agric. Food Sci. 2016, 44, 219–228. [Google Scholar]
  59. Fang, J.; Fang, J.; Fang, J.; Fang, J.; Wood, A.M.; Chen, Y.; Chen, Y.; Yue, J.; Ming, R.; Ming, R. Genomic variation between PRSV resistant transgenic SunUp and its progenitor cultivar Sunset. BMC Genom. 2020, 21. [Google Scholar] [CrossRef]
  60. Nadeem, M.A.; Nawaz, M.A.; Shahid, M.Q.; Doğan, Y.; Comertpay, G.; Yıldız, M.; Hatipoğlu, R.; Ahmad, F.; Alsaleh, A.; Labhane, N.; et al. DNA molecular markers in plant breeding: Current status and recent advancements in genomic selection and genome editing. Biotechnol. Biotechnol. Equip. 2018, 32, 261–285. [Google Scholar] [CrossRef] [Green Version]
  61. Vidal, N.M.; Grazziotin, A.L.; Ramos, H.C.C.; Pereira, M.G.; Venancio, T.M. Development of a gene-centered SSR atlas as a resource for papaya (Carica papaya) marker-assisted selection and population genetic studies. PLoS ONE 2014, 9, e112654. [Google Scholar] [CrossRef] [Green Version]
  62. Nantawan, U.; Kanchana-udomkan, C.; Drew, R.; Ford, R. Development of polymorphic simple sequence repeat (SSR) markers from genome re-sequencing of Carica papaya L. ‘Sunrise Solo’ and ‘RB2′ for marker-assisted breeding. Acta Hortic. 2018, 1205, 687–695. [Google Scholar] [CrossRef] [Green Version]
  63. Saxena, R.K.; Edwards, D.; Varshney, R.K. Structural variations in plant genomes. Brief. Funct. Genom. Proteom. 2014, 13, 296–307. [Google Scholar] [CrossRef] [Green Version]
  64. Bohry, D.; Ramos, H.C.C.; dos Santos, P.H.D.; Boechat, M.S.B.; Arêdes, F.A.S.; Pirovani, A.A.V.; Pereira, M.G. Discovery of SNPs and InDels in papaya genotypes and its potential for marker assisted selection of fruit quality traits. Sci. Rep. 2021, 11, 292. [Google Scholar] [CrossRef]
  65. Liao, Z.; Zhang, X.; Zhang, S.; Lin, Z.; Zhang, X.; Ming, R. Structural variations in papaya genomes. BMC Genom. 2021, 22, 335. [Google Scholar] [CrossRef]
  66. Collard, B.C.Y.; Jahufer, M.Z.Z.; Brouwer, J.B.; Pang, E.C.K. An introduction to markers, quantitative trait loci (QTL) mapping and marker-assisted selection for crop improvement: The basic concepts. Euphytica 2005, 142, 169–196. [Google Scholar] [CrossRef]
  67. Deputy, J.C.; Ming, R.; Ma, H.; Liu, Z.; Fitch, M.M.M.; Wang, M.; Manshardt, R.; Stiles, J.I. Molecular markers for sex determination in papaya (Carica papaya L.). Theor. Appl. Genet. 2002, 106, 107–111. [Google Scholar] [CrossRef]
  68. Blas, A.L.; Yu, Q.; Veatch, O.J.; Paull, R.E.; Moore, P.H.; Ming, R. Genetic mapping of quantitative trait loci controlling fruit size and shape in papaya. Mol. Breed. 2012, 29, 457–466. [Google Scholar] [CrossRef]
  69. Nantawan, U.; Kanchana-Udomkan, C.; Bar, I.; Ford, R. Linkage mapping and quantitative trait loci analysis of sweetness and other fruit quality traits in papaya. BMC Plant Biol. 2019, 19, 449. [Google Scholar] [CrossRef] [PubMed]
  70. Tao, Y.; Zhao, X.; Mace, E.; Henry, R.; Jordan, D. Exploring and Exploiting Pan-genomics for Crop Improvement. Mol. Plant 2019, 12, 156–169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Neik, T.X.; Amas, J.; Barbetti, M.; Edwards, D.; Batley, J. Understanding host–pathogen interactions in brassica napus in the omics era. Plants 2020, 9, 1336. [Google Scholar] [CrossRef] [PubMed]
  72. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009, 10, 57–63. [Google Scholar] [CrossRef]
  73. Gamboa-Tuz, S.D.; Pereira-Santana, A.; Zamora-Briseño, J.A.; Castano, E.; Espadas-Gil, F.; Ayala-Sumuano, J.T.; Keb-Llanes, M.Á.; Sanchez-Teyer, F.; Rodríguez-Zapata, L.C. Transcriptomics and co-expression networks reveal tissue-specific responses and regulatory hubs under mild and severe drought in papaya (Carica papaya L.). Sci. Rep. 2018, 8, 14539. [Google Scholar] [CrossRef]
  74. Zhu, X.; Ye, L.; Ding, X.; Gao, Q.; Xiao, S.; Tan, Q.; Huang, J.; Chen, W.; Li, X. Transcriptomic analysis reveals key factors in fruit ripening and rubbery texture caused by 1-MCP in papaya. BMC Plant Biol. 2019, 19, 1–23. [Google Scholar] [CrossRef]
  75. Shen, Y.H.; Lu, B.G.; Feng, L.; Yang, F.Y.; Geng, J.J.; Ming, R.; Chen, X.J. Isolation of ripening-related genes from ethylene/1-MCP treated papaya through RNA-seq. BMC Genom. 2017, 18, 1–13. [Google Scholar] [CrossRef]
  76. Lin, H.; Liao, Z.; Zhang, L.; Yu, Q. Transcriptome analysis of the male-to-hermaphrodite sex reversal induced by low temperature in papaya. Tree Genet. Genomes 2016, 12, 94. [Google Scholar] [CrossRef]
  77. Zhou, P.; Zhang, X.; Fatima, M.; Ma, X.; Fang, H.; Yan, H.; Ming, R. DNA methylome and transcriptome landscapes revealed differential characteristics of dioecious flowers in papaya. Hortic. Res. 2020, 7, 81. [Google Scholar] [CrossRef]
  78. Fang, J.; Lin, A.; Qiu, W.; Cai, H.; Umar, M.; Chen, R.; Ming, R. Transcriptome profiling revealed stress-induced and disease resistance genes up-regulated in PRSV resistant transgenic papaya. Front. Plant Sci. 2016, 7, 855. [Google Scholar] [CrossRef] [Green Version]
  79. Madroñero, J.; Rodrigues, S.P.; Antunes, T.F.S.; Abreu, P.M.V.; Ventura, J.A.; Fernandes, A.A.R.; Fernandes, P.M.B. Transcriptome analysis provides insights into the delayed sticky disease symptoms in Carica papaya. Plant Cell Rep. 2018, 37, 967–980. [Google Scholar] [CrossRef]
  80. Tian, T.; Liu, Y.; Yan, H.; You, Q.; Yi, X.; Du, Z.; Xu, W.; Su, Z. AgriGO v2.0: A GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 2017, 45, W122–W129. [Google Scholar] [CrossRef]
  81. Ge, S.X.; Jung, D.; Yao, R. ShinyGO: A graphical gene-set enrichment tool for animals and plants. Bioinformatics 2020, 36, 2628–2629. [Google Scholar] [CrossRef]
  82. Aslam, B.; Basit, M.; Nisar, M.A.; Khurshid, M.; Rasool, M.H. Proteomics: Technologies and their applications. J. Chromatogr. Sci. 2017, 55, 182–196. [Google Scholar] [CrossRef] [Green Version]
  83. Nogueira, S.B.; Labate, C.A.; Gozzo, F.C.; Pilau, E.J.; Lajolo, F.M.; Oliveira do Nascimento, J.R. Proteomic analysis of papaya fruit ripening using 2DE-DIGE. J. Proteom. 2011, 75, 1428–1439. [Google Scholar] [CrossRef]
  84. Huerta-Ocampo, J.Á.; Osuna-Castro, J.A.; Lino-López, G.J.; Barrera-Pacheco, A.; Mendoza-Hernández, G.; De León-Rodríguez, A.; Barba de la Rosa, A.P. Proteomic analysis of differentially accumulated proteins during ripening and in response to 1-MCP in papaya fruit. J. Proteom. 2012, 75, 2160–2169. [Google Scholar] [CrossRef]
  85. Jiang, B.; Ou, S.; Xu, L.; Mai, W.; Ye, M.; Gu, H.; Zhang, T.; Yuan, C.; Shen, C.; Wang, J.; et al. Comparative proteomic analysis provides novel insights into the regulation mechanism underlying papaya (Carica papaya L.) exocarp during fruit ripening process. BMC Plant Biol. 2019, 19. [Google Scholar] [CrossRef]
  86. Rodrigues, S.P.; Ventura, J.A.; Aguilar, C.; Nakayasu, E.S.; Almeida, I.C.; Fernandes, P.M.B.; Zingali, R.B. Proteomic analysis of papaya (Carica papaya L.) displaying typical sticky disease symptoms. Proteomics 2011, 11, 2592–2602. [Google Scholar] [CrossRef]
  87. Abu Bakar, N.; Badrun, R.; Shaharuddin, N.A.; Labuh, R. iTRAQ Proteins Analysis of Early Infected Papaya Plants with Papaya Dieback Pathogen. Asian J. Plant Biol. 2015, 3, 1–7. [Google Scholar]
  88. Abu Bakar, N.; Badrun, R.; Rozano, L.; Ahmad, L.; Raih, M.M.; Tarmizi, A.A. Identification and validation of putative Erwinia mallotivora effectors via quantitative proteomics and Real Time Analysis. J. Agric. Food. Technol. 2017, 7, 10–21. [Google Scholar]
  89. Soares, E.A.; Werth, E.G.; Madroñero, L.J.; Ventura, J.A.; Rodrigues, S.P.; Hicks, L.M.; Fernandes, P.M.B. Label-free quantitative proteomic analysis of pre-flowering PMeV-infected Carica papaya L. J. Proteom. 2016, 151, 275–283. [Google Scholar] [CrossRef] [PubMed]
  90. Vizcaíno, J.A.; Deutsch, E.W.; Wang, R.; Csordas, A.; Reisinger, F.; Ríos, D.; Dianes, J.A.; Sun, Z.; Farrah, T.; Bandeira, N.; et al. ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat. Biotechnol. 2014, 32, 223–226. [Google Scholar] [CrossRef] [PubMed]
  91. Kusano, M.; Yang, Z.; Okazaki, Y.; Nakabayashi, R.; Fukushima, A.; Saito, K. Using metabolomic approaches to explore chemical diversity in rice. Mol. Plant 2015, 8, 58–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Okazaki, Y.; Saito, K. Integrated metabolomics and phytochemical genomics approaches for studies on rice. Gigascience 2016, 5, 1–8. [Google Scholar] [CrossRef] [Green Version]
  93. Sanimah, S.; Sarip, J. Metabolomic analysis of Carica papaya variety Eksotika and Sekaki. J. Trop. Agric. Food Sci. 2015, 43, 103–117. [Google Scholar]
  94. Kelebek, H.; Selli, S.; Gubbuk, H.; Gunes, E. Comparative evaluation of volatiles, phenolics, sugars, organic acids and antioxidant properties of Sel-42 and Tainung papaya varieties. Food Chem. 2015, 173, 912–919. [Google Scholar] [CrossRef]
  95. Jing, G.; Li, T.; Qu, H.; Yun, Z.; Jia, Y.; Zheng, X.; Jiang, Y. Carotenoids and volatile profiles of yellow- and red-fleshed papaya fruit in relation to the expression of carotenoid cleavage dioxygenase genes. Postharvest Biol. Technol. 2015, 109, 114–119. [Google Scholar] [CrossRef]
  96. Santana, L. Nutraceutical Potential of Carica papaya in Metabolic syndrome. Nutrients 2019, 5, 1608. [Google Scholar] [CrossRef] [Green Version]
  97. Wu, Q.; Li, Z.; Chen, X.; Yun, Z.; Li, T.; Jiang, Y. Comparative metabolites profiling of harvested papaya (Carica papaya L.) peel in response to chilling stress. J. Sci. Food Agric. 2019, 99, 6868–6881. [Google Scholar] [CrossRef]
  98. Gogna, N.; Hamid, N.; Dorai, K. Metabolomic profiling of the phytomedicinal constituents of Carica papaya L. leaves and seeds by 1H NMR spectroscopy and multivariate statistical analysis. J. Pharm. Biomed. Anal. 2015, 115, 74–85. [Google Scholar] [CrossRef]
  99. Harini, R.S.; Saivikashini, A.; Keerthana, G. Profiling metabolites of Carica papaya Linn. variety CO7 through GC-MS analysis. J. Pharmacogn. Phytochem. 2016, 5, 200–203. [Google Scholar]
  100. Vuong, Q.V.; Hirun, S.; Roach, P.D.; Bowyer, M.C.; Phillips, P.A.; Scarlett, C.J. Effect of extraction conditions on total phenolic compounds and antioxidant activities of Carica papaya leaf aqueous extracts. J. Herb. Med. 2013, 3, 104–111. [Google Scholar] [CrossRef]
  101. Oikawa, A.; Otsuka, T.; Nakabayashi, R.; Jikumaru, Y.; Isuzugawa, K.; Murayama, H.; Saito, K.; Shiratake, K. Metabolic profiling of developing pear fruits reveals dynamic variation in primary and secondary metabolites, including plant hormones. PLoS ONE 2015, 10, 1–18. [Google Scholar] [CrossRef]
  102. Yun, Z.; Li, T.; Gao, H.; Zhu, H.; Gupta, V.K.; Jiang, Y. Integrated transcriptomic, proteomic and metabolomics analysis reveals peel ripening of harvested banana under natural condition. Biomolecules 2019, 9, 167. [Google Scholar] [CrossRef] [Green Version]
  103. Zheng, S.; Hao, Y.; Fan, S.; Cai, J.; Chen, W.; Li, X.; Zhu, X. Metabolomic and transcriptomic profiling provide novel insights into fruit ripening and ripening disorder caused by 1-MCP treatments in papaya. Int. J. Mol. Sci. 2021, 22, 916. [Google Scholar] [CrossRef]
  104. Alcalá-Briseño, R.I.; Casarrubias-Castillo, K.; López-Ley, D.; Garrett, K.A.; Silva-Rosales, L. Network analysis of the papaya orchard virome from two agroecological regions of Chiapas, Mexico. mSystems 2020, 5, 1–17. [Google Scholar] [CrossRef] [Green Version]
  105. Benes, B.; Guan, K.; Lang, M.; Long, S.P.; Lynch, J.P.; Marshall-Colón, A.; Peng, B.; Schnable, J.; Sweetlove, L.J.; Turk, M.J. Multiscale computational models can guide experimentation and targeted measurements for crop improvement. Plant J. 2020, 103, 21–31. [Google Scholar] [CrossRef]
  106. Li, Y.H.; Zhou, G.; Ma, J.; Jiang, W.; Jin, L.G.; Zhang, Z.; Guo, Y.; Zhang, J.; Sui, Y.; Zheng, L.; et al. De novo assembly of soybean wild relatives for pan-genome analysis of diversity and agronomic traits. Nat. Biotechnol. 2014, 32, 1045–1052. [Google Scholar] [CrossRef] [Green Version]
  107. Golicz, A.A.; Bayer, P.E.; Barker, G.C.; Edger, P.P.; Kim, H.R.; Martinez, P.A.; Chan, C.K.K.; Severn-Ellis, A.; McCombie, W.R.; Parkin, I.A.P.; et al. The pangenome of an agronomically important crop plant Brassica oleracea. Nat. Commun. 2016, 7, 1–8. [Google Scholar] [CrossRef] [PubMed]
  108. Kaur, N.; Alok, A.; Shivani; Kaur, N.; Pandey, P.; Awasthi, P.; Tiwari, S. CRISPR/Cas9-mediated efficient editing in phytoene desaturase (PDS) demonstrates precise manipulation in banana cv. Rasthali genome. Funct. Integr. Genom. 2018, 18, 89–99. [Google Scholar] [CrossRef]
  109. Charrier, A.; Vergne, E.; Dousset, N.; Richer, A.; Petiteau, A.; Chevreau, E. Efficient targeted mutagenesis in apple and first time edition of pear using the CRISPR-Cas9 system. Front. Plant Sci. 2019, 10, 40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  110. Wang, Z.; Wang, S.; Li, D.; Zhang, Q.; Li, L.; Zhong, C.; Liu, Y.; Huang, H. Optimized paired-sgRNA/Cas9 cloning and expression cassette triggers high-efficiency multiplex genome editing in kiwifruit. Plant Biotechnol. J. 2018, 16, 1424–1433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The breeding cycle of papaya. DAG denotes day after germination.
Figure 1. The breeding cycle of papaya. DAG denotes day after germination.
Agronomy 11 01356 g001
Figure 2. A summary of the omics approaches that are used in the improvement of papaya traits. Four traits have been studied using omics and bioinformatic approaches, including ripening, abiotic stress, disease resistance, and fruit quality.
Figure 2. A summary of the omics approaches that are used in the improvement of papaya traits. Four traits have been studied using omics and bioinformatic approaches, including ripening, abiotic stress, disease resistance, and fruit quality.
Agronomy 11 01356 g002
Figure 3. Workflow for comparative genomics analysis of the papaya gene families.
Figure 3. Workflow for comparative genomics analysis of the papaya gene families.
Agronomy 11 01356 g003
Figure 4. Schematic representation of the papaya RNA-seq analysis using a bioinformatics approach. The potential genes and transcription factor can be validated using qRT-PCR detection, which can be applied in precision breeding through gene editing or marker-assisted selection (MAS).
Figure 4. Schematic representation of the papaya RNA-seq analysis using a bioinformatics approach. The potential genes and transcription factor can be validated using qRT-PCR detection, which can be applied in precision breeding through gene editing or marker-assisted selection (MAS).
Agronomy 11 01356 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zainal-Abidin, R.-A.; Ruhaizat-Ooi, I.-H.; Harun, S. A Review of Omics Technologies and Bioinformatics to Accelerate Improvement of Papaya Traits. Agronomy 2021, 11, 1356. https://doi.org/10.3390/agronomy11071356

AMA Style

Zainal-Abidin R-A, Ruhaizat-Ooi I-H, Harun S. A Review of Omics Technologies and Bioinformatics to Accelerate Improvement of Papaya Traits. Agronomy. 2021; 11(7):1356. https://doi.org/10.3390/agronomy11071356

Chicago/Turabian Style

Zainal-Abidin, Rabiatul-Adawiah, Insyirah-Hannah Ruhaizat-Ooi, and Sarahani Harun. 2021. "A Review of Omics Technologies and Bioinformatics to Accelerate Improvement of Papaya Traits" Agronomy 11, no. 7: 1356. https://doi.org/10.3390/agronomy11071356

APA Style

Zainal-Abidin, R. -A., Ruhaizat-Ooi, I. -H., & Harun, S. (2021). A Review of Omics Technologies and Bioinformatics to Accelerate Improvement of Papaya Traits. Agronomy, 11(7), 1356. https://doi.org/10.3390/agronomy11071356

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop