Gibberellins regulate masculinization through the SpGAI-SpSTM module in dioecious spinach.

The sex of dioecious plants is mainly determined by genetic factors, but it can also be converted by environmental cues such as exogenous phytohormones. Gibberellic acids (GAs) are well-known inducers of flowering and sexual development, yet the pathway of gibberellin-induced sex conversion in dioecious spinach (Spinacia oleracea L.) remains elusive. Based on sex detection before and after GA3 application using T11A and SSR19 molecular markers, we confirmed and elevated the masculinization effect of GA on a single female plant through exogenous applications of GA3, showing complete conversion and functional stamens. Silencing of GIBBERELLIC ACID INSENSITIVE (SpGAI), a single DELLA family protein that is a central GA signaling repressor, results in similar masculinization. We also show that SpGAI can physically interact with the spinach KNOX transcription factor SHOOT MERISTEMLESS (SpSTM), which is a homolog of the flower meristem identity regulator STM in Arabidopsis. The silencing of SpSTM also masculinized female flowers in spinach. Furthermore, SpSTM could directly bind the intron of SpPI to repress SpPI expression in developing female flowers. Overall, our results suggest that GA induces a female masculinization process through the SpGAI-SpSTM-SpPI regulatory module in spinach. These insights may help to clarify the molecular mechanism underlying the sex conversion system in dioecious plants while also elucidating the physiological basis for the generation of unisexual flowers so as to establish dioecy in plants.

Genome-wide characterization, functional analysis, and expression profiling of the Aux/IAA gene family in spinach.

BACKGROUND: The auxin/indole-3-acetic acid (Aux/IAA) gene family is a crucial element of the auxin signaling pathway, significantly influencing plant growth and development. Hence, we conducted a comprehensive investigation of Aux/IAAs gene family using the Sp75 and Monoe-Viroflay genomes in spinach. RESULTS: A total of 24 definitive Aux/IAA genes were identified, exhibiting diverse attributes in terms of amino acid length, molecular weight, and isoelectric points. This diversity underscores potential specific roles within the family, such as growth regulation and stress response. Structural analysis revealed significant variations in gene length and molecular weight. These variations indicate distinct roles within the Aux/IAA gene family. Chromosomal distribution analysis exhibited a dispersed pattern, with chromosomes 4 and 1 hosting the highest and lowest numbers of Aux/IAA genes, respectively. Phylogenetic analysis grouped the identified genes into distinct clades, revealing potential evolutionary relationships. Notably, the phylogenetic tree highlighted specific gene clusters suggesting shared genetic ancestry and potential functional synergies within spinach. Expression analysis under NAA treatment unveiled gene-specific and time-dependent responses, with certain genes exhibiting distinct temporal expression patterns. Specifically, SpoIAA5 displayed a substantial increase at 2 h post-NAA treatment, while SpoIAA7 and SpoIAA9 demonstrated continuous rises, peaking at the 4-hour time point. CONCLUSIONS: These observations indicate a complex interplay of gene-specific and temporal regulation in response to auxin. Moreover, the comparison with other plant species emphasized both shared characteristics and unique features in Aux/IAA gene numbers, providing insights into the evolutionary dynamics of this gene family. This comprehensive characterization of Aux/IAA genes in spinach not only establishes the foundation for understanding their specific functions in spinach development but also provides a valuable resource for experimental validation and further exploration of their roles in the intricate network of auxin signaling pathways.

Dual transcriptional characterization of spinach and Peronospora effusa during resistant and susceptible race-cultivar interactions.

