A chromosome-level genome assembly provides insights into ascorbic acid accumulation and fruit softening in guava (Psidium guajava)

Development of genome-wide SSR markers through in silico mining of guava (Psidium guajava L.) genome for genetic diversity analysis and transferability studies across species and genera

Guava (Psidium guajava L.) is one of the economically major fruit crops, abundant in nutrients and found growing in tropical-subtropical regions around the world. Ensuring sufficient genomic resources is crucial for crop species to enhance breeding efficiency and facilitate molecular breeding. However, genomic resources, especially microsatellite or simple sequence repeat (SSR) markers, are limited in guava. Therefore, novel genome-wide SSR markers were developed by utilizing chromosome assembly (GCA_016432845.1) of the "New Age" cultivar through GMATA, a comprehensive software. The software evaluated about 397.8 million base pairs (Mbp) of the guava genome sequence, where 87,372 SSR loci were utilized to design primers, ultimately creating 75,084 new SSR markers. After in silico analysis, a total of 75 g-SSR markers were chosen to screen 35 guava genotypes, encompassing wild Psidium species and five jamun genotypes. Of the 72 amplified novel g-SSR markers (FHTGSSRs), 53 showed polymorphism, suggesting significant genetic variation among the guava genotypes, including wild species. The 53 polymorphic g-SSR markers had an average of 3.04 alleles per locus for 35 selected guava genotypes. Besides, in this study, the mean values recorded for major allele frequency, gene diversity, observed heterozygosity, and polymorphism information content were 0.73, 0.38, 0.13, and 0.33, respectively. Among the wild Psidium species studied, the transferability of these novel g-SSR loci across different species was found to be 45.83% to 90.28%. Furthermore, 17 novel g-SSR markers were successfully amplified in all the selected Syzygium genotypes, of which only four markers could differentiate between two Syzygium species. A neighbour-joining (N-J) tree was constructed using 53 polymorphic g-SSR markers and classified 35 guava genotypes into four clades and one outlier, emphasizing the genetic uniqueness of wild Psidium species compared to cultivated genotypes. Model-based structure analysis divided the guava genotypes into two distinct genetic groups, a classification that was strongly supported by Principal Coordinate Analysis (PCoA). In addition, the AMOVA and PCoA analyses also indicated substantial genetic diversity among the selected guava genotypes, including wild Psidium species. Hence, the developed novel genome-wide genomic SSRs could enhance the availability of genomic resources and assist in the molecular breeding of guava.

Comprehensive insights on association mapping in perennial fruit crops breeding - Its implications, current status and future perspectives

In order to provide food and nutritional security for the world's rapidly expanding population, fruit crop researchers have identified two critical priorities: increasing production and preserving fruit quality during the pre- and post-harvest periods. The genetic basis of these complex, commercially important fruit traits which are uniquely regulated by polygenes or multi-allelic genes that interact with one another and the environment can be analyzed with the aid of trait mapping tools. The most interesting trait mapping approach that offers the genetic level investigation for marker-trait associations (MTAs) for these complex fruit traits, without the development of mapping population, is association mapping. This approach was used during the genetic improvement program, emphasizing the obstacles (breeding strategies adopted, generation interval, and their genomic status) pertaining to perennial fruit crops. This method of studying population diversity and linkage disequilibrium in perennial fruit crops has been made possible by recent developments in genotyping, phenotyping, and statistical analysis. Thus, the purpose of this review is to provide an overview of different trait mapping techniques, with a focus on association mapping (method, essential components, viability, constraints, and future perspective) and its advantages, disadvantages, and possibilities for breeding perennial fruit crops.

Characterization of the aroma and flavor profiles of guava fruit (Psidium guajava) during developing by HS-SPME-GC/MS and RNA sequencing

The flavor of guava, an important tropical fruit, is influenced by secondary metabolites. However, the mechanisms and processes underlying flavor formation in guava remain unclear. In this study, dynamic changes in volatile organic compounds (VOCs), sugars, and organic acids in guava peel and flesh across different developmental stages were investigated using headspace solid-phase microextraction (HS-SPME) combined with gas chromatography-mass spectrometry (GC-MS) and high-performance liquid chromatography (HPLC). Here, we identified 90 VOCs, three sugars and eight organic acids. The dynamics of VOCs differ between the flesh and peel. The early developmental stages are more critical in influencing the variation of VOCs in the flesh, while VOC changes in peel occur more progressively across the developmental stages. By integrating transcriptomic and metabolomic analyses, we identified several key genes involved in VOC, sugar, and acid metabolism. This is the first study to describe the expression patterns of these genes throughout guava development, providing new insights into guava flavor development.

The Telomere-to-Telomere Genome of Jaboticaba Reveals the Genetic Basis of Fruit Color and Citric Acid Content

Jaboticaba is a typical tropical plant that blossoms and bears fruit on the tree trunks and branches. The fruits resemble grapes in appearance and texture and are also known as "treegrapes". Currently, research on the genomics of jaboticaba is lacking. In this study, we constructed an integrated, telomere-to-telomere (T2T) gap-free reference genome and two nearly complete haploid genomes, thereby providing a high-quality genomic resource. Furthermore, we unveiled the evolutionary history of several species within the Myrtaceae family, highlighting significant expansions in metabolic pathways such as the citric acid cycle, glycolysis/gluconeogenesis, and phenylpropanoid biosynthesis throughout their evolutionary process. Transcriptome analysis of jaboticaba fruits of different colors revealed that the development of fruit skin color in jaboticaba is associated with the phenylpropanoid and flavonoid biosynthesis pathways, with the flavanone 3-hydroxylase (F3H) gene potentially regulating fruit skin color. Additionally, by constructing the regulatory pathway of the citric acid cycle, we found that low citric acid content is correlated with high expression levels of genes such as thiamin diphosphate (ThDP) and low expression of phosphoenolpyruvate carboxykinase (PEPCK), indicating that PEPCK positively regulates citric acid content. These T2T genomic resources will accelerate jaboticaba pepper genetic improvement and help to understand jaboticaba genome evolution.

FEAtl: a comprehensive web-based expression atlas for functional genomics in tropical and subtropical fruit crops

BACKGROUND

Fruit crops, including tropical and subtropical fruits like Avocado (Persea americana), Fig (Ficus carica), Date Palm (Phoenix dactylifera), Mango (Mangifera indica), Guava (Psidium guajava), Papaya (Carica papaya), Pineapple (Ananas comosus), and Banana (Musa acuminata) are economically vital, contributing significantly to global agricultural output, as classified by the FAO's World Programme for the Census of Agriculture. Advancements in next-generation sequencing, have transformed fruit crop breeding by providing in-depth genomic and transcriptomic data. RNA sequencing enables high-throughput analysis of gene expression, and functional genomics, crucial for addressing horticultural challenges and enhancing fruit production. The genomic and expression data for key tropical and sub-tropical fruit crops is currently lacking a comprehensive expression atlas, revealing a significant gap in resources for horticulturists who require a unified platform with diverse datasets across various conditions and cultivars.

RESULTS

