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.

Lineage-specific amplification and epigenetic regulation of LTR-retrotransposons contribute to the structure, evolution, and function of Fabaceae species.

BACKGROUND: Long terminal repeat (LTR)-retrotransposons (LTR-RTs) are ubiquitous and make up the majority of nearly all sequenced plant genomes, whereas their pivotal roles in genome evolution, gene expression regulation as well as their epigenetic regulation are still not well understood, especially in a large number of closely related species. RESULTS: Here, we analyzed the abundance and dynamic evolution of LTR-RTs in 54 species from an economically and agronomically important family, Fabaceae, and also selected two representative species for further analysis in expression of associated genes, transcriptional activity and DNA methylation patterns of LTR-RTs. Annotation results revealed highly varied proportions of LTR-RTs in these genomes (5.1%~68.4%) and their correlation with genome size was highly positive, and they were significantly contributed to the variance in genome size through species-specific unique amplifications. Almost all of the intact LTR-RTs were inserted into the genomes 4 Mya (million years ago), and more than 50% of them were inserted in the last 0.5 million years, suggesting that recent amplifications of LTR-RTs were an important force driving genome evolution. In addition, expression levels of genes with intronic, promoter, and downstream LTR-RT insertions of Glycine max and Vigna radiata, two agronomically important crops in Fabaceae, showed that the LTR-RTs located in promoter or downstream regions suppressed associated gene expression. However, the LTR-RTs within introns promoted gene expression or had no contribution to gene expression. Additionally, shorter and younger LTR-RTs maintained higher mobility and transpositional potential. Compared with the transcriptionally silent LTR-RTs, the active elements showed significantly lower DNA methylation levels in all three contexts. The distributions of transcriptionally active and silent LTR-RT methylation varied across different lineages due to the position of LTR-RTs located or potentially epigenetic regulation. CONCLUSION: Lineage-specific amplification patterns were observed and higher methylation level may repress the activity of LTR-RTs, further influence evolution in Fabaceae species. This study offers valuable clues into the evolution, function, transcriptional activity and epigenetic regulation of LTR-RTs in Fabaceae genomes.

A rapid method for assembly of single chromosome and identification of sex determination region based on single-chromosome sequencing.

The sex-determining-region (SDR) may offer the best prospects for studying sex-determining gene, recombination suppression, and chromosome heteromorphism. However, current progress of SDR identification and cloning showed following shortcomings: large near-isogenic lines need to be constructed, and a relatively large population is needed; the cost of whole-genome sequencing and assembly is high. Herein, the X/Y chromosomes of Spinacia oleracea L. subsp. turkestanica were successfully microdissected and assembled using single-chromosome sequencing. The assembly length of X and Y chromosome is c. 192.1 and 195.2 Mb, respectively. Three large inversions existed between X and Y chromosome. The SDR size of X and Y chromosome is c. 13.2 and 24.1 Mb, respectively. MSY region and six male-biased genes were identified. A Y-chromosome-specific marker in SDR was constructed and used to verify the chromosome assembly quality at cytological level via fluorescence in situ hybridization. Meanwhile, it was observed that the SDR located on long arm of Y chromosome and near the centromere. Overall, a technical system was successfully established for rapid cloning the SDR and it is also applicable to rapid assembly of specific chromosome in other plants. Furthermore, this study laid a foundation for studying the molecular mechanism of sex chromosome evolution in spinach.

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.

Impact of LTR-Retrotransposons on Genome Structure, Evolution, and Function in Curcurbitaceae Species.