BACKGROUND: Spinach downy mildew, caused by the obligate oomycete pathogen, Peronospora effusa remains a major concern for spinach production. Disease control is predominantly based on development of resistant spinach cultivars. However, new races and novel isolates of the pathogen continue to emerge and overcome cultivar resistance. Currently there are 20 known races of P. effusa. Here we characterized the transcriptomes of spinach, Spinacia oleracea, and P. effusa during disease progression using the spinach cultivar Viroflay, the near isogenic lines NIL1 and NIL3, and P. effusa races, R13 and R19, at 24 h post inoculation and 6 days post inoculation. A total of 54 samples were collected and subjected to sequencing and transcriptomic analysis. RESULTS: Differentially expressed gene (DEG) analysis in resistant spinach interactions of R13-NIL1 and R19-NIL3 revealed spinach DEGs from protein kinase-like and P-loop containing families, which have roles in plant defense. The homologous plant defense genes included but were not limited to, receptor-like protein kinases (Spiol0281C06495, Spiol06Chr21559 and Spiol06Chr24027), a BAK1 homolog (Spiol0223C05961), genes with leucine rich repeat motifs (Spiol04Chr08771, Spiol04Chr01972, Spiol05Chr26812, Spiol04Chr11049, Spiol0084S08137, Spiol03Chr20299) and ABC-transporters (Spiol02Chr28975, Spiol06Chr22112, Spiol06Chr03998 and Spiol04Chr09723). Additionally, analysis of the expression of eight homologous to previously reported downy mildew resistance genes revealed that some are differentially expressed during resistant reactions but not during susceptible reactions. Examination of P. effusa gene expression during infection of susceptible cultivars identified expressed genes present in R19 or R13 including predicted RxLR and Crinkler effector genes that may be responsible for race-specific virulence on NIL1 or NIL3 spinach hosts, respectively. CONCLUSIONS: These findings deliver foundational insight to gene expression in both spinach and P. effusa during susceptible and resistant interactions and provide a library of candidate genes for further exploration and functional analysis. Such resources will be beneficial to spinach breeding efforts for disease resistance in addition to better understanding the virulence mechanisms of this obligate pathogen.

Insights into spinach domestication from genome sequences of two wild spinach progenitors, Spinacia turkestanica and Spinacia tetrandra.

Cultivated spinach (Spinacia oleracea) is a dioecious species. We report high-quality genome sequences for its two closest wild relatives, Spinacia turkestanica and Spinacia tetrandra, which are also dioecious, and are used to study the genetics of spinach domestication. Using a combination of genomic approaches, we assembled genomes of both these species and analyzed them in comparison with the previously assembled S. oleracea genome. These species diverged c. 6.3 million years ago (Ma), while cultivated spinach split from S. turkestanica 0.8 Ma. In all three species, all six chromosomes include very large gene-poor, repeat-rich regions, which, in S. oleracea, are pericentromeric regions with very low recombination rates in both male and female genetic maps. We describe population genomic evidence that the similar regions in the wild species also recombine rarely. We characterized 282 structural variants (SVs) that have been selected during domestication. These regions include genes associated with leaf margin type and flowering time. We also describe evidence that the downy mildew resistance loci of cultivated spinach are derived from introgression from both wild spinach species. Collectively, this study reveals the genome architecture of spinach assemblies and highlights the importance of SVs during the domestication of cultivated spinach.

Evolution of the spinach sex-linked region within a rarely recombining pericentromeric region.

Sex chromosomes have evolved independently in many different plant lineages. Here, we describe reference genomes for spinach (Spinacia oleracea) X and Y haplotypes by sequencing homozygous XX females and YY males. The long arm of 185-Mb chromosome 4 carries a 13-Mb X-linked region (XLR) and 24.1-Mb Y-linked region (YLR), of which 10 Mb is Y specific. We describe evidence that this reflects insertions of autosomal sequences creating a "Y duplication region" or "YDR" whose presence probably directly reduces genetic recombination in the immediately flanking regions, although both the X and Y sex-linked regions are within a large pericentromeric region of chromosome 4 that recombines rarely in meiosis of both sexes. Sequence divergence estimates using synonymous sites indicate that YDR genes started diverging from their likely autosomal progenitors about 3 MYA, around the time when the flanking YLR stopped recombining with the XLR. These flanking regions have a higher density of repetitive sequences in the YY than the XX assembly and include slightly more pseudogenes compared with the XLR, and the YLR has lost about 11% of the ancestral genes, suggesting some degeneration. Insertion of a male-determining factor would have caused Y linkage across the entire pericentromeric region, creating physically small, highly recombining, terminal pseudoautosomal regions. These findings provide a broader understanding of the origin of sex chromosomes in spinach.

Genome-Wide Identification and Characterization of MYB Gene Family and Analysis of Its Sex-Biased Expression Pattern in Spinacia oleracea L.