The Fruit Expression Atlas (FEAtl), available at http://backlin.cabgrid.res.in/FEAtl/ , is a first-ever extensive and unified expression atlas for tropical and subtropical fruit crops developed using 3-tier architecture. The expressivity of coding and non-coding genes, encompassing 2,060 RNA-Seq samples across 91 tissue types and 177 BioProjects, it provides a comprehensive view of gene expression patterns for different tissues under various conditions. FEAtl features multiple tabs that cater to different aspects of the dataset, namely, Home, About, Analyze, Statistics, and Team and contains seven central functional modules: Transcript Information,Sample Information, Expression Profiles in FPKM and TPM, Functional Analysis, Genes Based on Tau Score, and Search for Specific Gene. The expression of a transcript of interest can be easily queried by searching by tissue ID and transcript type. Expression data can be displayed as a heat map, along with functional descriptions as well as Gene Ontology and Kyoto Encyclopedia of Genes and Genomes.

CONCLUSIONS

This atlas represents a groundbreaking compilation of a wide array of information pertaining to eight distinct fruit crops and serves as a fundamental resource for comparative analysis among different fruit species and is a catalyst for functional genomic studies. Database availability: http://backlin.cabgrid.res.in/FEAtl/ .

A first de novo transcriptome assembly of feijoa (Acca sellowiana [Berg] Burret) reveals key genes involved in flavonoid biosynthesis

Acca sellowiana [Berg] Burret, a cultivated fruit tree originating from South America, is gaining the attention of the nutraceutical and pharmaceutical industries due to their high content of flavonoids and other phenolic compounds in fruits, leaves, and flowers. Flavonoids are a diverse group of secondary metabolites with antioxidant, anti-inflammatory, and antimicrobial properties. They also play a crucial role in plant immune response. Despite their importance, the lack of research on A. sellowiana genomics and transcriptomics hinders a deeper understanding of the molecular mechanisms behind flavonoid biosynthesis and its regulation. Here, we de novo assembled and benchmarked 11 A. sellowiana transcriptomes from leaves and floral tissues at three developmental stages using high-throughput sequencing. We selected and annotated the best assembly according to commonly used metrics and databases. This reference transcriptome consisted of 221,649 nonredundant transcripts, of which 107,612 were functionally annotated. We then used this reference transcriptome to explore the expression profiling of key secondary metabolite genes. Transcripts from genes involved in the flavonoid and anthocyanin biosynthesis pathways were identified. We also identified 4068 putative transcription factors, with the most abundant families being bHLH, C2H2, NAC, MYB, and MYB-related. Transcript expression profiling revealed distinct patterns of gene expression during flower development. Particularly, we found 71 differentially expressed transcripts representing 14 enzymes of the flavonoid pathway, suggesting major changes in flavonoid accumulation across floral stages. Our findings will contribute to understanding the genetic basis of flavonoids and provide a foundation for further research and exploitation of the economic potential of this species.

Comparative genomic analysis reveals expansion of the DnaJ gene family in Lagerstroemia indica and its members response to salt stress

DnaJs/Hsp40s/JPDs are obligate co-chaperones of heat shock proteins (Hsp70), performing crucial biological functions within organisms. A comparative genome analysis of four genomes (Vitis vinifera, Eucalyptus grandis, Lagerstroemia indica, and Punica granatum) revealed that the DnaJ gene family in L. indica has undergone expansion, although not to the extent observed in P. granatum. Inter-genome collinearity analysis of four plants indicates that members belonging to Class A and B are more conserved during evolution. In L. indica, the expanded members primarily belong to Class-C. Tissue expression patterns and the biochemical characterization of LiDnaJs further suggested that DnaJs may be involved in numerous biological processes in L. indica. Transcriptome and qPCR analyses of salt stressed leaves identified at least ten LiDnaJs that responded to salt stress. In summary, we have elucidated the expansion mechanism of the LiDnaJs, which is attributed to a recent whole-genome triplication. This research laid the foundation for functional analysis of LiDnaJs and provides gene resources for breeding salt-tolerant varieties of L. indica.

Effects of trimer repeats on Psidium guajava L. gene expression and prospection of functional microsatellite markers

Most research on trinucleotide repeats (TRs) focuses on human diseases, with few on the impact of TR expansions on plant gene expression. This work investigates TRs' effect on global gene expression in Psidium guajava L., a plant species with widespread distribution and significant relevance in the food, pharmacology, and economics sectors. We analyzed TR-containing coding sequences in 1,107 transcripts from 2,256 genes across root, shoot, young leaf, old leaf, and flower bud tissues of the Brazilian guava cultivars Cortibel RM and Paluma. Structural analysis revealed TR sequences with small repeat numbers (5-9) starting with cytosine or guanine or containing these bases. Functional annotation indicated TR-containing genes' involvement in cellular structures and processes (especially cell membranes and signal recognition), stress response, and resistance. Gene expression analysis showed significant variation, with a subset of highly expressed genes in both cultivars. Differential expression highlighted numerous down-regulated genes in Cortibel RM tissues, but not in Paluma, suggesting interplay between tissues and cultivars. Among 72 differentially expressed genes with TRs, 24 form miRNAs, 13 encode transcription factors, and 11 are associated with transposable elements. In addition, a set of 20 SSR-annotated, transcribed, and differentially expressed genes with TRs was selected as phenotypic markers for Psidium guajava and, potentially for closely related species as well.

Postharvest Physiology and Handling of Guava Fruit

Guavas are typical tropical fruit with high nutritional and commercial value. Because of their thin skin and high metabolic rate, guavas are highly susceptible to water loss, physical damage, and spoilage, severely limiting their shelf-life. Guavas can typically only be stored for approximately one week at room temperature, making transportation, storage, and handling difficult, resulting in low profit margins in the industry. This review focuses on the physiological and biochemical changes and their molecular mechanisms which occur in postharvest guavas, and summarizes the various management strategies for extending the shelf-life of these sensitive fruits by means of physical and chemical preservation and their combinations. This review also suggests future directions and reference ideas for the development of safe and efficient shelf-life extension techniques.

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica

BACKGROUND

Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds.

RESULTS

In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay.

CONCLUSIONS

This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

A high-density linkage map construction in guava (Psidium guajava L.) using genotyping by sequencing and identification of QTLs for leaf, peel, and pulp color in an intervarietal mapping population

Psidium guajava L. is an important fruit crop in the tropical and subtropical regions of the world. The advanced breeding methods are not employed for important commercial traits like peel and pulp color, seed hardiness, fruit size, etc., due to the scarcity of genome-wide molecular markers and high-density linkage maps. In this study, we employed single-nucleotide polymorphism (SNP) markers and identified quantitative trait loci (QTL) regions that are associated with color traits of leaf, peel, and pulp in the guava intervarietal mapping population. The mapping population was developed from the contrasting genotypes of fruit and leaf color. Variations in color among the segregating hybrids were recorded both visually and using a Color reader. A high-density linkage map of guava was constructed using the SNP markers from genotyping by sequencing (GBS) of 150 hybrid individuals of the cross 'Arka Poorna' (green) x 'Purple Local' (purple). The integrated linkage map consisted of 1426 SNPs mapped on 11 linkage groups (LG), spanning a total distance of around 730 cM with an average of 129.6 markers per LG. Through QTL analysis for color traits, a minor QTL region was identified for visually scored leaf color and peel color on LG1, whereas a major QTL was detected for pulp color in LG4. The Hunter color values (L* and, a*) also had major QTLs with overlapping marker intervals for leaf and peel colors, establishing the association of SNP markers to the trait. The QTLs harbored genes and transcription factors involved in lycopene and anthocyanin pigment biosynthesis. This is the first report of a high-density linkage map based on SNP markers in guava and QTL mapping for color characters in leaf, fruit peel and pulp. The genotyping information generated in this study can aid in genetic engineering and marker-assisted breeding in guava.

