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

Complete and highly accurate reference genomes and gene annotations are indispensable for basic biological research and trait improvement of woody tree species. In this study, we integrated single-molecule sequencing and high-throughput chromosome conformation capture techniques to produce a high-quality and long-range contiguity chromosome-scale genome assembly of the soft-seeded pomegranate cultivar 'Tunisia'. The genome covers 320.31 Mb (scaffold N50 = 39.96 Mb; contig N50 = 4.49 Mb) and includes 33 594 protein-coding genes. We also resequenced 26 pomegranate varieties that varied regarding seed hardness. Comparative genomic analyses revealed many genetic differences between soft- and hard-seeded pomegranate varieties. A set of selective loci containing SUC8-like, SUC6, FoxO and MAPK were identified by the selective sweep analysis between hard- and soft-seeded populations. An exceptionally large selective region (26.2 Mb) was identified on chromosome 1. Our assembled pomegranate genome is more complete than other currently available genome assemblies. Our results indicate that genomic variations and selective genes may have contributed to the genetic divergence between soft- and hard-seeded pomegranate varieties.

Identification of the SUT Gene Family in Pomegranate (Punica granatum L.) and Functional Analysis of PgL0145810.1.

Sucrose, an important sugar, is transported from source to sink tissues through the phloem, and plays important role in the development of important traits in plants. However, the SUT gene family is still not well characterized in pomegranate. In this study, we first identified the pomegranate sucrose transporter (SUT) gene family from the whole genome. Then, the phylogenetic relationship of SUT genes, gene structure and their promoters were analyzed. Additionally, their expression patterns were detected during the development of the seed. Lastly, genetic transformation and cytological observation were used to study the function of PgL0145810.1. A total of ten pomegranate SUT genes were identified from the whole genome of pomegranate 'Tunisia'. The promoter region of all the pomegranate SUT genes contained myeloblastosis (MYB) elements. Four of the SUT genes, PgL0328370.1, PgL0099690.1, PgL0145810.1 and PgL0145770.1, were differentially expressed during seed development. We further noticed that PgL0145810.1 was expressed most prominently in the stem parts in transgenic plants compared to other tissue parts (leaves, flowers and silique). The cells in the xylem vessels were small and lignin content was lower in the transgenic plants as compared to wild Arabidopsis plants. In general, our result suggests that the MYB cis-elements in the promoter region might regulate PgL0145810.1 expression to control the structure of xylem, thereby affecting seed hardness in pomegranate.

Small RNA and mRNA Sequencing Reveal the Roles of microRNAs Involved in Pomegranate Female Sterility.

Female sterility is a key factor restricting plant reproduction. Our previous studies have revealed that pomegranate female sterility mainly arose from the abnormality of ovule development. MicroRNAs (miRNAs) play important roles in ovule development. However, little is known about the roles of miRNAs in female sterility. In this study, a combined high-throughput sequencing approach was used to investigate the miRNAs and their targeted transcripts involved in female development. A total of 103 conserved and 58 novel miRNAs were identified. Comparative profiling indicated that the expression of 43 known miRNAs and 14 novel miRNAs were differentially expressed between functional male flowers (FMFs) and bisexual flowers (BFs), 30 known miRNAs and nine novel miRNAs showed significant differences among different stages of BFs, and 20 known miRNAs and 18 novel miRNAs exhibited remarkable expression differences among different stages of FMFs. Gene ontology (GO) analyses of 144 predicted targets of differentially expressed miRNAs indicated that the "reproduction process" and "floral whorl development" processes were significantly enriched. The miRNA-mRNA interaction analyses revealed six pairs of candidate miRNAs and their targets associated with female sterility. Interestingly, pg-miR166a-3p was accumulated, whereas its predicted targets (Gglean012177.1 and Gglean013966.1) were repressed in functional male flowers (FMFs), and the interaction between pg-miR166a-3p and its targets (Gglean012177.1 and Gglean013966.1) were confirmed by transient assay. A. thaliana transformed with 35S-pre-pg-miR166a-3p verified the role of pg-miR166a-3p in ovule development, which indicated pg-miR166a-3p's potential role in pomegranate female sterility. The results provide new insights into molecular mechanisms underlying the female sterility at the miRNA level.

Comparative Transcriptome Analysis of Genes Involved in Anthocyanin Biosynthesis in Red and Green Walnut (Juglans regia L.).