The members of the myeloblastosis (MYB) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, such as circadian rhythm, metabolism, and flower development. However, the characterization of MYB genes across the genomes of spinach Spinacia oleracea L. has not been reported. Here, we identified 140 MYB genes in spinach and described their characteristics using bioinformatics approaches. Among the MYB genes, 54 were 1R-MYB, 80 were 2R-MYB, 5 were 3R-MYB, and 1 was 4R-MYB. Almost all MYB genes were located in the 0-30 Mb region of autosomes; however, the 20 MYB genes were enriched at both ends of the sex chromosome (chromosome 4). Based on phylogeny, conserved motifs, and the structure of genes, 2R-MYB exhibited higher conservation relative to 1R-MYB genes. Tandem duplication and collinearity of spinach MYB genes drive their evolution, enabling the functional diversification of spinach genes. Subcellular localization prediction indicated that spinach MYB genes were mainly located in the nucleus. Cis-acting element analysis confirmed that MYB genes were involved in various processes of spinach growth and development, such as circadian rhythm, cell differentiation, and reproduction through hormone synthesis. Furthermore, through the transcriptome data analysis of male and female flower organs at five different periods, ten candidate genes showed biased expression in spinach males, suggesting that these genes might be related to the development of spinach anthers. Collectively, this study provides useful information for further investigating the function of MYB TFs and novel insights into the regulation of sex determination in spinach.

Gene Regulatory Network Controlling Flower Development in Spinach (Spinacia oleracea L.).

Spinach (Spinacia oleracea L.) is a dioecious, diploid, wind-pollinated crop cultivated worldwide. Sex determination plays an important role in spinach breeding. Hence, this study aimed to understand the differences in sexual differentiation and floral organ development of dioecious flowers, as well as the differences in the regulatory mechanisms of floral organ development of dioecious and monoecious flowers. We compared transcriptional-level differences between different genders and identified differentially expressed genes (DEGs) related to spinach floral development, as well as sex-biased genes to investigate the flower development mechanisms in spinach. In this study, 9189 DEGs were identified among the different genders. DEG analysis showed the participation of four main transcription factor families, MIKC_MADS, MYB, NAC, and bHLH, in spinach flower development. In our key findings, abscisic acid (ABA) and gibberellic acid (GA) signal transduction pathways play major roles in male flower development, while auxin regulates both male and female flower development. By constructing a gene regulatory network (GRN) for floral organ development, core transcription factors (TFs) controlling organ initiation and growth were discovered. This analysis of the development of female, male, and monoecious flowers in spinach provides new insights into the molecular mechanisms of floral organ development and sexual differentiation in dioecious and monoecious plants in spinach.

Real-Time PCR Assays for Races of the Spinach Fusarium Wilt Pathogen, Fusarium oxysporum f. sp. spinaciae.

Fusarium wilt of spinach, caused by Fusarium oxysporum f. sp. spinaciae, is a significant limitation for producers of vegetative spinach and spinach seed crops during warm temperatures and/or on acid soils. Identification of isolates of F. oxysporum f. sp. spinaciae, and distinction of isolates of the two known races, entails time-intensive pathogenicity tests. In this study, two real-time PCR assays were developed: one for a candidate effector gene common to both races of F. oxysporum f. sp. spinaciae, and another for a candidate effector gene unique to isolates of race 2. The assays were specific to isolates of F. oxysporum f. sp. spinaciae (n = 44) and isolates of race 2 (n = 23), respectively. Neither assay amplified DNA from 10 avirulent isolates of F. oxysporum associated with spinach, 57 isolates of other formae speciales and Fusarium spp., or 7 isolates of other spinach pathogens. When the assays were used to detect DNA extracted from spinach plants infected with an isolate of race 1, race 2, or a 1:1 mixture of both races, the amount of target DNA detected increased with increasing severity of wilt. Plants infected with one or both isolates could be distinguished based on the ratio in copy number for each target locus. The real-time PCR assays enable rapid diagnosis of Fusarium wilt of spinach and will facilitate research on the epidemiology and management of this disease, as well as surveys on the prevalence of this understudied pathogen in regions of spinach and/or spinach seed production.