Long terminal repeat (LTR)-retrotransposons (LTR-RTs) comprise a major portion of many plant genomes and may exert a profound impact on genome structure, function, and evolution. Although many studies have focused on these elements in an individual species, their dynamics on a family level remains elusive. Here, we investigated the abundance, evolutionary dynamics, and impact on associated genes of LTR-RTs in 16 species in an economically important plant family, Cucurbitaceae. Results showed that full-length LTR-RT numbers and LTR-RT content varied greatly among different species, and they were highly correlated with genome size. Most of the full-length LTR-RTs were amplified after the speciation event, reflecting the ongoing rapid evolution of these genomes. LTR-RTs highly contributed to genome size variation via species-specific distinct proliferations. The Angela and Tekay lineages with a greater evolutionary age were amplified in Trichosanthes anguina, whereas a recent activity burst of Reina and another ancient round of Tekay activity burst were examined in Sechium edule. In addition, Tekay and Retand lineages belonging to the Gypsy superfamily underwent a recent burst in Gynostemma pentaphyllum. Detailed investigation of genes with intronic and promoter LTR-RT insertion showed diverse functions, but the term of metabolism was enriched in most species. Further gene expression analysis in G.pentaphyllum revealed that the LTR-RTs within introns suppress the corresponding gene expression, whereas the LTR-RTs within promoters exert a complex influence on the downstream gene expression, with the main function of promoting gene expression. This study provides novel insights into the organization, evolution, and function of LTR-RTs in Cucurbitaceae genomes.

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.

Cytogenetic and genomic organization analyses of chloroplast DNA invasions in the nuclear genome of Asparagus officinalis L. provides signatures of evolutionary complexity and informativity in sex chromosome evolution.

BACKGROUND: The transfer of chloroplast DNA into nuclear genome is a common process in plants. These transfers form nuclear integrants of plastid DNAs (NUPTs), which are thought to be driving forces in genome evolution, including sex chromosome evolution. In this study, NUPTs in the genome of a dioecious plant Asparagus officinalis L. were systematically analyzed, in order to investigate the characteristics of NUPTs in the nuclear genome and the relationship between NUPTs and sex chromosome evolution in this species. RESULTS: A total of 3155 NUPT insertions were detected, and they represented approximated 0.06% of the nuclear genome. About 45% of the NUPTs were organized in clusters. These clusters were derived from various evolutionary events. The Y chromosome contained the highest number and largest proportion of NUPTs, suggesting more accumulation of NUPTs on sex chromosomes. NUPTs were distributed widely in all of the chromosomes, and some regions preferred these insertions. The highest density of NUPTs was found in a 47 kb region in the Y chromosome; more than 75% of this region was occupied by NUPTs. Further cytogenetic and sequence alignment analysis revealed that this region was likely the centromeric region of the sex chromosomes. On the other hand, the male-specific region of the Y chromosome (MSY) and the adjacent regions did not have NUPT insertions. CONCLUSIONS: These results indicated that NUPTs were involved in shaping the genome of A. officinalis through complicated process. NUPTs may play important roles in the centromere shaping of the sex chromosomes of A. officinalis, but were not implicated in MSY formation.

Genome-wide characterization of microsatellites and genetic diversity assessment of spinach in the Chinese germplasm collection.

Spinach is a nutritional leafy green vegetable, and it also serves as a model species for studying sex chromosome evolution. Genetic marker development and genome structure analysis are important in breeding practice and theoretical evolution studies of spinach. In this study, the frequency and distribution of different microsatellites in the recently released draft spinach genome were characterized. A total of 261,002 perfect microsatellites were identified (estimated frequency: ~262.1 loci/Mbp). The most abundant microsatellites were tetranucleotide and trinucleotide, accounting for 33.2% and 27.7% of the total number of microsatellites, respectively. A total of 105 primer pairs were designed and screened, and 34 were polymorphic among the detected spinach cultivars. Combined with seven primer sets developed previously, 41 primer pairs were used to investigate genetic diversity among 43 spinach cultivars in China. The average polymorphism information content value of the 41 markers was 0.43, representing an intermediate level. The spinach cultivars had a low genetic diversity, and no detectable common factors were shared by each group in the UPGMA dendrogram. This study's findings facilitate further investigations on the organization of the microsatellites in spinach genome and provide clues for future breeding applications of spinach in China.

Comparative transcriptome analysis reveals differentially expressed genes associated with sex expression in garden asparagus (Asparagus officinalis).

