Molecular insights into the non-recombining nature of the spinach male-determining region

QTL analysis of femaleness in monoecious spinach and fine mapping of a major QTL using an updated version of chromosome-scale pseudomolecules

Although spinach is predominantly dioecious, monoecious plants with varying proportions of female and male flowers are also present. Recently, monoecious inbred lines with highly female and male conditions have been preferentially used as parents for F1-hybrids, rather than dioecious lines. Accordingly, identifying the loci for monoecism is an important issue for spinach breeding. We here used long-read sequencing and Hi-C technology to construct SOL_r2.0_pseudomolecule, a set of six pseudomolecules of spinach chromosomes (total length: 879.2 Mb; BUSCO complete 97.0%) that are longer and more genetically complete than our previous version of pseudomolecules (688.0 Mb; 81.5%). Three QTLs, qFem2.1, qFem3.1, and qFem6.1, responsible for monoecism were mapped to SOL_r2.0_pseudomolecule. qFem3.1 had the highest LOD score and corresponded to the M locus, which was previously identified as a determinant of monoecious expression, by genetic analysis of progeny from female and monoecious plants. The other QTLs were shown to modulate the ratio of female to male flowers in monoecious plants harboring a dominant allele of the M gene. Our findings will enable breeders to efficiently produce highly female- and male-monoecious parental lines for F1-hybrids by pyramiding the three QTLs. Through fine-mapping, we narrowed the candidate region for the M locus to a 19.5 kb interval containing three protein-coding genes and one long non-coding RNA gene. Among them, only RADIALIS-like-2a showed a higher expression in the reproductive organs, suggesting that it might play a role in reproductive organogenesis. However, there is no evidence that it is involved in the regulation of stamen and pistil initiation, which are directly related to the floral sex differentiation system in spinach. Given that auxin is involved in reproductive organ formation in many plant species, genes related to auxin transport/response, in addition to floral organ formation, were identified as candidates for regulators of floral sex-differentiation from qFem2.1 and qFem6.1.

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.

Accumulation of retrotransposons contributes to W chromosome differentiation in the willow beauty Peribatodes rhomboidaria (Lepidoptera: Geometridae)

The W chromosome of Lepidoptera is typically gene-poor, repeat-rich and composed of heterochromatin. Pioneering studies investigating this chromosome reported an abundance of mobile elements. However, the actual composition of the W chromosome varies greatly between species, as repeatedly demonstrated by comparative genomic hybridization (CGH) or fluorescence in situ hybridization (FISH). Here we present an analysis of repeats on the W chromosome in the willow beauty, Peribatodes rhomboidaria (Geometridae), a species in which CGH predicted an abundance of W-enriched or W-specific sequences. Indeed, comparative analysis of male and female genomes using RepeatExplorer identified ten putative W chromosome-enriched repeats, most of which are LTR or LINE mobile elements. We analysed the two most abundant: PRW LINE-like and PRW Bel-Pao. The results of FISH mapping and bioinformatic analysis confirmed their enrichment on the W chromosome, supporting the hypothesis that mobile elements are the driving force of W chromosome differentiation in Lepidoptera. As the W chromosome is highly underrepresented in chromosome-level genome assemblies of Lepidoptera, this recently introduced approach, combining bioinformatic comparative genome analysis with molecular cytogenetics, provides an elegant tool for studying this elusive and rapidly evolving part of the genome.

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.

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.

Do Ty3/Gypsy Transposable Elements Play Preferential Roles in Sex Chromosome Differentiation?