Fruit color is an important economic trait. The color of red walnut cultivars is mainly attributed to anthocyanins. The aim of this study was to explore the differences in the molecular mechanism of leaf and peel color change between red and green walnut. A reference transcriptome of walnut was sequenced and annotated to identify genes related to fruit color at the ripening stage. More than 290 million high-quality reads were assembled into 39,411 genes using a combined assembly strategy. Using Illumina digital gene expression profiling, we identified 4568 differentially expressed genes (DEGs) between red and green walnut leaf and 3038 DEGs between red and green walnut peel at the ripening stage. We also identified some transcription factor families (MYB, bHLH, and WD40) involved in the control of anthocyanin biosynthesis. The trends in the expression levels of several genes encoding anthocyanin biosynthetic enzymes and transcription factors in the leaf and peel of red and green walnut were verified by quantitative real-time PCR. Together, our results identified the genes involved in anthocyanin accumulation in red walnut. These data provide a valuable resource for understanding the coloration of red walnut.

Comparative transcriptome analysis reveals a global insight into molecular processes regulating citrate accumulation in sweet orange (Citrus sinensis).

Citrate, the predominant organic acid in citrus, determines the taste of these fruits. However, little is known about the synergic molecular processes regulating citrate accumulation. Using 'Dahongtiancheng' (Citrus sinensis) and 'Bingtangcheng' (C. sinensis) with significant difference in citrate, the objectives of this study were to understand the global mechanisms of high-citrate accumulation in sweet orange. 'Dahongtiancheng' and 'Bingtangcheng' exhibit significantly different patterns in citrate accumulation throughout fruit development, with the largest differences observed at 50-70 days after full bloom (DAFB). Comparative transcriptome profiling was performed for the endocarps of both cultivars at 50 and 70 DAFB. Over 34.5 million clean reads per library were successfully mapped to the reference database and 670-2630 differentially expressed genes (DEGs) were found in four libraries. Among the genes, five transcription factors were ascertained to be the candidates regulating citrate accumulation. Functional assignments of the DEGs indicated that photosynthesis, the citrate cycle and amino acid metabolism were significantly altered in 'Dahongtiancheng'. Physiological and molecular analyses suggested that high photosynthetic efficiency and partial impairment of citrate catabolism were crucial for the high-citrate trait, and amino acid biosynthesis was one of the important directions for citrate flux. The results reveal a global insight into the gene expression changes in a high-citrate compared with a low-citrate sweet orange. High accumulating efficiency and impaired degradation of citrate may be associated with the high-citrate trait of 'Dahongtiancheng'. Findings in this study increase understanding of the molecular processes regulating citrate accumulation in sweet orange.

Promoter Variation of the Key Apple Fruit Texture Related Gene MdPG1 and the Upstream Regulation Analysis.

MdPG1 encoding polygalacturonase in apple (Malus × domestica) is a key gene associated with fruit firmness and texture variations among apple cultivars. However, the causative variants of MdPG1 are still not known. In this study, we identified a SNPA/C variant within an ERF-binding element located in the promoter region of MdPG1. The promoter containing the ERF-binding element with SNPA, rather than the SNPC, could be strongly bound and activated by MdCBF2, a member of the AP2/ERF transcription factor family, as determined by yeast-one-hybrid and dual-luciferase reporter assays. We also demonstrated that the presence of a novel long non-coding RNA, lncRNAPG1, in the promoter of MdPG1 was a causative variant. lncRNAPG1 was specifically expressed in fruit tissues postharvest. lncRNAPG1 could reduce promoter activity when it was fused to the promoter of MdPG1 and a tobacco gene encoding Mg-chelatase H subunit (NtCHLH) in transgenic tobacco cells but could not reduce promoter activity when it was supplied in a separate gene construct, indicating a cis-regulatory effect. Our results provide new insights into genetic regulation of MdPG1 allele expression and are also useful for the development of elite apple cultivars.

Characterization of a NAC transcription factor involved in the regulation of pomegranate seed hardness (Punica granatum L.).

