Genome-wide cis-decoding for expression design in tomato using cistrome data and explainable deep learning.

In the evolutionary history of plants, variation in cis-regulatory elements (CREs) resulting in diversification of gene expression has played a central role in driving the evolution of lineage-specific traits. However, it is difficult to predict expression behaviors from CRE patterns to properly harness them, mainly because the biological processes are complex. In this study, we used cistrome datasets and explainable convolutional neural network (CNN) frameworks to predict genome-wide expression patterns in tomato (Solanum lycopersicum) fruit from the DNA sequences in gene regulatory regions. By fixing the effects of trans-acting factors using single cell-type spatiotemporal transcriptome data for the response variables, we developed a prediction model for crucial expression patterns in the initiation of tomato fruit ripening. Feature visualization of the CNNs identified nucleotide residues critical to the objective expression pattern in each gene, and their effects were validated experimentally in ripening tomato fruit. This cis-decoding framework will not only contribute to the understanding of the regulatory networks derived from CREs and transcription factor interactions, but also provides a flexible means of designing alleles for optimized expression.

Distinct Functions of Ethylene and ACC in the Basal Land Plant Marchantia polymorpha.

Ethylene is a gaseous phytohormone involved in various physiological processes, including fruit ripening, senescence, root hair development and stress responses. Recent genomics studies have suggested that most homologous genes of ethylene biosynthesis and signaling are conserved from algae to angiosperms, whereas the function and biosynthesis of ethylene remain unknown in basal plants. Here, we examined the physiological effects of ethylene, an ethylene precursor, 1-aminocyclopropane-1-carboxylic acid (ACC) and an inhibitor of ethylene perception, silver thiosulfate (STS), in a basal land plant, Marchantia polymorpha. M. polymorpha plants biosynthesized ethylene, and treatment with high concentrations of ACC slightly promoted ethylene production. ACC remarkably suppressed the growth of thalli (vegetative organs) and rhizoids (root-hair-like cells), whereas exogenous ethylene slightly promoted thallus growth. STS suppressed thallus growth and induced ectopic rhizoid formation on the dorsal surface of thalli. Thus, ACC and ethylene have different effects on the vegetative growth of M. polymorpha. We generated single and double mutants of ACC synthase-like (ACSL) genes, MpACSL1 and MpACSL2. The mutants did not show obvious defects in thallus growth, ACC content and ethylene production, indicating that MpACSL genes are not essential for the vegetative growth and biosynthesis of ACC and ethylene. Gene expression analysis suggested the involvement of MpACSL1 and MpACSL2 in stress responses. Collectively, our results imply ethylene-independent function of ACC and the absence of ACC-mediated ethylene biosynthesis in M. polymorpha.

Physiological Role of Aerobic Fermentation Constitutively Expressed in an Aluminum-Tolerant Cell Line of Tobacco (Nicotiana tabacum).

Aluminum (Al)-tolerant tobacco cell line ALT301 derived from SL (wild-type) hardly exhibits Al-triggered reactive oxygen species (ROS) compared with SL. Molecular mechanism leading to this phenotype was investigated comparatively with SL. Under normal growth condition, metabolome data suggested the activation of glycolysis and lactate fermentation but the repression of the tricarboxylic acid (TCA) cycle in ALT301, namely aerobic fermentation, which seemed to be transcriptionally controlled partly by higher expression of genes encoding lactate dehydrogenase and pyruvate dehydrogenase kinase. Microarray and gene ontology analyses revealed the upregulation of the gene encoding related to APETALA2.3 (RAP2.3)-like protein, one of the group VII ethylene response factors (ERFVIIs), in ALT301. ERFVII transcription factors are known to be key regulators for hypoxia response that promotes substrate-level ATP production by glycolysis and fermentation. ERFVIIs are degraded under normoxia by the N-end rule pathway of proteolysis depending on both oxygen and nitric oxide (NO), and NO is produced mainly by nitrate reductase (NR) in plants. In ALT301, levels of the NR gene expression (NIA2), NR activity and NO production were all lower compared with SL. Consistently, the known effects of NO on respiratory pathways were also repressed in ALT301. Under Al-treatment condition, NO level increased in both lines but was lower in ALT301. These results suggest that the upregulation of the RAP2.3-like gene and the downregulation of the NIA2 gene and resultant NO depletion in ALT301 coordinately enhance aerobic fermentation, which seems to be related to a higher capacity to prevent ROS production in mitochondria under Al stress.

Low temperature modulates natural peel degreening in lemon fruit independently of endogenous ethylene.

