Comparative transcriptome and functional analyses provide insights into the key factors regulating shoot regeneration in highbush blueberry.

Establishing an efficient plant regeneration system is a crucial prerequisite for genetic engineering technology in plants. However, the regeneration rate exhibits considerable variability among genotypes, and the key factors underlying shoot regeneration capacity remain largely elusive. Blueberry leaf explants cultured on a medium rich in cytokinins exhibit direct shoot organogenesis without prominent callus formation, which holds promise for expediting genetic transformation while minimizing somatic mutations during culture. The objective of this study is to unravel the molecular and genetic determinants that govern cultivar-specific shoot regeneration potential in highbush blueberry (Vaccinium corymbosum L.). We conducted comparative transcriptome analysis using two highbush blueberry genotypes: 'Blue Muffin' ('BM') displaying a high regeneration rate (>80%) and 'O'Neal' ('ON') exhibiting a low regeneration rate (<10%). The findings revealed differential expression of numerous auxin-related genes; notably, 'BM' exhibited higher expression of auxin signaling genes compared to 'ON'. Among blueberry orthologs of transcription factors involved in meristem formation in Arabidopsis, expression of VcENHANCER OF SHOOT REGENERATION (VcESR), VcWUSCHEL (VcWUS), and VcCUP-SHAPED COTYLEDON 2.1 were significantly higher in 'BM' relative to 'ON'. Exogenous application of auxin promoted regeneration, as well as VcESR and VcWUS expression, whereas inhibition of auxin biosynthesis yielded the opposite effects. Overexpression of VcESR in 'BM' promoted shoot regeneration under phytohormone-free conditions by activating the expression of cytokinin- and auxin-related genes. These findings provide new insights into the molecular mechanisms underlying blueberry regeneration and have practical implications for enhancing plant regeneration and transformation techniques.

Regulatory role of Prunus mume DAM6 on lipid body accumulation and phytohormone metabolism in the dormant vegetative meristem.

Bud dormancy is a crucial process in the annual growth cycle of woody perennials. In Rosaceae fruit tree species, DORMANCY-ASSOCIATED MADS-box (DAM) transcription factor genes regulating bud dormancy have been identified, but their molecular roles in meristematic tissues have not been thoroughly characterized. In this study, molecular and physiological analyses of transgenic apple plants overexpressing the Japanese apricot DAM6 gene (PmDAM6) and Japanese apricot cultivars and F1 individuals with contrasting dormancy characteristics revealed the metabolic pathways controlled by PmDAM6. Our transcriptome analysis and transmission electron microscopy examination demonstrated that PmDAM6 promotes the accumulation of lipid bodies and inhibits cell division in the dormant vegetative meristem by down-regulating the expression of lipid catabolism genes (GDSL ESTERASE/LIPASE and OIL BODY LIPASE) and CYCLIN genes, respectively. Our findings also indicate PmDAM6 promotes abscisic acid (ABA) accumulation and decreases cytokinin (CTK) accumulation in vegetative buds by up-regulating the expression of the ABA biosynthesis gene ARABIDOPSIS ALDEHYDE OXIDASE and the CTK catabolism gene CYTOKININ DEHYDROGENASE, while also down-regulating the expression of the CTK biosynthesis genes ISOPENTENYL TRANSFERASE (IPT) and CYP735A. Additionally, PmDAM6 modulates gibberellin (GA) metabolism by up-regulating GA2-OXIDASE expression and down-regulating GA3-OXIDASE expression. Furthermore, PmDAM6 may also indirectly promote lipid accumulation and restrict cell division by limiting the accumulation of CTK and GA in buds. In conclusion, using our valuable genetic platform, we clarified how PmDAM6 modifies diverse cellular processes, including lipid catabolism, phytohormone (ABA, CTK, and GA) biosynthesis and catabolism, and cell division, in the dormant vegetative meristem.

Three-dimensional fruit growth analysis clarifies developmental mechanisms underlying complex shape diversity in persimmon fruit.