Multiple Physiological and Biochemical Functions of Ascorbic Acid in Plant Growth, Development, and Abiotic Stress Response

Ascorbic acid (AsA) is an important nutrient for human health and disease cures, and it is also a crucial indicator for the quality of fruit and vegetables. As a reductant, AsA plays a pivotal role in maintaining the intracellular redox balance throughout all the stages of plant growth and development, fruit ripening, and abiotic stress responses. In recent years, the de novo synthesis and regulation at the transcriptional level and post-transcriptional level of AsA in plants have been studied relatively thoroughly. However, a comprehensive and systematic summary about AsA-involved biochemical pathways, as well as AsA's physiological functions in plants, is still lacking. In this review, we summarize and discuss the multiple physiological and biochemical functions of AsA in plants, including its involvement as a cofactor, substrate, antioxidant, and pro-oxidant. This review will help to facilitate a better understanding of the multiple functions of AsA in plant cells, as well as provide information on how to utilize AsA more efficiently by using modern molecular biology methods.

Comprehensive genomic exploration of class III peroxidase genes in guava unravels physiology, evolution, and postharvest storage responses

Peroxidases (PRXs) play multifaceted roles in plant growth, development, and stress responses. Here, we present a comprehensive analysis of the PRX gene family in guava, a globally significant fruit. In the guava genome, we identified 37 PRX genes, a number lower than that of Arabidopsis, suggesting a distinctive gene family expansion pattern. Phylogenetic analysis unveiled close relationships with Arabidopsis PRXs, with 12 PgPRX genes forming ortholog pairs, indicating a specific expansion pattern. Predictions placed most PRX proteins in the chloroplast and extracellular regions. Structural analysis of PgPRX proteins revealed commonalities in domain structures and motif organization. Synteny analysis underscored the dynamic role of segmental duplication in the evolution of guava's PRX genes. We explored the dynamic expression of PgPRX genes across guava tissues, exposing functional diversity. Furthermore, we examined changes in peroxidase levels and gene expressions during postharvest fruit storage, providing insights for preserving fruit quality. This study offers an initial genome-wide identification and characterization of Class III peroxidases in guava, laying the foundation for future functional analyses.

Polyploid Genome Assembly Provides Insights into Morphological Development and Ascorbic Acid Accumulation of Sauropus androgynus

Sauropus androgynus (S. androgynus) (2n = 4x = 52) is one of the most popular functional leafy vegetables in South and Southeast Asia. With its rich nutritional and pharmaceutical values, it has traditionally had widespread use for dietary and herbal purposes. Here, the genome of S. androgynus was sequenced and assembled, revealing a genome size of 1.55 Gb with 26 pseudo-chromosomes. Phylogenetic analysis traced back the divergence of Sauropus from Phyllanthus to approximately 29.67 million years ago (Mya). Genome analysis revealed that S. androgynus polyploidized around 20.51 Mya and shared a γ event about 132.95 Mya. Gene function analysis suggested that the expansion of pathways related to phloem development, lignin biosynthesis, and photosynthesis tended to result in the morphological differences among species within the Phyllanthaceae family, characterized by varying ploidy levels. The high accumulation of ascorbic acid in S. androgynus was attributed to the high expression of genes associated with the L-galactose pathway and recycling pathway. Moreover, the expanded gene families of S. androgynus exhibited multiple biochemical pathways associated with its comprehensive pharmacological activity, geographic adaptation and distinctive pleasurable flavor. Altogether, our findings represent a crucial genomic asset for S. androgynus, casting light on the intricate ploidy within the Phyllanthaceae family.

The haplotype-resolved autotetraploid genome assembly provides insights into the genomic evolution and fruit divergence in wax apple (Syzygium samarangense (Blume) Merr. and Perry)

Wax apple (Syzygium samarangense) is an economically important fruit crop with great potential value to human health because of its richness in antioxidant substances. Here, we present a haplotype-resolved autotetraploid genome assembly of the wax apple with a size of 1.59 Gb. Comparative genomic analysis revealed three rounds of whole-genome duplication (WGD) events, including two independent WGDs after WGT-γ. Resequencing analysis of 35 accessions partitioned these individuals into two distinct groups, including 28 landraces and seven cultivated species, and several genes subject to selective sweeps possibly contributed to fruit growth, including the KRP1-like, IAA17-like, GME-like, and FLACCA-like genes. Transcriptome analysis of three different varieties during flower and fruit development identified key genes related to fruit size, sugar content, and male sterility. We found that AP2 also affected fruit size by regulating sepal development in wax apples. The expression of sugar transport-related genes (SWEETs and SUTs) was high in 'ZY', likely contributing to its high sugar content. Male sterility in 'Tub' was associated with tapetal abnormalities due to the decreased expression of DYT1, TDF1, and AMS, which affected early tapetum development. The chromosome-scale genome and large-scale transcriptome data presented in this study offer new valuable resources for biological research on S. samarangense and shed new light on fruit size control, sugar metabolism, and male sterility regulatory metabolism in wax apple.

Chromosome-scale assemblies of S. malaccense, S. aqueum, S. jambos, and S. syzygioides provide insights into the evolution of Syzygium genomes

Syzygium is a large and diverse tree genus in the Myrtaceae family. Genome assemblies for clove (Syzygium aromaticum, 370 Mb) and sea apple (Syzygium grande, 405 Mb) provided the first insights into the genomic features and evolution of the Syzygium genus. Here, we present additional de novo chromosome-scale genome assemblies for Syzygium malaccense, Syzygium aqueum, Syzygium jambos, and Syzygium syzygioides. Genome profiling analyses show that S. malaccense, like S. aromaticum and S. grande, is diploid (2n = 2x = 22), while the S. aqueum, S. jambos, and S. syzygioides specimens are autotetraploid (2n = 4x = 44). The genome assemblies of S. malaccense (430 Mb), S. aqueum (392 Mb), S. jambos (426 Mb), and S. syzygioides (431 Mb) are highly complete (BUSCO scores of 98%). Comparative genomics analyses showed conserved organization of the 11 chromosomes with S. aromaticum and S. grande, and revealed species-specific evolutionary dynamics of the long terminal repeat retrotransposon elements belonging to the Gypsy and Copia lineages. This set of Syzygium genomes is a valuable resource for future structural and functional comparative genomic studies on Myrtaceae species.

Genome of Phyllanthus emblica: the medicinal plant Amla with super antioxidant properties

