The pomegranate (Punica granatum L.) genome and the genomics of punicalagin biosynthesis

Integrated transcriptomic, metabolomic, and functional analyses unravel the mechanism of bagging delaying fruit cracking of pomegranate (Punica granatum L.)

The economic impact of fruit cracking in pomegranate products is substantial. In this study, we present the inaugural comprehensive analysis of transcriptome and metabolome in the outermost pericarp of pomegranate fruit in bagging conditions. Our investigation revealed a notable upregulation of differentially expressed genes (DEGs) associated with the calcium signaling pathway (76.92%) and xyloglucan endotransglucosylase/hydrolase (XTH) genes (87.50%) in the fruit peel of non-cracking fruit under bagging. Metabolomic analysis revealed that multiple phenolics, flavonoids, and tannins were identified in pomegranate. Among these, calmodulin-like 23 (PgCML23) exhibited a significant correlation with triterpenoids and demonstrated a marked upregulation under bagging treatment. The transgenic tomatoes overexpressing PgCML23 exhibited significantly higher cellulose content and xyloglucan endotransglucosylase (XET) enzyme activity in the pericarp at the red ripening stage compared to the wild type. Conversely, water-soluble pectin content, polygalacturonase (PG), and β-galactosidase (β-GAL) enzyme activities were significantly lower in the transgenic tomatoes. Importantly, the heterologous expression of PgCML23 led to a substantial reduction in the fruit cracking rate in tomatoes. Our findings highlight the reduction of fruit cracking in bagging conditions through the manipulation of PgCML23 expression.

A Genome-Wide Analysis of the BAM Gene Family and Identification of the Cold-Responsive Genes in Pomegranate (Punica granatum L.)

Beta-amylases (BAMs, EC 3.2.1.2), belonging to a multigene family, play a pivotal role in starch breakdown and are also involved in hormonal and stress responses, notably to cold stress. Pomegranate trees (Punica granatum L.) are adapted to warm climates and are sensitive to cold temperatures. In this study, we analyzed eight PgBAM genes from the pomegranate genome dataset. These members unevenly distributed across chromosomes and were categorized into four groups based on their orthologous members. The motif composition was highly consistent among most members. In contrast, exon numbers and arrangements were conserved within groups or subgroups, whereas significant diversity was observed between different groups. A syntenic analysis revealed that three PgBAM members (PgBAM1/4/5) showed a total of 11 syntenic relationships with the BAM members from Arabidopsis, kiwifruit, and Chinese white pear, respectively. Promoter binding motif prediction suggested potential roles for PgBAMs' genes in light, stress, hormones, and development signaling. Gene expression indicated that PgBAM4 was predominantly expressed in leaves, PgBAM7 in flowers, and PgBAM8 in roots and leaves and during fruit ripening, particularly in pericarp development. A transcriptome analysis identified the starch and sucrose metabolism pathway (map00500) as a key factor in the cold stress response of cold-sensitive cultivar 'Tunisia' seedlings. PgBAM4 exhibited remarkable expression and was closely associated with the cold-responsive BAM genes, characterized by a closer phylogenetic relationship, conserved catalytic residues, and similar secondary and tertiary structures. Moreover, the differences in soluble sugar levels and PgBAM4 expression were closely associated with the varying cold stress resistance observed between 'Tunisia' and 'Sanbai' seedlings. Furthermore, yeast one-hybrid assays confirmed that PgCBF7, a critical transcription factor for enhancing freezing tolerance, binds to the promoter region of PgBAM4. Our findings provide a systematic overview of the PgBAM gene family and shed new light on the regulatory mechanisms underlying cold stress tolerance in pomegranate.

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.

Identification of Candidate Expansin Genes Associated with Seed Weight in Pomegranate (Punica granatum L.)

Seed weight is an important target trait in pomegranate breeding and culture. Expansins act by loosening plant cell walls and cellulosic materials, permitting turgor-driven cell enlargement. However, the role of expansin genes (EXPs) in pomegranate seed weight remains elusive. A total of 29 PgrEXPs were identified in the 'Dabenzi' genome. These genes were classified into four subfamilies and 14 subgroups, including 22 PgrEXPAs, 5 PgrEXPBs, 1 PgrEXPLA, and 1 PgrEXPLB. Transcriptome analysis of PgrEXPs in different tissues (root, leaf, flower, peel, and seed testa) in 'Dabenzi', and the seed testa of the hard-seeded pomegranate cultivar 'Dabenzi' and soft-seeded cultivar 'Tunisia' at three development stages showed that three PgrEXPs (PgrEXPA11, PgrEXPA22, PgrEXPA6) were highly expressed throughout seed development, especially in the sarcotesta. SNP/Indel markers of these PgrEXPs were developed and used to genotype 101 pomegranate accessions. The association of polymorphic PgrEXPs with seed weight-related traits (100-seed weight, 100-kernel weight, 100-sarcotesta weight, and the percentage of 100-sarcotesta to 100-seed weight) were analyzed. PgrEXP22 was significantly associated with 100-seed weight and 100-sarcotesta weight and is a likely candidate for regulating seed weight and sarcotesta development in particular. This study provides an effective tool for the genetic improvement of seed weight in pomegranate breeding programs.

Genome-Wide Analysis and Expression Profiling of Trehalose-6-Phosphate Phosphatase (TPP) in Punica granatum in Response to Abscisic-Acid-Mediated Drought Stress

Trehalose, a nonreducing disaccharide, has been linked to plant growth and development as well as stress response. The enzyme trehalose-6-phosphate phosphatase (TPP) plays a crucial role in the production of trehalose in higher plants. This study identified a total of seven TPP family genes within the pomegranate species (PgTPP1-PgTPP7). Three subgroups of the seven PgTPPs were identified through phylogenetic analysis. The gene length, coding sequence (CD) length, and chromosomal location of the PgTPP genes were studied. In addition, the PgTPP proteins' length, isoelectric point (Ip), grand average of hydropathicity (GRAVY), conserved domains, conserved motifs, synteny, and phylogenetic relationships with Arabidopsis and tomato TPP proteins were examined. The cis-acting elements in the promoter region and the expression of the PgTPP genes under abscisic acid (ABA)-mediated drought stress as well as the differences in expression in the root, flower, and leaf tissues were also assessed. The PgTPP2 and PgTPP5 genes are involved in the response to abscisic-acid-mediated drought stress, as shown by drought-mediated stress transcriptomes. The PgTPP1 and PgTPP2 genes were expressed only in floral tissue and roots, respectively. The remaining PgTPPs did not exhibit any significant alterations in gene expression in roots, flowers, or leaves. The current study has the potential to provide a comprehensive understanding of the biological characteristics of PgTPP proteins in various developmental processes and their role in the pomegranate plant's response to different stressors. However, further research is required to explore their precise biological role. Hence, conducting a comprehensive functional validation study on PgTPPs could contribute to the development of stress-resistant agricultural cultivars.

Assembly and analysis of the first complete mitochondrial genome of Punica granatum and the gene transfer from chloroplast genome

Pomegranate (Punica granatum L.) is one of the oldest fruits with edible, medicinal and ornamental values. However, there is no report on the mitochondrial genome of pomegranate. In this study, the mitochondrial genome of P. granatum was sequenced, assembled and analyzed in detail, while the chloroplast genome was assembled using the same set of data. The results showed that the P. granatum mitogenome had a multi branched structure, using BGI + Nanopore mixed assembly strategy. The total genome length was 404,807 bp, with the GC content of 46.09%, and there were 37 protein coding genes, 20 tRNA genes and three rRNA genes. In the whole genome, 146 SSRs were identified. Besides, 400 pairs of dispersed repeats were detected, including 179 palindromic, 220 forward and one reverse. In the P. granatum mitochondrial genome, 14 homologous fragments of chloroplast genome were found, accounting for 0.54% of the total length. Phylogenetic analysis showed that among the published mitochondrial genomes of related genera, P. granatum had the closest genetic relationship with Lagerstroemia indica of Lythraceae. The 580 and 432 RNA editing sites were predicted on 37 protein coding genes of mitochondrial genome using BEDTools software and online website PREPACT respectively, but all were from C to U, of which ccmB and nad4 gene were most frequently edited, with 47 sites. This study provides a theoretical basis for understanding the evolution of higher plants, species classification and identification, and will also be useful for further utilization of pomegranate germplasm resources.

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.