Genome-wide identification and expression analysis reveals spinach brassinosteroid-signaling kinase (BSK) gene family functions in temperature stress response.

BACKGROUND: Brassinosteroid (BR)- signaling kinase (BSK) is a critical family of receptor-like cytoplasmic kinase for BR signal transduction, which plays important roles in plant development, immunity, and abiotic stress responses. Spinach (Spinacia oleracea) is cold- tolerant but heat- sensitive green leafy vegetable. A study on BSK family members and BSKs- mediated metabolic processes in spinach has not been performed. RESULTS: We identified and cloned seven SoBSKs in spinach. Phylogenetic and collinearity analyses suggested that SoBSKs had close relationship with dicotyledonous sugar beet (Beta vulgaris) rather than monocotyledons. The analyses of gene structure and conserved protein domain/ motif indicated that most SoBSKs were relative conserved, while SoBSK6 could be a truncated member. The prediction of post-translation modification (PTM) sites in SoBSKs implied their possible roles in signal transduction, redox regulation, and protein turnover of SoBSKs, especially the N-terminal myristoylation site was critical for BSK localization to cell periphery. Cis-acting elements for their responses to light, drought, temperature (heat and cold), and hormone distributed widely in the promoters of SoBSKs, implying the pivotal roles of SoBSKs in response to diverse abiotic stresses and phytohormone stimuli. Most SoBSKs were highly expressed in leaves, except for SoBSK7 in roots. Many SoBSKs were differentially regulated in spinach heat- sensitive variety Sp73 and heat- tolerant variety Sp75 under the treatments of heat, cold, as well as exogenous brassinolide (BL) and abscisic acid (ABA). The bsk134678 mutant Arabidopsis seedlings exhibited more heat tolerance than wild- type and SoBSK1- overexpressed seedlings. CONCLUSIONS: A comprehensive genome- wide analysis of the BSK gene family in spinach presented a global identification and functional prediction of SoBSKs. Seven SoBSKs had relatively- conserved gene structure and protein function domains. Except for SoBSK6, all the other SoBSKs had similar motifs and conserved PTM sites. Most SoBSKs participated in the responses to heat, cold, BR, and ABA. These findings paved the way for further functional analysis on BSK- mediated regulatory mechanisms in spinach development and stress response.

Spinach-based RNA mimicking GFP in plant cells.

Spinach RNA-mimicking GFP (S-RMG) has been successfully used to monitor cellular RNAs including microRNAs in bacterium, yeast, and human cells. However, S-RMG has not been established in plants. In this study, we found that like bacterial, yeast, and human cellular tRNAs, plant tRNAs such as tRNALys can protect and/or stabilize the Spinach RNA aptamer interaction with the fluorophore DFHBI enabling detectable levels of green fluorescence to be emitted. The tRNALys-Spinach-tRNALys, once delivered into "chloroplast-free" onion epidermal cells can emit strong green fluorescence in the presence of DFHBI. Our results demonstrate for the first time that Spinach-based RNA visualization has the potential for in vivo monitoring of RNAs in plant cells.

Sexual reproduction contributes to the evolution of resistance-breaking isolates of the spinach pathogen Peronospora effusa.

Peronospora effusa causes downy mildew, the economically most important disease of cultivated spinach worldwide. To date, 19 P. effusa races have been denominated based on their capacity to break spinach resistances, but their genetic diversity and the evolutionary processes that contribute to race emergence are unknown. Here, we performed the first systematic analysis of P. effusa races showing that those emerge by both asexual and sexual reproduction. Specifically, we studied the diversity of 26 P. effusa isolates from 16 denominated races based on mitochondrial and nuclear comparative genomics. Mitochondrial genomes based on long-read sequencing coupled with diversity assessment based on short-read sequencing uncovered two mitochondrial haplogroups, each with distinct genome organization. Nuclear genome-wide comparisons of the 26 isolates revealed that 10 isolates from six races could clearly be divided into three asexually evolving groups, in concordance with their mitochondrial phylogeny. The remaining isolates showed signals of reticulated evolution and discordance between nuclear and mitochondrial phylogenies, suggesting that these evolved through sexual reproduction. Increased understanding of this pathogen's reproductive modes will provide the framework for future studies into the molecular mechanisms underlying race emergence and into the P. effusa-spinach interaction, thus assisting in sustainable production of spinach through knowledge-driven resistance breeding.

