The pomegranate (Punica granatum L.) draft genome dissects genetic divergence between soft- and hard-seeded cultivars

A telomere-to-telomere gap-free assembly integrating multi-omics uncovers the genetic mechanism of fruit quality and important agronomic trait associations in pomegranate

Pomegranate is an important perennial fruit tree distributed worldwide. Reference genomes with gaps and limit gene identification controlling important agronomic traits hinder its functional genomics and genetic improvements. Here, we reported a telomere-to-telomere (T2T) gap-free genome assembly of the distinctive cultivar 'Moshiliu'. The Moshiliu reference genome was assembled into eight chromosomes without gaps, totalling ~366.71 Mb, with 32 158 predicted protein-coding genes. All 16 telomeres and eight centromeres were characterized; combined with FISH analysis, we revealed the atypical telomere units in pomegranate as TTTTAGGG. Furthermore, a total of 16 loci associated with 15 important agronomic traits were identified based on GWAS of 146 accessions. Gene editing and biochemical experiments demonstrated that a 37.2-Kb unique chromosome translocation disrupting the coding domain sequence of PgANS was responsible for anthocyanin-less, knockout of PgANS in pomegranate exhibited a defect in anthocyanin production; a unique repeat expansion in the promoter of PgANR may affected its expression, resulting in black peel; notably, the G → A transversion located at the 166-bp coding domain of PgNST3, which caused a E56K mutation in the PgNST3 protein, closely linked with soft-seed trait. Overexpression of PgNST3A in tomato presented smaller and softer seed coats. The E56K mutation in PgNST3 protein, eliminated the binding ability of PgNST3 to the PgMYB46 promoter, which subsequently affected the thickness of the inner seed coat of soft-seeded pomegranates. Collectively, the validated gap-free genome, the identified genes controlling important traits and the CRISPR-Cas9-mediated gene knockout system all provided invaluable resources for pomegranate precise breeding.

Genome assembly of pomegranate highlights structural variations driving population differentiation and key loci underpinning cold adaption

Cold damage poses a significant challenge to the cultivation of soft-seeded pomegranate varieties, hindering the growth of the pomegranate industry. The genetic basis of cold tolerance in pomegranates has remained elusive, largely due to the lack of high-quality genome assemblies for cold-tolerant varieties and comprehensive population-scale genomic studies. In this study, we addressed these challenges by assembling a high-quality chromosome-level reference genome for 'Sanbai', a pomegranate variety renowned for its freezing resistance, achieving an impressive contig N50 of 15.93 Mb. This robust assembly, enhanced by long-read sequencing of 38 pomegranate accessions, facilitated the identification of 14 239 polymorphic structural variants, revealing their critical roles in genomic diversity and population differentiation related to cold tolerance. Of particular significance was the discovery of a ~ 5.4-Mb inversion on chromosome 1, which emerged as an important factor affecting cold tolerance in pomegranate. Moreover, through the integration of bulked segregant analysis, differential selection analysis, and genetic transformation techniques, we identified and validated the interaction between the PgNAC12 transcription factor and PgCBF1, disclosing their pivotal roles in response to cold stress. These findings mark a significant advancement in pomegranate genomics, offering novel insights into the genetic mechanisms of cold tolerance and providing valuable resources for the genetic improvement of soft-seeded pomegranate varieties.

Genome assembly, annotation and evolutionary insights from the draft genome of wild pomegranate