The determination of fruit size and shape are of considerable interest in horticulture and developmental biology. Fruit typically exhibits three-dimensional structures characterized by geometric features that are dependent on the genotype. Although minor developmental variations have been recognized, few studies have fully visualized and measured these variations throughout fruit growth. Here, a high-resolution 3D scanner was used to investigate the fruit development of 51 persimmon (Diospyros kaki) cultivars with various complex shapes. We obtained 2380 3D models that fully represented fruit appearance, and enabled precise and automated measurements of shape features throughout fruit development, including horizontal and vertical grooves, length-to-width ratio, and roundness. The 3D fruit model analysis identified key stages that determined the shape attributes at maturity. Typically, genetic diversity was found in vertical groove development, and these grooves could be filled by tissue expansion in the carpel fusion zone during fruit development. In addition, transcriptome analysis of fruit tissues from groove and non-groove tissues revealed gene co-expression networks that were highly associated with groove depth variation. The presence of YABBY homologs was most closely associated with groove depth and indicated the possibility that this pathway is a key molecular contributor to vertical groove depth variation. Overall, our results revealed deterministic patterns of complex shape traits in persimmon fruit and showed that different growth patterns among tissues are the main factor contributing to the shape of both vertical and horizontal grooves.

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.

Genetic basis of lineage-specific evolution of fruit traits in hexaploid persimmon.

Frequent polyploidization events in plants have led to the establishment of many lineage-specific traits representing each species. Little is known about the genetic bases for these specific traits in polyploids, presumably due to plant genomic complexity and their difficulties in applying genetic approaches. Hexaploid Oriental persimmon (Diospyros kaki) has evolved specific fruit characteristics, including wide variations in fruit shapes and astringency. In this study, using whole-genome diploidized/quantitative genotypes from ddRAD-Seq data of 173 persimmon cultivars, we examined their population structures and potential correlations between their structural transitions and variations in nine fruit traits. The population structures of persimmon cultivars were highly randomized and not substantially correlated with the representative fruit traits focused on in this study, except for fruit astringency. With genome-wide association analytic tools considering polyploid alleles, we identified the loci associated with the nine fruit traits; we mainly focused on fruit-shape variations, which have been numerically characterized by principal component analysis of elliptic Fourier descriptors. The genomic regions that putatively underwent selective sweep exhibited no overlap with the loci associated with these persimmon-specific fruit traits. These insights will contribute to understanding the genetic mechanisms by which fruit traits are independently established, possibly due to polyploidization events.

Integer programming for selecting set of informative markers in paternity inference.

BACKGROUND: Parentage information is fundamental to various life sciences. Recent advances in sequencing technologies have made it possible to accurately infer parentage even in non-model species. The optimization of sets of genome-wide markers is valuable for cost-effective applications but requires extremely large amounts of computation, which presses for the development of new efficient algorithms. RESULTS: Here, for a closed half-sib population, we generalized the process of marker loci selection as a binary integer programming problem. The proposed systematic formulation considered marker localization and the family structure of the potential parental population, resulting in an accurate assignment with a small set of markers. We also proposed an efficient heuristic approach, which effectively improved the number of markers, localization, and tolerance to missing data of the set. Applying this method to the actual genotypes of apple (Malus × domestica) germplasm, we identified a set of 34 SNP markers that distinguished 300 potential parents crossed to a particular cultivar with a greater than 99% accuracy. CONCLUSIONS: We present a novel approach for selecting informative markers based on binary integer programming. Since the data generated by high-throughput sequencing technology far exceeds the requirement for parentage assignment, a combination of the systematic marker selection with targeted SNP genotyping, such as KASP, allows flexibly enlarging the analysis up to a scale that has been unrealistic in various species. The method developed in this study can be directly applied to unsolved large-scale problems in breeding, reproduction, and ecological research, and is expected to lead to novel knowledge in various biological fields. The implementation is available at https://github.com/SoNishiyama/IP-SIMPAT .

H3K4me3 plays a key role in establishing permissive chromatin states during bud dormancy and bud break in apple.