BACKGROUND: Garden asparagus (Asparagus officinalis) is a highly valuable vegetable crop of commercial and nutritional interest. It is also commonly used to investigate the mechanisms of sex determination and differentiation in plants. However, the sex expression mechanisms in asparagus remain poorly understood. RESULTS: De novo transcriptome sequencing via Illumina paired-end sequencing revealed more than 26 billion bases of high-quality sequence data from male and female asparagus flower buds. A total of 72,626 unigenes with an average length of 979 bp were assembled. In comparative transcriptome analysis, 4876 differentially expressed genes (DEGs) were identified in the possible sex-determining stage of female and male/supermale flower buds. Of these DEGs, 433, including 285 male/supermale-biased and 149 female-biased genes, were annotated as flower related. Of the male/supermale-biased flower-related genes, 102 were probably involved in anther development. In addition, 43 DEGs implicated in hormone response and biosynthesis putatively associated with sex expression and reproduction were discovered. Moreover, 128 transcription factor (TF)-related genes belonging to various families were found to be differentially expressed, and this finding implied the essential roles of TF in sex determination or differentiation in asparagus. Correlation analysis indicated that miRNA-DEG pairs were also implicated in asparagus sexual development. CONCLUSIONS: Our study identified a large number of DEGs involved in the sex expression and reproduction of asparagus, including known genes participating in plant reproduction, plant hormone signaling, TF encoding, and genes with unclear functions. We also found that miRNAs might be involved in the sex differentiation process. Our study could provide a valuable basis for further investigations on the regulatory networks of sex determination and differentiation in asparagus and facilitate further genetic and genomic studies on this dioecious species.

Identification of male-specific AFLP and SCAR markers in the dioecious plant Humulus scandens.

In this study, 17 male-specific amplified fragment length polymorphism (AFLP) markers were identified between male and female Humulus scandens plants. BLAST analysis revealed that 7 of the 17 sex-linked sequences were highly similar to retrotransposons. Two stable male-specific sequence-characterized amplified regions (SCAR) markers were developed. These AFLP and SCAR markers are novel molecular probes that can be used efficiently to identify the genetic gender of H. scandens and may provide a basis for further investigations on the evolution of sex chromosomes.

Repetitive sequences and epigenetic modification: inseparable partners play important roles in the evolution of plant sex chromosomes.

MAIN CONCLUSION: The present review discusses the roles of repetitive sequences played in plant sex chromosome evolution, and highlights epigenetic modification as potential mechanism of repetitive sequences involved in sex chromosome evolution. Sex determination in plants is mostly based on sex chromosomes. Classic theory proposes that sex chromosomes evolve from a specific pair of autosomes with emergence of a sex-determining gene(s). Subsequently, the newly formed sex chromosomes stop recombination in a small region around the sex-determining locus, and over time, the non-recombining region expands to almost all parts of the sex chromosomes. Accumulation of repetitive sequences, mostly transposable elements and tandem repeats, is a conspicuous feature of the non-recombining region of the Y chromosome, even in primitive one. Repetitive sequences may play multiple roles in sex chromosome evolution, such as triggering heterochromatization and causing recombination suppression, leading to structural and morphological differentiation of sex chromosomes, and promoting Y chromosome degeneration and X chromosome dosage compensation. In this article, we review the current status of this field, and based on preliminary evidence, we posit that repetitive sequences are involved in sex chromosome evolution probably via epigenetic modification, such as DNA and histone methylation, with small interfering RNAs as the mediator.

Genome-wide identification and validation of simple sequence repeats (SSRs) from Asparagus officinalis.