The pomegranate, Punica granatum L., which has been cultivated since antiquity, is known to be a superfruit, possessing an array of functional anti-oxidants and various other health benefits. The hardness of pomegranate seeds is an important indicator of fruit quality, which in turn affects economic value and market demand. However, the molecular mechanism underlying pomegranate seed hardness remains to be fully understood. In this study, we found a positive correlation between seed hardness and lignin content in two pomegranate varieties: "Tunisia" and "Sanbai". Specifically, genes associated with lignin biosynthesis were differentially expressed in soft-seed and hard-seed pomegranate varieties. Among these differential genes, we cloned and characterized the NAC transcription factor PgSND1-like. Sequence alignment found a single base replacement at the 166-bp position of CDS in the PgSND1-like gene from "Tunisia" and "Sanbai". Both PgSND1-like (Sanbai) and PgSND1-like (Tunisia) proteins are localized in the cell nucleus and have a transcription activation domain in the C-terminus. Yeast two-hybrid analysis indicated that PgSND1-like protein interacts with itself to form a homodimer. Overexpression of PgSND1-like (Sanbai) in Arabidopsis showed a higher lignin content in inflorescence stem and mature seed compared with wild-type Arabidopsis. Accordingly, the expression levels of several lignin biosynthesis-associated genes were upregulated in stem cells and mature seeds of transgenic plants. However, PgSND1-like (Tunisia) transgenic Arabidopsis showed no phenotypic differences with wild-type Arabidopsis. Taken together, we suggest that PgSND1-like may regulate at least two different functions in two pomegranate varieties, promoting lignin biosynthesis and seed hardness of pomegranate.

Complementary iTRAQ-based proteomic and RNA sequencing-based transcriptomic analyses reveal a complex network regulating pomegranate (Punica granatum L.) fruit peel colour.

Peel colour is an important factor affecting the marketability of pomegranate fruits. Therefore, elucidating the genetic mechanism of fruit peel colour development may be useful for breeding pomegranate cultivars with enhanced fruit peel colours. In this study, we combined an iTRAQ-based proteome-level analysis with an RNA sequencing-based transcriptome-level analysis to detect the proteins and genes related to fruit peel colour development in pomegranate. We analysed the 'Tunisia' (red fruit) and 'White' (white fruit) pomegranate cultivars at two stages of fruit development. A total of 27 differentially abundant proteins (increased abundance) and 54 differentially expressed genes (16 up-regulated and 38 down-regulated) were identified from our proteomics and transcriptomics data. The identified proteins and genes contribute to pomegranate fruit peel colour by participating in the biosynthesis of anthocyanins, stilbenoids, diarylheptanoids, gingerols, flavonoids, and phenylpropanoids. Several candidate proteins and genes corresponded to enzymes related to general reactions (PAL, 4CL, DFR, LDOX/ANS, CHS, and F3'5'H) and glycosylation (GT1 and UGAT) of compounds and pigments related to the colour of pomegranate fruit peel. Complementary proteome- and transcriptome-level analyses revealed a complex molecular network controlling fruit peel colour. The candidate genes identified in this study may be useful for the marker-based breeding of new pomegranate cultivars.

Integrated microRNA and mRNA expression profiling reveals a complex network regulating pomegranate (Punica granatum L.) seed hardness.

The breeding of new soft-seeded pomegranate cultivars provides new products for the market and increases farmers' incomes, yet the genetic architecture mediating seed hardness is largely unknown. Here, the seed hardness and hundred-seed weights of 26 cultivars were determined in 2 successive years. We conducted miRNA and mRNA sequencing to analyse the seeds of two varieties of Punica granatum: soft-seeded Tunisia and hard-seeded Sanbai, at 60 and 120 d after flowering. Seed hardness was strongly positively correlated with hundred-seed weight. We detected 25 and 12 differentially expressed miRNA-mRNA pairs with negative regulatory relationships between the two genotypes at 60 and 120 d after flowering, respectively. These miRNA-mRNA pairs mainly regulated seed hardness by altering cell wall structure. Transcription factors including NAC1, WRKY and MYC, which are involved in seed hardness, were targeted by differentially expressed mdm-miR164e and mdm-miR172b. Thus, seed hardness is the result of a complex biological process regulated by a miRNA-mRNA network in pomegranate. These results will help us understand the complexity of seed hardness and help to elucidate the miRNA-mediated molecular mechanisms that contribute to seed hardness in pomegranate.

De novo transcriptome assembly and quantification reveal differentially expressed genes between soft-seed and hard-seed pomegranate (Punica granatum L.).