Bud dormancy helps woody perennials survive winter and activate robust plant development in the spring. For apple (Malus × domestica), short-term chilling induces bud dormancy in autumn, then prolonged chilling leads to dormancy release and a shift to a quiescent state in winter, with subsequent warm periods promoting bud break in spring. Epigenetic regulation contributes to seasonal responses such as vernalization. However, how histone modifications integrate seasonal cues and internal signals during bud dormancy in woody perennials remains largely unknown. Here, we show that H3K4me3 plays a key role in establishing permissive chromatin states during bud dormancy and bud break in apple. The global changes in gene expression strongly correlated with changes in H3K4me3, but not H3K27me3. High expression of DORMANCY-ASSOCIATED MADS-box (DAM) genes, key regulators of dormancy, in autumn was associated with high H3K4me3 levels. In addition, known DAM/SHORT VEGETATIVE PHASE (SVP) target genes significantly overlapped with H3K4me3-modified genes as bud dormancy progressed. These data suggest that H3K4me3 contributes to the central dormancy circuit, consisting of DAM/SVP and abscisic acid (ABA), in autumn. In winter, the lower expression and H3K4me3 levels at DAMs and gibberellin metabolism genes control chilling-induced release of dormancy. Warming conditions in spring facilitate the expression of genes related to phytohormones, the cell cycle, and cell wall modification by increasing H3K4me3 toward bud break. Our study also revealed that activation of auxin and repression of ABA sensitivity in spring are conditioned at least partly through temperature-mediated epigenetic regulation in winter.

Novel insights into the dissemination route of Japanese apricot (Prunus mume Sieb. et Zucc.) based on genomics.

Japanese apricot (Prunus mume) is an attractive fruit tree originating from China, and its cultivation history dates back 7000 years. In this study, we investigated the genetic diversity, population structure, and genetic relationship of Japanese apricots in different regions of China and Japan. The analyses of the genetic variation between wild and cultivated populations improved our understanding of the general mechanisms of domestication and improvement. A total of 146 accessions of Japanese apricot from different geographic locations were sequenced. The genetic diversity of wild and domesticated accessions (3.60 × 10-3 and 3.51 × 10-3 , respectively) from China was high, and the effect of artificial selection pressure on domesticated accessions was small; however, the genetic diversity of artificially bred accessions decreased significantly (2.68 × 10-3 ) compared to domesticated accessions, which had an obvious improvement bottleneck effect. The chloroplast genome results showed that 41 haplotypes were detected, and Japanese apricots from the Yunnan region had the most haplotypes and the highest genetic diversity. The results revealed the dissemination route of Japanese apricot, not only along the Yangtze River basin system (from southwest China to Hunan, Jiangxi, and Anhui, and finally to the Jiangsu, Zhejiang, and Shanghai areas). Additionally, we discovered a second route for Japanese apricot dispersion, which was mostly in the Pearl River basin system, from southwest China to Libo of Guizhou and then to the Guangdong, Fujian, and Taiwan areas. This also showed that Japanese-bred accessions originated from Zhejiang, China. In addition, selective sweep analysis showed that most of the high-impact single nucleotide polymorphisms were identified in genes related to glucose metabolism, aromatic compound metabolism, flowering time, dormancy, and resistance to abiotic stress during the domestication and improvement of Japanese apricot.

Reinvention of hermaphroditism via activation of a RADIALIS-like gene in hexaploid persimmon.

In flowering plants, different lineages have independently transitioned from the ancestral hermaphroditic state into and out of various sexual systems1. Polyploidizations are often associated with this plasticity in sexual systems2,3. Persimmons (the genus Diospyros) have evolved dioecy via lineage-specific palaeoploidizations. More recently, hexaploid D. kaki has established monoecy and also exhibits reversions from male to hermaphrodite flowers in response to natural environmental signals (natural hermaphroditism, NH), or to artificial cytokinin treatment (artificial hermaphroditism, AH). We sought to identify the molecular pathways underlying these polyploid-specific reversions to hermaphroditism. Co-expression network analyses identified regulatory pathways specific to NH or AH transitions. Surprisingly, the two pathways appeared to be antagonistic, with abscisic acid and cytokinin signalling for NH and AH, respectively. Among the genes common to both pathways leading to hermaphroditic flowers, we identified a small-Myb RADIALIS-like gene, named DkRAD, which is specifically activated in hexaploid D. kaki. Consistently, ectopic overexpression of DkRAD in two model plants resulted in hypergrowth of the gynoecium. These results suggest that production of hermaphrodite flowers via polyploidization depends on DkRAD activation, which is not associated with a loss-of-function within the existing sex determination pathway, but rather represents a new path to (or reinvention of) hermaphroditism.

Genomics-based discrimination of 2n gamete formation mechanisms in polyploids: a case study in nonaploid Diospyros kaki 'Akiou'.