Wild pomegranate is a potent medicinal plant known for its medicinal and nutritional attributes. Despite its healing and curative properties, the genome of this wild species remains elusive, thus limiting our understanding on the genetic processes involved in the biosynthesis of functional molecules. This study presents the annotation of a de novo genome assembly of wild pomegranate, with a genome size of 279.0 Mb. From the assembly, 34.8 GB of the data was retained, encompassing 72,055 scaffolds. A total of 49,178 genes were predicted, with an average of 5.36 exons per gene and a GC content of 49%. About 14,400 genes were annotated in biological, cellular and molecular processes related mostly to carbohydrate metabolism, intracellular signal transduction, mRNA binding and DNA helicase activity. KEGG enrichment analysis revealed maximum number of genes associated with biosynthesis of secondary metabolites mainly phenypropanoid pathway, followed by ribosome and plant hormone signal transduction. From the identified functional genes, 230 genes scaffolds encoded for transcription factors belonging to 25 families with highest recorded for MYB gene family. Study of annotated transposable elements unveiled the existence of long terminal repeats and retrotransposons. Additionally, our investigation involves the comparative analysis and identification of orthologous genes among the genomes of wild and cultivated species of Punica granatum and also across selected five plant species Eucalyptus grandis, Vitis vinifera, Jatropha curcas, Theobroma cacao and Gossypium raimondii, revealing the functional and evolutionary dynamics across species. To the best of our knowledge, this is the first report on the genome assembly, annotation and gene prediction in wild pomegranate. Also, information regarding the terpenoid pathway genes has been unravelled for the first time in the present study. Inclusively, the current study offers thorough details on important aspects of the wild pomegranate genome that would be useful in comprehending its genetics and will facilitate discovery of genes against various biotic and abiotic stresses.

Integrated genomic and transcriptomic analyses reveal the genetic and molecular mechanisms underlying hawthorn peel color and seed hardness diversity

Hawthorn (Crataegus pinnatifida) fruit peel color and seed hardness are key traits that significantly impact economic value. We present here the high-quality chromosome-scale genomes of two cultivars, including the hard-seed, yellow-peel C. pinnatifida "Jinruyi" (JRY) and the soft-seed, red-peel C. pinnatifida "Ruanzi" (RZ). The assembled genomes comprising 17 chromosomes are 809.1 Mb and 760.5 Mb in size, achieving scaffold N50 values of 48.5 Mb and 46.8 Mb for JRY and RZ, respectively. Comparative genomic analysis identifies 3.6-3.8 million single nucleotide polymorphisms, 8.5-9.3 million insertions/deletions, and approximately 30 Mb of presence/absence variations across different hawthorn genomes. Through integrating differentially expressed genes and accumulated metabolites, we filter candidate genes CpMYB114 and CpMYB44 associated with differences in hawthorn fruit peel color and seed hardness, respectively. Functional validation confirms that the CpMYB114-CpANS regulates anthocyanin biosynthesis in hawthorn peels, contributing to the observed variation in peel color. CpMYB44-CpCOMT is significantly upregulated in JRY and is verified to promote lignin biosynthesis, resulting in the distinction in seed hardness. Overall, this study reveals the new insights into understanding of distinct peel pigmentation and seed hardness in hawthorn and provides an abundant resource for molecular breeding.

Genome assembly and multiomic analyses reveal insights into flower and bark colors of Lagerstroemia excelsa

Lagerstroemia excelsa is a unique plant species from China, holds a significant aesthetic and economic value, and plays a crucial role in landscape architecture and horticulture. Thus far, there is little genetic and genomic information available about this species, which limits its use in development of new cultivars. In this study, a high-quality genome map of L. excelsa was obtained via whole-genome sequencing. Results showed that its genome size is about 330.4 Mb and a scaffold mapping rate is approximately 97.20%, resulting in 24 pseudochromosomes. L. excelsa might have undergone a recent whole-genome triplication event and diverged from the pomegranate about 32.3 million years ago (MYA). Subsequently, the divergence time between L. indica and L. excelsa was around 5.9 MYA. The transcriptomic and metabolomic analyses of L. excelsa and L. indica indicated that the chalcone synthase pathway may play a key role in regulating flower color differentiation between the two species. Additionally, a transcription factor LeMYB103 may be involved in regulating anthocyanin synthesis by interacting with LeMYB66, resulting in the accumulation of anthocyanins in the stem bark. This study is the first step toward genomic analysis of L. excelsa, which may provide a foundation for further molecular investigation of this species and offer valuable insights into the molecular mechanisms underlying the flower and stem bark colors in L. excelsa, two important ornamental traits in Lagerstroemia breeding.

Pomegranate ATP-binding cassette transporter PgABCG9 plays a negative regulatory role in lignin accumulation