Garden asparagus (Asparagus officinalis), an important vegetable cultivated worldwide, can also serve as a model dioecious plant species in the study of sex determination and sex chromosome evolution. However, limited DNA marker resources have been developed and used for this species. To expand these resources, we examined the DNA sequences for simple sequence repeats (SSRs) in 163,406 scaffolds representing approximately 400 Mbp of the A. officinalis genome. A total of 87,576 SSRs were identified in 59,565 scaffolds. The most abundant SSR repeats were trinucleotide and tetranucleotide, accounting for 29.2 and 29.1% of the total SSRs, respectively, followed by di-, penta-, hexa-, hepta-, and octanucleotides. The AG motif was most common among dinucleotides and was also the most frequent motif in the entire A. officinalis genome, representing 14.7% of all SSRs. A total of 41,917 SSR primers pairs were designed to amplify SSRs. Twenty-two genomic SSR markers were tested in 39 asparagus accessions belonging to ten cultivars and one accession of Asparagus setaceus for determination of genetic diversity. The intra-species polymorphism information content (PIC) values of the 22 genomic SSR markers were intermediate, with an average of 0.41. The genetic diversity between the ten A. officinalis cultivars was low, and the UPGMA dendrogram was largely unrelated to cultivars. It is here suggested that the sex of individuals is an important factor influencing the clustering results. The information reported here provides new information about the organization of the microsatellites in A. officinalis genome and lays a foundation for further genetic studies and breeding applications of A. officinalis and related species.

Genes encoding Δ(8)-sphingolipid desaturase from various plants: identification, biochemical functions, and evolution.

∆(8)-sphingolipid desaturase catalyzes the C8 desaturation of a long chain base, which is the characteristic structure of various complex sphingolipids. The genes of 20 ∆(8)-sphingolipid desaturases from 12 plants were identified and functionally detected by using Saccharomyces cerevisiae system to elucidate the relationship between the biochemical function and evolution of this enzyme. Results showed that the 20 genes all can encode a functional ∆(8)-sphingolipid desaturase, which catalyzes different ratios of two products, namely, 8(Z) and 8(E)-C18-phytosphingenine. The coded enzymes could be divided into two groups on the basis of biochemical functions: ∆(8)-sphingolipid desaturase with a preference for an E-isomer product and ∆(8)-sphingolipid desaturase with a preference for a Z-isomer product. The conversion rate of the latter was generally lower than that of the former. Phylogenetic analysis revealed that the 20 desaturases could also be clustered into two groups, and this grouping is consistent with that of the biochemical functions. Thus, the biochemical function of ∆(8)-sphingolipid desaturase is correlated with its evolution. The two groups of ∆(8)-sphingolipid desaturases could arise from distinct ancestors in higher plants. However, they might have initially evolved from ∆(8)-sphingolipid desaturases in lower organisms, such as yeasts, which can produce E-isomer products only. Furthermore, almost all of the transgenic yeasts harboring ∆(8)-sphingolipid desaturase genes exhibit an improvement in aluminum tolerance. Our study provided new insights into the biochemical function and evolution of ∆(8)-sphingolipid desaturases in plants.

Isolation of differentially expressed sex genes in garden asparagus using suppression subtractive hybridization.

Garden asparagus (Asparagus officinalis L.) is a dioecious species whose male and female flowers are found in separate unisexual individuals. A region called the M-locus, located on a pair of homomorphic sex chromosomes, controls sexual dimorphism in asparagus. To date, no sex determining gene has been isolated from asparagus. To identify more genes involved in flower development in asparagus, subtractive hybridization library of male flowers in asparagus was constructed by suppression subtraction hybridization. A total of 107 expressed sequence tags (ESTs) were identified. BLASTX analysis showed that the library contained several genes that could be related to flower development. The expression patterns of seven selected genes believed to be involved in the development of asparagus male flower were further analyzed by semi-quantitative or real-time reverse-transcription polymerase chain reaction (RT-PCR). Results showed that AOEST4-5, AOEST12-40, and AOEST13-38 were strongly expressed in the male flower stage, whereas no transcript level of AOEST13-38 was detected in the female flower stage. The expression levels of AOEST13-87, AOEST13-92, AOEST13-40, and AOEST18-87 in the male flower stage were also higher than those in the female flower stage, although these transcripts were also expressed in other tissues. The identified genes can provide a strong starting point for further studies on the underlying molecular differences between the male and female flowers of asparagus.