Peel degreening is an important aspect of fruit ripening in many citrus fruit, and previous studies have shown that it can be advanced by ethylene treatment or by low-temperature storage. However, the important regulators and pathways involved in natural peel degreening remain largely unknown. To determine how natural peel degreening is regulated in lemon fruit (Citrus limon), we studied transcriptome and physiochemical changes in the flavedo in response to ethylene treatment and low temperatures. Treatment with ethylene induced rapid peel degreening, which was strongly inhibited by the ethylene antagonist, 1-methylcyclopropene (1-MCP). Compared with 25 ºC, moderately low storage temperatures of 5-20 °C also triggered peel degreening. Surprisingly, repeated 1-MCP treatments failed to inhibit the peel degreening induced by low temperature. Transcriptome analysis revealed that low temperature and ethylene independently regulated genes associated with chlorophyll degradation, carotenoid metabolism, photosystem proteins, phytohormone biosynthesis and signalling, and transcription factors. Peel degreening of fruit on trees occurred in association with drops in ambient temperature, and it coincided with the differential expression of low temperature-regulated genes. In contrast, genes that were uniquely regulated by ethylene showed no significant expression changes during on-tree peel degreening. Based on these findings, we hypothesize that low temperature plays a prominent role in regulating natural peel degreening independently of ethylene in citrus fruit.

Aberrant Stamen Development is Associated with Parthenocarpic Fruit Set Through Up-Regulation of Gibberellin Biosynthesis in Tomato.

Parthenocarpy, a process in which fruit set occurs without fertilization, leads to the production of seedless fruit. A number of floral homeotic mutants with abnormal stamen development exhibit parthenocarpic fruit set. Flower development is thought to repress ovary growth before anthesis. However, the mechanism of parthenocarpic fruit development caused by aberrant flower formation is poorly understood. To investigate the molecular mechanism of parthenocarpic fruit development in floral homeotic mutants, we performed functional analysis of Tomato APETALA3 (TAP3) by loss-of-function approaches. Organ-specific promoter was used to induce organ-specific loss of function in stamen and ovary/fruit. We observed increased cell expansion in tap3 mutants and TAP3-RNAi lines during parthenocarpic fruit growth. These were predominantly accompanied by the up-regulation of GA biosynthesis genes, including SlGA20ox1, SlGA20ox2, and SlGA20ox3, as well as reduced expression of the GA-inactivating gene SlGA2ox1 and the auxin signaling gene SlARF7 involved in a crosstalk between GA and auxin. These transcriptional profiles are in agreement with the GA levels in these lines. These results suggest that stamen development negatively regulates fruit set by repressing the GA biosynthesis.

[Toxoplasmosis-associated central nervous system vasculitis accompanied by multiple cerebral hemorrhages developing subsequent to cord blood transplantation].

A 57-year-old man with high-risk myelodysplastic syndrome underwent umbilical cord blood transplantation. He began receiving steroids on day 14 for acute graft-versus-host disease, and experienced dizziness on day 75 during gradual dose reduction. Multiple hemorrhages were observed in the cerebrum, cerebellum, and brainstem. His bleeding increased, and he underwent a brain biopsy on day 91. Subsequently, he was diagnosed with central nervous system vasculitis (CNSV) on the basis of the observed aggregation of mature CD3+ lymphocytes around small vessels and vascular wall invasion by lymphocytes and macrophages. After receiving high-dose steroid therapy, cerebral hemorrhage stopped; however, dysphasia occurred on day 113 and the patient died of cerebral edema on day 128. Toxoplasma DNA and tachyzoites were detected in the brain biopsy specimen during additional examinations; therefore, we suspected that the toxoplasmosis was related to the onset of CNSV. CNSV is a rare, rapidly progressing disease that may present as a fatal post-transplantation central nervous system complication. Investigating the causes of CNSV, including CNSV associated with toxoplasmosis, is critically important for improving the prognosis of patients with CNSV.

Comparative transcriptome analysis reveals distinct ethylene-independent regulation of ripening in response to low temperature in kiwifruit.

