Cytological, Phytohormone, and Transcriptome Analyses Provide Insights into Persimmon Fruit Shape Formation (Diospyros kaki Thunb.).

Fruit shape is an important external feature when consumers choose their preferred fruit varieties. Studying persimmon (Diospyros kaki Thunb.) fruit shape is beneficial to increasing its commodity value. However, research on persimmon fruit shape is still in the initial stage. In this study, the mechanism of fruit shape formation was studied by cytological observations, phytohormone assays, and transcriptome analysis using the long fruit and flat fruit produced by 'Yaoxianwuhua' hermaphroditic flowers. The results showed that stage 2-3 (June 11-June 25) was the critical period for persimmon fruit shape formation. Persimmon fruit shape is determined by cell number in the transverse direction and cell length in the longitudinal direction. High IAA, GA4, ZT, and BR levels may promote long fruit formation by promoting cell elongation in the longitudinal direction, and high GA3 and ABA levels may be more conducive to flat fruit formation by increasing the cell number in the transverse direction and inhibiting cell elongation in the longitudinal direction, respectively. Thirty-two DEGs related to phytohormone biosynthesis and signaling pathways and nine DEGs related to cell division and cell expansion may be involved in the persimmon fruit shape formation process. These results provide valuable information for regulatory mechanism research on persimmon fruit formation.

Genetic insights into the dissolution of dioecy in diploid persimmon Diospyros oleifera Cheng.

BACKGROUND: Dioecy, a sexual system of single-sexual (gynoecious/androecious) individuals, is rare in flowering plants. This rarity may be a result of the frequent transition from dioecy into systems with co-sexual individuals. RESULTS: In this study, co-sexual expression (monoecy and hermaphroditic development), previously thought to be polyploid-specific in Diospyros species, was identified in the diploid D. oleifeara historically. We characterized potential genetic mechanisms that underlie the dissolution of dioecy to monoecy and andro(gyno)monoecy, based on multiscale genome-wide investigations of 150 accessions of Diospyros oleifera. We found all co-sexual plants, including monoecious and andro(gyno)monoecious individuals, possessed the male determinant gene OGI, implying the presence of genetic factors controlling gynoecia development in genetically male D. oleifera. Importantly, discrepancies in the OGI/MeGI module were found in diploid monoecious D. oleifera compared with polyploid monoecious D. kaki, including no Kali insertion on the promoter of OGI, no different abundance of smRNAs targeting MeGI (a counterpart of OGI), and no different expression of MeGI between female and male floral buds. On the contrary, in both single- and co-sexual plants, female function was expressed in the presence of a genome-wide decrease in methylation levels, along with sexually distinct regulatory networks of smRNAs and their targets. Furthermore, a genome-wide association study (GWAS) identified a genomic region and a DUF247 gene cluster strongly associated with the monoecious phenotype and several regions that may contribute to andromonoecy. CONCLUSIONS: Collectively, our findings demonstrate stable breakdown of the dioecious system in D. oleifera, presumably also a result of genomic features of the Y-linked region.

Comparative Metabolomic and Transcriptomic Analyses Reveal Distinct Ascorbic Acid (AsA) Accumulation Patterns between PCA and PCNA Persimmon Developing Fruit.