Seed hardness is an important quality characteristic of pomegranate fruit. The development of seed hardness relies on the deposition of lignin in the inner seed coat, but the underlying molecular mechanisms remain unclear. In this study, we identified a member of ABCG transporters, PgABCG9, which may function in seed hardening by negatively regulating lignin biosynthesis. PgABCG9 was expressed at high levels in the inner seed coats of pomegranate fruit, and its transcript level was negatively correlated with seed hardness. PgABCG9-transgenic Arabidopsis plants exhibited weaker growth and thinner stems than the wild-type. The number of xylem cells, xylem cell wall thickness, and lignin deposition in the PgABCG9 transgenic plants were significantly reduced. In addition, overexpression of PgABCG9 in Arabidopsis enhanced plant tolerance to exogenous monolignols. Targeted metabolite profiling revealed that the contents of metabolites involved in lignin biosynthesis, including monolignols and monolignol precursors, were also reduced in PgABCG9- transgenic plants. We found that PgABCG9 is localized to the Golgi. These findings indicate that PgABCG9 plays a negative regulatory role in lignin biosynthesis and potentially contributes to soft-seed development in pomegranate through a mechanism that includes the reduction of lignin content in the seed coat by sequestration of monolignols in intracellular compartments.

A Chromosome level assembly of pomegranate (Punica granatum L.) variety grown in arid environment

The pomegranate (Punica granatum L.) is an ancient fruit-bearing tree known for its nutritional and antioxidant properties. They originated from the Middle East in regions having large farms including mountainous regions of Al-Baha in Saudi Arabia. Pomegranates can tolerate arid climates and are considered among the fruits that will play a major role in food security. However, the genomics resources of pomegranate growing in arid regions are scarce. Here, we present a high-quality chromosome-level reference genome using PacBio HiFi long reads. The final assembly was 384.65 Mb with N50 contig size of 43.11 Mb, with 353.42 Mb being anchored on the eight pseudo chromosomes. Annotation revealed that 48.79% of the genome comprises repetitive elements and contains 21,620 protein-coding genes. The new reference genome will contribute to identifying stress resistance traits in pomegranates thriving in arid environments as well as new dietary antioxidants and antimicrobial peptides with pharmaceutical and therapeutic applications.

Herbal Multiomics Provide Insights into Gene Discovery and Bioproduction of Triterpenoids by Engineered Microbes

Triterpenoids are natural products found in plants that exhibit industrial and agricultural importance. Triterpenoids are typically synthesized through two main pathways: the mevalonate (MVA) and methylerythritol 4-phosphate (MEP) pathways. They then undergo structural diversification with the help of squalene cyclases (OSCs), cytochrome P450 monooxygenases (P450s), UDP glycosyltransferases (UGTs), and acyltransferases (ATs). Advances in multiomics technologies for herbal plants have led to the identification of novel triterpenoid biosynthetic pathways. The application of various analytical techniques facilitates the qualitative and quantitative analysis of triterpenoids. Progress in synthetic biology and metabolic engineering has also facilitated the heterologous production of triterpenoids in microorganisms, such as Escherichia coli and Saccharomyces cerevisiae. This review summarizes recent advances in biotechnological approaches aimed at elucidating the complex pathway of triterpenoid biosynthesis. It also discusses the metabolic engineering strategies employed to increase the level of triterpenoid production in chassis cells.

Engineering Saccharomyces cerevisiae for the Production of Punicic Acid-Rich Yeast Biomass

Punicic acid (PuA), an unusual conjugated linolenic acid found in pomegranate, offers diverse health benefits and has potential applications in the food industry. Due to the limited availability of PuA from natural plant sources, there is growing interest in producing it through microbial fermentation. In this study, the yeast Saccharomyces cerevisiae, which is classified as "generally recognized as safe", was engineered to produce PuA using a results-driven approach. Genes potentially involved in PuA synthesis were integrated directly into the yeast genome, targeting Ty retrotransposon sites. Screening of the yeast transformants, followed by optimization of culture conditions, resulted in the production of 26.7% PuA within the yeast's total fatty acids. Further analysis revealed that the strain's triacylglycerol fraction contained over 22% PuA. By incorporating this health-promoting lipid into the nutritional profile of S. cerevisiae, the engineered strain could serve as a sustainable source of yeast biomass with enhanced nutritional value.

