Molecular evidence for recent divergence of X- and Y-linked gene pairs in Spinacia oleracea L

Genome-wide characterization of the xyloglucan endotransglucosylase/hydrolase family genes and their response to plant hormone in sugar beet

Xyloglucan endotransglucosylase/hydrolases (XTHs) play a crucial role in plant growth and development. However, their functional response to phytohormone in sugar beet still remains obscure. In this study, we identified 30 putative BvXTH genes in the sugar beet genome. Phylogenetic and evolutionary relationship analysis revealed that they were clustered into three groups and have gone through eight tandem duplication events under purifying selection. Gene structure and motif composition analysis demonstrated that they were highly conserved and all contained one conserved glycoside hydrolase family 16 domain (Glyco_hydro_16) and one xyloglucan endotransglycosylase C-terminus (XET_C) domain. Transcriptional expression analysis exhibited that all BvXTHs were ubiquitously expressed in leaves, root hairs and tuberous roots, and most of them were up-regulated by brassinolide (BR), jasmonic acid (JA), abscisic acid (ABA) and gibberellic acid (GA3). Further mutant complementary experiment demonstrated that expression of BvXTH17 rescued the retarded growth phenotype of xth22, an Arabidopsis knock out mutant of AtXTH22. The findings in our work provide fundamental information on the structure and evolutionary relationship of the XTH family genes in sugar beet, and reveal the potential function of BvXTH17 in plant growth and hormone response.

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.

Revisiting the amalgaviral landscapes in plant transcriptomes expands the host range of plant amalgaviruses

Plant amalgaviruses are monopartite, double-stranded RNA viruses, capable of vertical transmission through seeds. An attempt to revisit plant transcriptome-assembled contigs for amalgaviral sequences identified 40 putative novel amalgaviruses in 35 plant species, nearly doubling the number of plant amalgaviruses. Of the 35 plant species, 33 are reported to host amalgaviruses for the first time, including a pteridophytic and two gymnospermic species. Coding-complete genomes of all identified viruses were recovered and the putative +1 programmed ribosomal frameshift (PRF) sites were determined. Genomes of 35 identified amalgaviruses contained the conserved +1 PRF motif 'UUU_CGN', while variant versions were predicted in five genomes. Phylogenetic analysis grouped pteridophyte- and gymnosperm-infecting amalgaviruses together in divergent sub-clades while few of the related angiosperm-infecting amalgaviruses infect members of the same plant family, reiterating the co-evolution of plant amalgaviruses and their hosts. The current study paves way for further studies on understanding biological properties of identified viruses.

Remarkable Divergence of the Sex-Linked Region between Two Wild Spinach Progenitors, Spinacia turkestanica and Spinacia tetrandra

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 contributions of Nettie Stevens to the field of sex chromosome biology

The early 1900s delivered many foundational discoveries in genetics, including re-discovery of Mendel's research and the chromosomal theory of inheritance. Following these insights, many focused their research on whether the development of separate sexes had a chromosomal basis or if instead it was caused by environmental factors. It is Dr Nettie M. Stevens' Studies in spermatogenesis (1905) that provided the unequivocal evidence that the inheritance of the Y chromosome initiated male development in mealworms. This result established that sex is indeed a Mendelian trait with a genetic basis and that the sex chromosomes play a critical role. In Part II of Studies in spermatogenesis (1906), an XY pair was identified in dozens of additional species, further validating the function of sex chromosomes. Since this formative work, a wealth of studies in animals and plants have examined the genetic basis of sex. The goal of this review is to shine a light again on Stevens’ Studies in spermatogenesis and the lasting impact of this work. We additionally focus on key findings in plant systems over the last century and open questions that are best answered, as in Stevens' work, by synthesizing across many systems.

This article is part of the theme issue ‘Sex determination and sex chromosome evolution in land plants’.

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.

The spinach YY genome reveals sex chromosome evolution, domestication, and introgression history of the species

Background

Spinach (Spinacia oleracea L.) is a dioecious species with an XY sex chromosome system, but its Y chromosome has not been fully characterized. Our knowledge about the history of its domestication and improvement remains limited.

Results