BACKGROUND: Kiwifruit are classified as climacteric since exogenous ethylene (or its analogue propylene) induces rapid ripening accompanied by ethylene production under positive feedback regulation. However, most of the ripening-associated changes (Phase 1 ripening) in kiwifruit during storage and on-vine occur largely in the absence of any detectable ethylene. This ripening behavior is often attributed to basal levels of system I ethylene, although it is suggested to be modulated by low temperature. RESULTS: To elucidate the mechanisms regulating Phase 1 ripening in kiwifruit, a comparative transcriptome analysis using fruit continuously exposed to propylene (at 20 °C), and during storage at 5 °C and 20 °C was conducted. Propylene exposure induced kiwifruit softening, reduction of titratable acidity (TA), increase in soluble solids content (SSC) and ethylene production within 5 days. During storage, softening and reduction of TA occurred faster in fruit at 5 °C compared to 20 °C although no endogenous ethylene production was detected. Transcriptome analysis revealed 3761 ripening-related differentially expressed genes (DEGs), of which 2742 were up-regulated by propylene while 1058 were up-regulated by low temperature. Propylene exclusively up-regulated 2112 DEGs including those associated with ethylene biosynthesis and ripening such as AcACS1, AcACO2, AcPL1, AcXET1, Acß-GAL, AcAAT, AcERF6 and AcNAC7. Similarly, low temperature exclusively up-regulated 467 DEGS including AcACO3, AcPL2, AcPMEi, AcADH, Acß-AMY2, AcGA2ox2, AcNAC5 and AcbZIP2 among others. A considerable number of DEGs such as AcPG, AcEXP1, AcXET2, Acß-AMY1, AcGA2ox1, AcNAC6, AcMADS1 and AcbZIP1 were up-regulated by either propylene or low temperature. Frequent 1-MCP treatments failed to inhibit the accelerated ripening and up-regulation of associated DEGs by low temperature indicating that the changes were independent of ethylene. On-vine kiwifruit ripening proceeded in the absence of any detectable endogenous ethylene production, and coincided with increased expression of low temperature-responsive DEGs as well as the decrease in environmental temperature. CONCLUSIONS: These results indicate that kiwifruit possess both ethylene-dependent and low temperature-modulated ripening mechanisms that are distinct and independent of each other. The current work provides a foundation for elaborating the control of these two ripening mechanisms in kiwifruit.

TOMATOMICS: A Web Database for Integrated Omics Information in Tomato.