Two major chromosome evolution events with unrivaled conserved gene content in pomegranate

Pomegranate has a unique evolutionary history given that different cultivars have eight or nine bivalent chromosomes with possible crossability between the two classes. Therefore, it is important to study chromosome evolution in pomegranate to understand the dynamics of its population. Here, we de novo assembled the Azerbaijani cultivar "Azerbaijan guloyshasi" (AG2017; 2n = 16) and re-sequenced six cultivars to track the evolution of pomegranate and to compare it with previously published de novo assembled and re-sequenced cultivars. High synteny was observed between AG2017, Bhagawa (2n = 16), Tunisia (2n = 16), and Dabenzi (2n = 18), but these four cultivars diverged from the cultivar Taishanhong (2n = 18) with several rearrangements indicating the presence of two major chromosome evolution events. Major presence/absence variations were not observed as >99% of the five genomes aligned across the cultivars, while >99% of the pan-genic content was represented by Tunisia and Taishanhong only. We also revisited the divergence between soft- and hard-seeded cultivars with less structured population genomic data, compared to previous studies, to refine the selected genomic regions and detect global migration routes for pomegranate. We reported a unique admixture between soft- and hard-seeded cultivars that can be exploited to improve the diversity, quality, and adaptability of local pomegranate varieties around the world. Our study adds body knowledge to understanding the evolution of the pomegranate genome and its implications for the population structure of global pomegranate diversity, as well as planning breeding programs aiming to develop improved cultivars.

Flavan-3-ol Galloylation-Related Functional Gene Cluster and the Functional Diversification of SCPL Paralogs in Camellia sp

The high accumulation of galloylated flavan-3-ols in Camellia sp. is a noteworthy phenomenon. We identified a flavan-3-ol galloylation-related functional gene cluster in tannin-rich plant Camellia sp., which included UGT84A22 and SCPL-AT gene clusters. We investigated the possible correlation between the accumulation of metabolites and the expression of SCPL-ATs and UGT84A22. The results revealed that C. sinensis, C. ptilophylla, and C. oleifera accumulated galloylated cis-flavan-3-ols (EGCG), galloylated trans-flavan-3-ols (GCG), and hydrolyzed tannins, respectively; however, C. nitidissima did not accumulate any galloylated compounds. C. nitidissima exhibited no expression of SCPL-AT or UGT84A22, whereas the other three species of Camellia exhibited various expression patterns. This indicated that the functions of the paralogs of SCPL-AT vary. Enzymatic analysis revealed that SCPL5 was neofunctionalized as a noncatalytic chaperone paralog, a type of chaerone-like protein, associating with flavan-3-ol galloylation; moreover, CsSCPL4 was subfunctionalized in association with the galloylation of cis- and trans-flavan-3-ols. In C. nitidissima, an SCPL4 homolog was noted with mutations in two cysteine residues forming a disulfide bond, which suggested that this homolog was defunctionalized. The findings of this study improve our understanding of the functional diversification of SCPL paralogs in Camellia sp.

Identification and Functional Analysis of CAD Gene Family in Pomegranate (Punica granatum)