A high-quality YY genome of spinach is assembled into 952 Mb in six pseudo-chromosomes. By a combination of genetic mapping, Genome-Wide Association Studies, and genomic analysis, we characterize a 17.42-Mb sex determination region (SDR) on chromosome 1. The sex chromosomes of spinach evolved when an insertion containing sex determination genes occurred, followed by a large genomic inversion about 1.98 Mya. A subsequent burst of SDR-specific repeats (0.1–0.15 Mya) explains the large size of this SDR. We identify a Y-specific gene, NRT1/PTR 6.4 which resides in this insertion, as a strong candidate for the sex determination or differentiation factor. Resequencing of 112 spinach genomes reveals a severe domestication bottleneck approximately 10.87 Kya, which dates the domestication of spinach 7000 years earlier than the archeological record. We demonstrate that a strong selection signal associated with internode elongation and leaf area expansion is associated with domestication of edibility traits in spinach. We find that several strong genomic introgressions from the wild species Spinacia turkestanica and Spinacia tetrandra harbor desirable alleles of genes related to downy mildew resistance, frost resistance, leaf morphology, and flowering-time shift, which likely contribute to spinach improvement.

Conclusions

Analysis of the YY genome uncovers evolutionary forces shaping nascent sex chromosome evolution in spinach. Our findings provide novel insights about the domestication and improvement of spinach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02633-x.

A spinach genome assembly with remarkable completeness, and its use for rapid identification of candidate genes for agronomic traits

Spinach (Spinacia oleracea) is grown as a nutritious leafy vegetable worldwide. To accelerate spinach breeding efficiency, a high-quality reference genome sequence with great completeness and continuity is needed as a basic infrastructure. Here, we used long-read and linked-read technologies to construct a de novo spinach genome assembly, designated SOL_r1.1, which was comprised of 287 scaffolds (total size: 935.7 Mb; N50 = 11.3 Mb) with a low proportion of undetermined nucleotides (Ns = 0.34%) and with high gene completeness (BUSCO complete 96.9%). A genome-wide survey of resistance gene analogues identified 695 genes encoding nucleotide-binding site domains, receptor-like protein kinases, receptor-like proteins and transmembrane-coiled coil domains. Based on a high-density double-digest restriction-site associated DNA sequencing-based linkage map, the genome assembly was anchored to six pseudomolecules representing ∼73.5% of the whole genome assembly. In addition, we used SOL_r1.1 to identify quantitative trait loci for bolting timing and fruit/seed shape, which harbour biologically plausible candidate genes, such as homologues of the FLOWERING LOCUS T and EPIDERMAL PATTERNING FACTOR-LIKE genes. The new genome assembly, SOL_r1.1, will serve as a useful resource for identifying loci associated with important agronomic traits and for developing molecular markers for spinach breeding/selection programs.

The timing of genetic degeneration of sex chromosomes

Genetic degeneration is an extraordinary feature of sex chromosomes, with the loss of functions of Y-linked genes in species with XY systems, and W-linked genes in ZW systems, eventually affecting almost all genes. Although degeneration is familiar to most biologists, important aspects are not yet well understood, including how quickly a Y or W chromosome can become completely degenerated. I review the current understanding of the time-course of degeneration. Degeneration starts after crossing over between the sex chromosome pair stops, and theoretical models predict an initially fast degeneration rate and a later much slower one. It has become possible to estimate the two quantities that the models suggest are the most important in determining degeneration rates-the size of the sex-linked region, and the time when recombination became suppressed (which can be estimated using Y-X or W-Z sequence divergence). However, quantifying degeneration is still difficult. I review evidence on gene losses (based on coverage analysis) or loss of function (by classifying coding sequences into functional alleles and pseudogenes). I also review evidence about whether small genome regions degenerate, or only large ones, whether selective constraints on the genes in a sex-linked region also strongly affect degeneration rates, and about how long it takes before all (or almost all) genes are lost. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

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 BC₁ 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 BC₁ population in 2015. Based on the results of 2015, the specific linkage groups (LG) detected QTLs were generated using BC₁ 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.

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.

The Diversity of Plant Sex Chromosomes Highlighted through Advances in Genome Sequencing

For centuries, scientists have been intrigued by the origin of dioecy in plants, characterizing sex-specific development, uncovering cytological differences between the sexes, and developing theoretical models. Through the invention and continued improvements in genomic technologies, we have truly begun to unlock the genetic basis of dioecy in many species. Here we broadly review the advances in research on dioecy and sex chromosomes. We start by first discussing the early works that built the foundation for current studies and the advances in genome sequencing that have facilitated more-recent findings. We next discuss the analyses of sex chromosomes and sex-determination genes uncovered by genome sequencing. We synthesize these results to find some patterns are emerging, such as the role of duplications, the involvement of hormones in sex-determination, and support for the two-locus model for the origin of dioecy. Though across systems, there are also many novel insights into how sex chromosomes evolve, including different sex-determining genes and routes to suppressed recombination. We propose the future of research in plant sex chromosomes should involve interdisciplinary approaches, combining cutting-edge technologies with the classics to unravel the patterns that can be found across the hundreds of independent origins.

