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.

Spinach (Spinacia oleracea L.) is commonly considered a dioecious plant with heterogametic (XY) and homogametic (XX) sex chromosomes. The characteristic is also utilized for the production of spinach hybrid seeds. However, the molecular mechanisms of sex determination in spinach are still unclear because of a lack of genomic and transcriptomic information. In this study, RNA sequencing (RNA-seq) was performed in male and female inflorescences to provide insight into the molecular basis of sex determination in spinach. Comparative transcriptome analyses showed that 2278 differentially expressed genes (DEGs) were identified between male and female inflorescences. A high correlation between the RNA-Seq and qRT-PCR validation for DEGs was observed. Among these, 182 DEGs were annotated to transcription factors including the MYB family protein, bHLH family, and MADS family, suggesting these factors might play a vital role in sex determination. Moreover, 26 DEGs related to flower development, including nine ABCE class genes, were detected. Expression analyses of hormone pathways showed that brassinosteroids may be key hormones related to sex determination in spinach. Overall, this study provides a large amount of DEGs related to sexual expression and lays a foundation for unraveling the regulatory mechanism of sex determination in spinach.

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 BC₁ 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 BC₁ 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.

The sex-linked region (SLR) plays an important role in determining the sex of a plant. The SLR of the Y chromosome, composed of a 14.1-Mb inversion and a 10-Mb Y-duplication region (YDR), was deciphered in Spinacia oleracea previously. However, our understanding of the SLR in its wild relatives, S. turkestanica and S. tetrandra, remains limited. In this study, we used 63 resequencing data from the three Spinacia species to infer the evolution of the SLR among the Spinacia species. In the SLR, all the cultivated spinach and S. turkestanica accessions were clustered into two distinct categories with both sexes, while the S. tetrandra accessions of both sexes were grouped. This suggests that S. oleracea shared a similar SLR with S. turkestanica, but not with S. tetrandra, which was further confirmed based on the population structure and principal component analysis. Furthermore, we identified 3910 fully sex-linked SNPs in S. oleracea and 92.82% of them were available in S. turkestanica, while none of the SNPs were adopted in S. tetrandra. Genome coverage in males and females supported the hypothesis that the YDR increasingly expanded during its evolution. Otherwise, we identified 13 sex-linked transposable element insertion polymorphisms within the inversion in both S. oleracea and S. turkestanica, demonstrating that the transposable element insertions might have occurred before the recombination suppression event of the inversion. The SLR was conserved compared with the pseudoautosomal region given that the genetic hitchhiking process occurred in the SLR during its evolution. Our findings will significantly advance our understanding of the characteristics and evolution of the SLR in Spinacia species.

The species Brassica rapa includes enormous leafy vegetables with extreme leaf morphological diversity. Leaf traits such as size, shape, weight, and ratio of the leaf blade to the petiole contribute to yield, appearance, and desirability to consumers. These leaf-related traits are controlled by quantitative trait loci (QTLs). The construction of high-density bin maps using low-coverage sequencing is a powerful method for QTL fine-mapping and gene identification. In this study, we performed whole-genome re-sequencing of Wutacai ‘Zhongbaye’ and Chinese cabbage ‘HN53’ and 150 F₂ individuals to construct a high-density bin map for QTL mapping of 11 leaf-related traits. The parental lines and F₂ population were re-sequenced at 10x and 1x coverage, respectively. A map containing 565 bin markers was constructed based on parental single-nucleotide polymorphisms and a modified sliding window approach. The total map length was 944.6 cM and the average distance of the bins was 1.65 cM. In total, 60 significant QTLs controlling 11 leaf-related traits were detected. We further identified candidate genes responsible for these complex leaf-related traits. These findings suggest that this cost-effective bin-mapping approach is capable of rapid identification of QTLs and candidate genes, and will thus facilitate the dissection of the underlying molecular basis of leaf morphological variations and accelerate the improvement of B. rapa vegetable breeding.

Bolting is a symbol of the transition from vegetative to reproductive growth in plants. Late bolting can effectively prolong the commercial value of spinach and is of great importance for spinach breeding. Bolting has complex regulatory networks, and current research on spinach bolting is relatively weak, with specific regulatory pathways and genes unclear. To clarify the regulatory characteristics and key genes related to bolting in spinach, we conducted a comparative transcriptome analysis. In this study, 18 samples from three periods of bolting-tolerant spinach material 12S3 and bolting-susceptible material 12S4 were analyzed using RNA-seq on, resulting in 10,693 differentially expressed genes (DEGs). Functional enrichment and co-expression trend analysis indicated that most DEGs were enriched in the photoperiod pathway, the hormone signaling pathway, and the cutin, suberin, and wax biosynthetic pathways. According to the weighted gene co-expression network analysis (WGCNA), SpFT (SOV4g003400), SOV4g040250, and SpGASA1 (SOV6g017600) were likely to regulate bolting through the gibberellin and photoperiod pathways, and SpELF4 (SOV1g028600) and SpPAT1 (SOV4g058860) caused differences in early and late bolting among different cultivars. These results provide important insights into the genetic control of bolting in spinach and will help elucidate the molecular mechanisms of bolting in leafy vegetables.