Transcriptomic Analysis Reveals the Mechanism of Color Formation in the Peel of an Evergreen Pomegranate Cultivar 'Danruo No.1' During Fruit Development

Pomegranate (Punica granatum L.) is an ancient fruit crop that has been cultivated worldwide and is known for its attractive appearance and functional metabolites. Fruit color is an important index of fruit quality, but the color formation pattern in the peel of evergreen pomegranate and the relevant molecular mechanism is still unknown. In this study, the contents of pigments including anthocyanins, carotenoids, and chlorophyll in the peel of 'Danruo No. 1' pomegranate fruit during three developmental stages were measured, and RNA-seq was conducted to screen key genes regulating fruit color formation. The results show that pomegranate fruit turned from green to red during development, with a dramatic increase in a* value, indicating redness and anthocyanins concentration, and a decrease of chlorophyll content. Moreover, carotenoids exhibited a decrease-increase accumulation pattern. Through RNA-seq, totals of 30, 18, and 17 structural genes related to anthocyanin biosynthesis, carotenoid biosynthesis and chlorophyll metabolism were identified from differentially expressed genes (DEGs), respectively. Transcription factors (TFs) such as MYB, bHLH, WRKY and AP2/ERF were identified as key candidates regulating pigment metabolism by K-means analysis and weighted gene co-expression network analysis (WGCNA). The results provide an insight into the theory of peel color formation in evergreen pomegranate fruit.

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.

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.

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.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

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.

High-quality maple genome reveals duplication-facilitated leaf color diversity

Acer truncatum is a horticultural tree species with individuals that display either yellow or red leaves in autumn, giving it high ornamental and economic value. 'Lihong' of A. truncatum is an excellent cultivar due to its characteristic of having autumn leaves that turn a bright and beautiful shade of red, while its closely related cultivar 'Bunge' does not. However, the molecular mechanism underlying the color change in the cultivar 'Lihong' is still unclear. Here, we assembled a high-quality genome sequence of Acer truncatum 'Lihong' (genome size = 688 Mb, scaffold N50 = 9.14 Mb) with 28,438 protein-coding genes predicted. Through comparative genomic analysis, we found that 'Lihong' had experienced more tandem duplication events although it's a high degree of collinearity with 'Bunge'. Especially, the expansion of key enzymes in the anthocyanin synthesis pathway was significantly uneven between the two varieties, with 'Lihong' genome containing a significantly higher number of tandem/dispersed duplication key genes. Further transcriptomic, metabolomic, and molecular functional analyses demonstrated that several UFGT genes, mainly resulting from tandem/dispersed duplication, followed by the promoter sequence variation, may contribute greatly to the leaf color phenotype, which provides new insights into the mechanism of divergent anthocyanin accumulation process in the 'Lihong' and 'Bunge' with yellow leaves in autumn. Further, constitutive expression of two UFGT genes, which showed higher expression in 'Lihong', elevated the anthocyanin content. We proposed that the small-scale duplication events could contribute to phenotype innovation.

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.

High-quality genome of a pioneer mangrove Laguncularia racemosa explains its advantages for intertidal zone reforestation