Persimmon fruit has a high nutritional value and significantly varies between pollination-constant astringent (PCA) and pollination-constant non-astringent (PCNA) persimmons. The astringency type affects sugar, flavonoids, and tannin accumulation and is well known in persimmon fruit. However, the impact of the fruit astringency type on ascorbic acid (AsA) accumulation is limited. In this study, typical PCA varieties ('Huojing' and 'Zhongshi5') and PCNA varieties ('Yohou' and 'Jiro') of persimmon fruit were sampled at four developing stages (S1-S4) to provide valuable information on AsA content variation in PCA and PCNA persimmon. Persimmon fruit is rich in ascorbic acid; the AsA content of the four varieties 'Zhongshi5', 'Huojing', 'Jiro', and 'Youhou' mature fruit reached 104.49, 48.69, 69.69, and 47.48 mg/100 g. Fruit of the same astringency type persimmon showed a similar AsA accumulation pattern. AsA content was significantly higher in PCA than PCNA fruit at S1-S3. The initial KEGG analysis of metabolites showed that galactose metabolism is the major biosynthetic pathway of AsA in persimmon fruit. There were significant differences in galactose pathway-related metabolite content in developing PCA and PCNA fruit, such as Lactose, D-Tagatose, and D-Sorbitol content in PCA being higher than that of PCNA. Combined gene expression and WGCNA analyses showed that the expression of the GME (evm.TU.contig4144.37) gene was higher in PCA-type than in PCNA-type fruit in S1-S3 and exhibited the highest correlation with AsA content (r = 690 **, p < 0.01). Four hub genes, including the DNA methylation gene, methyltransferase gene, F-box, and Actin-like Protein, were identified as potential regulators of the GME gene. These results provide basic information on how astringency types affect AsA accumulation and will provide valuable information for further investigation on AsA content variation in persimmon fruit.

Transcriptome and Metabolome Reveal Distinct Sugar Accumulation Pattern between PCNA and PCA Mature Persimmon Fruit.

Persimmon (Diospyros kaki) fruit have significant variation between pollination-constant non-astringent (PCNA) and pollination-constant astringent (PCA) persimmons. The astringency type affects not only the soluble tannin concentration but also the accumulation of individual sugars. Thus, we comprehensively investigate the gene expression and metabolite profiles of individual sugars to resolve the formation of flavor differences in PCNA and PCA persimmon fruit. The results showed that soluble sugar, starch content, sucrose synthase, and sucrose invertase were significantly different between PCNA and PCA persimmon fruit. The sucrose and starch metabolism pathway was considerably enriched, and six sugar metabolites involving this pathway were significantly differentially accumulated. In addition, the expression patterns of diferentially expressed genes (such as bglX, eglC, Cel, TPS, SUS, and TREH genes) were significantly correlated with the content of deferentially accumulated metabolites (such as starch, sucrose, and trehalose) in the sucrose and starch metabolism pathway. These results indicated that the sucrose and starch metabolism pathway maintained a central position of sugar metabolism between PCNA and PCA persimmon fruit. Our results provide a theoretical basis for exploring functional genes related to sugar metabolism and provide useful resources for future studies on the flavor differences between PCNA and PCA persimmon fruit.

Allele-aware chromosome-level genome assembly of the autohexaploid Diospyros kaki Thunb.

Artificially improving persimmon (Diospyros kaki Thunb.), one of the most important fruit trees, remains challenging owing to the lack of reference genomes. In this study, we generated an allele-aware chromosome-level genome assembly for the autohexaploid persimmon 'Xiaoguotianshi' (Chinese-PCNA type) using PacBio CCS and Hi-C technology. The final assembly contained 4.52 Gb, with a contig N50 value of 5.28 Mb and scaffold N50 value of 44.01 Mb, of which 4.06 Gb (89.87%) of the assembly were anchored onto 90 chromosome-level pseudomolecules comprising 15 homologous groups with 6 allelic chromosomes in each. A total of 153,288 protein-coding genes were predicted, of which 98.60% were functionally annotated. Repetitive sequences accounted for 64.02% of the genome; and 110,480 rRNAs, 12,297 tRNAs, 1,483 miRNAs, and 3,510 snRNA genes were also identified. This genome assembly fills the knowledge gap in the autohexaploid persimmon genome, which is conducive in the study on the regulatory mechanisms underlying the major economically advantageous traits of persimmons and promoting breeding programs.

Regulatory mechanism of MeGI on sexuality in Diospyros oleifera.

