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

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

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

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

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

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

Chromosome-level assembly of the Neolamarckia cadamba genome provides insights into the evolution of cadambine biosynthesis.

Neolamarckia cadamba (Roxb.), a close relative of Coffea canephora and Ophiorrhiza pumila, is an important traditional medicine in Southeast Asia. Three major glycosidic monoterpenoid indole alkaloids (MIAs), cadambine and its derivatives 3ß-isodihydrocadambine and 3ß-dihydrocadambine, accumulate in the bark and leaves, and exhibit antimalarial, antiproliferative, antioxidant, anticancer and anti-inflammatory activities. Here, we report a chromosome-scale N. cadamba genome, with 744.5 Mb assembled into 22 pseudochromosomes with contig N50 and scaffold N50 of 824.14 Kb and 29.20 Mb, respectively. Comparative genomic analysis of N. cadamba with Co. canephora revealed that N. cadamba underwent a relatively recent whole-genome duplication (WGD) event after diverging from Co. canephora, which contributed to the evolution of the MIA biosynthetic pathway. We determined the key intermediates of the cadambine biosynthetic pathway and further showed that NcSTR1 catalyzed the synthesis of strictosidine in N. cadamba. A new component, epoxystrictosidine (C27H34N2O10, m/z 547.2285), was identified in the cadambine biosynthetic pathway. Combining genome-wide association study (GWAS), population analysis, multi-omics analysis and metabolic gene cluster prediction, this study will shed light on the evolution of MIA biosynthetic pathway genes. This N. cadamba reference sequence will accelerate the understanding of the evolutionary history of specific metabolic pathways and facilitate the development of tools for enhancing bioactive productivity by metabolic engineering in microbes or by molecular breeding in plants.

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.

Nuclear Integrants of Organellar DNA Contribute to Genome Structure and Evolution in Plants.

The transfer of genetic material from the mitochondria and plastid to the nucleus gives rise to nuclear integrants of mitochondrial DNA (NUMTs) and nuclear integrants of plastid DNA (NUPTs). This frequently occurring DNA transfer is ongoing and has important evolutionary implications. In this review, based on previous studies and the analysis of NUMT/NUPT insertions of more than 200 sequenced plant genomes, we analyzed and summarized the general features of NUMTs/NUPTs and highlighted the genetic consequence of organellar DNA insertions. The statistics of organellar DNA integrants among various plant genomes revealed that organellar DNA-derived sequence content is positively correlated with the nuclear genome size. After integration, the nuclear organellar DNA could undergo different fates, including elimination, mutation, rearrangement, fragmentation, and proliferation. The integrated organellar DNAs play important roles in increasing genetic diversity, promoting gene and genome evolution, and are involved in sex chromosome evolution in dioecious plants. The integrating mechanisms, involving non-homologous end joining at double-strand breaks were also discussed.

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.

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.