Ecological restoration of mangrove ecosystems that became susceptible to recent habitat perturbations is crucial for tropical coast conservation. The white mangrove Laguncularia racemosa, a pioneer species inhabiting intertidal environments of the Atlantic East Pacific (AEP) region, has been used for reforestation in China for decades. However, the molecular mechanisms underlying its fast growth and high adaptive potential remain unknown. Using PacBio single-molecule real-time sequencing, we completed a high-quality L. racemosa genome assembly covering 1105 Mb with scaffold N50 of 3.46 Mb. Genomic phylogeny shows that L. racemosa invaded intertidal zones during a period of global warming. Multi-level genomic convergence analyses between L. racemosa and three native dominant mangrove clades show that they experienced convergent changes in genes involved in nutrient absorption and high salinity tolerance. This may explain successful L. racemosa adaptation to stressful intertidal environments after introduction. Without recent whole-genome duplications or activated transposable elements, L. racemosa has retained many tandem gene duplications. Some of them are involved in auxin biosynthesis, intense light stress and cold stress response pathways, associated with L. racemosa's ability to grow fast under high light or cold conditions when used for reforestation. In summary, our study identifies shared mechanisms of intertidal environmental adaptation and unique genetic changes underlying fast growth in mangrove-unfavourable conditions and sheds light on the molecular mechanisms of the white mangrove utility in ecological restoration.

Genome assembly and resequencing analyses provide new insights into the evolution, domestication and ornamental traits of crape myrtle

Crape myrtle (Lagerstroemia indica) is a globally used ornamental woody plant and is the representative species of Lagerstroemia. However, studies on the evolution and genomic breeding of L. indica have been hindered by the lack of a reference genome. Here we assembled the first high-quality genome of L. indica using PacBio combined with Hi-C scaffolding to anchor the 329.14-Mb genome assembly into 24 pseudochromosomes. We detected a previously undescribed independent whole-genome triplication event occurring 35.5 million years ago in L. indica following its divergence from Punica granatum. After resequencing 73 accessions of Lagerstroemia, the main parents of modern crape myrtle cultivars were found to be L. indica and L. fauriei. During the process of domestication, genetic diversity tended to decrease in many plants, but this was not observed in L. indica. We constructed a high-density genetic linkage map with an average map distance of 0.33 cM. Furthermore, we integrated the results of quantitative trait locus (QTL) using genetic mapping and bulk segregant analysis (BSA), revealing that the major-effect interval controlling internode length (IL) is located on chr1, which contains CDL15, CRG98, and GID1b1 associated with the phytohormone pathways. Analysis of gene expression of the red, purple, and white flower-colour flavonoid pathways revealed that differential expression of multiple genes determined the flower colour of L. indica, with white flowers having the lowest gene expression. In addition, BSA of purple- and green-leaved individuals of populations of L. indica was performed, and the leaf colour loci were mapped to chr12 and chr17. Within these intervals, we identified MYB35, NCED, and KAS1. Our genome assembly provided a foundation for investigating the evolution, population structure, and differentiation of Myrtaceae species and accelerating the molecular breeding of L. indica.

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.

Chromosome-level assemblies of cultivated water chestnut Trapa bicornis and its wild relative Trapa incisa

Water chestnut (Trapa L.) is a floating-leaved aquatic plant with high edible and medicinal value. In this study, we presented chromosome-level genome assemblies of cultivated large-seed species Trapa bicornis and its wild small-seed relative Trapa incisa by using PacBio HiFi long reads and Hi-C technology. The T. bicornis and T. incisa assemblies consisted of 479.90 Mb and 463.97 Mb contigs with N50 values of 13.52 Mb and 13.77 Mb, respectively, and repeat contents of 62.88% and 62.49%, respectively. A total of 33,306 and 33,315 protein-coding genes were predicted in T. bicornis and T. incisa assemblies, respectively. There were 159,232 structural variants affecting more than 11 thousand genes detected between the two genomes. The phylogenetic analysis indicated that the lineage leading to Trapa was diverged from the lineage to Sonneratia approximately 23 million years ago. These two assemblies provide valuable resources for future evolutionary and functional genomic research and molecular breeding of water chestnut.

Genome-wide alternative polyadenylation dynamics underlying plant growth retardant-induced dwarfing of pomegranate