Unreduced gametes (2n gametes), possessing double the haploid genome, whatever ploidy that happens to be, are a common source of ploidy variation in plant populations. First and second division restitution (FDR and SDR) are the dominant mechanisms of 2n gamete production; all else being equal, FDR gametes have a higher degree of heterozygosity, thus they are advantageous in breeding. The discrimination of these mechanisms from the consequence of hybridization is challenging, especially in higher polyploids, and usually requires information on centromere location. In this study, we propose a genotyping-based strategy to uncover the mechanisms of 2n gamete formation in progeny that has a higher ploidy than its parents. Simulation of 2n gamete production revealed that FDR and SDR pathways can be discriminated based on allele transmission patterns alone without information on centromere location. We applied this strategy to study the formation mechanism of a nonaploid Diospyros kaki 'Akiou', which was bred via hybridization between D. kaki hexaploid cultivars. The result demonstrated that 'Akiou' was derived from the fertilization of a normal female gamete by a 2n male gamete and that this 2n gamete was produced through FDR. Consequently, the distinct duplex transmission pattern in the FDR gamete enabled us to infer the genomic characteristics of polyploid persimmon. The method could be tested only for the plant being polypoid, which allows for the ability to discriminate causes of 2n gamete formation using allele dosage in progeny, and will be useful in future studies of polyploid genomics.

Preharvest long-term exposure to UV-B radiation promotes fruit ripening and modifies stage-specific anthocyanin metabolism in highbush blueberry.

Ultraviolet-B (UV-B) light (280-315 nm) is an important environmental signal that regulates plant development and photomorphogenesis, while also affecting the flavonoid pathway, including anthocyanin biosynthesis. Regarding the effects of UV-B radiation on fruits, the effects of a short-term or postharvest irradiation on fruit quality have been well-documented, but the effects of a long-term preharvest UV-B irradiation on fruit growth and coloration remain unclear. Thus, in this study, we investigated the effects of a long-term treatment involving an environmentally relevant UV-B dose on highbush blueberry (Vaccinium corymbosum) fruit. The preharvest UV-B treatment quickly promoted fruit growth and sugar accumulation, which is not commonly observed in other fruit tree species. The UV-B exposure also accelerated fruit ripening and coloration. The dual-luciferase assay proved that in blueberries, expression of VcUFGT encoding anthocyanin biosynthesis key enzyme, is positively and negatively regulated by VcMYBA1 and VcMYBC2, respectively. Throughout the fruit development stage, the UV-B treatment up-regulated VcMYBPA1 expression, which increased VcUFGT expression via VcMYBA1. In the green fruit stage, the UV-B treatment increased HY5 encoding UV receptor, which up-regulates VcMYBPA1 and down-regulates VcMYBC2, thereby promotes the accumulation of anthocyanins. On the other hand, excessive anthocyanin synthesis was inhibited by increased VcMYBC2 levels in mature fruits when exposed to UV-B light through HY5-independent pathway. In conclusion, anthocyanin-related MYB activators and repressor may coordinately balance the accumulation of anthocyanins in blueberry fruits, with UV-B treatments possibly influencing their effects in a stage-specific manner. The potential utility of preharvest UV-B treatments for improving blueberry fruit quality is discussed herein.

Quantitative analysis of auxin metabolites in lychee flowers.

To investigate the modulation of endogenous indole-3-acetic acid (IAA) level by biosynthesis and inactivation during floral development, IAA and its metabolites were analyzed by LC-ESI/MS/MS in Lychee (Litchi chinensis Sonn.) flowers. In the bloomed flowers, the level of free IAA was higher in males than in females. In contrast, the total sum level of IAA metabolites was higher in females than in males, suggesting a higher biosynthetic activity of IAA in the females before the bloom. A detailed time-course analysis from the bud stage to the developing flower stage showed higher levels of IAA in females than males. The major metabolites were oxidized IAA in both sexes. The results suggest that IAA is involved in the maturation of female floral tissues in lychee, and oxidative metabolism plays an essential role in controlling the free IAA levels therein.

Functional and expressional analyses of apple FLC-like in relation to dormancy progress and flower bud development.