[Role of transposons in origin and evolution of plant XY sex chromosomes].

The XY sex-determination system is crucial for plant reproduction. However, little is known about the mechanism of the origin and evolution of the XY sex chromosomes. It has been believed that a pair of autosomes is evolved to produce young sex chromosomes (neo-X chromosome and neo-Y chromosome) by loss of function or gain of function mutation, which influences the development of pistil or stamen. With the aggravation of the recombination suppression between neo-X and neo-Y and consequent expanding of the non-recombination region, the proto-sex chromosomes were finally developed to heteromorphic sex chromosomes. Accumulation of repetitive sequences and DNA methylation were probably involved in this process. Transposons, as the most abundant repetitive sequences in the genome, might be the initial motivation factors for the evolution of sex chromosome. Moreover, transposons may also increase heterochromatin expansion and recombination suppression of sex chromosome by local epigenetics modification. In this review, we summarize the function of transposon accumulation and the relationship between transposon and heterochromatization in the evolution of plant sex chromosome.

DOP-PCR based painting of rye chromosomes in a wheat background.

To determine the appropriateness of chromosome painting for identifying genomic elements in rye, we microdissected the 1R and 1RS chromosomes from rye (Secale cereale L. var. King II) and wheat-rye addition line 1RS, respectively. Degenerate oligonucleotide primed - polymerase chain reaction (DOP-PCR) amplification of 1R and 1RS products from dissected chromosomes were used as probes to hybridize to metaphase chromosomes of rye, wheat-rye addition lines 1R and 1RS, translocation line 1RS.1BL, and allohexaploid triticale. The results showed that (i) the hybridization signal distribution patterns on rye chromosomes using 1R-derived DOP-PCR products as the probe were similar to those using 1RS-derived DOP-PCR products as the probe; (ii) 1R and (or) 1RS could not be distinguished from other rye chromosomes solely by the hybridization patterns using 1R- and (or) 1RS-derived DOP-PCR products as the probe; (iii) rye chromosomes and (or) rye chromosome fragments could be clearly identified in wheat-rye hybrids using either 1R- or 1RS-derived DOP-PCR products as the probe and could be more accurate in the nontelomeric region than using genomic in situ hybridization (GISH). Our results suggested that 1R- and (or) 1RS-derived DOP-PCR products contain many repetitive DNA sequences, are similar on different rye chromosomes, are R-genome specific, and can be used to identify rye chromosomes and chromosome fragments in wheat-rye hybrids. Our research widens the application range of chromosome painting in plants.

Detection of genome DNA methylation change in spinach induced by 5-azaC.

DNA methylation has been implicated in the regulation of gene expression, genome imprinting, and chromatin remodeling in eukaryotes. In this study, we analyzed possible alterations in levels and patterns of cytosine methylation in male and female spinach plants after treatment with demethylation agent 5-azacytidine (5-azaC) using two methods: (1) direct determination of 5-methylcytidine (5 mC) amounts in genomic DNA by high-performance liquid chromatography (HPLC) separation and quantification of nucleosides and (2) methylation-sensitive inter-simple sequence repeat (MS-ISSR) technique. HPLC analysis revealed that the DNA methylation events in male and female spinach leaves markedly decreased upon 30 µM 5-azaC treatment, and the methylation level gradually decreased with the increase in 5-azaC concentration. To study the altered DNA methylation patterns in spinach after 5-azaC treatment, untreated and 500 µM 5-azaC-treated samples were analyzed by MS-ISSR assay. A total of 385 informative profiles were resolved using 35 ISSR primer sets. MS-ISSR analysis showed various altered methylation patterns between untreated and 5-azaC-treated spinach plants. These alterations were mainly demethylation events, which were largely consistent with the HPLC results. Both HPLC and MS-ISSR analyses showed that the changes in DNA methylation levels and patterns were similar in male and female spinach leaves, which implies that sex was not the main factor influencing DNA methylation levels and patterns in the vegetative organs of spinach. This study could provide a molecular basis of the altered DNA methylation induced by 5-azaC, and lay a foundation for further investigation of the relationship between methylation and sex determination and development in this dioecious plant spinach.

