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.

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.

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.

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.

The application of MutMap in forward genetic studies based on whole-genome sequencing.

Classical forward genetic analysis relies on construction of complicated progeny populations and development of many molecular markers for linkage analysis in genetic mapping, which is both time- and cost-consuming. The recently developed MutMap is a new forward genetic approach based on high-throughput next-generation sequencing technologies. It is more efficient and affordable than traditional methods. Moreover, new extended methods based on MutMap have been developed: MutMap+, which is based on self-crossing; MutMap-Gap, which is used to recognize the causative variations occurring in genome gap regions; QTL-seq, a method similar to MutMap for mapping quantitative trait loci. These methods are free from constructing complicated mapping population, genetic hybridization and linkage information. They have greatly accelerated the identification of genetic elements associated with interested phenotypic variation. Here, we review the basic principles of MutMap, and discuss their future applications in next generation sequencing-based forward genetic mapping and crop improvement.

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.

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.

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.

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.

[Preliminary study on molecular mechanism of lotus mutants induced by ion implantation].

Ion implantation, as a new biophysically mutagenic technique, has shown a great potential for horticultural plant breeding. Up to date, little is known about the mutation mechanism of ion implantation at the DNA level. To reveal the mutation effect of Fe+ ion implantation on Baiyangdian red lotus, the random amplified polymorphic DNA (RAPD) was used, and then the bands of mutants and the control in the radiation-sensitive sites were cloned to be sequenced for comparing their DNA sequences. The results indicated that the total base mutation rate of mutants was 0.87%, and there was different in the six mutants. The types of base changes included base transition, transversion, deletion, and insertion. Among the 159 base changes detected, the frequency of single base substitutions (61.01%) was higher than that of base deletions and insertions (38.99%), and the frequency of base transitions (44.65%) was 2.7 times of that of the base transversions (16.35%). The transitions between C and T accounted for largest proportion, A→G transitions and A→T transversions were also present at high frequency. Adenine, thymine, guanine or cytosine could be replaced by any of other three bases, except that there was no C → G substitution. However, thymine was more sensitive to the irradiation than other bases. In our study, we found many purine bases around the purine mutational sites, and many pyrimidine bases around the pyrimidine mutational sites. These will further help us to understand the mechanism of mutagenesis by ion implantation.

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.

[Role of repetitive sequence and heterochromatize in recombination suppression of plant sex chromosomes].

Suppression of recombination is the prerequisite for plant sex chromosome evolution from a pair of autosomes. Recombination suppression around the locus controlling sex determination results in sex chromosome degeneration and differentiation. Important events such as repetitive sequence accumulation, heterochromatize, and DNA methylation have relation to recombination suppression. Accumulation of repetitive DNA sequence, including transposable elements and satellite DNA, leads to primitive sex chromosome differentiated on morphological and molecular structure, and also gives rise to chromosome heterochromatize, and thus recombination between sex chromosomes was suppressed. Here, we re-viewed the advances in this field, meanwhile, the function of DNA methylation in recombination suppression was analyzed.

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.