Dioecy system is an important strategy for maintaining genetic diversity. The transcription factor MeGI, contributes to dioecy by promoting gynoecium development in Diospyros lotus and D. kaki. However, the function of MeGI in D. oleifera has not been identified. In this study, we confirmed that MeGI, cloned from D. oleifera, repressed the androecium development in Arabidopsis thaliana. Subsequently, chromatin immunoprecipitation-sequencing (ChIP-seq), DNA affinity purification-sequencing (DAP-seq), and RNA-seq were used to uncover the gene expression response to MeGI. The results showed that the genes upregulated and downregulated in response to MeGI were mainly enriched in the circadian rhythm-related and flavonoid biosynthetic pathways, respectively. Additionally, the WRKY DNA-binding protein 28 (WRKY28) gene, which was detected by ChIP-seq, DAP-seq, and RNA-seq, was emphasized. WRKY28 has been reported to inhibit salicylic acid (SA) biosynthesis and was upregulated in MeGI-overexpressing A. thaliana flowers, suggesting that MeGI represses the SA level by increasing the expression level of WRKY28. This was confirmed that SA level was lower in D. oleifera female floral buds than male. Overall, our findings indicate that the MeGI mediates its sex control function in D. oleifera mainly by regulating genes in the circadian rhythm, SA biosynthetic, and flavonoid biosynthetic pathways.

Comparative transcriptomic and metabolomic analyses reveal differences in flavonoid biosynthesis between PCNA and PCA persimmon fruit.

The fruit of the persimmon (Diospyros kaki.) has high economic and nutritional value and is rich in flavonoids. Flavonoids are essential secondary metabolisms in plants. The association between persimmon astringency and changes in the proanthocyanidins (a flavonoid subclass) content is well-known. However, information on the relationships between different astringency types and other flavonoid subclasses and biosynthetic genes is more limited. In this study, an initial correlation analysis between total flavonoids and fruit astringency type, and KEGG analysis of metabolites showed that flavonoid-related pathways were linked to differences between mature pollination-constant non-astringent (PCNA) varieties ('Jiro' and 'Yohou') and pollination-constant astringent (PCA) fruit varieties ('Zhongshi5' and 'Huojing'). Based on these findings, variations in the expression of genes and metabolites associated with flavonoid biosynthesis were investigated between typical PCNA ('Jiro') and PCA ('Huojing') persimmons during fruit development. The flavonoid concentration in 'Huojing' fruit was significantly higher than that of 'Jiro' fruit, especially, in levels of proanthocyanin precursor epicatechin and anthocyanin cyanidin derivatives. Combined WGCNA and KEGG analyses showed that genes such as PAL, C4H, CHI, CHS, F3H, F3'5'H, FLS, DFR, ANR, ANS, and UF3GT in the phenylpropanoid and flavonoid biosynthesis pathways may be significant factors impacting the proanthocyanin precursor and anthocyanin contents. Moreover, interactions between the R2R3MYB (evm.TU.contig7272.598) and WD40 (evm.TU.contig3208.5) transcription factors were found to be associated with the above structural genes. These findings provide essential information on flavonoid biosynthesis and its regulation in the persimmon and lay a foundation for further investigation into how astringency types affect flavor components in PCNA and PCA persimmons.

Transcriptomic profiling analysis to identify genes associated with PA biosynthesis and insolubilization in the late stage of fruit development in C-PCNA persimmon.