Phyllanthus emblica or Indian gooseberry, commonly known as amla, is an important medicinal horticultural plant used in traditional and modern medicines. It bears stone fruits with immense antioxidant properties due to being one of the richest natural sources of vitamin C and numerous flavonoids. This study presents the first genome sequencing of this species performed using 10x Genomics and Oxford Nanopore Technology. The draft genome assembly was 519 Mbp in size and consisted of 4,384 contigs, N50 of 597 Kbp, 98.4% BUSCO score, and 37,858 coding sequences. This study also reports the genome-wide phylogeny of this species with 26 other plant species that resolved the phylogenetic position of P. emblica. The presence of three ascorbate biosynthesis pathways including L-galactose, galacturonate, and myo-inositol pathways was confirmed in this genome. A comprehensive comparative evolutionary genomic analysis including gene family expansion/contraction and identification of multiple signatures of adaptive evolution provided evolutionary insights into ascorbate and flavonoid biosynthesis pathways and stone fruit formation through lignin biosynthesis. The availability of this genome will be beneficial for its horticultural, medicinal, dietary, and cosmetic applications and will also help in comparative genomics analysis studies.

E. urophylla × E. grandis high-quality genome and comparative genomics provide insights on evolution and diversification of eucalyptus

BACKGROUND

Eucalyptus urophylla × Eucalyptus grandis, an economically important forest tree, provides important raw material for energy and reduces damage to native forests. However, the absence of a high-quality E. urophylla × E. grandis reference genome has significantly hindered its evolution and genetic analysis.

RESULTS

We successfully presented a high-quality reference genome of E. urophylla × E. grandis (545.75 Mb; scaffold N50, 51.62 Mb) using a combination of the Illumina, PacBio HiFi, and Hi-C sequencing platforms. A total of 34,502 genes and 58.56% of the repetitive sequences in this genome were annotated. Using genome evolution analyses, we identified a recent whole-genome duplication (WGD) event in E. urophylla × E. grandis. We further found that gene families associated with starch and sucrose metabolism, flavonoid biosynthesis, and plant-pathogen interaction were significantly expanded in E. urophylla × E. grandis. Moreover, comparative genomic and evolutionary analyses showed large structural variations among the different chromosomes of the 34 Eucalyptus accessions, which were divided into six clades.

CONCLUSIONS

Overall, our findings provide a valuable resource for expanding our understanding of the E. urophylla × E. grandis genome evolution, genetic improvement, and its comparative biology.

Complete Genome of Rose Myrtle, Rhodomyrtus tomentosa, and Its Population Genetics in Thai Peninsula

Several parts of rose myrtle, Rhodomyrtus tomentosa, exhibited profound antibacterial and anti-inflammatory activities, suggesting its potential in healthcare and cosmetics applications. During the past few years, the demand for biologically active compounds in the industrial sectors increased. Therefore, gathering comprehensive information on all aspects of this plant species is essential. Here, the genome sequencing using short and long reads was used to understand the genome biology of R. tomentosa. Inter-simple sequence repeats (ISSR) and simple sequence repeats (SSR) markers, and geometric morphometrics of the leaves of R. tomentosa collected across Thai Peninsula, were determined for population differentiation analysis. The genome size of R. tomentosa was 442 Mb, and the divergence time between R. tomentosa and Rhodamnia argentea, the white myrtle of eastern Australia, was around 15 million years. No population structure was observed between R. tomentosa on the eastern and western sides of the Thai Peninsula using the ISSR and SSR markers. However, significant differences in leaf size and shape of R. tomentosa were observed in all locations.

Gap-free genome assembly and comparative analysis reveal the evolution and anthocyanin accumulation mechanism of Rhodomyrtus tomentosa

Rhodomyrtus tomentosa is an important fleshy-fruited tree and a well-known medicinal plant of the Myrtaceae family that is widely cultivated in tropical and subtropical areas of the world. However, studies on the evolution and genomic breeding of R. tomentosa were hindered by the lack of a reference genome. Here, we presented a chromosome-level gap-free T2T genome assembly of R. tomentosa using PacBio and ONT long read sequencing. We assembled the genome with size of 470.35 Mb and contig N50 of ~43.80 Mb with 11 pseudochromosomes. A total of 33 382 genes and 239.31 Mb of repetitive sequences were annotated in this genome. Phylogenetic analysis elucidated the independent evolution of R. tomentosa starting from 14.37MYA and shared a recent WGD event with other Myrtaceae species. We identified four major compounds of anthocyanins and their synthetic pathways in R. tomentosa. Comparative genomic and gene expression analysis suggested the coloring and high anthocyanin accumulation in R. tomentosa tends to be determined by the activation of anthocyanin synthesis pathway. The positive selection and up-regulation of MYB transcription factors were the implicit factors in this process. The copy number increase of downstream anthocyanin transport-related OMT and GST gene were also detected in R. tomentosa. Expression analysis and pathway identification enriched the importance of starch degradation, response to stimuli, effect of hormones, and cell wall metabolism during the fleshy fruit development in Myrtaceae. Our genome assembly provided a foundation for investigating the origins and differentiation of Myrtaceae species and accelerated the genetic improvement of R. tomentosa.

Genomic and epigenomic variation in Psidium species and their outcome under the yield and composition of essential oils

Diploid and polyploid species derived from the euploid series x = 11 occur in the genus Psidium, as well as intraspecific cytotypes. Euploidy in the genus can alter the gene copy number, resulting in several "omics" variations. We revisited the euploidy, reported genomic (nuclear 2C value, GC%, and copy number of secondary metabolism genes) and epigenomic (5-mC%) differences in Psidium, and related them to essential oil yield and composition. Mean 2C values ranged from 0.90 pg (P. guajava) to 7.40 pg (P. gaudichaudianum). 2C value is intraspecifically varied in P. cattleyanum and P. gaudichaudianum, evidencing cytotypes that can be formed from euploid (non-reduced) and/or aneuploid reproductive cells. GC% ranged from 34.33% (P. guineense) to 48.95% (P. myrtoides), and intraspecific variations occurred even for species without 2C value intraspecific variation. Essential oil yield increased in relation to 2C value and to GC%. We showed that P. guajava (diploid) possesses two and P. guineense (tetraploid) four copies of the one specific TPS gene, as well as eight and sixteen copies respectively of the conserved regions that occur in eight TPS genes. We provide a wide "omics'' characterization of Psidium and show the outcome of the genome and epigenome variation in secondary metabolism.

Gene expression profiling in Rosa roxburghii fruit and overexpressing RrGGP2 in tobacco and tomato indicates the key control point of AsA biosynthesis

Rosa roxburghii Tratt. is an important commercial horticultural crop endemic to China, which is recognized for its extremely high content of L-ascorbic acid (AsA). To understand the mechanisms underlying AsA overproduction in fruit of R. roxburghii, content levels, accumulation rate, and the expression of genes putatively in the biosynthesis of AsA during fruit development have been characterized. The content of AsA increased with fruit weight during development, and AsA accumulation rate was found to be highest between 60 and 90 days after anthesis (DAA), with approximately 60% of the total amount being accumulated during this period. In vitro incubating analysis of 70DAA fruit flesh tissues confirmed that AsA was synthesized mainly via the L-galactose pathway although L-Gulono-1, 4-lactone was also an effective precursor elevating AsA biosynthesis. Furthermore, in transcript level, AsA content was significantly associated with GDP-L-galactose phosphorylase (RrGGP2) gene expression. Virus-induced RrGGP2 silencing reduced the AsA content in R. roxburghii fruit by 28.9%. Overexpressing RrGGP2 increased AsA content by 8-12-fold in tobacco leaves and 2.33-3.11-fold in tomato fruit, respectively, and it showed enhanced resistance to oxidative stress caused by paraquat in transformed tobacco. These results further justified the importance of RrGGP2 as a major control step to AsA biosynthesis in R. roxburghii fruit.