Dwarfed stature is a desired agronomic trait for pomegranate (Punica granatum L.), with its advantages such as lower cost and increased yield. A comprehensive understanding of regulatory mechanisms underlying the growth repression would provide a genetic foundation to molecular-assisted dwarfing cultivation of pomegranate. Our previous study induced dwarfed pomegranate seedlings via exogenous application of plant growth retardants (PGRs) and highlighted the important roles of differential expression of plant growth-related genes in eliciting the dwarfed phenotype of pomegranate. Alternative polyadenylation (APA) is an important post-transcriptional mechanism and has been demonstrated to act as a key regulator in plant growth and development. However, no attention has been paid to the role of APA in PGR-induced dwarfing in pomegranate. In this study, we characterized and compared APA-mediated regulation events underlying PGR-induced treatments and normal growth condition. Genome-wide alterations in the usage of poly(A) sites were elicited by PGR treatments, and these changes were involved in modulating the growth and development of pomegranate seedlings. Importantly, ample specificities were observed in APA dynamics among the different PGR treatments, which mirrors their distinct nature. Despite the asynchrony between APA events and differential gene expression, APA was found to regulate transcriptome via influencing microRNA (miRNA)-mediated mRNA cleavage or translation inhibition. A global preference for lengthening of 3' untranslated regions (3' UTRs) was observed under PGR treatments, which was likely to host more miRNA target sites in 3' UTRs and thus suppress the expression of the corresponding genes, especially those associated with developmental growth, lateral root branching, and maintenance of shoot apical meristem. Together, these results highlighted the key role of APA-mediated regulations in fine-tuning the PGR-induced dwarfed stature of pomegranate, which provides new insights into the genetic basis underlying the growth and development of pomegranate.

Development of Molecular Marker Linked to Seed Hardness in Pomegranate Using Bulked Segregant Analysis

The pomegranate (Punica granatum L.) is one of the fruit species with the oldest cultural history. There are many traits to determine the quality of pomegranate fruits. Among them, soft-seeded feature of pomegranate fruit is important trait for the market value of the fruit. For this reason, the demand for pomegranate varieties with soft seeds has been increasing, especially in recent years. In this study, molecular markers associated with seed hardness were developed to distinguish pomegranate cultivars with soft-seeded feature based on genomic DNA at the early stages of the pomegranate breeding process. For this purpose, pomegranate genotypes and/or cultivars from the population involved in reciprocal crosses of hard-seeded Ernar, medium-hard-seeded Hicaznar, and soft-seeded Fellahyemez cultivars were grouped as soft-seeded or hard-seeded. Further, leaf samples were collected from individuals belonging to each group. Then, the genomic DNA was isolated from each plant separately, and equal amount of genomic DNA from individuals with the similar seed hardness were mixed for bulked segregant analysis (BSA). The bulked genomic DNAs of opposite characters were analyzed by polymerase chain reaction (PCR) using random decamer primers to develop random amplified polymorphic DNA (RAPD) markers associated with soft-seeded or hard-seeded pomegranates. A total of three RAPD markers were determined to distinguish the individuals having soft- or hard-seeded pomegranate genotypes and/or cultivars. As a result of the comparison of the DNA sequences of these RAPD markers, insertion-deletions (inDels) primers were designed to developed and validate a PCR assay to distinguish the soft- and hard-seeded pomegranate genotypes/cultivars from each other. The molecular markers developed in this study will enable us to distinguish soft-seeded pomegranate types easily in a short time at the early stages of the pomegranate breeding programs.

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.

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.

Biotechnological interventions in reducing losses of tropical fruits and vegetables

Tropical fruits and vegetables are predominantly cultivated in warm climate zones, resulting in cultivar diversity in terms of structure, features, and physiology. These constitute a variety of bioactive ingredients such as vitamins, minerals, phenolic acids, anthocyanins, flavonoids, fatty acids, fiber, and their distinctive appearances attract customers across the world. The global production of fruit and vegetables has been attained a tremendous increase for the last few decades. However, huge losses at pre- and postharvest levels are major constraints in their judicious use. Traditional breeding strategies were used to minimize these losses, but their functionality is limited due to their time and labor intensiveness. Recent biotechnological, computational, and multiomics approaches not only address the losses concern but also aid in boosting crop productivity and nutritional values. This article emphasizes molecular tools that have been used to reduce losses of tropical fruits and vegetables at pre- and postharvest levels.