PA-enhanced content causes astringency in persimmon fruit. PCNA persimmons can lose their astringency naturally and they become edible when still on the tree, which allows for conserves of physical and financial resources. C-PCNA persimmon originates in China. Its deastringency trait primarily depends on decreased PA biosynthesis and PA insolubilization at the late stage of fruit development. Although some genes and transcription factors that may be involved in the deastringency of C-PCNA persimmon have been reported, the expression patterns of these genes during the key deastringency stage are reported less. To investigate the variation in PA contents and the expression patterns of deastringency-related genes during typical C-PCNA persimmon 'Xiaoguo-tianshi' fruit development and ripening, PA content and transcriptional profiling were carried out at five late stages from 70 to 160 DAF. The combinational analysis phenotype, PA content, and DEG enrichment revealed that 120-140 DAF and 140-160 DAF were the critical phases for PA biosynthesis reduction and PA insolubilization, respectively. The expression of PA biosynthesis-associated genes indicated that the downregulation of the ANR gene at 140-160 DAF may be associated with PA biosynthesis and is decreased by inhibiting its precursor cis-flavan-3-ols. We also found that a decrease in acetaldehyde metabolism-associated ALDH genes and an increase in ADH and PDC genes might result in C-PCNA persimmon PA insolubilization. In addition, a few MYB-bHLH-WD40 (MBW) homologous transcription factors in persimmon might play important roles in persimmon PA accumulation. Furthermore, combined coexpression network analysis and phylogenetic analysis of MBW suggested that three putative transcription factors WD40 (evm.TU.contig1.155), MYB (evm.TU.contig8910.486) and bHLH (evm.TU.contig1398.203), might connect and co-regulate both PA biosynthesis and its insolubilization in C-PCNA persimmon. The present study elucidated transcriptional insights into PA biosynthesis and insolubilization during the late development stages based on the C-PCNA D. kaki genome (unpublished). Thus, we focused on PA content variation and the expression patterns of genes involved in PA biosynthesis and insolubilization. Our work has provided additional evidence on previous knowledge and a basis for further exploration of the natural deastringency of C-PCNA persimmon.

Effects of Different Chemicals on Sexual Regulation in Persimmon (Diospyros kaki Thunb.) Flowers.

Research on crop sexuality is important for establishing systems for germplasm innovation and cultivating improved varieties. In this study, androecious persimmon trees were treated with various concentrations of ethrel (100, 500, and 1,000 mg/L) and zeatin (1, 5, and 10 mg/L) to investigate the morphological, physiological, and molecular characteristics of persimmon. Ethrel at 1,000 mg/L and zeatin at 10 mg/L both significantly reduced the stamen length and pollen grain diameter in androecious trees. Ethrel treatment also led to reduced stamen development with degenerated cellular contents; zeatin treatment promoted the development of arrested pistils via maintaining relatively normal mitochondrial morphology. Both treatments altered carbohydrate, amino acid, and endogenous phytohormone contents, as well as genes associated with hormone production and floral organ development. Thereafter, we explored the combined effects of four chemicals, including ethrel and zeatin, as well as zebularine and 5-azacytidine, both of which are DNA methylation inhibitors, on androecious persimmon flower development. Morphological comparisons showed that stamen length, pollen viability, and pollen grain diameter were significantly inhibited after combined treatment. Large numbers of genes involving in carbohydrate metabolic, mitogen-activated protein kinase (MAPK) signaling, and ribosome pathways, and metabolites including uridine monophosphate (UMP) and cyclamic acid were identified in response to the treatment, indicating complex regulatory mechanisms. An association analysis of transcriptomic and metabolomic data indicated that ribosomal genes have distinct effects on UMP and cyclamic acid metabolites, explaining how male floral buds of androecious persimmon trees respond to these exogenous chemicals. These findings extend the knowledge concerning sexual differentiation in persimmon; they also provide a theoretical basis for molecular breeding, high-yield cultivation, and quality improvement in persimmon.

Phytohormone and integrated mRNA and miRNA transcriptome analyses and differentiation of male between hermaphroditic floral buds of andromonoecious Diospyros kaki Thunb.