Dye-Conjugated Spinach RNA by Genetic Alphabet Expansion.

Light-emitting systems using an RNA aptamer-dye pair, such as Spinach RNA, are an attractive method for imaging and tracing RNA expression in vitro and in vivo. We present an alternative Spinach method by genetic alphabet expansion using an unnatural base pair system, in which a dye-conjugated unnatural base substrate is site-specifically incorporated at a specific position in Spinach RNA by transcription involving the third base pair. The incorporation position was predicted by molecular dynamics simulations. This dye-conjugated Spinach RNA increased the thermal stability of the fluorescence, the robustness against ion sensitivity, and the resistance against photobleaching. Furthermore, we applied our method to Baby Spinach, a shorter version of Spinach, for dye conjugation toward the visible detection of transcripts. This is the first demonstration of an alternative RNA imaging method for a detection system using genetic alphabet expansion.

Plant terpenoid metabolism co-opts a component of the cell wall biosynthesis machinery.

Glycosylation is one of the most prevalent molecular modifications in nature. Single or multiple sugars can decorate a wide range of acceptors from proteins to lipids, cell wall glycans and small molecules, dramatically affecting their activity. Here, we discovered that by 'hijacking' an enzyme of the cellulose synthesis machinery involved in cell wall assembly, plants evolved cellulose synthase-like enzymes (Csls) and acquired the capacity to glucuronidate specialized metabolites, that is, triterpenoid saponins. Apparently, endoplasmic reticulum-membrane localization of Csls and of other pathway proteins was part of evolving a new glycosyltransferase function, as plant metabolite glycosyltransferases typically act in the cytosol. Discovery of glucuronic acid transferases across several plant orders uncovered the long-pursued enzymatic reaction in the production of a low-calorie sweetener from licorice roots. Our work opens the way for engineering potent saponins through microbial fermentation and plant-based systems.

Early Detection of the Spinach Downy Mildew Pathogen in Leaves by Recombinase Polymerase Amplification.

Downy mildew of spinach, caused by Peronospora effusa, is a major economic threat to both organic and conventional spinach production. Symptomatic spinach leaves are unmarketable and spinach with latent infections are problematic because symptoms can develop postharvest. Therefore, early detection methods for P. effusa could help producers identify infection before visible symptoms appear. Recombinase polymerase amplification (RPA) provides sensitive and specific detection of pathogen DNA and is a rapid, field-applicable method that does not require advanced technical knowledge or equipment-heavy DNA extraction. Here, we used comparative genomics to identify a unique region of the P. effusa mitochondrial genome to develop an RPA assay for the early detection of P. effusa in spinach leaves. In tandem, we established a TaqMan quantitative PCR (qPCR) assay and used this assay to validate the P. effusa specificity of the locus across Peronospora spp. and to compare assay performance. Neither the TaqMan qPCR nor the RPA showed cross reactivity with the closely related beet downy mildew pathogen, P. schachtii. TaqMan qPCR and RPA have detection thresholds of 100 and 900 fg of DNA, respectively. Both assays could detect P. effusa in presymptomatic leaves, with RPA-based detection occurring as early as 5 days before the appearance of symptoms and TaqMan qPCR-based detection occurring after 24 h of plant exposure to airborne spores. Implementation of the RPA detection method could provide real-time information for point-of-care management strategies at field sites.

Fine Mapping and Identification of a Candidate Gene of Downy Mildew Resistance, RPF2, in Spinach (Spinacia oleracea L.).

Downy mildew is a major threat to the economic value of spinach. The most effective approach to managing spinach downy mildew is breeding cultivars with resistance genes. The resistance allele RPF2 is effective against races 1-10 and 15 of Peronospora farinosa f. sp. Spinaciae (P. effusa) and is widely used as a resistance gene. However, the gene and the linked marker of RPF2 remain unclear, which limit its utilization. Herein, we located the RPF2 gene in a 0.61 Mb region using a BC1 population derived from Sp39 (rr) and Sp62 (RR) cultivars via kompetitive allele specific PCR (KASP) markers. Within this region, only one R gene, Spo12821, was identified based on annotation information. The amino acid sequence analysis showed that there were large differences in the length of the LRR domain between the parents. Additionally, a molecular marker, RPF2-IN12821, was developed based on the sequence variation in the Spo12821, and the evaluation in the BC1 population produced a 100% match with resistance/susceptibility. The finding of the study could be valuable for improving our understanding of the genetic basis of resistance against the downy mildew pathogen and breeding resistance lines in the future.