Comparative analysis of gene expression by microarray analysis of male and female flowers of Asparagus officinalis.

To identify rapidly a number of genes probably involved in sex determination and differentiation of the dioecious plant Asparagus officinalis, gene expression profiles in early flower development for male and female plants were investigated by microarray assay with 8,665 probes. In total, 638 male-biased and 543 female-biased genes were identified. These genes with biased-expression for male and female were involved in a variety of processes associated with molecular functions, cellular components, and biological processes, suggesting that a complex mechanism underlies the sex development of asparagus. Among the differentially expressed genes involved in the reproductive process, a number of genes associated with floral development were identified. Reverse transcription-PCR was performed for validation, and the results were largely consistent with those obtained by microarray analysis. The findings of this study might contribute to understanding of the molecular mechanisms of sex determination and differentiation in dioecious asparagus and provide a foundation for further studies of this plant.

Microdissection and painting of the Y chromosome in spinach (Spinacia oleracea).

Spinach has long been used as a model for genetic and physiological studies of sex determination and expression. Although trisomic analysis from a cross between diploid and triploid plants identified the XY chromosome as the largest chromosome, no direct evidence has been provided to support this at the molecular level. In this study, the largest chromosomes of spinach from mitotic metaphase spreads were microdissected using glass needles. Degenerate oligonucleotide primed polymerase chain reaction was used to amplify the dissected chromosomes. The amplified products from the Y chromosome were identified using the male-specific marker T11A. For the first time, the largest spinach chromosome was confirmed to be a sex chromosome at the molecular level. PCR products from the isolated chromosomes were used in an in situ probe mixture for painting the Y chromosome. The fluorescence signals were mainly distributed on all chromosomes and four pair of weaker punctate fluorescence signal sites were observed on the terminal region of two pair of autosomes. These findings provide a foundation for the study of sex chromosome evolution in spinach.

Development and applications of a collection of single copy gene-based cytogenetic DNA markers in garden asparagus.

Garden asparagus (Asparagus officinalis, 2n = 2x = 20 chromosomes) is an important dioecious vegetable crop and a model species for studying sex chromosome formation and evolution. However, few molecular cytogenetic studies on garden asparagus have been reported because of its small metaphase chromosomes, the scarcity of distinguished cytogenetic markers, and the high content of repetitive sequences. In this study, a set of single copy genes free of repetitive sequences with sizes ranging from 4.3 kb to 8.2 kb were screened and used as probes for fluorescence in situ hybridization (FISH) to identify individual chromosomes of garden asparagus. The chromosome-specific signal distribution patterns of these probes enabled the distinguishment of each pair of chromosomes. The sequence assembly and cytogenetic map were successfully integrated, and the results confirmed that the chromosome 1 representing the sex chromosome in the genome assembly is chromosome 5 in the karyotype analysis. The cytogenetic identification of the male-specific region of the Y chromosome (MSY) was implemented using a mixed probe derived from a number of MSY-specific single copy sequences. In addition, the chromosome orthologous relationship between garden asparagus (A1-A10, karyotypic analysis) and its hermaphrodite close relative, A. setaceus (B1-B10, karyotypic analysis), was analyzed using this collection of chromosome-specific cytological markers. The results showed that B3 is the ortholog of sex chromosome A5 and thus may represent the ancestral autosome of the current sex chromosome in garden asparagus. Chromosomes B5, B4, B1, B8, B7, and B9 are the orthologs of A2, A3, A4, A7, A8, and A10, respectively. The chromosome identification, cytogenetic recognition of MSY, and the orthologous relationship analysis between garden asparagus and A. setaceus are valuable for the further investigation of the sex chromosome emergence and evolutionary mechanism of garden asparagus and genome structure evolution in the Asparagus genus.