Identification of the egusi seed trait locus (eg) and its suppressor gene associated with the thin seed coat trait in watermelon

Egusi watermelon has a unique egusi seed type, which could be useful for breeding both edible seeds and edible flesh in watermelon. However, the genetic basis of the unique egusi seed type is not clear. In the present study, we first reported that at least two genes with inhibitory epistasis were responsible for the thin seed coat (unique egusi seed type) in watermelon. Inheritance analysis of five populations, including F₂, BC, and BCF₂, suggested that the thin seed coat trait was controlled by a suppressor gene together with the egusi seed locus (eg) in egusi watermelon. Based on high-throughput sequencing technology, two quantitative trait loci located on chromosome 1 and chromosome 6 were identified for the thin seed coat trait in watermelon. One of the loci, the eg locus on chromosome 6, was finely mapped to a genomic region of 15.7 kb, which contained only one candidate gene. Comparative transcriptome analysis highlighted differentially expressed genes involved in cellulose and lignin synthesis between watermelon genotypes varying in the thickness of the seed coat and provided several potential candidate genes for the thin seed coat trait. Taken together, our data suggest that at least two genes are complementarily involved in the thin seed coat trait and will be useful for cloning novel genes. The results presented here provide a new reference for uncovering egusi seed genetic mechanisms and valuable information for marker-assisted selection in seed coat breeding.

Multiple introductions and recombination events underlie the emergence of a hyper-transmissible Cryptosporidium hominis subtype in the USA

The parasite Cryptosporidium hominis is a leading cause of the diarrheal disease cryptosporidiosis, whose incidence in the United States has increased since 2005. Here, we show that the newly emerged and hyper-transmissible subtype IfA12G1R5 is now dominant in the United States. In a comparative analysis of 127 newly sequenced and 95 published C. hominis genomes, IfA12G1R5 isolates from the United States place into three of the 14 clusters (Pop6, Pop13, and Pop14), indicating that this subtype has multiple ancestral origins. Pop6 (IfA12G1R5a) has an East Africa origin and has recombined with autochthonous subtypes after its arrival. Pop13 (IfA12G1R5b) is imported from Europe, where it has recombined with the prevalent local subtype, whereas Pop14 (IfA12G1R5c) is a progeny of secondary recombination between Pop6 and Pop13. Selective sweeps in invasion-associated genes have accompanied the emergence of the dominant Pop14. These observations offer insights into the emergence and evolution of hyper-transmissible pathogens.

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.

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.

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.

Overexpression of PgCBF3 and PgCBF7 Transcription Factors from Pomegranate Enhances Freezing Tolerance in Arabidopsis under the Promoter Activity Positively Regulated by PgICE1

Cold stress limits plant growth, development and yields, and the C-repeat binding factors (CBFs) function in the cold resistance in plants. However, how pomegranate CBF transcription factors respond to cold signal remains unclear. Considering the significantly up-regulated expression of PgCBF3 and PgCBF7 in cold-tolerant Punica granatum ‘Yudazi’ in comparison with cold-sensitive ‘Tunisia’ under 4 °C, the present study focused on the two CBF genes. PgCBF3 was localized in the nucleus, while PgCBF7 was localized in the cell membrane, cytoplasm, and nucleus, both owning transcriptional activation activity in yeast. Yeast one-hybrid and dual-luciferase reporter assay further confirmed that PgICE1 could specifically bind to and significantly enhance the activation activity of the promoters of PgCBF3 and PgCBF7. Compared with the wild-type plants, the PgCBF3 and PgCBF7 transgenic Arabidopsis thaliana lines had the higher survival rate after cold treatment; exhibited increased the contents of soluble sugar and proline, while lower electrolyte leakage, malondialdehyde content, and reactive oxygen species production, accompanying with elevated enzyme activity of catalase, peroxidase, and superoxide dismutase; and upregulated the expression of AtCOR15A, AtCOR47, AtRD29A, and AtKIN1. Collectively, PgCBFs were positively regulated by the upstream PgICE1 and mediated the downstream COR genes expression, thereby enhancing freezing tolerance.