BACKGROUND: Persimmon (Diospyros kaki Thunb.) has various labile sex types, and studying its sex differentiation can improve breeding efficiency. However, studies on sexual regulation patterns in persimmon have focused mainly on monoecy and dioecy, whereas little research has been published on andromonoecy. In order to reveal the sex differentiation regulation mechanism of andromonoecious persimmon, we performed histological and cytological observations, evaluated OGI and MeGI expression and conducted phytohormones assays and mRNA and small RNA transcriptome analyses of the male and hermaphroditic floral buds of the andromonoecious persimmon 'Longyanyeshi 1'. RESULTS: Stages 2 and 4 were identified as the critical morphological periods for sex differentiation of 'Longyanyeshi 1' by histological and cytological observation. At both stages, OGI was differentially expressed in male and hermaphroditic buds, but MeGI was not. This was different from their expressions in dioecious and monoecious persimmons. Meantime, the results of phytohormones assays showed that high IAA, ABA, GA3, and JA levels at stage 2 may have promoted male floral bud differentiation. However, high JA levels at stage 4 and high ZT levels at stages 2 and 4 may have promoted hermaphroditic floral bud differentiation. In these phytohormone biosynthesis and signaling pathways, 52 and 54 differential expression genes (including Aux/IAA, ARFs, DELLA, AHP, A-ARR, B-ARR, CYP735A, CRE1, PP2C, JAZ, MYC2, COI1, CTR1, SIMKK, ACO, and MPK6) were identified, respectively. During the development of male floral buds, five metacaspases genes may have been involved in pistil abortion. In addition, MYB, FAR1, bHLH, WRKY, and MADS transcription factors might play important roles in persimmon floral bud sex differentiation. Noteworthy, miR169v_1, miR169e_3, miR319_1, and miR319 were predicted to contribute to phytohormone biosynthesis and signaling pathways and floral organogenesis and may also regulate floral bud sex differentiation. CONCLUSION: The present study revealed the differences in morphology and phytohormones content between male and hermaphroditic floral buds of 'Longyanyeshi 1' during the process of sex differentiation, and identified a subset of candidate genes and miRNAs putatively associated with its sex differentiation. These findings can provide a foundation for molecular regulatory mechanism researching on andromonoecious persimmon.

Study of Sexual-Linked Genes (OGI and MeGI) on the Performance of Androecious Persimmons (Diospyros kaki Thunb.).

It is reported that the production of floral sexual phenotype in hexaploid monoecious persimmon (Diospyros kaki) is closely related to a pseudogene called OGI, and a short interspersed nuclear element (SINE)-like insertion (named Kali) in the OGI promoter leads to the gene silence. As a result, DNA methylation level of MeGI promoter determines the development of male or female flowers. However, the molecular mechanism in androecious D. kaki, which only bear male flowers, remains elusive. Here, real-time quantitative polymerase chain reaction (RT-qPCR), molecular cloning, and bisulfite PCR sequencing technique were carried out using 87 materials, including 56 androecious resources, 15 monoecious, and 16 gynoecious cultivars, to investigate the performance of OGI and MeGI on the specific androecious type of D. kaki in China. In conclusion, the Kali insertion was exactly located in the OGI promoter region, and the OGI gene and the Kali sequence were existing and conserved in androecious D. kaki. Meanwhile, we also demonstrated that the MeGI gene was widespread in our investigated samples. Ultimately, our result convincingly provided evidence that the low expression of OGI is probably ascribed to the presence of Kali displaying strong methylation in the OGI promoter, and low expression of MeGI, as well as high DNA methylation level, in the promoter was closely connected with the production of male flowers; this result was consistent with the monoecious persimmon model. Our findings provide predominant genetic aspects for investigation into androecious D. kaki, and future perfecting the sex-determining mechanisms in persimmon.

A high-quality chromosomal genome assembly of Diospyros oleiferaCheng.

