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.

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.

Genome-wide identification, gene cloning, subcellular location and expression analysis of SPL gene family in P. granatum L.

BACKGROUNDS: Pomegranate is an excellent tree species with nutritional, medicinal, ornamental and ecological values. Studies have confirmed that SPL factors play an important role in floral transition and flower development. RESULTS: Used bioinformatics methods, 15 SPL (SQUAMOSA promoter-binding protein-like) genes were identified and analyzed from the 'Taishanhong' pomegranate (P. granatum L.) genome. Phylogenetic analysis showed that PgSPLs were divided into six subfamilies (G1 ~ G6). PgSPL promoter sequences contained multiple cis-acting elements associated with abiotic stress or hormonal response. Based on the transcriptome data, expression profiles of different tissues and different developmental stages showed that PgSPL genes had distinct temporal and spatial expression characteristics. The expression analysis of miR156 in small RNA sequencing results showed that miR156 negatively regulated the expression of target genes. qRT-PCR analysis showed that the expression levels of PgSPL2, PgSPL3, PgSPL6, PgSPL11 and PgSPL14 in leaves were significantly higher than those in buds and stems (p < 0.05). The expression levels of PgSPL5, PgSPL12 and PgSPL13 in flower buds were significantly higher than that in leaves and stems (p < 0.05). The full-length of coding sequence of PgSPL5 and PgSPL13 were obtained by homologous cloning technology. The full length of PgSPL5 is 1020 bp, and PgSPL13 is 489 bp, which encodes 339 and 162 amino acids, respectively. Further investigation revealed that PgSPL5 and PgSPL13 proteins were located in the nucleus. Exogenous plant growth regulator induction experiments showed that PgSPL5 was up-regulated in leaves and stems. PgSPL13 was up-regulated in leaves and down-regulated in stems. When sprayed with 6-BA, IBA and PP333 respectively, PgSPL5 and PgSPL13 were up-regulated most significantly at P2 (bud vertical diameter was 5.1 ~ 12.0 mm) stage of bisexual and functional male flowers. CONCLUSIONS: Our findings suggested that PgSPL2, PgSPL3, PgSPL6, PgSPL11 and PgSPL14 played roles in leaves development of pomegranate. PgSPL5, PgSPL12 and PgSPL13 played roles in pomegranate flower development. PgSPL5 and PgSPL13 were involved in the response process of different plant hormone signal transduction in pomegranate development. This study provided a robust basis for further functional analyses of SPL genes in pomegranate.

Epigenetically modified N6-methyladenine inhibits DNA replication by human DNA polymerase iota.

Among human four Y-family DNA polymerases, hPol ι is exceptionally error-prone in DNA synthesis. 6 mA plays significant roles in epigenetic regulation of numerous biological processes. Nonetheless, its effects on DNA replication by hPol ι is still unclear. In this work, we found that 6 mA and Hyp, the intermediate of 6 mA, inhibited the replication of DNA by hPol ι. 6 mA lost priority in extension beyond 6 mA:T pair, partially reducing dTTP incorporation efficiency and inhibiting next-base extension. Hyp was prone to dCTP incorporation and extension beyond Hyp:C instead of Hyp:T pair. Statistically, 6 mA primarily reduced the burst incorporation rate (kpol) and slightly increased the dissociation constant (Kd,dTTP). However, Hyp mainly increased the Kd,dCTP yet did not affect the kpol, both reducing the burst incorporation efficiency (kpol/Kd,dCTP). 6 mA together with Hyp weakened the binding affinity of hPol ι to DNA in binary or ternary complex. The misincorporation opposite 6 mA or Hyp further weakened this binding affinity. The methyl group in 6 mA doesn't almost affect the H-bond formation with dTTP, therefore mildly inhibiting dTTP incorporation. As an analogue of G, Hyp can form only two H-bonds with dCTP, thus reducing dCTP incorporation. This work provides a new insight in how the epigenetically modified 6 mA and its intermediate Hyp affect replication of DNA by human DNA polymerase ι.

Ribonucleoside triphosphates promote T7 DNA replication and the lysis of T7-Infected Escherichia coli.

rNTPs are structurally similar to dNTPs, but their concentrations are much higher than those of dNTPs in cells. rNTPs in solutions or rNMP at the primer terminus or embedded in template always inhibit or block DNA replication, due to the reduced Mg2+ apparent concentration, competition of rNTPs with dNTPs, and the extra repulsive interaction of rNTP or rNMP with polymerase active site. In this work, unexpectedly, we found rNTPs can promote T7 DNA replication with the maximal promotion at rNTPs/dNTPs concentration ratio of 20. This promotion was not due to the optimized Mg2+ apparent concentration or the direct incorporation of extra rNMPs into DNA. This promotion was dependent on the concentrations and types of rNTPs. Kinetic analysis showed that this promotion was originated from the increased fraction of polymerase-DNA productive complex and the accelerated DNA polymerization. Further evidence showed that more polymerase-DNA complex was formed and their binding affinity was also enhanced in the presence of extra rNTPs. Moreover, this promotion in T7 DNA replication also accelerated the lysis of T7-infected host Escherichia coli. This work discovered that rNTPs could promote DNA replication, completely different from the traditional concept that rNTPs always inhibit DNA replication.