[Objective] Cinnamyl alcohol dehydrogenase (CAD) is a key enzyme in lignin biosynthesis. The aim of this study was to identify CAD gene family members in pomegranate and its expression correlation with seed hardness. [Methods] Based on the reported CAD sequence of Arabidopsis, the CAD gene family of pomegranate was identified by homologous comparison, and then phylogenetic, molecular characterization, and expression profile analysis were performed. [Results] Pomegranate CAD gene family has 25 members, distributed on seven chromosomes of pomegranate. All pomegranate CAD proteins have similar physical and chemical properties. We divide the family into four groups based on evolutionary relationships. The member of group I, called bona fide CAD, was involved in lignin synthesis. Most of the members of group II were involved in stress resistance. The functions of groups III and IV need to be explored. We found four duplicated modes (whole genome duplication or segmental (WGD), tandem duplication (TD), dispersed duplication (DSD), proximal duplication (PD) in this family; TD (36%) had the largest number of them. We predicted that 20 cis-acting elements were involved in lignin synthesis, stress resistance, and response to various hormones. Gene expression profiles further demonstrated that the PgCAD gene family had multiple functions. [Conclusions] Pomegranate CAD gene family is involved in lignin synthesis of hard-seeded cultivar Hongyushizi and Baiyushizi, but its role in seed hardness of soft-seeded cultivar Tunisia needs to be further studied.

Transcriptional profiling of long non-coding RNAs regulating fruit cracking in Punica granatum L. under bagging

Fruit cracking tremendously damages the appearance of fruit, easily leads to pathogen invasion, greatly reduces the marketability and causes immense economic losses. The pivotal role of long non-coding RNAs (lncRNAs) in diverse biological processes has been confirmed, while the roles of lncRNAs underlying fruit cracking remain poorly understood. In this study, the incidence of fruit cracking was 7.26% under the bagging treatment, the control group was 38.11%, indicating that bagging considerably diminished the fruit cracking rate. LncRNA libraries for fruit cracking (FC), fruit non-cracking (FNC) and fruit non-cracking under bagging (FB) in pomegranate (Punica granatum L.) were performed and analysed via high-throughput transcriptome sequencing. A total of 3194 lncRNAs were obtained with a total length of 4898846 nt and an average length of 1533.77 nt in pomegranate. We identified 42 differentially expressed lncRNAs (DELs) and 137 differentially expressed mRNAs (DEGs) in FC vs FNC and 35 DELs and 160 DEGs in FB vs FC that formed co-expression networks respectively, suggesting that there are involved in phytohormone signaling pathway, lignin catabolic process, lipid transport/binding, cutin biosynthetic process and cell wall organization. We also found that 18 cis-acting DELs regulated 18 target genes, and 10 trans-acting DELs regulated 24 target genes in FC vs FNC, 23 DELs regulate 23 target genes for the cis-acting lncRNAs and 12 DELs regulated 36 target genes in FB vs FC, which provides an understanding for the regulation of the fruit cracking. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis results demonstrated that DELs participated in calcium ion binding, glycerophospholipid metabolism, flavonoid biosynthetic process, cell wall biogenesis, xyloglucan metabolic process, hormone signal transduction and starch and sucrose metabolism. Our findings provide new insights into the roles of lncRNAs in regulating the fruit cracking and lay the foundation for further improvement of pomegranate quality.

Characterization of the ABC Transporter G Subfamily in Pomegranate and Function Analysis of PgrABCG14

ATP-binding cassette subfamily G (ABCG) proteins play important roles in plant growth and development by transporting metabolites across cell membranes. To date, the genetic characteristics and potential functions of pomegranate ABCG proteins (PgrABCGs) have remained largely unknown. In this study, we found that 47 PgrABCGs were divided into five groups according to a phylogenetic analysis; groups I, II, III, and IV members are half-size proteins, and group V members are full-size proteins. PgrABCG14, PgrABCG21, and PgrABCG47 were highly expressed in the inner seed coat but had very low expression levels in the outer seed coat, and the expression levels of these three PgrABCG genes in the inner seed coats of hard-seeded pomegranate ‘Dabenzi’ were higher than those of soft-seeded pomegranate ‘Tunisia’. In addition, the expression of these three PgrABCG genes was highly correlated with the expression of genes involved in lignin biosynthesis and hormone signaling pathways. The evolution of PgrABCG14 presents a highly similar trend to the origin and evolution of lignin biosynthesis during land plant evolution. Ectopic expression of PgrABCG14 in Arabidopsis promoted plant growth and lignin accumulation compared to wild type plants; meanwhile, the expression levels of lignin biosynthesis-related genes (CAD5, C4H, and Prx71) and cytokinin response marker genes (ARR5 and ARR15) were significantly upregulated in transgenic plants, which suggests the potential role of PgrABCG14 in promoting plant growth and lignin accumulation. Taken together, these findings not only provide insight into the characteristics and evolution of PgrABCGs, but also shed a light on the potential functions of PgrABCGs in seed hardness development.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

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.

Reference quality genome sequence of Indian pomegranate cv. 'Bhagawa' (Punica granatum L.)

Pomegranate is an important fruit crop for ensuring livelihood and nutrition security in fragile semi-arid regions of the globe having limited irrigation resources. This is a high-value, nutritionally rich, and export-oriented agri-commodity that ensures high returns on investment to growers across the world. Although it is a valuable fruit crop, it has received only a limited genomics research outcome. To fast-track the pomegranate improvement program, de novo whole-genome sequencing of the main Indian cultivar 'Bhagawa' was initiated by the Indian Council of Agricultural Research-National Research Center on Pomegranate (ICAR-NRCP). We have demonstrated that a combination of commercially available technologies from Illumina, PacBio, 10X Genomics, and BioNano Genomics could be used efficiently for sequencing and reference-grade de novo assembly of the pomegranate genome. The research led to a final reference-quality genome assembly for 'Bhagawa' of 346.08 Mb in 342 scaffolds and an average N50 of 16.12 Mb and N90 of 1088.62 Kb. This assembly covered more than 98% of the estimated pomegranate genome size, 352.54 Mb. The LTR assembly index (LAI) value of 10 and 93.68% Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness score over the 1,440 ortholog genes of the completed pomegranate genome indicates the quality of the assembled pomegranate genome. Furthermore, 29,435 gene models were discovered with a mean transcript length of 2,954 bp and a mean coding sequence length 1,090 bp. Four transcript data samples of pomegranate tissues were mapped over the assembled 'Bhagawa' genome up to 95% significant matches, indicating the high quality of the assembled genome. We have compared the 'Bhagawa' genome with the genomes of the pomegranate cultivars 'Dabenzi' and 'Taishanhong.' We have also performed whole-genome phylogenetic analysis using Computational Analysis of Gene Family Evolution (CAFE) and found that Eucalyptus grandis and pomegranate diverged 64 (60-70) million years ago. About 1,573 protein-coding resistance genes identified in the 'Bhagawa' genome were classified into 32 domains. In all, 314 copies of miRNA belonging to 26 different families were identified in the 'Bhagawa' genome. The reference-quality genome assembly of 'Bhagawa' is certainly a significant genomic resource for accelerated pomegranate improvement.

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.

Identification of sugar transporter (SWEET) genes involved in pomegranate seed coat sugar accumulation

Sugar content of the outer seed coat and hardness of the inner seed coat are important traits of the pomegranate fruit. The translocation of sugars across biological membranes, mediated by SWEET transporters, is critical to seed development. In this study, we identified 16 PgrSWEET genes distributed on six chromosomes in the pomegranate genome. According to the phylogenetic analysis, PgrSWEET proteins were divided into four groups. Tandem and segmental duplications contributed to the expansion of the PgrSWEET family, while functional redundancy and diversification may have occurred among SWEET members according to analyses of evolution and gene expression. RNA-seq and qRT-PCR analyses revealed that PgrSWEET1a and PgrSWEET9 were highly expressed in the inner seed coat, and the expression levels gradually increased during seed development. Moreover, the relative expression levels of PgrSWEET1a and PgrSWEET9 in a hard-seeded cultivar were higher than those in a soft-seeded cultivar, indicating that PgrSWEET1a and PgrSWEET9 might function in the inner seed coat development by accumulating sugar metabolites. We also found that PgrSWEET2 was highly expressed in the outer seed coat during seed development, and the protein was localized to the tonoplast, indicating that PgrSWEET2 is likely a candidate regulating sugar accumulation or reutilization in the vacuoles of the outer seed coat. Genes encoding transcription factors probably regulating the candidate PgrSWEET genes were chosen by co-expression analysis. These results not only helped to characterize PgrSWEET genes but also provided an insight into their functions in relation to seed coat development.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03248-6.

Chromosome-specific potential intron polymorphism markers for large-scale genotyping applications in pomegranate

Despite the availability of whole genome assemblies, the identification and utilization of gene-based marker systems has been limited in pomegranate. In the present study, we performed a genome-wide survey of intron length (IL) markers in the 36,524 annotated genes of the Tunisia genome. We identified and designed a total of 8,812 potential intron polymorphism (PIP) markers specific to 3,445 (13.40%) gene models that span 8 Tunisia chromosomes. The ePCR validation of all these PIP markers on the Tunisia genome revealed single-locus amplification for 1,233 (14%) markers corresponding to 958 (27.80%) genes. The markers yielding single amplicons were then mapped onto Tunisia chromosomes to develop a saturated linkage map. The functional categorization of 958 genes revealed them to be a part of the nucleus and the cytoplasm having protein binding and catalytic activity, and these genes are mainly involved in the metabolic process, including photosynthesis. Further, through ePCR, 1,233 PIP markers were assayed on multiple genomes, which resulted in the identification of 886 polymorphic markers with an average PIC value of 0.62. In silico comparative mapping based on physically mapped PIP markers indicates a higher synteny of Tunisia with the Dabenzi and Taishanhong genomes (>98%) in comparison with the AG2017 genome (95%). We then performed experimental validation of a subset of 100 PIP primers on eight pomegranate genotypes and identified 76 polymorphic markers, with 15 having PIC values ≥0.50. We demonstrated the potential utility of the developed markers by analyzing the genetic diversity of 31 pomegranate genotypes using 24 PIP markers. This study reports for the first time large-scale development of gene-based and chromosome-specific PIP markers, which would serve as a rich marker resource for genetic variation studies, functional gene discovery, and genomics-assisted breeding of pomegranate.

Ripening of Pomegranate Skin as Revealed by Developmental Transcriptomics

The appearance of pomegranate (Punica granatum L.) fruit is highly important for its marketing. The primary concerns are obtaining sufficient red pigment accumulation and minimal cracking of the fruit skin (the outer red layer of the peel). We analyzed the skin transcriptome of pomegranate cv. Wonderful at distinct time points of fruit development to characterize the processes that occur in the skin during fruit ripening and which may reflect on processes in the whole fruit, such as the non-climacteric nature of pomegranate. The data suggested a ripening mechanism in pomegranate skin that differs from that in strawberry—the model plant for non-climacteric fruit where abscisic acid is the growth regulator that drives ripening—involving ethylene, polyamine, and jasmonic acid pathways. The biosynthetic pathways of important metabolites in pomegranate—hydrolyzable tannins and anthocyanins—were co-upregulated at the ripening stage, in line with the visual enhancement of red coloration. Interestingly, cuticle- and cell-wall-related genes that showed differential expression between the developmental stages were mainly upregulated in the skin of early fruit, with lower expression at mid-growth and ripening stages. Nevertheless, lignification may be involved in skin hardening in the mature fruit.

Physical Characteristics, Mineral Content, and Antioxidant and Antibacterial Activities of Punica granatum or Citrus sinensis Peel Extracts and Their Applications to Improve Cake Quality

One-third of all food produced for human use is discarded as waste, resulting in environmental pollution and impaired food security. Fruit peels have bioactive compounds that may be used as antimicrobials and antioxidants, and the use of fruit peels is considered an alternative way to reduce environmental problems and agro-industrial waste. The aim of this study was to evaluate the phytochemical, mineral, extraction yield, total phenolic, total flavonoids, antioxidant, and antibacterial activity of several peel fruits, including Citrus sinensis (orange) and Punica granatum (pomegranate). The results revealed that pomegranate peel powder contains the highest amounts of ash, fiber, total carbohydrates, Ca, Fe, Mg, and Cu, while orange peel contains the highest amounts of moisture, protein, crude fat, P, and K. Furthermore, the aqueous and methanolic pomegranate peel extracts yielded higher total phenolic and total flavonoids than the orange peel extract. The identification and quantification of polyphenol compounds belonging to different classes, such as tannins, phenolic acids, and flavonoids in pomegranate peel and flavonoid compounds in orange peel were performed using UPLC-MS/MS. In addition, GC-MS analysis of orange peel essential oil discovered that the predominant compound is D-Limonene (95.7%). The aqueous and methanolic extracts of pomegranate peel were proven to be efficient against both gram-positive and gram-negative bacteria linked to human infections. Sponge cake substituting wheat flour with 3% pomegranate peel and 10% orange peel powder had the highest total phenolic, flavonoid compounds, and antioxidant activity as compared to the control cake. Our results concluded that pomegranate and orange peel flour can be used in cake preparation and natural food preservers.

Transcriptomics and metabolomics analyses identified key genes associated with sugar and acid metabolism in sweet and sour pomegranate cultivars during the developmental period

Pomegranate (Punica granatum), an important fruit tree in the world, is rich in bioactive substances and has broad prospects for development. In this study, gene expression levels and the concentrations of metabolites involved in the metabolism of soluble sugars and organic acids were investigated in sweet and sour pomegranate cultivars at the S1 (July 25) stage, S2 (August 26) stage, and S3 (September 24) stage. The results showed that glucose, fructose, citric acid, and malic acid were predominantly present in pomegranate. The expression of invertase 2 (INV2), INV1, FRK2, FRK7, PFK2, PFK7, and HK1 was closely correlated with the fructose and glucose contents during different developmental stages, whereas the expression of sucrose synthase 3 (SUS3) and INV1 was negatively correlated with the sucrose content. The expression of MDH (c28468_g3) and WRKY42 (c20711_g1) genes were closely related to the content of sucrose, malic acid, citric acid, and succinic acid during different developmental stages. Gene expression and metabolite concentrations varied between the two cultivars. The results provide valuable information for gene discovery, marker-assisted selection, and investigation of metabolism mechanisms in pomegranate fruits.

Hybrid Assembly and Annotation of the Genome of the Indian Punica granatum, a Superfood

The wonder fruit pomegranate (Punica granatum, family Lythraceae) is one of India’s economically important fruit crops that can grow in different agro-climatic conditions ranging from tropical to temperate regions. This study reports high-quality de novo draft hybrid genome assembly of diploid Punica cultivar “Bhagwa” and identifies its genomic features. This cultivar is most common among the farmers due to its high sustainability, glossy red color, soft seed, and nutraceutical properties with high market value. The draft genome assembly is about 361.76 Mb (N50 = 40 Mb), ∼9.0 Mb more than the genome size estimated by flow cytometry. The genome is 90.9% complete, and only 26.68% of the genome is occupied by transposable elements and has a relative abundance of 369.93 SSRs/Mb of the genome. A total of 30,803 proteins and their putative functions were predicted. Comparative whole-genome analysis revealed Eucalyptus grandis as the nearest neighbor. KEGG-KASS annotations indicated an abundance of genes involved in the biosynthesis of flavonoids, phenylpropanoids, and secondary metabolites, which are responsible for various medicinal properties of pomegranate, including anticancer, antihyperglycemic, antioxidant, and anti-inflammatory activities. The genome and gene annotations provide new insights into the pharmacological properties of the secondary metabolites synthesized in pomegranate. They will also serve as a valuable resource in mining biosynthetic pathways for key metabolites, novel genes, and variations associated with disease resistance, which can facilitate the breeding of new varieties with high yield and superior quality.

The Pomegranate Deciduous Trait Is Genetically Controlled by a PgPolyQ-MADS Gene

The pomegranate (Punica granatum L.) is a deciduous fruit tree that grows worldwide. However, there are variants, which stay green in mild winter conditions and are determined evergreen. The evergreen trait is of commercial and scientific importance as it extends the period of fruit production and provides opportunity to identify genetic functions that are involved in sensing environmental cues. Several different evergreen pomegranate accessions from different genetic sources grow in the Israeli pomegranate collection. The leaves of deciduous pomegranates begin to lose chlorophyll during mid of September, while evergreen accessions continue to generate new buds. When winter temperature decreases 10°C, evergreen variants cease growing, but as soon as temperatures arise budding starts, weeks before the response of the deciduous varieties. In order to understand the genetic components that control the evergreen/deciduous phenotype, several segregating populations were constructed, and high-resolution genetic maps were assembled. Analysis of three segregating populations showed that the evergreen/deciduous trait in pomegranate is controlled by one major gene that mapped to linkage group 3. Fine mapping with advanced F3 and F4 populations and data from the pomegranate genome sequences revealed that a gene encoding for a putative and unique MADS transcription factor (PgPolyQ-MADS) is responsible for the evergreen trait. Ectopic expression of PgPolyQ-MADS in Arabidopsis generated small plants and early flowering. The deduced protein of PgPolyQ-MADS includes eight glutamines (polyQ) at the N-terminus. Three-dimensional protein model suggests that the polyQ domain structure might be involved in DNA binding of PgMADS. Interestingly, all the evergreen pomegranate varieties contain a mutation within the polyQ that cause a stop codon at the N terminal. The polyQ domain of PgPolyQ-MADS resembles that of the ELF3 prion-like domain recently reported to act as a thermo-sensor in Arabidopsis, suggesting that similar function could be attributed to PgPolyQ-MADS protein in control of dormancy. The study of the evergreen trait broadens our understanding of the molecular mechanism related to response to environmental cues. This enables the development of new cultivars that are better adapted to a wide range of climatic conditions.

Novel miRNA-SSRs for Improving Seed Hardness Trait of Pomegranate (Punica granatum L.)

Present research discovered novel miRNA-SSRs for seed type trait from 761 potential precursor miRNA sequences of pomegranate. SSR mining and BLASTx of the unique sequences identified 69 non-coding pre-miRNA sequences, which were then searched for BLASTn homology against Dabenzi genome. Sixty three true pri-miRNA contigs encoding 213 pre-miRNAs were predicted. Analysis of the resulting sequences enabled discovery of SSRs within pri-miRNA (227) and pre-miRNA sequences (79). A total of 132 miRNA-SSRs were developed for seed type trait from 63 true pri-miRNAs, of which 46 were specific to pre-miRNAs. Through ePCR, 123 primers were validated and mapped on eight Tunisia chromosomes. Further, 80 SSRs producing specific amplicons were ePCR-confirmed on multiple genomes i.e. Dabenzi, Taishanhong, AG2017 and Tunisia, yielding a set of 63 polymorphic SSRs (polymorphism information content ≥0.5). Of these, 32 miRNA-SSRs revealed higher polymorphism level (89.29%) when assayed on six pomegranate genotypes. Furthermore, target prediction and network analysis suggested a possible association of miRNA-SSRs i.e. miRNA_SH_SSR69, miRNA_SH_SSR36, miRNA_SH_SSR103, miRNA_SH_SSR35 and miRNA_SH_SSR53 with seed type trait. These miRNA-SSRs would serve as important genomic resource for rapid and targeted improvement of seed type trait of pomegranate.

Genome-wide identification and characterization of bZIP gene family and cloning of candidate genes for anthocyanin biosynthesis in pomegranate (Punica granatum)

BACKGROUND

The basic leucine zipper (bZIP) transcription factor is one of the most abundant and conserved gene families in eukaryotes. In addition to participating in plant development and growth, bZIP transcription factors play crucial roles in various abiotic stress responses and anthocyanin accumulation. Up to now, analysis of bZIP gene family members in pomegranate (Punica granatum) has not been reported. Three published pomegranate genome sequences provide valuable resources for further gene function analysis.

RESULTS

Using bioinformatics analysis, 65 PgbZIPs were identified and analyzed from the 'Taishanhong' pomegranate genome. We divided them into 13 groups (A, B, C, D, E, F, G, H, I, J, K, M, and S) according to the phylogenetic relationship with those of Arabidopsis, each containing a different number of genes. The regularity of exon/intron number and distribution was consistent with the classification of groups in the evolutionary tree. Transcriptome analysis of different tissues showed that members of the PgbZIP gene family were differentially expressed in different developmental stages and tissues of pomegranate. Among them, we selected PgbZIP16 and PgbZIP34 as candidate genes which affect anthocyanin accumulation. The full-length CDS region of PgbZIP16 and PgbZIP34 were cloned from pomegranate petals by homologous cloning technique, encoding 170 and 174 amino acids, which were 510 bp and 522 bp, respectively. Subcellular localization assays suggested that both PgbZIP16 and PgbZIP34 were nucleus-localized. Real-time quantitative PCR (qPCR) was used to explore the expression of PgbZIP16 and PgbZIP34 in the petals of three kinds of ornamental pomegranates at the full flowering stage. The results demonstrated that the expression of PgbZIP16 in red petals was 5.83 times of that in white petals, while PgbZIP34 was 3.9 times. The results of transient expression in tobacco showed that consistent trends were observed in anthocyanin concentration and expression levels of related genes, which both increased and then decreased. Both PgbZIP16 and PgbZIP34 could promote anthocyanin accumulation in tobacco leaves. We obtained transgenic strains overexpressing PgbZIP16, and the histochemical staining for GUS activity showed that overexpressed PgbZIP16 seedlings were expressed in the stem. Transgenic experiments indicated that overexpression of PgbZIP16 significantly upregulated UF3GT, ANS and DFR genes in Arabidopsis and enhanced anthocyanin accumulation.

CONCLUSIONS

The whole genome identification, gene structure, phylogeny, gene cloning, subcellular location and functional verification of the pomegranate bZIP gene family provide a theoretical foundation for the functional study of the PgbZIP gene family and candidate genes for anthocyanin biosynthesis.

Systematic Analysis and Expression Profiles of the 4-Coumarate: CoA Ligase (4CL) Gene Family in Pomegranate (Punica granatum L.)

4-Coumarate:CoA ligase (4CL, EC6.2.1.12), located at the end of the phenylpropanoid metabolic pathway, regulates the metabolic direction of phenylpropanoid derivatives and plays a pivotal role in the biosynthesis of flavonoids, lignin, and other secondary metabolites. In order to understand the molecular characteristics and potential biological functions of the 4CL gene family in the pomegranate, a bioinformatics analysis was carried out on the identified 4CLs. In this study, 12 Pg4CLs were identified in the pomegranate genome, which contained two conserved amino acid domains: AMP-binding domain Box I (SSGTTGLPKGV) and Box II (GEICIRG). During the identification, it was found that Pg4CL2 was missing Box II. The gene cloning and sequencing verified that this partial amino acid deletion was caused by genome sequencing and splicing errors, and the gene cloning results corrected the Pg4CL2 sequence information in the ‘Taishanhong’ genome. According to the phylogenetic tree, Pg4CLs were divided into three subfamilies, and each subfamily had 1, 1, and 10 members, respectively. Analysis of cis-acting elements found that all the upstream sequences of Pg4CLs contained at least one phytohormone response element. An RNA-seq and protein interaction network analysis suggested that Pg4CL5 was highly expressed in different tissues and may participate in lignin synthesis of pomegranate. The expression of Pg4CL in developing pomegranate fruits was analyzed by quantitative real-time PCR (qRT-PCR), and the expression level of Pg4CL2 demonstrated a decreasing trend, similar to the trend of flavonoid content, indicating Pg4CL2 may involve in flavonoid synthesis and pigment accumulation. Pg4CL3, Pg4CL7, Pg4CL8, and Pg4CL10 were almost not expressed or lowly expressed, the expression level of Pg4CL4 was higher in the later stage of fruit development, suggesting that Pg4CL4 played a crucial role in fruit ripening. The expression levels of 4CL genes were significantly different in various fruit development stages. The results laid the foundation for an in-depth analysis of pomegranate 4CL gene functions.

Identification and Characterization of Two Regiospecific Tricetin UDP-Dependent Glycosyltransferases from Pomegranate (Punica granatum L.)

Tricetin (5,7,3′,4′,5′-pentahydroxyflavone) is a dietary flavone from flowers of Myrtales plants with demonstrated functions in promoting human health. By contrast, the bioactivity of its glucosylated derivative tricetin 4′-O-glucoside has not been extensively explored. We conducted metabolite profiling analysis of pomegranate (a Myrtales plant) floral tissues and revealed that tricetin and tricetin 4′-O-glucoside accumulate in anthers, but not petals. In addition, the comparative analysis of anther and petal transcriptomes identified 10 UGTs that are more highly expressed in anthers than petals. Of the 10 UGTs, PgUGT76Z1 and PgUGT73AL1 glucosylated specifically at the 4′-O position of tricetin to form tricetin 4′-O-glucoside. The phylogenetic analysis indicated that PgUGT76Z1 and PgUGT73AL1 belong to different plant UGT groups, suggesting a convergent evolution of these tricetin UGTs. Overall, identification and characterization of PgUGT76Z1 and PgUGT73AL1 not only provides evolutionary insights into tricetin glucosylation, but also offers an opportunity to produce tricetin 4′-O-glucoside in large quantities through microbial biotransformation or plant metabolic engineering, thus facilitating the investigation of tricetin 4′-O-glucoside bioactivities.

ARF6s Identification and Function Analysis Provide Insights Into Flower Development of Punica granatum L

Based on the genome and small-RNA sequencing of pomegranate, miRNA167 and three target genes PgARF6 were identified in "Taishanhong" genome. Three PgARF6 genes and their corresponding protein sequences, expression patterns in pomegranate flower development and under exogenous hormones treatments were systematically analyzed in this paper. We found that PgARF6s are nuclear proteins with conserved structures. However, PgARF6s had different protein structures and expression profiles in pomegranate flower development. At the critical stages of pomegranate ovule sterility (8.1-14.0 mm), the expression levels of PgARF6s in bisexual flowers were lower than those in functional male flowers. Interestingly, PgARF6c expression level was significantly higher than PgARF6a and PgARF6b. Under the treatment of exogenous IBA and 6-BA, PgARF6s were down-regulated, and the expression of PgARF6c was significantly inhibited. PgmiR167a and PgmiR167d had the binding site on PgARF6 genes sequences, and PgARF6a has the directly targeted regulatory relationship with PgmiR167a in pomegranate. At the critical stage of ovule development (8.1-12.0 mm), exogenous IBA and 6-BA promoted the content of GA and ZR accumulation, inhibited BR accumulation. There was a strong correlation between the expression of PgARF6a and PgARF6b. Under exogenous hormone treatment, the content of ZR, BR, GA, and ABA were negatively correlated with the expressions of PgARF6 genes. However, JA was positively correlated with PgARF6a and PgARF6c under IBA treatment. Thus, our results provide new evidence for PgARF6 genes involving in ovule sterility in pomegranate flowers.

Identification, Analysis and Gene Cloning of the SWEET Gene Family Provide Insights into Sugar Transport in Pomegranate (Punica granatum)

Members of the sugars will eventually be exported transporter (SWEET) family regulate the transport of different sugars through the cell membrane and control the distribution of sugars inside and outside the cell. The SWEET gene family also plays important roles in plant growth and development and physiological processes. So far, there are no reports on the SWEET family in pomegranate. Meanwhile, pomegranate is rich in sugar, and three published pomegranate genome sequences provide resources for the study of the SWEET gene family. 20 PgSWEETs from pomegranate and the known Arabidopsis and grape SWEETs were divided into four clades (Ⅰ, Ⅱ, Ⅲ and Ⅳ) according to the phylogenetic relationships. PgSWEETs of the same clade share similar gene structures, predicting their similar biological functions. RNA-Seq data suggested that PgSWEET genes have a tissue-specific expression pattern. Foliar application of tripotassium phosphate significantly increased the total soluble sugar content of pomegranate fruits and leaves and significantly affected the expression levels of PgSWEETs. The plant growth hormone regulator assay also significantly affected the PgSWEETs expression both in buds of bisexual and functional male flowers. Among them, we selected PgSWEET17a as a candidate gene that plays a role in fructose transport in leaves. The 798 bp CDS sequence of PgSWEET17a was cloned, which encodes 265 amino acids. The subcellular localization of PgSWEET17a showed that it was localized to the cell membrane, indicating its involvement in sugar transport. Transient expression results showed that tobacco fructose content was significantly increased with the up-regulation of PgSWEET17a, while both sucrose and glucose contents were significantly down-regulated. The integration of the PgSWEET phylogenetic tree, gene structure and RNA-Seq data provide a genome-wide trait and expression pattern. Our findings suggest that tripotassium phosphate and plant exogenous hormone treatments could alter PgSWEET expression patterns. These provide a reference for further functional verification and sugar metabolism pathway regulation of PgSWEETs.

Understanding the role of SWEET genes in fruit development and abiotic stress in pomegranate (Punica granatum L.)

BACKGROUND

The Sugar Will Eventually Be Exported Transporters (SWEET), consisting of the MtN3 and salvia domain, are sugar transporters having an active role in diverse activities in plants such as pollen nutrition, phloem loading, nectar secretion, reproductive tissue development, and plant-pathogen interaction. The SWEET genes have been characterized only in a few fruit crop species.

METHODS AND RESULTS

In this study, a total of 15 SWEET genes were identified in the pomegranate (Punica granatum) genome. The gene structure, transmembrane (TM) helices, domain architecture, and phylogenetic relationships of these genes were evaluated using computational approaches. Genes were further classified as Semi-SWEETs or SWEETs based on the TM domains. Similarly, pomegranate, Arabidopsis, rice, and soybean SWEETs were studied together to classify into major groups. In addition, analysis of RNAseq transcriptome data was performed to study SWEEET gene expression dynamics in different tissue. The expression suggests that SWEETs are mostly expressed in pomegranate peel. In addition, PgSWEET13 was found to be differentially expressed under high salinity stress in pomegranate. Further, quantitative PCR analysis confirmed the expression of four candidate genes in leaf and stem tissues.

CONCLUSION

The information provided here will help to understand the role of SWEET genes in fruit development and under abiotic stress conditions in pomegranate.

Profiling Phenolic Composition in Pomegranate Peel From Nine Selected Cultivars Using UHPLC-QTOF-MS and UPLC-QQQ-MS

Pomegranate is widely cultivated across China, and the phenolics in its peel are principal components associated with health benefits. Ultra-high performance liquid chromatography coupled to a quadrupole time-of-flight mass spectrometer (UHPLC-QTOF-MS) and ultra-performance liquid chromatography coupled to a triple quadrupole mass spectrometer (UPLC-QQQ-MS) were used in this study, aiming at profiling the total phenolic composition in pomegranate peel from nine selected cultivars in 7 production areas. Sixty-four phenolic compounds were identified or annotated, and 23 of them were firstly reported in pomegranate peel. Principal component analysis (PCA) plots show differences and similarities of phenolics among nine cultivars. Furthermore, 15 phenolic compounds were quantified with the standards, and punicalagin, ellagic acid, gallocatechin, punicalin, catechin, and corilagin were found to be dominant. Punicalagin weighed the highest content (28.03–104.14 mg/g). This study can provide a deeper and more detailed insight into the phenolic composition in pomegranate peel and facilitate the health-promoting utilization of phenolics.

Identification and expression analysis of the sucrose synthase gene family in pomegranate (Punica granatum L.)

BACKGROUND

Sucrose synthase (SUS, EC 2.4.1.13) is one of the major enzymes of sucrose metabolism in higher plants. It has been associated with C allocation, biomass accumulation, and sink strength. The SUS gene families have been broadly explored and characterized in a number of plants. The pomegranate (Punica granatum) genome is known, however, it lacks a comprehensive study on its SUS genes family.

METHODS

PgSUS genes were identified from the pomegranate genome using a genome-wide search method. The PgSUS gene family was comprehensively analyzed by physicochemical properties, evolutionary relationship, gene structure, conserved motifs and domains, protein structure, syntenic relationships, and cis-acting elements using bioinformatics methods. The expression pattern of the PgSUS gene in different organs and fruit development stages were assayed with RNA-seq obtained from the NCBI SRA database as well as real-time quantitative polymerase chain reaction (qPCR).

RESULTS

Five pomegranate SUS genes, located on four different chromosomes, were divided into three subgroupsaccording to the classification of other seven species. The PgSUS family was found to be highly conserved during evolution after studying the gene structure, motifs, and domain analysis. Furthermore, the predicted PgSUS proteins showed similar secondary and tertiary structures. Syntenic analysis demonstrated that four PgSUS genes showed syntenic relationships with four species, with the exception of PgSUS2. Predictive promoter analysis indicated that PgSUS genes may be responsive to light, hormone signaling, and stress stimulation. RNA-seq analysis revealed that PgSUS1/3/4 were highly expressed in sink organs, including the root, flower, and fruit, and particularly in the outer seed coats. qPCR analysis showed also that PgSUS1, PgSUS3, and PgSUS4 were remarkably expressed during fruit seed coat development. Our results provide a systematic overview of the PgSUS gene family in pomegranate, developing the framework for further research and use of functional PgSUS genes.

Synchronized flowering in pomegranate, following pruning, is associated with expression of the FLOWERING LOCUS T homolog, PgFT1

Flowering in angiosperms is a crucial event that marks the transition from the vegetative to the reproductive phase. In many perennials, pruning is an important horticultural practice that induces synchronized and profuse flowering. In pomegranate, vegetative growth immediately after pruning is associated with activation of PgCENa, a flowering suppressor of the phosphatidyl ethanolamine binding protein (PEBP) family, while a reduction is associated with synchronous flowering. We show that flowering in pomegranate is activated by expression of another PEBP family member, PgFT1, a homolog of the FLOWERING LOCUS T (FT) gene that promotes flowering. PgFT1 shows a rapid reduction in expression during the extensive vegetative growth immediately after pruning but shows robust expression during synchronous flowering post-pruning, in flower-bearing shoots but not in branches that do not bear flowers. A continuous low-level flowering in the absence of pruning is associated with continuous but reduced expression of PgFT1. Flowering by heterologous expression of PgFT1 in Arabidopsis is affected by a single amino acid change in the C-terminal region of PgFT1, which upon correction, promotes flowering in Arabidopsis. Our study provides insights into the molecular mechanisms by which pruning affects flowering pathways in tropical perennial fruit plants such as pomegranate.

The Euscaphis japonica genome and the evolution of malvids

Malvids is one of the largest clades of rosids, includes 58 families and exhibits remarkable morphological and ecological diversity. Here, we report a high-quality chromosome-level genome assembly for Euscaphis japonica, an early-diverging species within malvids. Genome-based phylogenetic analysis suggests that the unstable phylogenetic position of E. japonica may result from incomplete lineage sorting and hybridization event during the diversification of the ancestral population of malvids. Euscaphis japonica experienced two polyploidization events: the ancient whole genome triplication event shared with most eudicots (commonly known as the γ event) and a more recent whole genome duplication event, unique to E. japonica. By resequencing 101 samples from 11 populations, we speculate that the temperature has led to the differentiation of the evergreen and deciduous of E. japonica and the completely different population histories of these two groups. In total, 1012 candidate positively selected genes in the evergreen were detected, some of which are involved in flower and fruit development. We found that reddening and dehiscence of the E. japonica pericarp and long fruit-hanging time promoted the reproduction of E. japonica populations, and revealed the expression patterns of genes related to fruit reddening, dehiscence and abscission. The key genes involved in pentacyclic triterpene synthesis in E. japonica were identified, and different expression patterns of these genes may contribute to pentacyclic triterpene diversification. Our work sheds light on the evolution of E. japonica and malvids, particularly on the diversification of E. japonica and the genetic basis for their fruit dehiscence and abscission.

Rambutan genome revealed gene networks for spine formation and aril development

Rambutan is a popular tropical fruit known for its exotic appearance, has long flexible spines on shells, extraordinary aril growth, desirable nutrition, and a favorable taste. The genome of an elite rambutan cultivar Baoyan 7 was assembled into 328 Mb in 16 pseudo-chromosomes. Comparative genomics analysis between rambutan and lychee revealed that rambutan chromosomes 8 and 12 are collinear with lychee chromosome 1, which resulted in a chromosome fission event in rambutan (n = 16) or a fusion event in lychee (n = 15) after their divergence from a common ancestor 15.7 million years ago. Root development genes played a crucial role in spine development, such as endoplasmic reticulum pathway genes, jasmonic acid response genes, vascular bundle development genes, and K⁺ transport genes. Aril development was regulated by D-class genes (STK and SHP1), plant hormone and phenylpropanoid biosynthesis genes, and sugar metabolism genes. The lower rate of male sterility of hermaphroditic flowers appears to be regulated by MYB24. Population genomic analyses revealed genes in selective sweeps during domestication that are related to fruit morphology and environment stress response. These findings enhance our understanding of spine and aril development and provide genomic resources for rambutan improvement.

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.

Antioxidant and Antibacterial Profiling of Pomegranate-pericarp Extract Functionalized-zinc Oxide Nanocomposite

With the advancement in green nanotechnology, considerable attention is being given to the synthesis of different kinds of nanomaterials for biological applications. In this study, zinc oxide nanocomposites (ZnO NPs) were synthesized using Punica granatum L. (Pomegranate) pericarp ethanolic extract (PE) by the chemical precipitation method. The prepared ZnO NPs showed a characteristic peak at 270 nm in the UV-Vis spectrophotometer and chemical bond stretching in the Fourier transforms infrared spectroscopy (FT-IR) spectra, indicated the formation of PE-functionalized zinc oxide nanocomposite (PE-ZnO NPs). The SEM results showed agglomerated PE-ZnO NPs of a spherical shape with an average size of 80-100 nm. Moreover, biological assessment of the PE-ZnO NPs revealed significant scavenging activity in DPPH (116.5%) and ABTS^·+ (95.2%) radical assay methods, and substantial antibacterial activity against Bacillus cereus, Bacillus licheniformis, and Escherichia coli. Furthermore, PE-ZnO NPs showed about 96.3% of cell viability for human HaCaT cells at the maximum concentration (100 µg/mL), marked as a reliable bioactive agent. Therefore, the developed PE-ZnO NPs were elucidated with substantial ROS scavenger and non-antibiotic antibacterial agent and hence, can be applied in respective biological applications.

Methyl jasmonate elicits distinctive hydrolyzable tannin, flavonoid, and phyto-oxylipin responses in pomegranate (Punica granatum L.) leaves

Transcriptome and biochemical analyses suggested that, while suppression of multiple flavonoids and anthocyanins occurs at least partially at the transcriptional level, increased biosynthesis of non-jasmonate phyto-oxylipins is likely controlled non-transcriptionally. Methyl jasmonate (MeJA) produced in plants can mediate their response to environmental stresses. Exogenous application of MeJA has also shown to activate signaling pathways and induce phytoalexin accumulation in many plant species. To understand how pomegranate plants respond biochemically to environmental stresses, metabolite analysis was conducted in pomegranate leaves subjected to MeJA application and revealed unique changes in hydrolyzable tannins, flavonoids, and phyto-oxylipins. Additionally, transcriptome and real-time qPCR analyses of mock- and MeJA-treated pomegranate leaves identified differentially expressed metabolic genes and transcription factors that are potentially involved in the control of hydrolyzable tannin, flavonoid, and phyto-oxylipin pathways. Molecular, biochemical, and bioinformatic characterization of the only lipoxygenase with sustained, MeJA-induced expression showed that it is capable of oxidizing polyunsaturated fatty acids, though not located in the subcellular compartment where non-jasmonate (non-JA) phyto-oxylipins were produced. These results collectively suggested that while the broad suppression of flavonoids and anthocyanins is at least partially controlled at the transcriptional level, the induced biosynthesis of non-JA phyto-oxylipins is likely not regulated transcriptionally. Overall, a better understanding of how pomegranate leaves respond to environmental stresses will not only promote plant health and productivity, but also have an impact on human health as fruits produced by pomegranate plants are a rich source of nutritional compounds.

Transcriptome-based analysis reveals that the biosynthesis of anthocyanins is more active than that of flavonols and proanthocyanins in the colorful flowers of Lagerstroemia indica

Mechanisms associated with the control of flower color in crape myrtle varieties have yet to be sufficiently elucidated, which has tended to hamper the use of modern molecular and genetic strategies in the breeding programs for this plant. The whole transcriptome of four L. indica varieties characterized by different flower colors (white, light purple, deep purplish pink, and strong red) was sequenced, and we performed bioinformatic, quantitative PCR, and co-expression analyses of R2R3 MYB transcription factor and anthocyanin/flavonol pathway genes. We obtained a total of 49,980 transcripts with full-length coding sequences. Both transcriptome and qPCR analyses revealed that anthocyanin/flavonol pathway genes were differentially expressed among the four different flowers types, with the expression of LiPAL, LiCHS, LiCHI, LiDFR, LiANS/LDOX, and LiUFGT being induced in colorful flowers, whereas that of LiF3´5´H, LiFLS, and LiLAR was found to be inhibited. Base on phylogenetic analysis, seven R2R3 MYB transcriptional factors were identified as putative regulators of flower color. The molecular characteristics and co-expression patterns indicated that these MYBs differentially modulate their target genes, with two probably acting as activators, three as repressors, and one contributing to the regulation of vacuolar pH. The findings of this study indicate that the anthocyanin biosynthesis is more active than the flavonol and proanthocyanin in the colorful flowers. These observations provide new genomic information on L. indica and contribute gene resources for the flower color-targeted breeding of crape myrtle.

Chloroplast Genome Sequencing, Comparative Analysis, and Discovery of Unique Cytoplasmic Variants in Pomegranate (Punica granatum L.)

Here we report on comprehensive chloroplast (cp) genome analysis of 16 pomegranate (Punica granatum L.) genotypes representing commercial cultivars, ornamental and wild types, through large-scale sequencing and assembling using next-generation sequencing (NGS) technology. Comparative genome analysis revealed that the size of cp genomes varied from 158,593 bp (in wild, “1201” and “1181”) to 158,662 bp (cultivar, “Gul-e-Shah Red”) among the genotypes, with characteristic quadripartite structures separated by a pair of inverted repeats (IRs). The higher conservation for the total number of coding and non-coding genes (rRNA and tRNA) and their sizes, and IRs (IR-A and IR-B) were observed across all the cp genomes. Interestingly, high variations were observed in sizes of large single copy (LSC, 88,976 to 89,044 bp) and small single copy (SSC, 18,682 to 18,684 bp) regions. Although, the structural organization of newly assembled cp genomes were comparable to that of previously reported cp genomes of pomegranate (“Helow,” “Tunisia,” and “Bhagawa”), the striking differences were observed with the Lagerstroemia lines, viz., Lagerstroemia intermedia (NC_0346620) and Lagerstroemia speciosa (NC_031414), which clearly confirmed previous findings. Furthermore, phylogenetic analysis also revealed that members outside the genus Punica were clubbed into a separate clade. The contraction and expansion analysis revealed that the structural variations in IRs, LSC, and SSC have significantly accounted for the evolution of cp genomes of Punica and L. intermedia over the periods. Microsatellite survey across cp genomes resulted in the identification of a total of 233 to 234 SSRs, with majority of them being mono- (A/T or C/G, 164–165 numbers), followed by di- (AT/AT or AG/CT, 54), tri- (6), tetra- (8), and pentanucleotides (1). Furthermore, the comparative structural variant analyses across cp genomes resulted in the identification of many varietal specific SNP/indel markers. In summary, our study has offered a successful development of large-scale cp genomics resources to leverage future genetic, taxonomical, and phylogenetic studies in pomegranate.

TCM-Blast for traditional Chinese medicine genome alignment with integrated resources

The traditional Chinese medicine (TCM) genome project aims to reveal the genetic information and regulatory network of herbal medicines, and to clarify their molecular mechanisms in the prevention and treatment of human diseases. Moreover, the TCM genome could provide the basis for the discovery of the functional genes of active ingredients in TCM, and for the breeding and improvement of TCM. The traditional Chinese Medicine Basic Local Alignment Search Tool (TCM-Blast) is a web interface for TCM protein and DNA sequence similarity searches. It contains approximately 40G of genome data on TCMs, including protein and DNA sequence for 36 TCMs with high medical value.The development of a publicly accessible TCM genome alignment database hosted on the TCM-Blast website ( http://viroblast.pungentdb.org.cn/TCM-Blast/viroblast.php ) has expanded to query multiple sequence databases to obtain TCM genome data, and provide user-friendly output for easy analysis and browsing of BLAST results. The genome sequencing of TCMs helps to elucidate the biosynthetic pathways of important secondary metabolites and provides an essential resource for gene discovery studies and molecular breeding. The TCMs genome provides a valuable resource for the investigation of novel bioactive compounds and drugs from these TCMs under the guidance of TCM clinical practice. Our database could be expanded to other TCMs after the determination of their genome data.

Genome-wide identification of candidate aquaporins involved in water accumulation of pomegranate outer seed coat

Aquaporins (AQPs) are a class of highly conserved integral membrane proteins that facilitate the uptake and transport of water and other small molecules across cell membranes. However, little is known about AQP genes in pomegranate (Punica granatum L.) and their potential role in water accumulation of the outer seed coat. We identified 38 PgrAQP genes in the pomegranate genome and divided them into five subfamilies based on a comparative analysis. Purifying selection played a role in the evolution of PgrAQP genes and a whole-genome duplication event in Myrtales may have contributed to the expansion of PgrTIP, PgrSIP, and PgrXIP genes. Transcriptome data analysis revealed that the PgrAQP genes exhibited different tissue-specific expression patterns. Among them, the transcript abundance of PgrPIPs were significantly higher than that of other subfamilies. The mRNA transcription levels of PgrPIP1.3, PgrPIP2.8, and PgrSIP1.2 showed a significant linear relationship with water accumulation in seed coats, indicating that PgrPIP1.3/PgrPIP2.8 located in the plasma membrane and PgrSIP1.2 proteins located on the tonoplast may be involved in water accumulation and contribute to the cell expansion of the outer seed coat, which then develops into juicy edible flesh. Overall, our results provided not only information on the characteristics and evolution of PgrAQPs, but also insights on the genetic improvement of outer seed coats.

Physiological and Immune Functions of Punicalagin

The aim of this publication is to compile a summary of the findings regarding punicalagin in various tissues described thus far in the literature, with an emphasis on the effect of this substance on immune reactions. Punicalagin (PUN) is an ellagitannin found in the peel of pomegranate (Punica granatum). It is a polyphenol with proven antioxidant, hepatoprotective, anti-atherosclerotic and chemopreventive activities, antiproliferative activity against tumor cells; it inhibits inflammatory pathways and the action of toxic substances, and is highly tolerated. This work describes the source, metabolism, functions and effects of punicalagin, its derivatives and metabolites. Furthermore, its anti-inflammatory and antioxidant effects are described.

Identification of suitable reference genes for expression studies in pomegranate under different biotic and abiotic stress conditions

Pomegranate (Punica granatum L.) is an important economic fruit crop, facing many biotic and abiotic challenges during cultivation. Several research programs are in progress to understand both biotic and abiotic stress factors and mitigate these challenges using gene expression studies based on the qPCR approach. However, research publications are not available yet to select the standard reference gene for normalizing target gene expression values in pomegranate. The most suitable candidate reference gene is required to ensure precise and reliable results for qPCR analysis. Eight candidate reference genes' stability was evaluated under different stress conditions using different algorithms such as ∆Ct, geNorm, BestKeeper, NormFinder, and RefFinder. The various algorithms revealed that EFA1 and 18S rRNA were common and most stable reference genes (RGs) under abiotic and wilt stress. Whereas comprehensive ranking by RefFinder showed GAPDH and CYPF were the most stable RGs under combined biotic (pooled samples of all biotic stress) and bacterial blight samples. For normalizing target gene expression under wilt, nematode, bacterial blight, and abiotic stress conditions both GAPDH and CYPFreference genes are adequate for qPCR. The above data provide comprehensive details for the selection of a candidate reference gene in various stresses in pomegranate.

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

Guava (Psidium guajava) is an important fleshy-fruited tree of the Myrtaceae family that is widely cultivated in tropical and subtropical areas of the world and has attracted considerable attention for the richness of ascorbic acid in its fruits. However, studies on the evolution and genetic breeding potential of guava are hindered by the lack of a reference genome. Here, we present a chromosome-level genomic assembly of guava using PacBio sequencing and Hi-C technology. We found that the genome assembly size was 443.8 Mb with a contig N50 of ~15.8 Mb. We annotated a total of 25 601 genes and 193.2 Mb of repetitive sequences for this genome. Comparative genomic analysis revealed that guava has undergone a recent whole-genome duplication (WGD) event shared by all species in Myrtaceae. In addition, through metabolic analysis, we determined that the L-galactose pathway plays a major role in ascorbic acid biosynthesis in guava fruits. Moreover, the softening of fruits of guava may result from both starch and cell wall degradation according to analyses of gene expression profiles and positively selected genes. Our data provide a foundational resource to support molecular breeding of guava and represent new insights into the evolution of soft, fleshy fruits in Myrtaceae.

Fine Mapping of the "black" Peel Color in Pomegranate (Punica granatum L.) Strongly Suggests That a Mutation in the Anthocyanidin Reductase (ANR) Gene Is Responsible for the Trait

Anthocyanins are important dietary and health-promoting substances present in high quantities in the peel and arils of the pomegranate (Punica granatum L.) fruit. Yet, there is a high variation in the content of anthocyanin among different pomegranate varieties. The 'Black' pomegranate variety (P.G.127-28) found in Israel contains exceptionally high levels of anthocyanins in its fruit peel which can reach up to two orders of magnitude higher content as compared to that of other pomegranate varieties' peel anthocyanins. Biochemical analysis reveals that delphinidin is highly abundant in the peel of 'Black' variety. The pattern of anthocyanin accumulation in the fruit peel during fruit development of 'Black' variety differs from that of other pomegranates. High anthocyanin levels are maintained during all developmental stages. Moreover, the accumulation of anthocyanin in the fruit peel of 'Black' variety is not dependent on light. Genetic analysis of an F₂ population segregating for the "black" phenotype reveals that it is determined by a single recessive gene. Genetic mapping of the F₂ population using single nucleotide polymorphism (SNP) markers identified few markers tightly linked to the "black" phenotype. Recombination analysis of the F₂ population and F₃ populations narrowed the "black" trait to an area of 178.5 kb on the draft genome sequence of pomegranate cv. 'Dabenzi.' A putative anthocyanidin reductase (ANR) gene is located in this area. Only pomegranate varieties displaying the "black" trait carry a base pair deletion toward the end of the gene, causing a frame shift resulting in a shorter protein. We propose that this mutation in the ANR gene is responsible for the different anthocyanin composition and high anthocyanin levels of the "black" trait in pomegranate.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.

Functioning of PPR Proteins in Organelle RNA Metabolism and Chloroplast Biogenesis

The pentatricopeptide repeat (PPR) proteins constitute one of the largest nuclear-encoded protein families in higher plants, with over 400 members in most sequenced plant species. The molecular functions of these proteins and their physiological roles during plant growth and development have been widely studied. Generally, there is mounting evidence that PPR proteins are involved in the post-transcriptional regulation of chloroplast and/or mitochondrial genes, including RNA maturation, editing, intron splicing, transcripts' stabilization, and translation initiation. The cooperative action of RNA metabolism has profound effects on the biogenesis and functioning of both chloroplasts and mitochondria and, consequently, on the photosynthesis, respiration, and development of plants and their environmental responses. In this review, we summarize the latest research on PPR proteins, specifically how they might function in the chloroplast, by documenting their mechanism of molecular function, their corresponding RNA targets, and their specific effects upon chloroplast biogenesis and host organisms.