BACKGROUND: Diospyros oleifera Cheng, of the family Ebenaceae, is an economically important tree. Phylogenetic analyses indicate that D. oleifera is closely related to Diospyros kaki Thunb. and could be used as a model plant for studies of D. kaki. Therefore, development of genomic resources of D. oleifera will facilitate auxiliary assembly of the hexaploid persimmon genome and elucidate the molecular mechanisms of important traits. FINDINGS: The D. oleifera genome was assembled with 443.6 Gb of raw reads using the Pacific Bioscience Sequel and Illumina HiSeq X Ten platforms. The final draft genome was ∼812.3 Mb and had a high level of continuity with N50 of 3.36 Mb. Fifteen scaffolds corresponding to the 15 chromosomes were assembled to a final size of 721.5 Mb using 332 scaffolds, accounting for 88.81% of the genome. Repeat sequences accounted for 54.8% of the genome. By de novo sequencing and analysis of homology with other plant species, 30,530 protein-coding genes with an average transcript size of 7,105.40 bp were annotated; of these, 28,580 protein-coding genes (93.61%) had conserved functional motifs or terms. In addition, 171 candidate genes involved in tannin synthesis and deastringency in persimmon were identified; of these chalcone synthase (CHS) genes were expanded in the D. oleifera genome compared with Diospyros lotus, Camellia sinensis, and Vitis vinifera. Moreover, 186 positively selected genes were identified, including chalcone isomerase (CHI) gene, a key enzyme in the flavonoid-anthocyanin pathway. Phylogenetic tree analysis indicated that the split of D. oleifera and D. lotus likely occurred 9.0 million years ago. In addition to the ancient γ event, a second whole-genome duplication event occurred in D. oleifera and D. lotus. CONCLUSIONS: We generated a high-quality chromosome-level draft genome for D. oleifera, which will facilitate assembly of the hexaploid persimmon genome and further studies of major economic traits in the genus Diospyros.

Selection and validation of reference genes for quantitative gene expression analyses in persimmon (Diospyros kaki Thunb.) using real-time quantitative PCR.

BACKGROUND AND AIMS: Persimmon (Diospyros kaki) is an economically important fruit tree species with complex flowering characteristics. To obtain accurate expression pattern analysis results, it is vital to select a reliable gene for the normalization of real-time quantitative polymerase chain reaction data. The aim of this study was to identify the optimal internal control gene among six candidate genes for gene expression analysis in different persimmon organs and developmental stages. MATERIALS AND METHODS: This analysis was conducted using geNorm and NormFinder software to show differences in the stability of the six reference genes among tissues and floral developmental stages of the same plant. RESULTS: Although genes that exhibited moderate expression in NormFinder revealed slightly different expression stabilities than those obtained by geNorm, both sets of results showed that GAPDH was the best reference gene in different organs and floral buds at different developmental stages, whereas 18SrRNA was the least stable gene. CONCLUSIONS: Based on the overall ranking, GAPDH is the most suitable reference gene and is highly recommended for gene expression studies in different organs and different developmental stages of persimmon. This study provides useful reference data for future gene expression studies and will contribute to improving the accuracy of gene expression results in persimmon.

Proteomic Changes Between Populus Allotriploids and Diploids Revealed Using an iTRAQ-based Quantitative Approach.

BACKGROUND: Polyploid breeding is a powerful approach for Populus genetic improve-ment because polyploid trees have valuable characteristics, including better timber quality and a higher degree of stress resistance compared with their full-sib diploids. However, the genetic mech-anism underlying this phenomenon remains unknown. OBJECTIVE: To better understand the proteomic changes between Populus allotriploids and diploids, we examined the proteomic profiles of allotriploid and diploid Populus by iTRAQ labeling coupled with two-dimensional liquid chromatography and MALDI-TOF/TOF mass spectrometry (MS). METHOD: iTRAQ labeling coupled with two-dimensional liquid chromatography and MALDI-TOF/TOF mass spectrometry (MS). RESULTS: Between the Populus allotriploid and the full-sib diploid, 932 differentially expressed proteins (DEPs) were identified. These DEPs were primarily involved in stress, defense, transportation, transcriptional and/or translational modification, and energy production. The pathway analysis indi-cated that most of the DEPs were implicated in carbohydrate transport and metabolism, nitrogen me-tabolism and glycolysis, and the ribosome assembly pathway. These data suggest high protein di-vergence between Populus allotriploids and diploids, and rapid changes during hybridization. CONCLUSION: The results provide new data for further understanding of the mechanisms of polyploid trees that generally display increased height growth compared with their full-sib diploids.

MicroRNA expression changes following synthesis of three full-sib Populus triploid populations with different heterozygosities.