Genomic and Bioinformatic Resources for Perennial Fruit Species

In the post-genomic era, data management and development of bioinformatic tools are critical for the adequate exploitation of genomics data. In this review, we address the actual situation for the subset of crops represented by the perennial fruit species. The agronomical singularity of these species compared to plant and crop model species provides significant challenges on the implementation of good practices generally not addressed in other species. Studies are usually performed over several years in non-controlled environments, usage of rootstock is common, and breeders heavily rely on vegetative propagation. A reference genome is now available for all the major species as well as many members of the economically important genera for breeding purposes. Development of pangenome for these species is beginning to gain momentum which will require a substantial effort in term of bioinformatic tool development. The available tools for genome annotation and functional analysis will also be presented.

Genome-Wide Identification and Comprehensive Analysis of the AP2/ERF Gene Family in Pomegranate Fruit Development and Postharvest Preservation

Pomegranate (Punica granatum L.) is a kind of fruit with significant economic, ecological and health values. AP2/ERF transcription factors belong to a large group of factors mainly found in plants and play key roles in plant growth and development. However, AP2/ERF genes in pomegranate and their implication in development and postharvest preservation have been little described. In this study, 116 PgAP2/ERF genes in pomegranate were identified and renamed based on their chromosomal distributions. Phylogenetic relationship with genes from other species, structures, duplications, annotations, cis-elements in promoter sequences, and protein-protein interaction networks among PgAP2/ERF proteins were comprehensively explored. Expression analysis revealed several PgAP2/ERFs associated with the phenotypes of pomegranate seed hardness, including PgAP2/ERF5, PgAP2/ERF36, PgAP2/ERF58, and PgAP2/ERF86. Subsequent analysis indicated that many differentially expressed PgAP2/ERF genes are potentially important regulators of pomegranate fruit development. Furthermore, expression of more than one-half of PgAP2/ERFs was repressed in 'Tunisian soft seed' pomegranate fruit under low-temperature cold storage. The results showed that 1-MCP implicated in promoting postharvest preservation of 'Tunisian soft seed' pomegranate upregulated the PgAP2/ERF4, PgAP2/ERF15, PgAP2/ERF26, PgAP2/ERF30, PgAP2/ERF35 and PgAP2/ERF45 genes compared to those under low-temperature cold storage. This indicates that these genes are important candidate genes involved in pomegranate postharvest preservation. In summary, the findings of the present study provide an important basis for characterizing the PgAP2/ERF family genes and provide information on the candidate genes involved in pomegranate fruit development and postharvest preservation.

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.

Molecular Dissection Unveiling Dwarfing Effects of Plant Growth Retardants on Pomegranate

Dwarfed stature is a desired trait for modern orchard production systems. One effective strategy for dwarfing cultivation is exogenously applying plant growth retardants (PGRs) to plants. However, for many economic fruit trees, the current knowledge on the regulatory mechanisms underlying the dwarfing effect of PGRs is limited, which largely restricts the agricultural application of PGRs. In this study, we exogenously applied three kinds of PGRs [paclobutrazol, daminozide (B9), and mannitol] to the seedlings of pomegranate (Punica granatum L.) and performed comparative transcriptome analysis to elucidate the molecular features of PGR-induced dwarfing in pomegranates. Our results showed that all the three PGRs could significantly suppress plant growth of pomegranate. The inhibition of auxin biosynthetic processes, as well as auxin-mediated shoot development, may be considered as the main reason for the dwarfing. Besides that, different PGRs were also found to induce dwarfing via specific mechanisms, for example, cellular response to strigolactone was particularly suppressed by the application of paclobutrazol, while the level of carbohydrate homeostasis and metabolism were downregulated in conditions of either B9 or mannitol treatments. Furthermore, exogenous PGR application was supposed to cause adverse impacts on the normal physiological process of pomegranate seedlings, which may bring extra burden to pomegranate plants. These novel findings unveiled the genetic basis underlying the dwarfing in pomegranates, which provides deeper insights into PGR-mediated dwarfing cultivation of pomegranates.

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.

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.

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.

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.

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.