Identification of Sex Differentiation-Related microRNAs in Spinach Female and Male Flower.

Sex determination and differentiation is an important biological process for unisexual flower development. Spinach is a model plant to study the mechanism of sex determination and differentiation of dioecious plant. Till now, little is known about spinach sex determination and differentiation mechanism. MicroRNAs are key factors in flower development. Herein, small RNA sequencing was performed to explore the roles of microRNAs in spinach sex determination and differentiation. As a result, 92 known and 3402 novel microRNAs were identified in 18 spinach female and male flower samples. 74 differentially expressed microRNAs were identified between female and male flowers, including 20 female-biased and 48 male-biased expression microRNAs. Target prediction identified 22 sex-biased microRNA-target pairs, which may be involved in spinach sex determination or differentiation. Among the differentially expressed microRNAs between FNS and M03, 55 microRNAs were found to reside in sex chromosome; one of them, sol-miR2550n, was functionally studied via genetic transformation. Silencing of sol-miR2550n resulted in abnormal anther while overexpression of sol-miR2550n induced early flowering, indicating sol-miR2550n was a male-promoting factor and validating the reliability of our small RNA sequencing data. Conclusively, this work can supply valuable information for exploring spinach sex determination and differentiation and provide a new insight in studying unisexual flower development.

High resolution mapping and candidate gene identification of downy mildew race 16 resistance in spinach.

BACKGROUND: Downy mildew, the most devastating disease of spinach (Spinacia oleracea L.), is caused by the oomycete Peronospora effusa [=P. farinosa f. sp. spinaciae]. The P. effusa shows race specificities to the resistant host and comprises 19 reported races and many novel isolates. Sixteen new P. effusa races were identified during the past three decades, and the new pathogen races are continually overcoming the genetic resistances used in commercial cultivars. A spinach breeding population derived from the cross between cultivars Whale and Lazio was inoculated with P. effusa race 16 in an environment-controlled facility; disease response was recorded and genotyped using genotyping by sequencing (GBS). The main objective of this study was to identify resistance-associated single nucleotide polymorphism (SNP) markers from the cultivar Whale against the P. effusa race 16. RESULTS: Association analysis conducted using GBS markers identified six significant SNPs (S3_658,306, S3_692697, S3_1050601, S3_1227787, S3_1227802, S3_1231197). The downy mildew resistance locus from cultivar Whale was mapped to a 0.57 Mb region on chromosome 3, including four disease resistance candidate genes (Spo12736, Spo12784, Spo12908, and Spo12821) within 2.69-11.28 Kb of the peak SNP. CONCLUSIONS: Genomewide association analysis approach was used to map the P. effusa race 16 resistance loci and identify associated SNP markers and the candidate genes. The results from this study could be valuable in understanding the genetic basis of downy mildew resistance, and the SNP marker will be useful in spinach breeding to select resistant lines.

Quantitative trait loci (QTL) analysis of leaf related traits in spinach (Spinacia oleracea L.).