KEY MESSAGE: Through high-throughput sequencing, we compared the relative expression levels of miRNA in three full-sib Populus triploid populations with that in their parents and one diploid hybrid population. We found similar numbers of miRNAs differentially expressed between the parents and the four progeny hybrid populations. In addition, unbalanced parental expression level dominance of miRNAs were found in the three allotriploid and interspecific hybrid populations, which may reprogram gene expression networks and contribute to the growth of Populus hybrids. These results indicated that hybridization has a great impact on the miRNA expression variation in the newly synthesized Populus triploid and diploid hybrid populations. However, we also found no significant differences in miRNA expression among one diploid and three triploid hybrid populations, hinting that miRNA abundances do not increase with the genome content. No dosage effect of miRNA expression could lead to dosage-dependent negative effects on target genes and their downstream pathway in polyploids. We speculate that polyploids may gain advantages from the slight decrease in miRNA regulation, suggesting an important molecular mechanism of polyploid advantage. Hybridization with three types of induced 2n gametes transmitted different parental heterozygosities has been proven as an efficient method for Populus triploid production. Several researches have shown that miRNA could be non-additively expressed in allopolyploids. However, it is still unclear whether the non-additively expressed miRNAs result from the effect of hybridization or polyploidization, and whether a dose response to the additional genomic content exists for the expression of miRNA. Toward this end, through high-throughput sequencing, we compared the expression levels of miRNA in three full-sib Populus triploid populations with that in their parents and one interspecific hybrid population. We found similar numbers of miRNAs differentially expressed between the parents and the four progeny hybrid populations. Unbalanced parental expression level dominance of miRNAs were found in the three triploid and diploid hybrid populations, which may reprogram gene expression networks and affect the growth of Populus hybrids. These results indicated that hybridization has a great impact on the miRNA expression variation in the newly synthesized Populus triploid and diploid hybrid populations. However, we also found no significant differences in miRNA expression among the three triploid populations and the diploid hybrid population. No dosage effect of miRNA expression could lead to dosage-dependent negative effects on target genes and their downstream pathway in polyploids. We speculate that polyploids may gain advantages from the decrease in miRNA negative regulation, suggesting an important molecular mechanism of polyploid advantage.

Transcriptomic changes following synthesis of a Populus full-sib diploid and allotriploid population with different heterozygosities driven by three types of 2n female gamete.

Diploid gametes are usually applied to produce triploids of Populus [originating from first-division restitution (FDR), second-division restitution (SDR), and postmeiotic restitution (PMR) 2n eggs]. Three types of 2n gametes transmitted different parental heterozygosities in Populus. Failed spindle formation and no chromosomal separation to opposite poles during meiosis I mean that FDR 2n gametes carry nonsister chromatids that are potentially heterozygous. By contrast, SDR 2n gametes result from failed sister chromatid separation in meiosis II, and therefore, they carry sister chromatid that are potentially homozygous. Completely homozygous 2n gametes can arise from the PMR mechanism. The alteration of gene expression resulting from allopolyploidization is a prominent feature in plants. We compared gene expression in the full-sib progeny of three allotriploid Populus populations (triploid-F, triploid-S, and triploid-P) with that in its parent species, and their full-sib diploid F1 hybrid. Genome-wide expression level dominance was biased toward the maternal in the diploid F1 hybrid and three allotriploid populations, whereas our data indicated important, but different, effects of the transmission of different heterozygosity by 2n female gametes in the expression patterns of allopolyploids. Because of the higher level of heterozygosity, the triploids had higher rates of non-additive and transgressive expression patterns in the triploid-F than in triploid-S and triploid-P. Compared with diploid F1, about 30-fold more genes (251) were differently expressed in the triploid-F than in the triploid-S (9) and triploid-P (8), respectively. These findings indicate that hybridization and polyploidization have immediate and distinct effects on the large-scale patterns of gene expression, and different effects on the transmission of heterozygosity by three 2n female gametes.

Inheritance and variation of Cytosine methylation in three populus allotriploid populations with different heterozygosity.