We previously identified the FLOWERING LOCUS C (FLC)-like gene, a MADS-box transcription factor gene that belongs to Arabidopsis thaliana L. FLC clade, in apple (Malus $\times$  domestica Borkh.), and its expression in dormant flower buds is positively correlated with cumulative cold exposure. To elucidate the role of the MdFLC-like in the dormancy process and flower development, we first characterized the phenotypes of MdFLC-like overexpressing lines with the Arabidopsis Columbia-0 background. The overexpression of MdFLC-like significantly delayed the bolting date and reduced the plant size, but it did not significantly affect the number of rosette leaves or flower organ formation. Thus, MdFLC-like may affect vegetative growth and development rather than flowering when expressed in Arabidopsis, which is not like Arabidopsis FLC that affects development of flowering. We compared seasonal expression patterns of MdFLC-like in low-chill 'Anna' and high-chill 'Fuji' and 'Tsugaru' apples collected from trees grown in a cold winter region in temperate zone and found an earlier upregulation in 'Anna' compared with 'Fuji' and 'Tsugaru'. Expression patterns were also compared in relation to developmental changes in the flower primordia during the chilling accumulation period. Overall, MdFLC-like was progressively upregulated during flower primordia differentiation and development in autumn to early winter and reached a maximum expression level at around the same time as the genotype-dependent chilling requirements were fulfilled in high-chill cultivars. Thus, we hypothesize MdFLC-like may be upregulated in response to cold exposure and flower primordia development during the progress of endodormancy. Our study also suggests MdFLC-like may have a growth-inhibiting function during the end of endodormancy and ecodormancy when the temperature is low and unfavorable for rapid bud outgrowth.

Genomic insight into the developmental history of southern highbush blueberry populations.

Interspecific hybridization is a common breeding approach for introducing novel traits and genetic diversity to breeding populations. Southern highbush blueberry (SHB) is a blueberry cultivar group that has been intensively bred over the last 60 years. Specifically, it was developed by multiple interspecific crosses between northern highbush blueberry [NHB, Vaccinium corymbosum L. (2n = 4x = 48)] and low-chill Vaccinium species to expand the geographic limits of highbush blueberry production. In this study, we genotyped polyploid blueberries, including 105 SHB, 17 NHB, and 10 rabbiteye blueberry (RE) (Vaccinium virgatum Aiton), from the accessions planted at Poplarville, Mississippi, and accessions distributed in Japan, based on the double-digest restriction site-associated DNA sequencing. The genome-wide SNP data clearly indicated that RE cultivars were genetically distinct from SHB and NHB cultivars, whereas NHB and SHB were genetically indistinguishable. The population structure results appeared to reflect the differences in the allele selection strategies that breeders used for developing germplasm adapted to local climates. The genotype data implied that there are no or very few genomic segments that were commonly introgressed from low-chill Vaccinium species to the SHB genome. Principal component analysis-based outlier detection analysis found a few loci associated with a variable that could partially differentiate NHB and SHB. These SNP loci were detected in Mb-scale haplotype blocks and may be close to the functional genes related to SHB development. Collectively, the data generated in this study suggest a polygenic adaptation of SHB to the southern climate, and may be relevant for future population-scale genome-wide analyses of blueberry.

Genome-Wide Identification of Loci Associated With Phenology-Related Traits and Their Adaptive Variations in a Highbush Blueberry Collection.

Genetic variation in phenological traits is the key in expanding production areas of crops. Southern highbush blueberry (SHB) is a blueberry cultivar group adapted to warmer climates and has been developed by multiple interspecific hybridizations between elite northern highbush blueberry (NHB) (Vaccinium corymbosum L.) and low-chill Vaccinium species native to the southern United States. In this study, we employed a collection of diverse SHB accessions and performed a genome-wide association study (GWAS) for five phenology-related traits [chilling requirement (CR), flowering date, ripening date, fruit development period, and continuous flowering] using polyploid GWAS models. Phenology-related traits showed higher heritability and larger correlation coefficients between year replications, which resulted in the detection of robust phenotype-genotype association peaks. Notably, a single association peak for the CR was detected on Chromosome 4. Comparison of genotypes at the GWAS peaks between NHB and SHB revealed the putative introgression of low-chill and late-flowering alleles into the highbush genetic pool. Our results provide basic insights into the diversity of phenological traits in blueberry and the genetic establishment of current highbush cultivar groups.

Genome-wide study on the polysomic genetic factors conferring plasticity of flower sexuality in hexaploid persimmon.