BACKGROUND: Spinach (Spinacia oleracea L.) is an important leafy vegetable crop, and leaf-related traits including leaf length, leaf width, and petiole length, are important commercial traits. However, the underlying genes remain unclear. The objective of the study was to conduct QTL mapping of leaf-related traits in spinach. RESULTS: A BC1 population was used to construct the linkage map and for QTL mapping of leaf length, leaf width, petiole length, and the ratio of leaf length to width in 2015 and 2019. Two genetic linkage maps were constructed by specific locus amplified fragment sequencing (SLAF-seq), and kompetitive allele specific PCR (KASP) technology, respectively using BC1 population in 2015. Based on the results of 2015, the specific linkage groups (LG) detected QTLs were generated using BC1 population in 2019. A total of 13 QTLs were detected for leaf-related traits, only five QTLs being repeatedly detected in multiple years or linkage maps. Interestingly, the major QTLs of leaf length, petiole length, and the ratio of leaf length to width were highly associated with the same SNP markers (KM3102838, KM1360385 and KM2191098). A major QTL of leaf width was mapped on chromosome 1 from 41.470-42.045 Mb. And 44 genes were identified within the region. Based on the GO analysis, these genes were significantly enriched on ribonuclease, lyase activity, phosphodiester bond hydrolysis process, and cell wall component, thus it might change cell size to determine leaves shape. CONCLUSIONS: Five QTLs for leaf-related traits were repeatedly detected at least two years or linkage maps. The major QTLs of leaf length, petiole length, and the ratio of leaf length to width were mapped on the same loci. And three genes (Spo10792, Spo21018, and Spo21019) were identified as important candidate genes for leaf width.

Transcriptome architecture reveals genetic networks of bolting regulation in spinach.

BACKGROUND: Bolting refers to the early flowering stem production on agricultural and horticultural crops before harvesting. Indeed, bolting is an event induced by the coordinated effects of various environmental factors and endogenous genetic components, which cause a large reduction in the quality and productivity of vegetable crops like spinach. However, little is known about the signaling pathways and molecular functions involved in bolting mechanisms in spinach. The genetic information regarding the transition from vegetative growth to the reproductive stage in spinach would represent an advantage to regulate bolting time and improvement of resistant cultivars to minimize performance loss. RESULTS: To investigate the key genes and their genetic networks controlling spinach bolting, we performed RNA-seq analysis on early bolting accession Kashan and late-bolting accession Viroflay at both vegetative and reproductive stages and found a significant number of differentially expressed genes (DEGs) ranging from 195 to 1230 in different comparisons. These genes were mainly associated with the signaling pathways of vernalization, photoperiod/circadian clock, gibberellin, autonomous, and aging pathways. Gene ontology analysis uncovered terms associated with carbohydrate metabolism, and detailed analysis of expression patterns for genes of Fructose-1, 6-bisphosphate aldolase, TREHALOSE-6-PHOSPHATE SYNTHASE 1, FLOWERING PROMOTING FACTOR 1, EARLY FLOWERING, GIGANTEA, and MADS-box proteins revealed their potential roles in the initiating or delaying of bolting. CONCLUSION: This study is the first report on identifying bolting and flowering-related genes based on transcriptome sequencing in spinach, which provides insight into bolting control and can be useful for molecular breeding programs and further study in the regulation of the genetic mechanisms related to bolting in other vegetable crops.

Auxin regulated metabolic changes underlying sepal retention and development after pollination in spinach.

BACKGROUND: Pollination accelerate sepal development that enhances plant fitness by protecting seeds in female spinach. This response requires pollination signals that result in the remodeling within the sepal cells for retention and development, but the regulatory mechanism for this response is still unclear. To investigate the early pollination-induced metabolic changes in sepal, we utilize the high-throughput RNA-seq approach. RESULTS: Spinach variety 'Cornel 9' was used for differentially expressed gene analysis followed by experiments of auxin analog and auxin inhibitor treatments. We first compared the candidate transcripts expressed differentially at different time points (12H, 48H, and 96H) after pollination and detected significant difference in Trp-dependent auxin biosynthesis and auxin modulation and transduction process. Furthermore, several auxin regulatory pathways i.e. cell division, cell wall expansion, and biogenesis were activated from pollination to early developmental symptoms in sepals following pollination. To further confirm the role auxin genes play in the sepal development, auxin analog (2, 4-D; IAA) and auxin transport inhibitor (NPA) with different concentrations gradient were sprayed to the spinach unpollinated and pollinated flowers, respectively. NPA treatment resulted in auxin transport weakening that led to inhibition of sepal development at concentration 0.1 and 1 mM after pollination. 2, 4-D and IAA treatment to unpollinated flowers resulted in sepal development at lower concentration but wilting at higher concentration. CONCLUSION: We hypothesized that sepal retention and development might have associated with auxin homeostasis that regulates the sepal size by modulating associated pathways. These findings advanced the understanding of this unusual phenomenon of sepal growth instead of abscission after pollination in spinach.