Genome-Wide Analysis of Transposable Elements and Satellite DNAs in Spinacia Species to Shed Light on Their Roles in Sex Chromosome Evolution

Sex chromosome evolution has mostly been studied in species with heteromorphic sex chromosomes. The Spinacia genus serves as an ideal model for investigating evolutionary mechanisms underlying the transition from homomorphic to heteromorphic sex chromosomes. Among evolutionary factors, repetitive sequences play multiple roles in sex chromosome evolution while their forces have not been fully explored in Spinacia species. Here, we identified major repetitive sequence classes in male and female genomes of Spinacia species and their ancestral relative sugar beet to elucidate the evolutionary processes of sex chromosome evolution using next-generation sequencing (NGS) data. Comparative analysis revealed that the repeat elements of Spinacia species are considerably higher than of sugar beet, especially the Ty3/Gypsy and Ty1/Copia retrotransposons. The long terminal repeat retroelements (LTR) Angela, Athila, and Ogre may be accounted for the higher proportion of repeats in the spinach genome. Comparison of the repeats proportion between female and male genomes of three Spinacia species indicated the different representation in Spinacia tetrandra samples but not in the S. oleracea or S. turkestanica samples. From these results, we speculated that emergence of repetitive DNA sequences may correlate the formation of sex chromosome and the transition from homomorphic sex chromosomes to heteromorphic sex chromosomes as heteromorphic sex chromosomes exclusively existed in Spinacia tetrandra. Three novel sugar beet-specific satellites were identified and confirmed by fluorescence in situ hybridization (FISH); six out of eight new spinach-specific satellites were mapped to the short arm of sex chromosomes. A total of 141 copies of SolSat01-171-s were found in the sex determination region (SDR). Thus, the accumulation of satellite DNA on the short arm of chromosome 1 may be involved in the sex chromosome evolution in Spinacia species. Our study provides a fundamental resource for understanding repeat sequences in Spinacia species and their roles in sex chromosome evolution.

Construction of high-density genetic maps defined sex determination region of the Y chromosome in spinach

Spinach (Spinacia olracea L.) is a dioecious leafy vegetable with a highly repetitive genome of around 990 Mb, which is challenging for de-novo genome assembly. In our study, a segregating F₁ (double pseudo-testcross) population from 'Viroflay' × 'Cornell-NO. 9' was used for genetic mapping by resequencing genotyping. In the paternal 'Cornell-NO. 9' map, 212,414 SNPs were mapped, and the total linkage distance was 476.83 cM; the maternal 'Viroflay' map included 29,282 SNPs with 401.28 cM total genetic distance. Both paternal and maternal maps have the expected number of six linkage groups (LGs). A non-recombining region with 5678 SNPs (39 bin markers) co-segregates with sex type which located at 45.2 cM of LG1 in the 'Cornell-NO. 9' map while indicates the sex determination region (SDR). Integration of two maps into a consensus map guided us to anchor additional 1242 contigs to six pseudomolecules from the published reference genome, which improved additional 233 Mb (23.4%) assembly based on spinach estimated genome size. Particularly, the X counterpart of SDR in our assembly is estimated around 18.4 Mb which locates at the largest chromosome, as consensus with sex-biased FISH signals from previous cytogenetics studies. The region is featured by reduced gene density, higher percentage of repetitive sequences, and no recombination. Our linkage maps provide the resource for improving spinach genome de-novo assembly and identification of sex-determining genes in spinach.

Comparative transcriptome analysis of male and female flowers in Spinacia oleracea L

Background

Dioecious spinach (Spinacia oleracea L.), a commercial and nutritional vegetable crop, serves as a model for studying the mechanisms of sex determination and differentiation in plants. However, this mechanism is still unclear. Herein, based on PacBio Iso-seq and Illumina RNA-seq data, comparative transcriptome analysis of male and female flowers were performed to explore the sex differentiation mechanism in spinach.

Results