Transposable elements (TEs) comprise a substantial portion of eukaryotic genomes. They have the unique ability to integrate into new locations and serve as the main source of genomic novelties by mediating chromosomal rearrangements and regulating portions of functional genes. Recent studies have revealed that TEs are abundant in sex chromosomes. In this review, we propose evolutionary relationships between specific TEs, such as Ty3/Gypsy, and sex chromosomes in different lineages based on the hypothesis that these elements contributed to sex chromosome differentiation processes. We highlight how TEs can drive the dynamics of sex-determining regions via suppression recombination under a selective force to affect the organization and structural evolution of sex chromosomes. The abundance of TEs in the sex-determining regions originates from TE-poor genomic regions, suggesting a link between TE accumulation and the emergence of the sex-determining regions. TEs are generally considered to be a hallmark of chromosome degeneration. Finally, we outline recent approaches to identify TEs and study their sex-related roles and effects in the differentiation and evolution of sex chromosomes.

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.

Identification of structural variation and polymorphisms of a sex co-segregating scaffold in spinach

Spinach is a common vegetable, and dioecy is maintained by a pair of XY sex chromosomes. Due to limited genomic resources and its highly repetitive genome, limited studies were conducted to investigate the genomic landscape of the region near sex-determining loci. In this study, we screened the structure variations (SVs) between Y-linked contigs and a 1.78-Mb X scaffold (Super_scaffold 66), which enabled the development of 12 sex co-segregating DNA markers. These markers were tested in one F₁ mapping population and 40 spinach accessions, which comprised 692 individual plants with the strong sex linkage pattern. In addition, we found that Super_scaffold 66 was highly repetitive along with the enriched LTR-RTs insertions and decreased microsatellite distribution compared with the rest genome, which matches extremely low gene density featured by only nine annotated genes. Synteny analysis between Y contigs and Superscaffold_66 revealed a 340-Kb accumulative Y contig (non-continuous) and a 500-Kb X counterpart along with SVs and wide-spread tandem duplications. Among the nine genes, one ABC transporter gene revealed noticeable SVs between Y contig and X counterpart, as an approximate 5-Kb recent Gypsy LTR-RT insertion in the Y-linked allele, but not the X allele. The gene paucity, SVs, and sex-linked polymorphisms attributed to the recombination suppression. We proposed that Super_scaffold 66 is part of the non-recombining region containing the sex determination genes. The spread of 12 sex co-segregating markers from this 1.78 Mb genomic region indicated the existence and expansion of sex determination region during progression of the Y chromosome.

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.

Cloning and physical localization of male-biased repetitive DNA sequences in Spinacia oleracea (Amaranthaceae)

Spinach (Spinacia oleracea Linnaeus, 1753) is an ideal material for studying molecular mechanisms of early-stage sex chromosome evolution in dioecious plants. Degenerate oligonucleotide-primed polymerase chain reaction (DOP-PCR) technique facilitates the retrotransposon-relevant studies by enriching specific repetitive DNA sequences from a micro-dissected single chromosome. We conducted genomic subtractive hybridization to screen sex-biased DNA sequences by using the DOP-PCR amplification products of micro-dissected spinach Y chromosome. The screening yielded 55 male-biased DNA sequences with 30 576 bp in length, of which, 32 DNA sequences (12 049 bp) contained repeat DNA sequences, including LTR/Copia, LTR/Gypsy, simple repeats, and DNA/CMC-EnSpm. Among these repetitive DNA sequences, four DNA sequences that contained a fragment of Ty3-gypsy retrotransposons (SP73, SP75, SP76, and SP77) were selected as fluorescence probes to hybridization on male and female spinach karyotypes. Fluorescence in situ hybridization (FISH) signals of SP73 and SP75 were captured mostly on the centromeres and their surrounding area for each homolog. Hybridization signals primarily appeared near the putative centromeres for each homologous chromosome pair by using SP76 and SP77 probes for FISH, and sporadic signals existed on the long arms. Results can be served as a basis to study the function of repetitive DNA sequences in sex chromosome evolution in spinach.

Transcriptome profiling of differentially expressed genes of male and female inflorescences in spinach (Spinacia oleracea L.)

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.

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.

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.

Young sex chromosomes in plants and animals