DNA methylation is an epigenetic mechanism with the potential to regulate gene expression and affect plant phenotypes. Both hybridization and genome doubling may affect the DNA methylation status of newly formed allopolyploid plants. Previous studies demonstrated that changes in cytosine methylation levels and patterns were different among individual hybrid plant, therefore, studies investigating the characteristics of variation in cytosine methylation status must be conducted at the population level to avoid sampling error. In the present study, an F1 hybrid diploid population and three allotriploid populations with different heterozygosity [originating from first-division restitution (FDR), second-division restitution (SDR), and post-meiotic restitution (PMR) 2n eggs of the same female parent] were used to investigate cytosine methylation inheritance and variation relative to their common parents using methylation-sensitive amplification polymorphism (MSAP). The variation in cytosine methylation in individuals in each population exhibited substantial differences, confirming the necessity of population epigenetics. The total methylation levels of the diploid population were significantly higher than in the parents, but those of the three allotriploid populations were significantly lower than in the parents, indicating that both hybridization and polyploidization contributed to cytosine methylation variation. The vast majority of methylated status could be inherited from the parents, and the average percentages of non-additive variation were 6.29, 3.27, 5.49 and 5.07% in the diploid, FDR, SDR and PMR progeny populations, respectively. This study lays a foundation for further research on population epigenetics in allopolyploids.

Differential transcriptome analysis between Populus and its synthesized allotriploids driven by second-division restitution.

In this report, we compared transcriptomic differences between a synthetic Populus section Tacamahaca triploid driven by second-division restitution and its parents using a high-throughput RNA-seq method. A total of 4,080 genes were differentially expressed between the high-growth vigor allotriploids (SDR-H) and their parents, and 719 genes were non-additively expressed in SDR-H. Differences in gene expression between the allotriploid and male parent were more significant than those between the allotriploid and female parent, which may be caused by maternal effects. We observed 3,559 differentially expressed genes (DEGs) between the SDR-H and male parent. Notably, the genes were mainly involved in metabolic process, cell proliferation, DNA methylation, cell division, and meristem and developmental growth. Among the 1,056 DEGs between SDR-H and female parent, many genes were associated with metabolic process and carbon utilization. In addition, 1,789 DEGs between high- and low-growth vigor allotriploid were mainly associated with metabolic process, auxin poplar transport, and regulation of meristem growth. Our results indicated that the higher poplar ploidy level can generate extensive transcriptomic diversity compared with its parents. Overall, these results increased our understanding of the driving force for phenotypic variation and adaptation in allopolyploids driven by second-division restitution.

A strategy for characterization of persistent heteroduplex DNA in higher plants.

Heteroduplex DNA (hDNA) generated during homologous recombination (HR) is an important component that shapes genetic diversity in sexually reproducing organisms. However, studies of this process in higher plants are limited. This is because hDNAs are difficult to capture in higher plants as their reproductive developmental model only produces normal gametes and does not preserve the mitotic products of the post-meiotic segregation (PMS) process which is crucial for studying hDNAs. In this study, using the model system for tree and woody perennial plant biology (Populus), we propose a strategy for characterizing hDNAs in higher plants. We captured hDNAs by constructing triploid hybrids originating from a cross between unreduced 2n eggs (containing hDNA information as a result of inhibition chromosome segregation at the PMS stage) with normal male gametes. These triploid hybrids allowed us to detect the frequency and location of persistent hDNAs resulting from HR at the molecular level. We found that the frequency of persistent hDNAs, which ranged from 5.3 to 76.6%, was related to locations of the simple sequence repeat markers at the chromosomes, such as the locus-centromere distance, the surrounding DNA sequence and epigenetic information, and the richness of protein-coding transcripts at these loci. In summary, this study provides a method for characterizing persistent hDNAs in higher plants. When high-throughput sequencing techniques can be incorporated, genome-wide persistent hDNA assays for higher plants can be easily carried out using the strategy presented in this study.