Sexuality is one of the fundamental mechanisms that work towards maintaining genetic diversity within a species. In diploid persimmons (Diospyros spp.), separated sexuality, the presence of separate male and female individuals (dioecy), is controlled by the Y chromosome-encoded small-RNA gene, OGI. On the other hand, sexuality in hexaploid Oriental persimmon (Diospyros kaki) is more plastic, with OGI-bearing genetically male individuals, able to produce both male and female flowers (monoecy). This is thought to be linked to the partial inactivation of OGI by a retrotransposon insertion, resulting in DNA methylation of the OGI promoter region. To identify the genetic factors regulating branch sexual conversion, genome-wide correlation/association analyses were conducted using ddRAD-Seq data from an F1 segregating population, and using both quantitative and diploidized genotypes, respectively. We found that allelic ratio at the Y-chromosomal region, including OGI, was correlated with male conversion based on quantitative genotypes, suggesting that OGI can be activated in cis in a dosage-dependent manner. Genome-wide association analysis based on diploidized genotypes, normalized for the effect of OGI allele dosage, detected three fundamental loci associated with male conversion. These loci underlie candidate genes, which could potentially act epigenetically for the activation of OGI expression.

Correction: The persimmon genome reveals clues to the evolution of a lineage-specific sex determination system in plants.

[This corrects the article DOI: 10.1371/journal.pgen.1008566.].

The persimmon genome reveals clues to the evolution of a lineage-specific sex determination system in plants.

Most angiosperms bear hermaphroditic flowers, but a few species have evolved outcrossing strategies, such as dioecy, the presence of separate male and female individuals. We previously investigated the mechanisms underlying dioecy in diploid persimmon (D. lotus) and found that male flowers are specified by repression of the autosomal gene MeGI by its paralog, the Y-encoded pseudo-gene OGI. This mechanism is thought to be lineage-specific, but its evolutionary path remains unknown. Here, we developed a full draft of the diploid persimmon genome (D. lotus), which revealed a lineage-specific whole-genome duplication event and provided information on the architecture of the Y chromosome. We also identified three paralogs, MeGI, OGI and newly identified Sister of MeGI (SiMeGI). Evolutionary analysis suggested that MeGI underwent adaptive evolution after the whole-genome duplication event. Transformation of tobacco plants with MeGI and SiMeGI revealed that MeGI specifically acquired a new function as a repressor of male organ development, while SiMeGI presumably maintained the original function. Later, a segmental duplication event spawned MeGI's regulator OGI on the Y-chromosome, completing the path leading to dioecy, and probably initiating the formation of the Y-chromosome. These findings exemplify how duplication events can provide flexible genetic material available to help respond to varying environments and provide interesting parallels for our understanding of the mechanisms underlying the transition into dieocy in plants.

Molecular Mechanism Underlying Derepressed Male Production in Hexaploid Persimmon.

Sex expression in plants is often flexible and contributes to the maintenance of genetic diversity within a species. In diploid persimmons (the genus Diospyros), the sexuality is controlled by the Y chromosome-encoded small-RNA gene, OGI, and its autosomal counterpart, MeGI. Hexaploid Oriental persimmon (Diospyros kaki) evolved more flexible sex expression, where genetically male individuals carrying OGI can produce both male and female flowers (monoecy). This is due to (semi-)inactivation of OGI by the Kali-SINE retrotransposon insertion on the promoter region and the resultant DNA methylations. Instead, flower sex determination in Oriental persimmon is also dependent on DNA methylation states of MeGI. Here, we focused on a cultivar, Kumemaru, which shows stable male flower production. Our results demonstrated that cv. Kumemaru carries OGI with Kali-SINE, which was highly methylated as well as in other monoecious cultivars; nevertheless, OGI gene could have a basal expression level. Transcriptomic analysis between cv. Kumemaru and 14 cultivars that predominantly produce female flowers showed differentially expressed genes (DEGs) specific to cv. Kumemaru, which is mainly involved in stress responses. Co-expression gene networks focusing on the DEGs also suggested the involvement of stress signals, mainly via gibberellin (GA), salicylic acid (SA), and especially jasmonic acid (JA) signal pathways. We also identified potential regulators of this co-expression module, represented by the TCP4 transcription factor. Furthermore, we attempted to identify cv. Kumemaru-specific transcript polymorphisms potentially contributing to derepressed OGI expression by cataloging subsequences (k-mers) in the transcriptomic reads from cv. Kumemaru and the other 14 female cultivars. Overall, although the direct genetic factor to activate OGI remains to be solved, our results implied the involvement of stress signals in the release of silenced OGI and the resultant continuous male production.