Epigenetically modified N6-methyladenine inhibits DNA replication by human DNA polymerase η.

N6-methyladenine (6mA), as a newly reported epigenetic marker, plays significant roles in regulation of various biological processes in eukaryotes. However, the effect of 6mA on human DNA replication remain elusive. In this work, we used Y-family human DNA polymerase η as a model to investigate the kinetics of bypass of 6mA by hPol η. We found 6mA and its intermediate hypoxanthine (I) on template partially inhibited DNA replication by hPol η. dTMP incorporation opposite 6mA and dCMP incorporation opposite I can be considered as correct incorporation. However, both 6mA and I reduced correct incorporation efficiency, next-base extension efficiency, and the priority in extension beyond correct base pair. Both dTMP incorporation opposite 6mA and dCTP opposite I showed fast burst phases. However, 6mA and I reduced the burst incorporation rates (kpol) and increased the dissociation constant (Kd,dNTP), compared with that of dTMP incorporation opposite unmodified A. Biophysical binding assays revealed that both 6mA and I on template reduced the binding affinity of hPol η to DNA in binary or ternary complex compared with unmodified A. All the results explain the inhibition effects of 6mA and I on DNA replication by hPol η, providing new insight in the effects of epigenetically modified 6mA on human DNA replication.

Epigenetic DNA Modification N6-Methyladenine Inhibits DNA Replication by DNA Polymerase of Pseudomonas aeruginosa Phage PaP1.

N6-methyladenine (6mA), a newly identified epigenetic modification, plays important roles in regulation of various biological processes. However, the effect of 6mA on DNA replication has been little addressed. In this work, we investigated how 6mA affected DNA replication by DNA polymerase of Pseudomonas aeruginosa Phage PaP1 (gp90 exo-). The presence of 6mA, as well as its intermediate hypoxanthine (Hyp), inhibited DNA replication by gp90 exo-. The 6mA reduced dTTP incorporation efficiency by 10-fold and inhibited next-base extension efficiency by 100-fold. Differently, dCTP was preferentially incorporated opposite Hyp among four dNTPs. Gp90 exo- reduced the extension priority beyond the 6mA:T pair rather than the 6mA:C mispair and preferred to extend beyond Hyp:C rather than the Hyp:T pair. Incorporation of dTTP opposite 6mA and dCTP opposite Hyp showed fast burst phases. The burst rate and burst amplitude were both reduced for 6mA compared with unmodified A. Moreover, the total incorporation efficiency ( kpol/ Kd,dNTP) was decreased for dTTP incorporation opposite 6mA and dCTP incorporation opposite Hyp compared with dTTP incorporation opposite A. 6mA reduced the incorporation rate ( kpol), and Hyp increased the dissociation constant ( Kd,dNTP). However, 6mA or Hyp on template did not affect the binding of DNA polymerase to DNA in binary or ternary complexes. This work provides new insight into the inhibited effects of epigenetic modification of 6mA on DNA replication in PaP1.

Consecutive ribonucleoside monophosphates on template inhibit DNA replication by T7 DNA polymerase or by T7 polymerase and helicase complex.

rNTPs are structurally similar to dNTPs but are largely molar excessive in cell than dNTPs. rNTPs are inevitably incorporated into DNA to form rNMPs. Long RNA primers can also be incorporated into lagging-strand DNA. However, the influence of these incorporated rNMPs on T7 DNA replication remains unknown. In this work, we investigated primer extension past consecutive rNMPs (rA, r(AC), r(ACC), or r(ACCA)) on template by T7 DNA polymerase or by its complex with helicase. Primer extension is gradually inhibited with increasing rNMP number. rNMPs decrease the dNTP incorporation efficiency, slightly weaken the binding affinity of polymerase to DNA in ternary complex, and reduce the protein interaction between polymerase and helicase at DNA fork, thereby decreasing the fraction of the productive enzyme·DNA complex and the average primer extension rate. Therefore, the consecutive rNMPs on template gradually inhibit T7 primer extension and strand-displacement DNA synthesis, providing the kinetic information for the inhibition of rNMPs on DNA replication.

General misincorporation frequency: Re-evaluation of the fidelity of DNA polymerases.

DNA replication in cells is performed in the presence of four dNTPs and four rNTPs. In this study, we re-evaluated the fidelity of DNA polymerases using the general misincorporation frequency consisting of three incorrect dNTPs and four rNTPs but not using the traditional special misincorporation frequency with only the three incorrect dNTPs. We analyzed both the general and special misincorporation frequencies of nucleotide incorporation opposite dG, rG, or 8-oxoG by Pseudomonas aeruginosa phage 1 (PaP1) DNA polymerase Gp90 or Sulfolobus solfataricus DNA polymerase Dpo4. Both misincorporation frequencies of other DNA polymerases published were also summarized and analyzed. The general misincorporation frequency is obviously higher than the special misincorporation frequency for many DNA polymerases, indicating the real fidelity of a DNA polymerase should be evaluated using the general misincorporation frequency.