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.

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.

Genome-wide analysis of transposable elements and satellite DNA in Humulus scandens, a dioecious plant with XX/XY1Y2 chromosomes.

Transposable elements (TEs) and satellite DNAs, two major categories of repetitive sequences, are expected to accumulate in non-recombining genome regions, including sex-linked regions, and contribute to sex chromosome evolution. The dioecious plant, Humulus scandens, can be used for studying the evolution of the XX/XY1Y2 sex chromosomes. In this study, we thoroughly examined the repetitive components of male and female H. scandens using next-generation sequencing data followed by bioinformatics analysis and florescence in situ hybridization (FISH). The H. scandens genome has a high overall repetitive sequence composition, 68.30% in the female and 66.78% in the male genome, with abundant long terminal repeat (LTR) retrotransposons (RTs), including more Ty3/Gypsy than Ty1/Copia elements, particularly two Ty3/Gypsy lineages, Tekay and Retand. Most LTR-RT lineages were found dispersed across the chromosomes, though CRM and Athila elements were predominately found within the centromeres and the pericentromeric regions. The Athila elements also showed clearly higher FISH signal intensities in the Y1 and Y2 chromosomes than in the X or autosomes. Three novel satellite DNAs were specifically distributed in the centromeric and/or telomeric regions, with markedly different distributions on the X, Y1, and Y2 chromosomes. Combined with FISH using satellite DNAs to stain chromosomes during meiotic diakinesis, we determined the synapsis pattern and distinguish pseudoautosomal regions (PARs). The results indicate that the XY1Y2 sex chromosomes of H. scandens might have originated from a centric fission event. This study improves our understanding of the repetitive sequence organization of H. scandens genome and provides a basis for further analysis of their chromosome evolution process.

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.

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.

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.

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.

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.

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.

Genomic landscape of mitochondrial DNA insertions in 23 bat genomes: characteristics, loci, phylogeny, and polymorphism.

The transfer of mitochondrial DNA to the nuclear genome gives rise to the nuclear DNA sequences of mitochondrial origin (NUMTs), considered as a driving force in genome evolution. In this study, NUMTs in 23 bat genomes were investigated and compared systematically. The results showed that NUMTs existed in 22 genomes except for Noctilio leporinus, suggesting that mitochondrial fragment insertion in the nuclear genome was a common event in bat genomes. However, remarkable variations in NUMTs number, cumulative length, and proportion in the nuclear genome were discovered across bat species. Further orthologous NUMT loci analysis of the Phyllostomidae family indicated that the NUMTs insertion events in bat genomes were homoplasy-free. The NUMTs were mainly inserted into the intergenic regions, particularly, co-localized with repetitive sequences (especially transposable elements). However, several NUMTs were inserted into genes, some of which were in the exon region of functional genes. One NUMT in the genome of Pteropus alecto surprisingly matched with cDNA of ATP8B3 that provided evidence of NUMTs with coding function. Phylogenic analysis on NUMTs originating from COXI and COXII loci highlighted 2 clusters of Yinpterochiroptera and Yangochiroptera for Chiroptera. Seven NUMTs from Rhinolophus ferrumequinum were amplified, and the sequencing results confirmed the reliability of the NUMT analysis. One of them was polymorphic for the presence or absence of the NUMT insertion, and each genotype of NUMT loci showed a distinct regional distribution pattern. The information obtained in this study provides novel insights into the NUMT organization and features in bat genomes and establishes a basis for further studying of the evolution of bat species.

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.

DNA methylation is involved in sexual differentiation and sex chromosome evolution in the dioecious plant garden asparagus.

DNA methylation is a crucial regulatory mechanism in many biological processes. However, limited studies have dissected the contribution of DNA methylation to sexual differentiation in dioecious plants. In this study, we investigated the variances in methylation and transcriptional patterns of male and female flowers of garden asparagus. Compared with male flowers, female flowers at the same stages showed higher levels of DNA methylation. Both male and female flowers gained DNA methylation globally from the premeiotic to meiotic stages. Detailed analysis revealed that the increased DNA methylation was largely due to increased CHH methylation. Correlation analysis of differentially expressed genes and differentially methylated regions suggested that DNA methylation might not have contributed to the expression variation of the sex-determining genes SOFF and TDF1 but probably played important roles in sexual differentiation and flower development of garden asparagus. The upregulated genes AoMS1, AoLAP3, AoAMS, and AoLAP5 with varied methylated CHH regions might have been involved in sexual differentiation and flower development of garden asparagus. Plant hormone signaling genes and transcription factor genes also participated in sexual differentiation and flower development with potential epigenetic regulation. In addition, the CG and CHG methylation levels in the Y chromosome were notably higher than those in the X chromosome, implying that DNA methylation might have been involved in Y chromosome evolution. These data provide insights into the epigenetic modification of sexual differentiation and flower development and improve our understanding of sex chromosome evolution in garden asparagus.

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.

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.