Effects of Pectinase Pre-Treatment on the Physicochemical Properties, Bioactive Compounds, and Volatile Components of Juices from Different Cultivars of Guava

In this study, the physicochemical properties, antioxidant capacity, and volatile compounds of the juices of different guava cultivars before and after pectinase treatment were evaluated. The results showed that the guava juice of the small fragrant (SF) cultivar exhibited the highest ascorbic acid concentration (1761.09 mg/L), and the highest contents of total phenolics (329.52 mg GAE/L) and total flavonoids (411.13 mg RE/L) were both found in the juice of the watermelon red (WR) cultivar. After pectinase treatment, the juice yield and the titratable acid, sugar components, total phenolics and total flavonoids, and antioxidant capacity levels of the guava juices were all higher than those of the non-pectinase group. However, lower sensory evaluation scores were obtained in the pectinase-treated guava juices. Aldehydes and terpenoids were the main flavor components in the guava juices, which were responsible for the aroma of the juice, while their relative contents were different in the four cultivar guava juices. Furthermore, pectinase treatment could change the amounts and relative contents of volatile compounds in the guava juice. During the pectinase treatment process, the relative contents of the main aroma constituents in the guava juices were significantly decreased. The findings of this research provide valuable information for the processing of guava juice.

A high-quality pseudo-phased genome for Melaleuca quinquenervia shows allelic diversity of NLR-type resistance genes

BACKGROUND

Melaleuca quinquenervia (broad-leaved paperbark) is a coastal wetland tree species that serves as a foundation species in eastern Australia, Indonesia, Papua New Guinea, and New Caledonia. While extensively cultivated for its ornamental value, it has also become invasive in regions like Florida, USA. Long-lived trees face diverse pest and pathogen pressures, and plant stress responses rely on immune receptors encoded by the nucleotide-binding leucine-rich repeat (NLR) gene family. However, the comprehensive annotation of NLR encoding genes has been challenging due to their clustering arrangement on chromosomes and highly repetitive domain structure; expansion of the NLR gene family is driven largely by tandem duplication. Additionally, the allelic diversity of the NLR gene family remains largely unexplored in outcrossing tree species, as many genomes are presented in their haploid, collapsed state.

RESULTS

We assembled a chromosome-level pseudo-phased genome for M. quinquenervia and described the allelic diversity of plant NLRs using the novel FindPlantNLRs pipeline. Analysis reveals variation in the number of NLR genes on each haplotype, distinct clustering patterns, and differences in the types and numbers of novel integrated domains.

CONCLUSIONS

The high-quality M. quinquenervia genome assembly establishes a new framework for functional and evolutionary studies of this significant tree species. Our findings suggest that maintaining allelic diversity within the NLR gene family is crucial for enabling responses to environmental stress, particularly in long-lived plants.

Chemical Constituents, Antioxidant, and Enzyme Inhibitory Activities Supported by In-Silico Study of n-Hexane Extract and Essential Oil of Guava Leaves

Psidium guajava (Guava tree) is one of the most widely known species in the family Myrtaceae. The Guava tree has been reported for its potential antioxidant, anti-inflammatory, antimicrobial, and cytotoxic activities. In the current study, the chemical compositions of the n-hexane extract and the essential oil of P. guajava were investigated using the GC/MS analysis, along with an evaluation of their antioxidant potential, and an investigation into the enzyme inhibition of acetylcholinesterase (AChE), butyrylcholinesterase (BchE), tyrosinase, α-amylase, and α-glucosidase. Moreover, molecular docking of the major identified active sites of the target enzymes were investigated. The chemical characterization of the n-hexane extract and essential oil revealed that squalene (9.76%), α-tocopherol (8.53%), and γ-sitosterol (3.90%) are the major compounds in the n-hexane extract. In contrast, the major constituents of the essential oil are D-limonene (36.68%) and viridiflorol (9.68%). The n-hexane extract showed more antioxidant potential in the cupric reducing antioxidant capacity (CUPRAC), the ferric reducing power (FRAP), and the metal chelating ability (MCA) assays, equivalent to 70.80 ± 1.46 mg TE/g, 26.01 ± 0.97 mg TE/g, and 24.83 ± 0.35 mg EDTAE/g, respectively. In the phosphomolybdenum (PM) assay, the essential oil showed more antioxidant activity equivalent to 2.58 ± 0.14 mmol TE/g. The essential oil demonstrated a potent BChE and tyrosinase inhibitory ability at 6.85 ± 0.03 mg GALAE/g and 61.70 ± 3.21 mg KAE/g, respectively. The α-amylase, and α-glucosidase inhibitory activity of the n-hexane extract and the essential oil varied from 0.52 to 1.49 mmol ACAE/g. Additionally, the molecular docking study revealed that the major compounds achieved acceptable binding scores upon docking with the tested enzymes. Consequently, the P. guajava n-hexane extract and oil can be used as a promising candidate for the development of novel treatment strategies for oxidative stress, neurodegeneration, and diabetes mellitus diseases.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Genetic associations with resistance to Meloidogyne enterolobii in guava (Psidium sp.) using cross-genera SNPs and comparative genomics to Eucalyptus highlight evolutionary conservation across the Myrtaceae

Tropical fruit tree species constitute a yet untapped supply of outstanding diversity of taste and nutritional value, barely developed from the genetics standpoint, with scarce or no genomic resources to tackle the challenges arising in modern breeding practice. We generated a de novo genome assembly of the Psidium guajava, the super fruit “apple of the tropics”, and successfully transferred 14,268 SNP probesets from Eucalyptus to Psidium at the nucleotide level, to detect genomic loci linked to resistance to the root knot nematode (RKN) Meloidogyne enterolobii derived from the wild relative P. guineense. Significantly associated loci with resistance across alternative analytical frameworks, were detected at two SNPs on chromosome 3 in a pseudo-assembly of Psidium guajava genome built using a syntenic path approach with the Eucalyptus grandis genome to determine the order and orientation of the contigs. The P. guineense-derived resistance response to RKN and disease onset is conceivably triggered by mineral nutrients and phytohormone homeostasis or signaling with the involvement of the miRNA pathway. Hotspots of mapped resistance quantitative trait loci and functional annotation in the same genomic region of Eucalyptus provide further indirect support to our results, highlighting the evolutionary conservation of genomes across genera of Myrtaceae in the adaptation to pathogens. Marker assisted introgression of the resistance loci mapped should accelerate the development of improved guava cultivars and hybrid rootstocks.

A high-quality Buxus austro-yunnanensis (Buxales) genome provides new insights into karyotype evolution in early eudicots

BACKGROUND