Compared with published genome of spinach, 10,800 transcripts were newly annotated; alternative splicing, alternative polyadenylation and lncRNA were analyzed for the first time, increasing the diversity of spinach transcriptome. A total of 2965 differentially expressed genes were identified between female and male flowers at three early development stages. The differential expression of RNA splicing-related genes, polyadenylation-related genes and lncRNAs suggested the involvement of alternative splicing, alternative polyadenylation and lncRNA in sex differentiation. Moreover, 1946 male-biased genes and 961 female-biased genes were found and several candidate genes related to gender development were identified, providing new clues to reveal the mechanism of sex differentiation. In addition, weighted gene co-expression network analysis showed that auxin and gibberellin were the common crucial factors in regulating female or male flower development; however, the closely co-expressed genes of these two factors were different between male and female flower, which may result in spinach sex differentiation.

Conclusions

In this study, 10,800 transcripts were newly annotated, and the alternative splicing, alternative polyadenylation and long-noncoding RNA were comprehensively analyzed for the first time in spinach, providing valuable information for functional genome study. Moreover, candidate genes related to gender development were identified, shedding new insight on studying the mechanism of sex determination and differentiation in plant.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07277-4.

Variability in somatic embryo-forming capacity of spinach

High variability in somatic embryo (SE)-forming capacity has previously been observed in several spinach cultivars. Such variability frequently accounted for more variation in embryogenic response of the explants than the factor being investigated. Hence, the variability in embryogenic capacity was examined in the present study at both the population and the single-seedling level, using seeds of spinach cultivar Matador obtained from nine European seed companies. Seed population obtained from Slovenia (Sl) was superior to others, with the highest regeneration frequency (100%) and the highest mean SE number (14.4). A total of 82% of these seedlings had 80–100% of regenerating explants, while in populations with intermediate embryogenic capacity approximately 40% of seedlings had 20–60% of regenerating explants. The explants from the majority of seedlings (52–100%) in the least responsive populations were irresponsive. Furthermore, the explants from Sl seedlings regenerated from 10–20 (43.5%) up to > 20 (27.6%) SEs on average, while the explants from the majority of seedlings belonging to other populations regenerated 1–10 SEs. The present study strongly indicates that the variability of plant material must not be overlooked, because choosing more responsive individuals for one treatment and less responsive ones for another may lead to misinterpretation of the data.

A single gene underlies the dynamic evolution of poplar sex determination

Although hundreds of plant lineages have independently evolved dioecy (that is, separation of the sexes), the underlying genetic basis remains largely elusive¹. Here we show that diverse poplar species carry partial duplicates of the ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17) orthologue in the male-specific region of the Y chromosome. These duplicates give rise to small RNAs apparently causing male-specific DNA methylation and silencing of the ARR17 gene. CRISPR-Cas9-induced mutations demonstrate that ARR17 functions as a sex switch, triggering female development when on and male development when off. Despite repeated turnover events, including a transition from the XY system to a ZW system, the sex-specific regulation of ARR17 is conserved across the poplar genus and probably beyond. Our data reveal how a single-gene-based mechanism of dioecy can enable highly dynamic sex-linked regions and contribute to maintaining recombination and integrity of sex chromosomes.

Application of cucumber mosaic virus to efficient induction and long-term maintenance of virus-induced gene silencing in spinach

Virus-induced gene silencing (VIGS) is a useful tool for functional genomics in plants. In this study, we tried to apply cucumber mosaic virus (CMV) to efficient induction of VIGS in spinach. Although VIGS for spinach had been previously developed based on two viruses (beet curly top virus and tobacco rattle virus), they still have some problems with systemic movement and long-term maintenance of VIGS in spinach. Although ordinary CMV strains infect spinach inducing distinct mosaic symptoms, using a CMV pseudorecombinant, we can modify the viral pathogenicity to attenuate viral symptoms that may mask the silencing phenotype. We here successfully demonstrated the viral ability to silence the phytoene desaturase (PDS) and the dihydroflavonol 4-reductase (DFR) genes in spinach. Because CMV could quickly induce VIGS even at 7-10 days postinoculation and the virus did not disappear even at the flowering stage, this CMV-based VIGS system would contribute to functional genomics in spinach and especially to the elucidation of molecular mechanisms for some properties unique to spinach such as plasticity of sex expression; the CMV-induced VIGS can last until the flowering stage after the virus was inoculated onto the seedling.