Pomegranate (Punica granatum L.) belongs to Punicaceae, and is valued for its social, ecological, economic, and aesthetic values, as well as more recently for its health benefits. The 'Tunisia' variety has softer seeds and big arils that are easily swallowed. It is a widely popular fruit; however, the molecular mechanisms of the formation of hard and soft seeds is not yet clear. We conducted a de novo assembly of the seed transcriptome in P. granatum L. and revealed differential gene expression between the soft-seed and hard-seed pomegranate varieties. A total of 35.1 Gb of data were acquired in this study, including 280,881,106 raw reads. Additionally, de novo transcriptome assembly generated 132,287 transcripts and 105,743 representative unigenes; approximately 13,805 unigenes (37.7%) were longer than 1,000 bp. Using bioinformatics annotation libraries, a total of 76,806 unigenes were annotated and, among the high-quality reads, 72.63% had at least one significant match to an existing gene model. Gene expression and differentially expressed genes were analyzed. The seed formation of the two pomegranate cultivars involves lignin biosynthesis and metabolism, including some genes encoding laccase and peroxidase, WRKY, MYB, and NAC transcription factors. In the hard-seed pomegranate, lignin-related genes and cellulose synthesis-related genes were highly expressed; in soft-seed pomegranates, expression of genes related to flavonoids and programmed cell death was slightly higher. We validated selection of the identified genes using qRT-PCR. This is the first transcriptome analysis of P. granatum L. This transcription sequencing greatly enriched the pomegranate molecular database, and the high-quality SSRs generated in this study will aid the gene cloning from pomegranate in the future. It provides important insights into the molecular mechanisms underlying the formation of soft seeds in pomegranate.

Transcriptome and Metabolome Analyses Provide Insights into the Occurrence of Peel Roughing Disorder on Satsuma Mandarin (Citrus unshiu Marc.) Fruit.

Roughing disorder (RD) is a significant quality barrier in citrus fruit, prevalent on easy-peeling mandarins. As RD is not yet well-understood, this study aimed to examine the changes and synergic molecular processes involved in peel RD. Peel with RD was induced by severely defruiting Satsuma mandarin trees. Morphology observations, RNA-sequencing, and targeted and untargeted metabolic analyses were conducted. The results showed that the primary metabolites of sugars, organic acids and amino acids are dramatically changed in RD peel. The RD peel was always characterized by higher magnesium content during development. Comparative transcriptome profiling was performed for CK and RD peels at 30, 80, and 170 days after full bloom (DAFB) which represented fruit at cell division stage, cell enlargement stage and fruit maturity stage, respectively. Physiological and molecular biological evidence suggested that the month after full bloom is a crucial stage for RD initiation. A total of 4,855 differentially expressed genes (DEGs) in RD peel, relative to CK peel were detected at cell division stage, about 2 to 4-fold more than other stages had. Among the differentially expressed transcription factors, the bHLH family were affected most by RD, and six bHLH transcription factors functionally involved in GA metabolism were assessed to associate with RD occurrence. Gene set enrichment analysis suggested that RD significantly altered starch and GA metabolism in peel. Higher starch content and hydrolysed chain status were found in RD peel at cell division stage. RD occurrence on the peel was influenced significantly by GA, especially abundant GA before July. These changes may mean a significant alteration in sink strength of RD peel. The findings of this study provide insights into the emergence, development and molecular mechanisms of RD.

Transcriptomic Analysis Reveals Candidate Genes for Female Sterility in Pomegranate Flowers.

Pomegranate has two types of flowers on the same plant: functional male flowers (FMF) and bisexual flowers (BF). BF are female-fertile flowers that can set fruits. FMF are female-sterile flowers that fail to set fruit and that eventually drop. The putative cause of pomegranate FMF female sterility is abnormal ovule development. However, the key stage at which the FMF pomegranate ovules become abnormal and the mechanism of regulation of pomegranate female sterility remain unknown. Here, we studied ovule development in FMF and BF, using scanning electron microscopy to explore the key stage at which ovule development was terminated and then analyzed genes differentially expressed (differentially expressed genes - DEGs) between FMF and BF to investigate the mechanism responsible for pomegranate female sterility. Ovule development in FMF ceased following the formation of the inner integument primordium. The key stage for the termination of FMF ovule development was when the bud vertical diameter was 5.0-13.0 mm. Candidate genes influencing ovule development may be crucial factors in pomegranate female sterility. INNER OUTER (INO/YABBY4) (Gglean016270) and AINTEGUMENTA (ANT) homolog genes (Gglean003340 and Gglean011480), which regulate the development of the integument, showed down-regulation in FMF at the key stage of ovule development cessation (ATNSII). Their upstream regulator genes, such as AGAMOUS-like (AG-like) (Gglean028014, Gglean026618, and Gglean028632) and SPOROCYTELESS (SPL) homolog genes (Gglean005812), also showed differential expression pattern between BF and FMF at this key stage. The differential expression of the ethylene response signal genes, ETR (ethylene-resistant) (Gglean022853) and ERF1/2 (ethylene-responsive factor) (Gglean022880), between FMF and BF indicated that ethylene signaling may also be an important factor in the development of pomegranate female sterility. The increase in BF observed after spraying with ethephon supported this interpretation. Results from qRT-PCR confirmed the findings of the transcriptomic analysis.