Eudicots are the most diverse group of flowering plants that compromise five well-defined lineages: core eudicots, Ranunculales, Proteales, Trochodendrales, and Buxales. However, the phylogenetic relationships between these five lineages and their chromosomal evolutions remain unclear, and a lack of high-quality genome analyses for Buxales has hindered many efforts to address this knowledge gap.

RESULTS

Here, we present a high-quality chromosome-level genome of Buxus austro-yunnanensis (Buxales). Our phylogenomic analyses revealed that Buxales and Trochodendrales are genetically similar and classified as sisters. Additionally, both are sisters to the core eudicots, while Ranunculales was found to be the first lineage to diverge from these groups. Incomplete lineage sorting and hybridization were identified as the main contributors to phylogenetic discordance (34.33%) between the lineages. In fact, B. austro-yunnanensis underwent only one whole-genome duplication event, and collinear gene phylogeny analyses suggested that separate independent polyploidizations occurred in the five eudicot lineages. Using representative genomes from these five lineages, we reconstructed the ancestral eudicot karyotype (AEK) and generated a nearly gapless karyotype projection for each eudicot species. Within core eudicots, we recovered one common chromosome fusion event in asterids and malvids, respectively. Further, we also found that the previously reported fused AEKs in Aquilegia (Ranunculales) and Vitis (core eudicots) have different fusion positions, which indicates that these two species have different karyotype evolution histories.

CONCLUSIONS

Based on our phylogenomic and karyotype evolution analyses, we revealed the likely relationships and evolutionary histories of early eudicots. Ultimately, our study expands genomic resources for early-diverging eudicots.

Adaptive radiation in Orinus, an endemic alpine grass of the Qinghai-Tibet Plateau, based on comparative transcriptomic analysis

The species of Orinus (Poaceae) are important alpine plants with a variety of phenotypic traits and potential usages in molecular breeding toward drought-tolerant forage crops. However, the genetic basis of evolutionary adaption and diversification in the genus is still unclear. In the present study, we obtained transcriptomes for the two most divergent species, O. thoroldii and O. kokonoricus, using the Illumina platform and de novo assembly. In total, we generated 23,029 and 24,086 unigenes with N50 values of 1188 and 1203 for O. thoroldii and O. kokonoricus respectively, and identified 19,005 pairs of putative orthologs between the two species of Orinus. For these orthologs, estimations of non-synonymous/synonymous substitution rate ratios indicated that 568 pairs may be under strongly positive selection (Ka/Ks > 1), and Gene Ontogeny (GO) enrichment analysis revealed that significantly enriched pathways were in DNA repair and resistance to abiotic stress. Meanwhile, the divergence times of species between O. thoroldii and O. kokonoricus occurred 3.2 million years ago (Mya), and the recent evolutionary branch is an allotetraploid species, Cleistogenes songorica. We also detected a Ks peak of ∼0.60 for Orinus. Additionally, we identified 188 pairs of differentially expressed genes (DEGs) between the two species of Orinus, which were significantly enrich in stress resistance and lateral root development. Thus, we considered that the species diversification and evolutionary adaption of this genus was initiated by environmental selection, followed by phenotypic differentiation, finally leading to niche separation in the Qinghai-Tibet Plateau.

Genome-wide identification of calcineurin B-like protein-interacting protein kinase gene family reveals members participating in abiotic stress in the ornamental woody plant Lagerstroemia indica

Calcineurin B-like protein-interacting protein kinases (CIPKs) play important roles in plant responses to stress. However, their function in the ornamental woody plant Lagerstroemia indica is remains unclear. In this study, the LiCIPK gene family was analyzed at the whole genome level. A total of 37 LiCIPKs, distributed across 17 chromosomes, were identified. Conserved motif analysis indicated that all LiCIPKs possess a protein kinase motif (S_TKc) and C-terminal regulatory motif (NAF), while seven LiCIPKs lack a protein phosphatase interaction (PPI) motif. 3D structure analysis further revealed that the N-terminal and C-terminal 3D-structure of 27 members are situated near to each other, while 4 members have a looser structure, and 6 members lack intact structures. The intra- and interspecies collinearity analysis, synonymous substitution rate (K _s ) peaks of duplicated LiCIPKs, revealed that ∼80% of LiCIPKs were retained by the two whole genome duplication (WGD) events that occurred approximately 56.12-61.16 million year ago (MYA) and 16.24-26.34 MYA ago. The promoter of each LiCIPK contains a number of auxin, abscisic acid, gibberellic acid, salicylic acid, and drought, anaerobic, defense, stress, and wound responsive cis-elements. Of the 21 members that were successfully amplified by qPCR, 18 LiCIPKs exhibited different expression patterns under NaCl, mannitol, PEG8000, and ABA treatments. Given that LiCIPK30, the AtSOS2 ortholog, responded to all four types of stress it was selected for functional verification. LiCIPK30 complements the atsos2 phenotype in vivo. 35S:LiCIPK-overexpressing lines exhibit increased leaf area increment, chlorophyll a and b content, reactive oxygen species scavenging enzyme activity, and expression of ABF3 and RD22, while the degree of membrane lipid oxidation decreases under NaCl treatment compared to WT. The evolutionary history, and potential mechanism by which LiCIPK30 may regulate plant tolerance to salt stress were also discussed. In summary, we identified LiCIPK members involved in abiotic stress and found that LiCIPK30 transgenic Arabidopsis exhibits more salt and osmotic stress tolerance than WT. This research provides a theoretical foundation for further investigation into the function of LiCIPKs, and for mining gene resources to facilitate the cultivation and breeding of new L. indica varieties in coastal saline-alkali soil.

Genomic insights into rapid speciation within the world's largest tree genus Syzygium

Species radiations, despite immense phenotypic variation, can be difficult to resolve phylogenetically when genetic change poorly matches the rapidity of diversification. Genomic potential furnished by palaeopolyploidy, and relative roles for adaptation, random drift and hybridisation in the apportionment of genetic variation, remain poorly understood factors. Here, we study these aspects in a model radiation, Syzygium, the most species-rich tree genus worldwide. Genomes of 182 distinct species and 58 unidentified taxa are compared against a chromosome-level reference genome of the sea apple, Syzygium grande. We show that while Syzygium shares an ancient genome doubling event with other Myrtales, little evidence exists for recent polyploidy events. Phylogenomics confirms that Syzygium originated in Australia-New Guinea and diversified in multiple migrations, eastward to the Pacific and westward to India and Africa, in bursts of speciation visible as poorly resolved branches on phylogenies. Furthermore, some sublineages demonstrate genomic clines that recapitulate cladogenetic events, suggesting that stepwise geographic speciation, a neutral process, has been important in Syzygium diversification.

Jackfruit genome and population genomics provide insights into fruit evolution and domestication history in China