A major reason for studying plant sex chromosomes is that they may often be 'young' systems. There is considerable evidence for the independent evolution of separate sexes within plant families or genera, in some cases showing that the maximum possible time during which their sex-determining genes have existed must be much shorter than those of several animal taxa. Consequently, their sex-linked regions could either have evolved soon after genetic sex determination arose or considerably later. Plants, therefore, include species with both young and old systems. I review several questions about the evolution of sex-determining systems and sex chromosomes that require studies of young systems, including: the kinds of mutations involved in the transition to unisexual reproduction from hermaphroditism or monoecy (a form of functional hermaphroditism); the times when they arose; and the extent to which the properties of sex-linked regions of genomes reflect responses to new selective situations created by the presence of a sex-determining locus. I also evaluate which questions are best studied in plants, vs other suitable candidate organisms. Studies of young plant systems can help understand general evolutionary processes that are shared with the sex chromosomes of other organisms.

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

Dioecy has evolved recently and independently from cosexual populations in many angiosperm lineages, providing opportunities to understand the evolutionary process underlying this transition. Spinach (Spinacia oleracea) is a dioecious plant with homomorphic sex chromosomes (XY). Although most of the spinach Y chromosome recombines with the X chromosome, a region around the male-determining locus on Y does not recombine with its X counterpart, suggesting that this region might be related to the evolution of dioecy in the species. To identify genes located in the non-recombining region (MSY, male-specific region of Y), RNA-seq analysis of male and female progeny plants (eight each) from a sib-cross of a dioecious line was performed. We discovered only 354 sex-chromosomal SNPs in 219 transcript sequences (genes). We randomly selected 39 sex-chromosomal genes to examine the reproducibility of the RNA-seq results and observed tight linkage to the male-determining locus in a spinach segregating population (140 individuals). Further analysis using a large-scale population (>1400) and over 100 spinach germplasm accessions and cultivars showed that SNPs in at least 12 genes are fully linked to the male-determining locus, suggesting that the genes reside in the spinach MSY. Synonymous substitution rates of the MSY genes and X homologues predict a recent divergence (0.40 ± 0.08 Mya). Furthermore, synonymous divergence between spinach and its wild relative (S. tetrandra), whose sex chromosomes (XY) originated from a common ancestral chromosome, predicted that the species diverged around 5.7 Mya. Assuming that dioecy in Spinacia evolved before speciation within the genus and has a monophyletic origin, our data suggest that recombination around the spinach sex-determining locus might have stopped significantly later than the evolution of dioecy in Spinacia.

Development of an X-specific marker and identification of YY individuals in spinach

Spinach is a popular vegetable native to central and western Asia. It is dioecious with a pair of nascent sex chromosomes. The difficulties of working with the non-recombining sex determination region of XY individuals have hindered the progress toward sequencing sex chromosomes of most dioecious species. Here we present important advances toward characterizing the non-recombining sex chromosomes in spinach. Of nearly 400 spinach accessions screened, we identified a single accession of spinach in which androdioecious XY individuals segregate YY spinach. The male and female genomes of the spinach cultivar Shami and USDA accession PI 664497 were sequenced at 12-17 × coverage. X-specific sequences were identified by comparing the depth of coverage differences between male and female alignments to a female draft genome. YY individuals were used as a negative control to validate X-specific markers found by depth of coverage analysis. Of 19 possible X chromosome sequences found by depth of coverage analysis, one was verified to be X-specific by a PCR-based marker, SpoX, which amplified genomic DNA from XX and XY, but not YY templates. Androdioecious XY individuals of accession PI 217425 (Cornell #9) were used to develop inbred lines, and at S7 generation, all XY individuals were androdioecious and all YY individuals were pure male. The sex reversal of the XY mutant to hermaphrodite is strong evidence that the sex chromosomes in spinach have a two-gene sex determination system. These results are crucial towards sequencing the X and Y chromosomes to advance sex chromosome research in spinach.