Solanum lycopersicum (tomato) is an important agronomic crop and a major model fruit-producing plant. To facilitate basic and applied research, comprehensive experimental resources and omics information on tomato are available following their development. Mutant lines and cDNA clones from a dwarf cultivar, Micro-Tom, are two of these genetic resources. Large-scale sequencing data for ESTs and full-length cDNAs from Micro-Tom continue to be gathered. In conjunction with information on the reference genome sequence of another cultivar, Heinz 1706, the Micro-Tom experimental resources have facilitated comprehensive functional analyses. To enhance the efficiency of acquiring omics information for tomato biology, we have integrated the information on the Micro-Tom experimental resources and the Heinz 1706 genome sequence. We have also inferred gene structure by comparison of sequences between the genome of Heinz 1706 and the transcriptome, which are comprised of Micro-Tom full-length cDNAs and Heinz 1706 RNA-seq data stored in the KaFTom and Sequence Read Archive databases. In order to provide large-scale omics information with streamlined connectivity we have developed and maintain a web database TOMATOMICS (http://bioinf.mind.meiji.ac.jp/tomatomics/). In TOMATOMICS, access to the information on the cDNA clone resources, full-length mRNA sequences, gene structures, expression profiles and functional annotations of genes is available through search functions and the genome browser, which has an intuitive graphical interface.

Two Members of the Aluminum-Activated Malate Transporter Family, SlALMT4 and SlALMT5, are Expressed during Fruit Development, and the Overexpression of SlALMT5 Alters Organic Acid Contents in Seeds in Tomato (Solanum lycopersicum).

The aluminum-activated malate transporter (ALMT) family of proteins transports malate and/or inorganic anions across plant membranes. To demonstrate the possible role of ALMT genes in tomato fruit development, we focused on SlALMT4 and SlALMT5, the two major genes expressed during fruit development. Predicted proteins were classified into clade 2 of the family, many members of which localize to endomembranes. Tissue-specific gene expression was determined using transgenic tomato expressing the ß-glucuronidase reporter gene controlled by their own promoters. Both the genes were expressed in vascular bundles connecting to developing seeds in fruit and in the embryo of mature seeds. Further, SlALMT5 was expressed in embryo in developing seeds in fruit. Subcellular localization of both proteins to the endoplasmic reticulum (ER) was established by transiently expressing the green fluorescent protein fusions in plant protoplasts. SlALMT5 probably localized to other endomembranes as well. Localization of SlALMT5 to the ER was also confirmed by immunoblot analysis. The transport function of both SlALMT proteins was investigated electrophysiologically in Xenopus oocytes. SlALMT5 transported malate and inorganic anions such as nitrate and chloride, but not citrate. SlALMT4 also transported malate, but the results were less consistent perhaps because it did not localize strongly to the plasma membrane. To elucidate the physiological role of SlALMT5 further, we overexpressed SlALMT5 in tomato. Compared with the wild type, overexpressors exhibited higher malate and citrate contents in mature seeds, but not in fruit. We conclude that the malate transport function of SlALMT5 expressed in developing fruit influences the organic acid contents in mature seeds.

Plant Omics Data Center: an integrated web repository for interspecies gene expression networks with NLP-based curation.

Comprehensive integration of large-scale omics resources such as genomes, transcriptomes and metabolomes will provide deeper insights into broader aspects of molecular biology. For better understanding of plant biology, we aim to construct a next-generation sequencing (NGS)-derived gene expression network (GEN) repository for a broad range of plant species. So far we have incorporated information about 745 high-quality mRNA sequencing (mRNA-Seq) samples from eight plant species (Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum, Sorghum bicolor, Vitis vinifera, Solanum tuberosum, Medicago truncatula and Glycine max) from the public short read archive, digitally profiled the entire set of gene expression profiles, and drawn GENs by using correspondence analysis (CA) to take advantage of gene expression similarities. In order to understand the evolutionary significance of the GENs from multiple species, they were linked according to the orthology of each node (gene) among species. In addition to other gene expression information, functional annotation of the genes will facilitate biological comprehension. Currently we are improving the given gene annotations with natural language processing (NLP) techniques and manual curation. Here we introduce the current status of our analyses and the web database, PODC (Plant Omics Data Center; http://bioinf.mind.meiji.ac.jp/podc/), now open to the public, providing GENs, functional annotations and additional comprehensive omics resources.

Corrigendum: Transcriptomic analysis in tomato fruit reveals divergences in genes involved in cold stress response and fruit ripening.

[This corrects the article DOI: 10.3389/fpls.2023.1227349.].

Isolation of the floral morph-related genes in heterostylous flax (Linum grandiflorum): the genetic polymorphism and the transcriptional and post-transcriptional regulations of the S locus.

Heterostylous species have two types of flowers, thrum and pin morphs, and these are controlled by a single diallelic locus designated the 'S locus'; fertilization between these two types of flowers is successful. The S gene and the molecular mechanism by which it operates remain to be uncovered, although heterostyly has been studied since the time of Darwin. We compared transcripts and proteins of the thrum and pin flowers of heterostylous flax (Linum grandiflorum) to characterize the molecular differences between them and to elucidate the molecular machinery of heterostyly. Twelve floral morph-related genes were eventually isolated by an integrated study of subtraction and 2D-PAGE analyses, and four genes, TSS1, LgAP1, LgMYB21 and LgSKS1, were predicted to be related to heterostyly. TSS1, a thrum style-specific gene, showed some features suitable for the S gene. Although its biological function is unclear, TSS1 was expressed only in the thrum style and is probably linked to the S locus. LgMYB21, another thrum style gene, would be involved in floral morphogenesis. LgMYB21 was highly expressed in the thrum style, which is shorter than the pin style, and its overexpression in Arabidopsis reduced pistil length. Furthermore, a comparison of transcript and protein accumulations showed no differences in the mRNA accumulation of some thrum-specific proteins, including LgSKS1, suggesting that these are regulated by floral morph-specific post-transcriptional regulation. The Linum S locus regulates not only S specificity but also many floral phenotypes. Dynamic regulation of transcripts and proteins would be necessary for the pleiotropic function of the Linum S locus.

Transcriptomic analysis in tomato fruit reveals divergences in genes involved in cold stress response and fruit ripening.

Cold storage is widely used to extend the postharvest life of most horticultural crops, including tomatoes, but this practice triggers cold stress and leads to the development of undesirable chilling injury (CI) symptoms. The underlying mechanisms of cold stress response and CI development in fruits remain unclear as they are often intermingled with fruit ripening changes. To gain insight into cold responses in fruits, we examined the effect of the potent ethylene signaling inhibitor 1-methylcyclopropene (1-MCP) on fruit ripening, CI occurrence and gene expression in mature green tomatoes during storage at 20°C and 5°C. 1-MCP treatments effectively inhibited ethylene production and peel color changes during storage at 20°C. Storage at 5°C also inhibited both ethylene production and peel color change; during rewarming at 20°C, 1-MCP treatments inhibited peel color change but failed to inhibit ethylene production. Furthermore, fruits stored at 5°C for 14 d developed CI symptoms (surface pitting and decay) during the rewarming period at 20°C regardless of 1-MCP treatment. Subsequent RNA-Seq analysis revealed that cold stress triggers a large-scale transcriptomic adjustment, as noticeably more genes were differentially expressed at 5°C (8,406) than at 20°C (4,814). More importantly, we have found some important divergences among genes involved in fruit ripening (up- or down-regulated at 20°C; inhibited by 1-MCP treatment) and those involved in cold stress (up- or down-regulated at 5°C; unaffected by 1-MCP treatment). Transcriptomic adjustments unique to cold stress response were associated with ribosome biogenesis, NcRNA metabolism, DNA methylation, chromatin formation/remodeling, and alternative splicing events. These data should foster further research into cold stress response mechanisms in fruits with the ultimate aim of improving tolerance to low temperature and reduction of CI symptoms during cold storage.

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.

Low-temperature-modulated fruit ripening is independent of ethylene in 'Sanuki Gold' kiwifruit.

Fruit ripening in response to treatments with propylene, 1-methycyclopropene (1-MCP), and low temperature was characterized in 'Sanuki Gold' kiwifruit, Actinidia chinensis Planch. Propylene treatment immediately induced rapid fruit softening, increased AC-PG (polygalacturonase) and AC-EXP (expansin) mRNA accumulation, and stimulated an increase in the soluble solid concentration (SSC) and a decrease in titratable acidity (TA). After 3 d exposure to propylene, ethylene production and AC-PL (pectate lyase) mRNA accumulation were observed. 1-MCP treatment after 24 h exposure to propylene eliminated AC-PG mRNA accumulation and suppressed continued changes in SSC and TA. Application of 1-MCP at the start of the treatment, followed by continuous propylene exposure, markedly delayed fruit softening, and the expression of the cell wall-modifying genes, and changes in the SSC and TA, indicating that kiwifruit become insensitive to ethylene at least for 3 d following 1-MCP exposure. Surprisingly, significant fruit softening, mRNA accumulation of AC-PG, AC-PL, and AC-EXP, and decreased TA were observed without ethylene production in intact fruit stored at low temperature for 1 month, but not in fruit stored at room temperature. Repeated 1-MCP treatments (twice a week) failed to inhibit the changes that occurred in low temperature storage. These observations indicate that low temperature modulates the ripening of kiwifruit in an ethylene-independent manner, suggesting that kiwifruit ripening is inducible by either ethylene or low temperature signals.

Examining the Role of Low Temperature in Satsuma Mandarin Fruit Peel Degreening via Comparative Physiological and Transcriptomic Analysis.

Peel degreening is the most conspicuous aspect of fruit ripening in many citrus fruits because of its importance for marketability. In this study, peel degreening in response to propylene (an ethylene analog) and at varying storage temperatures was characterized in Satsuma mandarin (Citrus unshiu Marc.) fruit. Propylene treatment triggered rapid peel degreening (within 4-6 days), indicated by an increase in the citrus color index (CCI) and chlorophyll loss. Peel degreening was also observed in fruit at 10°C and 15°C after 28-42 days, with gradual CCI increase and chlorophyll reduction. However, fruit at 5°C, 20°C, and 25°C remained green, and no substantial changes in peel CCI and chlorophyll content were recorded during the 42-day storage duration. The transcriptomes of peels of fruit treated with propylene for 4 days and those stored at varying temperatures for 28 days were then analyzed by RNA-Seq. We identified three categories of differentially expressed genes that were regulated by (i) propylene (and by analogy, ethylene) alone, (ii) low temperature (5°C, 10°C, or 15°C vs. 25°C) alone, and (iii) either propylene or low temperature. Gene-encoding proteins associated with chlorophyll degradation (such as CuSGR1, CuNOL, CuACD2, CuCAB2, and CuLHCB2) and a transcription factor (CuERF114) were differentially expressed by propylene or low temperature. To further examine temperature-induced pathways, we also monitored gene expression during on-tree fruit maturation vs. postharvest. The onset of on-tree peel degreening coincided with autumnal drops in field temperatures, and it was accompanied by differential expression of low temperature-regulated genes. On the contrary, genes that were exclusively regulated by propylene (such as CuCOPT1 and CuPOX-A2) displayed insignificant expression changes during on-tree peel degreening. These findings indicate that low temperatures could be involved in the fruit ripening-related peel degreening independently of ethylene.

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.

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.

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.