As the largest known tree-borne fruit in the world, jackfruit (Artocarpus heterophyllus) is an important cultivated crop in tropical regions of South and Southeast Asia. The species has been cultivated in China for more than 1000 years, but the history of its introduction to the country remains unclear. We assembled a high-quality chromosome-level genome of jackfruit into 985.63 Mb with scaffold N50 of 32.81 Mb. We analyzed whole-genome resequencing data of 295 landraces to investigate the domestication history in China and agronomic trait evolution of jackfruit. Population structure analysis revealed that jackfruits of China could be traced back to originate from Southeast Asia and South Asia independently. Selection signals between jackfruit and its edible congener, cempedak (Artocarpus integer), revealed several important candidate genes associated with fruit development and ripening. Moreover, analyses of selective sweeps and gene expression revealed that the AhePG1 gene may be the major factor in determining fruit texture. This study not only resolves the origins of jackfruit of China, but also provides valuable genomic resources for jackfruit breeding improvement and offers insights into fruit size evolution and fruit texture changes.

The clove (Syzygium aromaticum) genome provides insights into the eugenol biosynthesis pathway

The clove (Syzygium aromaticum) is an important tropical spice crop in global trade. Evolving environmental pressures necessitate modern characterization and selection techniques that are currently inaccessible to clove growers owing to the scarcity of genomic and genetic information. Here, we present a 370-Mb high-quality chromosome-scale genome assembly for clove. Comparative genomic analysis between S. aromaticum and Eucalyptus grandis-both species of the Myrtaceae family-reveals good genome structure conservation and intrachromosomal rearrangements on seven of the eleven chromosomes. We report genes that belong to families involved in the biosynthesis of eugenol, the major bioactive component of clove products. On the basis of our transcriptomic and metabolomic findings, we propose a hypothetical scenario in which eugenol acetate plays a key role in high eugenol accumulation in clove leaves and buds. The clove genome is a new contribution to omics resources for the Myrtaceae family and an important tool for clove research.

Genome-Wide Analysis of Ascorbic Acid Metabolism Related Genes in Fragaria × ananassa and Its Expression Pattern Analysis in Strawberry Fruits

Ascorbic acid (AsA) is an important antioxidant for scavenging reactive oxygen species and it is essential for human health. Strawberry (Fragaria × ananassa) fruits are rich in AsA. In recent years, strawberry has been regarded as a model for non-climacteric fruit ripening. However, in contrast to climacteric fruits, such as tomato, the regulatory mechanism of AsA accumulation in strawberry fruits remains largely unknown. In this study, we first identified 125 AsA metabolism-related genes from the cultivated strawberry "Camarosa" genome. The expression pattern analysis using an available RNA-seq data showed that the AsA biosynthetic-related genes in the D-mannose/L-galactose pathway were downregulated remarkably during fruit ripening which was opposite to the increasing AsA content in fruits. The D-galacturonate reductase gene (GalUR) in the D-Galacturonic acid pathway was extremely upregulated in strawberry receptacles during fruit ripening. The FaGalUR gene above belongs to the aldo-keto reductases (AKR) superfamily and has been proposed to participate in AsA biosynthesis in strawberry fruits. To explore whether there are other genes in the AKR superfamily involved in regulating AsA accumulation during strawberry fruit ripening, we further implemented a genome-wide analysis of the AKR superfamily using the octoploid strawberry genome. A total of 80 FaAKR genes were identified from the genome and divided into 20 subgroups based on phylogenetic analysis. These FaAKR genes were unevenly distributed on 23 chromosomes. Among them, nine genes showed increased expression in receptacles as the fruit ripened, and notably, FaAKR23 was the most dramatically upregulated FaAKR gene in receptacles. Compared with fruits at green stage, its expression level increased by 142-fold at red stage. The qRT-PCR results supported that the expression of FaAKR23 was increased significantly during fruit ripening. In particular, the FaAKR23 was the only FaAKR gene that was significantly upregulated by abscisic acid (ABA) and suppressed by nordihydroguaiaretic acid (NDGA, an ABA biosynthesis blocker), indicating FaAKR23 might play important roles in ABA-mediated strawberry fruit ripening. In a word, our study provides useful information on the AsA metabolism during strawberry fruit ripening and will help understand the mechanism of AsA accumulation in strawberry fruits.

Genome sequence and population genomics provide insights into chromosomal evolution and phytochemical innovation of Hippophae rhamnoides