Chromosome-level genome assembly, annotation and evolutionary analysis of the ornamental plant Asparagus setaceus.

Asparagus setaceus is a popular ornamental plant cultivated in tropical and subtropical regions globally. Here, we constructed a chromosome-scale reference genome of A. setaceus to facilitate the investigation of its genome characteristics and evolution. Using a combination of Nanopore long reads, Illumina short reads, 10× Genomics linked reads, and Hi-C data, we generated a high-quality genome assembly of A. setaceus covering 710.15 Mb, accounting for 98.63% of the estimated genome size. A total of 96.85% of the sequences were anchored to ten superscaffolds corresponding to the ten chromosomes. The genome of A. setaceus was predicted to contain 28,410 genes, 25,649 (90.28%) of which were functionally annotated. A total of 65.59% of the genome was occupied by repetitive sequences, among which long terminal repeats were predominant (42.51% of the whole genome). Evolutionary analysis revealed an estimated divergence time of A. setaceus from its close relative A. officinalis of ~9.66 million years ago, and A. setaceus underwent two rounds of whole-genome duplication. In addition, 762 specific gene families, 96 positively selected genes, and 76 resistance (R) genes were detected and functionally predicted in A. setaceus. These findings provide new knowledge about the characteristics and evolution of the A. setaceus genome, and will facilitate comparative genetic and genomic research on the genus Asparagus.

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.

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.

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.

The landscape of transposable elements and satellite DNAs in the genome of a dioecious plant spinach (Spinacia oleracea L.).

BACKGROUND: Repetitive sequences, including transposable elements (TEs) and satellite DNAs, occupy a considerable portion of plant genomes. Analysis of the repeat fraction benefits the understanding of genome structure and evolution. Spinach (Spinacia oleracea L.), an important vegetable crop, is also a model dioecious plant species for studying sex determination and sex chromosome evolution. However, the repetitive sequences of the spinach genome have not been fully investigated. RESULTS: We extensively analyzed the repetitive components of draft spinach genome, especially TEs and satellites, by different strategies. A total of 16,002 full-length TEs were identified. Among the most abundant long terminal repeat (LTR) retrotransposons (REs), Copia elements were overrepresented compared with Gypsy ones. Angela was the most dominating Copia lineage; Ogre/Tat was the most abundant Gypsy lineage. The mean insertion age of LTR-REs was 1.42 million years; approximately 83.7% of these elements were retrotransposed during the last two million years. RepeatMasker totally masked about 64.05% of the spinach genome, with LTR-REs, non-LTR-REs, and DNA transposons occupying 49.2, 2.4, and 5.6%, respectively. Fluorescence in situ hybridization (FISH) analysis showed that most LTR-REs dispersed all over the chromosomes, by contrast, elements of CRM lineage were distributed at the centromeric region of all chromosomes. In addition, Ogre/Tat lineage mainly accumulated on sex chromosomes, and satellites Spsat2 and Spsat3 were exclusively located at the telomeric region of the short arm of sex chromosomes. CONCLUSIONS: We reliably annotated the TE fraction of the draft genome of spinach. FISH analysis indicates that Ogre/Tat lineage and the sex chromosome-specific satellites DNAs might participate in sex chromosome formation and evolution. Based on FISH signals of microsatellites, together with 45S rDNA, a fine karyotype of spinach was established. This study improves our knowledge of repetitive sequence organization in spinach genome and aids in accurate spinach karyotype construction.

Characterization of the complete chloroplast genome of Asparagus setaceus.

Asparagus setaceus is an important ornamental plant and is important to study the sex chromosome evolution in genus Asparagus. Here the complete chloroplast (cp) genome sequence of A. setaceus was determined. The cpDNA is 156,978 bp in length, contains a large single copy region of 85,311 bp and a small single copy region of 18,641 bp, which was separated by a pair of inverted repeat regions of 26,513 bp each. The cp genome contains 132 genes, including 88 protein-coding genes, 6 ribosomal RNA genes, and 38 transfer RNA genes. Phylogenetic analysis indicates that A. setaceus is closely related to Asparagus officinalis and Asparagus schoberioides.