Plants of the Elaeagnaceae family are widely used to treat various health disorders owing to their natural phytochemicals. Seabuckthorn (Hippophae rhamnoides L.) is an economically and ecologically important species within the family with richness of biologically and pharmacologically active substances. Here, we present a chromosome-level genome assembly of seabuckthorn (http://hipp.shengxin.ren/), the first genome sequence of Elaeagnaceae, which has a total length of 849.04 Mb with scaffold N50 of 69.52 Mb and 30 864 annotated genes. Two sequential tetraploidizations with one occurring ~36-41 million years ago (Mya) and the last ~24-27 Mya were inferred, resulting in expansion of genes related to ascorbate and aldarate metabolism, lipid biosynthesis, and fatty acid elongation. Comparative genomic analysis reconstructed the evolutionary trajectories of the seabuckthorn genome with the predicted ancestral genome of 14 proto-chromosomes. Comparative transcriptomic and metabonomic analyses identified some key genes contributing to high content of polyunsaturated fatty acids and ascorbic acid (AsA). Additionally, we generated and analysed 55 whole-genome sequences of diverse accessions, and identified 9.80 million genetic variants in the seabuckthorn germplasms. Intriguingly, genes in selective sweep regions identified through population genomic analysis appeared to contribute to the richness of AsA and fatty acid in seabuckthorn fruits, among which GalLDH, GMPase and ACC, TER were the potentially major-effect causative genes controlling AsA and fatty acid content of the fruit, respectively. Our research offers novel insights into the molecular basis underlying phytochemical innovation of seabuckthorn, and provides valuable resources for exploring the evolution of the Elaeagnaceae family and molecular breeding.

Metabolic Profiling of Sugars and Organic Acids, and Expression Analyses of Metabolism-Associated Genes in Two Yellow-Peel Pitaya Species

Sugar and organic acids are important factors determining pitaya fruit quality. However, changes in sugars and acids, and expressions of metabolism-associated genes during fruit maturation of yellow-peel pitayas are not well-documented. In this study, metabolic and expression analyses in pulps of different fruit developmental stages of ‘Wucihuanglong’ (‘WCHL’, Hylocereus undatus) and ‘Youcihuanglong’ pitaya (‘YCHL’, Hylocereus megalanthus) were used to explore the sugar and organic acid metabolic process. Total phenols and flavonoids were mainly accumulated at S1 in pitaya pulps. Ascorbic acid contents of ‘WCHL’ pitaya were higher than that of ‘YCHL’ pitaya during fruit maturation. Starch was mainly accumulated at early fruit development stages while soluble sugars were rich in late stages. Sucrose, fructose, and glucose were the main sugar components of ‘YCHL’ pitaya while glucose was dominant in ‘WCHL’ pitaya. Malic and citric acids were the main organic acids in ‘WCHL’ and ‘YCHL’ pitayas, respectively. Based on the transcriptome analyses, 118 genes involved in pitaya sugar and organic acid metabolism were obtained. Results from the correlation analyses between the expression profiling of candidate genes and the contents of sugar and organic acid showed that 51 genes had a significant correlation relationship and probably perform key role in pitaya sugar and organic acid metabolism processes. The finding of the present study provides new information for quality regulation of pitayas.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

Identification in silico and expression analysis of a β-1-4-endoglucanase and β-galactosidase genes related to ripening in guava fruit

BACKGROUND

Guava fruit softening is a crucial process during ripening and this process involves a number of enzymes that modifies the cell wall. Two of the enzymes that regulate this process are (a) the β-1, 4-endoglucanase 17 (BEG) which hydrolyze β-1, 4 bonds from cellulose and hemicellulose, and (b) β-galactosidase (BGA) that hydrolyzes pectin chains. Bioinformatics and expression analysis information on these genes is limited in guava fruit.

RESULTS

A fragment of a β-1, 4-endoglucanase 17 (PgE17), and another of a β-galactosidase (PgGa1) were identified. These sequences have a similarity of more than 85% with those reported in the NCBI database. In the guava genome, one homologous sequence was found for PgE17 in Chr 4 and two homologous to PgGa1: one in Chr 3 and the other one in Chr 6. Putative protein PgE17 contains part of the glyco_hydro_9 domain. Putative protein PgGa1 has a part of the glyco_hydro_35 domain. Phylogenetic analysis of PgE17 and PgGa1 revealed that both are highly conserved inside the Myrtaceae family. In silico expression analysis showed that both PgE17 and PgGa1 work in a coordinated way with other cell wall modifier enzymes. Expression of these genes was found in all the guava samples analyzed. However, the highest expression was found in the fruit in the breaking and ripe states.

CONCLUSIONS

A β-1, 4-endoglucanase 17, and β-galactosidase 1 sequences were identified. PgE17 and PgGa1 are expressed in all the plant tissues, and fruit ripening states. Although, the highest expression was on breaker and ripe states.

Genomic Approaches for Improvement of Tropical Fruits: Fruit Quality, Shelf Life and Nutrient Content

The breeding of tropical fruit trees for improving fruit traits is complicated, due to the long juvenile phase, generation cycle, parthenocarpy, polyploidy, polyembryony, heterozygosity and biotic and abiotic factors, as well as a lack of good genomic resources. Many molecular techniques have recently evolved to assist and hasten conventional breeding efforts. Molecular markers linked to fruit development and fruit quality traits such as fruit shape, size, texture, aroma, peel and pulp colour were identified in tropical fruit crops, facilitating Marker-assisted breeding (MAB). An increase in the availability of genome sequences of tropical fruits further aided in the discovery of SNP variants/Indels, QTLs and genes that can ascertain the genetic determinants of fruit characters. Through multi-omics approaches such as genomics, transcriptomics, metabolomics and proteomics, the identification and quantification of transcripts, including non-coding RNAs, involved in sugar metabolism, fruit development and ripening, shelf life, and the biotic and abiotic stress that impacts fruit quality were made possible. Utilizing genomic assisted breeding methods such as genome wide association (GWAS), genomic selection (GS) and genetic modifications using CRISPR/Cas9 and transgenics has paved the way to studying gene function and developing cultivars with desirable fruit traits by overcoming long breeding cycles. Such comprehensive multi-omics approaches related to fruit characters in tropical fruits and their applications in breeding strategies and crop improvement are reviewed, discussed and presented here.

The Melastoma dodecandrum genome and the evolution of Myrtales

Melastomataceae have abundant morphological diversity with high economic and ornamental merit in Myrtales. The phylogenetic position of Myrtales is still contested. Here, we report the first chromosome-level genome assembly of Melastoma dodecandrum in Melastomataceae. The assembled genome size was 299.81 Mb with a contig N50 value of 3.00 Mb. Genome evolution analysis indicated that M. dodecandrum, Eucalyptus grandis and Punica granatum were clustered into a clade of Myrtales and formed a sister group with the ancestor of fabids and malvids. We found that M. dodecandrum experienced four whole-genome polyploidization events: the ancient event was shared with most eudicots, one event was shared with Myrtales, and the other two events were unique to M. dodecandrum. Moreover, we identified MADS-box genes and found that the AP1-like genes expanded, and AP3-like genes might have undergone subfunctionalization. We found that the SUAR63-like genes and AG-like genes showed different expression patterns in stamens, which may be associated with heteranthery. In addition, we found that LAZY1-like genes were involved in the negative regulation of stem branching development, which may be related to its creeping features. Our study sheds new light on the evolution of Melastomataceae and Myrtales, which provides a comprehensive genetic resource for future research.

Development of Genome-Wide Functional Markers Using Draft Genome Assembly of Guava (Psidium guajava L.) cv. Allahabad Safeda to Expedite Molecular Breeding

Guava (Psidium guajava L.), a rich source of nutrients, is an important tropical and subtropical fruit of the Myrtaceae family and exhibits magnificent diversity. Genetic diversity analysis is the first step toward the identification of parents for hybridization, genetic mapping, and molecular breeding in any crop species. A diversity analysis based on whole-genome functional markers increases the chances of identifying genetic associations with agronomically important traits. Therefore, here, we sequenced the genome of guava cv. Allahabad Safeda on an Illumina platform and generated a draft assembly of ~304 MB. The assembly of the Allahabad Safeda genome constituted >37.95% repeat sequences, gene prediction with RNA-seq data as evidence identified 14,115 genes, and BLAST n/r, Interproscan, PfamScan, BLAST2GO, and KEGG annotated 13,957 genes. A comparative protein transcript analysis of tree species revealed the close relatedness of guava with Eucalyptus. Comparative transcriptomics-based SSR/InDel/SNP-PCR ready genome-wide markers in greenish-yellow skinned and white fleshed-Allahabad Safeda to four contrasting cultivars viz apple-color-skinned and white-fleshed-Lalima, greenish-yellow-skinned and pink-fleshed-Punjab Pink, purple-black-skinned and purple-fleshed-Purple Local and widely used rootstock-Lucknow-49 were developed. The molecular markers developed here revealed a high level of individual heterozygosity within genotypes in 22 phenotypically diverse guava cultivars. Principal coordinate, STRUCTURE clustering, and neighbor-joining-based genetic diversity analysis identified distinct clusters associated with fruit skin and flesh color. The genome sequencing of guava, functional annotation, comparative transcriptomics-based genome-wide markers, and genetic diversity analysis will expand the knowledge of genomes of climacteric fruits, facilitating trait-based molecular breeding and diversifying the nutritional basket.

Advances in genomics and genome editing for breeding next generation of fruit and nut crops

Fruit tree crops are an essential part of the food production systems and are key to achieve food and nutrition security. Genetic improvement of fruit trees by conventional breeding has been slow due to the long juvenile phase. Advancements in genomics and molecular biology have paved the way for devising novel genetic improvement tools like genome editing, which can accelerate the breeding of these perennial crops to a great extent. In this article, advancements in genomics of fruit trees covering genome sequencing, transcriptome sequencing, genome editing technologies (GET), CRISPR-Cas system based genome editing, potential applications of CRISPR-Cas9 in fruit tree crops improvement, the factors influencing the CRISPR-Cas editing efficiency and the challenges for CRISPR-Cas9 applications in fruit tree crops improvement are reviewed. Besides, base editing, a recently emerging more precise editing system, and the future perspectives of genome editing in the improvement of fruit and nut crops are covered.