The asparagus genome sheds light on the origin and evolution of a young Y chromosome

Profiles of phenolics and their synthetic pathways in Asparagus officinalis L

The synthetic pathways of some phenolics compounds in asparagus have been reported, however, the diversified phenolics compounds including their modification and transcription regulation remains unknown. Thus, multi-omics strategies were applied to detect the phenolics profiles, contents, and screen the key genes for phenolics biosynthesis and regulation in asparagus. A total of 437 compounds, among which 204 phenolics including 105 flavonoids and 82 phenolic acids were detected with fluctuated concentrations in roots (Rs), spears (Ss) and flowering twigs (Fs) of the both green and purple cultivars. Based on the detected phenolics profiles and contents correlated to the gene expressions of screened synthetic enzymes and regulatory TFs, a full phenolics synthetic pathway of asparagus was proposed for the first time, essential for future breeding of asparagus and scaled healthy phenolics production using synthetic biological strategies.

Genome assembly and resequencing shed light on evolution, population selection, and sex identification in Vernicia montana

Vernicia montana is a dioecious plant widely cultivated for high-quality tung oil production and ornamental purposes in the Euphorbiaceae family. The lack of genomic information has severely hindered molecular breeding for genetic improvement and early sex identification in V. montana. Here, we present a chromosome-level reference genome of a male V. montana with a total size of 1.29 Gb and a contig N50 of 3.69 Mb. Genome analysis revealed that different repeat lineages drove the expansion of genome size. The model of chromosome evolution in the Euphorbiaceae family suggests that polyploidization-induced genomic structural variation reshaped the chromosome structure, giving rise to the diverse modern chromosomes. Based on whole-genome resequencing data and analyses of selective sweep and genetic diversity, several genes associated with stress resistance and flavonoid synthesis such as CYP450 genes and members of the LRR-RLK family, were identified and presumed to have been selected during the evolutionary process. Genome-wide association studies were conducted and a putative sex-linked insertion and deletion (InDel) (Chr 2: 102 799 917-102 799 933 bp) was identified and developed as a polymorphic molecular marker capable of effectively detecting the gender of V. montana. This InDel is located in the second intron of VmBASS4, suggesting a possible role of VmBASS4 in sex determination in V. montana. This study sheds light on the genome evolution and sex identification of V. montana, which will facilitate research on the development of agronomically important traits and genomics-assisted breeding.

Identification, Elucidation and Deployment of a Cytoplasmic Male Sterility System for Hybrid Potato

Recent advances in diploid F1 hybrid potato breeding rely on the production of inbred lines using the S-locus inhibitor (Sli) gene. As a result of this method, female parent lines are self-fertile and require emasculation before hybrid seed production. The resulting F1 hybrids are self-fertile as well and produce many undesirable berries in the field. Utilization of cytoplasmic male sterility would eliminate the need for emasculation, resulting in more efficient hybrid seed production and male sterile F1 hybrids. We observed plants that completely lacked anthers in an F2 population derived from an interspecific cross between diploid S. tuberosum and S. microdontum. We studied the antherless trait to determine its suitability for use in hybrid potato breeding. We mapped the causal locus to the short arm of Chromosome 6, developed KASP markers for the antherless (al) locus and introduced it into lines with T and A cytoplasm. We found that antherless type male sterility is not expressed in T and A cytoplasm, proving that it is a form of CMS. We hybridized male sterile al/al plants with P cytoplasm with pollen from al/al plants with T and A cytoplasm and we show that the resulting hybrids set significantly fewer berries in the field. Here, we show that the antherless CMS system can be readily deployed in diploid F1 hybrid potato breeding to improve hybridization efficiency and reduce berry set in the field.

Genome assembly of autotetraploid Actinidia arguta highlights adaptive evolution and enables dissection of important economic traits

Actinidia arguta, the most widely distributed Actinidia species and the second cultivated species in the genus, can be distinguished from the currently cultivated Actinidia chinensis on the basis of its small and smooth fruit, rapid softening, and excellent cold tolerance. Adaptive evolution of tetraploid Actinidia species and the genetic basis of their important agronomic traits are still unclear. Here, we generated a chromosome-scale genome assembly of an autotetraploid male A. arguta accession. The genome assembly was 2.77 Gb in length with a contig N50 of 9.97 Mb and was anchored onto 116 pseudo-chromosomes. Resequencing and clustering of 101 geographically representative accessions showed that they could be divided into two geographic groups, Southern and Northern, which first diverged 12.9 million years ago. A. arguta underwent two prominent expansions and one demographic bottleneck from the mid-Pleistocene climate transition to the late Pleistocene. Population genomics studies using paleoclimate data enabled us to discern the evolution of the species' adaptation to different historical environments. Three genes (AaCEL1, AaPME1, and AaDOF1) related to flesh softening were identified by multi-omics analysis, and their ability to accelerate flesh softening was verified through transient expression assays. A set of genes that characteristically regulate sexual dimorphism located on the sex chromosome (Chr3) or autosomal chromosomes showed biased expression during stamen or carpel development. This chromosome-level assembly of the autotetraploid A. arguta genome and the genes related to important agronomic traits will facilitate future functional genomics research and improvement of A. arguta.

Current status of community resources and priorities for weed genomics research

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

Crocus genome reveals the evolutionary origin of crocin biosynthesis

Crocus sativus (saffron) is a globally autumn-flowering plant, and its stigmas are the most expensive spice and valuable herb medicine. Crocus specialized metabolites, crocins, are biosynthesized in distant species, Gardenia (eudicot) and Crocus (monocot), and the evolution of crocin biosynthesis remains poorly understood. With the chromosome-level Crocus genome assembly, we revealed that two rounds of lineage-specific whole genome triplication occurred, contributing important roles in the production of carotenoids and apocarotenoids. According to the kingdom-wide identification, phylogenetic analysis, and functional assays of carotenoid cleavage dioxygenases (CCDs), we deduced that the duplication, site positive selection, and neofunctionalization of Crocus-specific CCD2 from CCD1 members are responsible for the crocin biosynthesis. In addition, site mutation of CsCCD2 revealed the key amino acids, including I143, L146, R161, E181, T259, and S292 related to the catalytic activity of zeaxanthin cleavage. Our study provides important insights into the origin and evolution of plant specialized metabolites, which are derived by duplication events of biosynthetic genes.

A chromosome-level genome assembly of Agave hybrid NO.11648 provides insights into the CAM photosynthesis

The subfamily Agavoideae comprises crassulacean acid metabolism (CAM), C3, and C4 plants with a young age of speciation and slower mutation accumulation, making it a model crop for studying CAM evolution. However, the genetic mechanism underlying CAM evolution remains unclear because of lacking genomic information. This study assembled the genome of Agave hybrid NO.11648, a constitutive CAM plant belonging to subfamily Agavoideae, at the chromosome level using data generated from high-throughput chromosome conformation capture, Nanopore, and Illumina techniques, resulting in 30 pseudo-chromosomes with a size of 4.87 Gb and scaffold N50 of 186.42 Mb. The genome annotation revealed 58 841 protein-coding genes and 76.91% repetitive sequences, with the dominant repetitive sequences being the I-type repeats (Copia and Gypsy accounting for 18.34% and 13.5% of the genome, respectively). Our findings also provide support for a whole genome duplication event in the lineage leading to A. hybrid, which occurred after its divergence from subfamily Asparagoideae. Moreover, we identified a gene duplication event in the phosphoenolpyruvate carboxylase kinase (PEPCK) gene family and revealed that three PEPCK genes (PEPCK3, PEPCK5, and PEPCK12) were involved in the CAM pathway. More importantly, we identified transcription factors enriched in the circadian rhythm, MAPK signaling, and plant hormone signal pathway that regulate the PEPCK3 expression by analysing the transcriptome and using yeast one-hybrid assays. Our results shed light on CAM evolution and offer an essential resource for the molecular breeding program of Agave spp.

Young evolutionary origins of dioecy in the genus Asparagus

PREMISE

Dioecy (separate sexes) has independently evolved numerous times across the angiosperm phylogeny and is recently derived in many lineages. However, our understanding is limited regarding the evolutionary mechanisms that drive the origins of dioecy in plants. The recent and repeated evolution of dioecy across angiosperms offers an opportunity to make strong inferences about the ecological, developmental, and molecular factors influencing the evolution of dioecy, and thus sex chromosomes. The genus Asparagus (Asparagaceae) is an emerging model taxon for studying dioecy and sex chromosome evolution, yet estimates for the age and origin of dioecy in the genus are lacking.

METHODS

We use plastome sequences and fossil time calibrations in phylogenetic analyses to investigate the age and origin of dioecy in the genus Asparagus. We also review the diversity of sexual systems present across the genus to address contradicting reports in the literature.

RESULTS

We estimate that dioecy evolved once or twice approximately 2.78-3.78 million years ago in Asparagus, of which roughly 27% of the species are dioecious and the remaining are hermaphroditic with monoclinous flowers.

CONCLUSIONS

Our findings support previous work implicating a young age and the possibility of two origins of dioecy in Asparagus, which appear to be associated with rapid radiations and range expansion out of Africa. Lastly, we speculate that paleoclimatic oscillations throughout northern Africa may have helped set the stage for the origin(s) of dioecy in Asparagus approximately 2.78-3.78 million years ago.

Genetic insights into the dissolution of dioecy in diploid persimmon Diospyros oleifera Cheng

BACKGROUND

Dioecy, a sexual system of single-sexual (gynoecious/androecious) individuals, is rare in flowering plants. This rarity may be a result of the frequent transition from dioecy into systems with co-sexual individuals.

RESULTS

In this study, co-sexual expression (monoecy and hermaphroditic development), previously thought to be polyploid-specific in Diospyros species, was identified in the diploid D. oleifeara historically. We characterized potential genetic mechanisms that underlie the dissolution of dioecy to monoecy and andro(gyno)monoecy, based on multiscale genome-wide investigations of 150 accessions of Diospyros oleifera. We found all co-sexual plants, including monoecious and andro(gyno)monoecious individuals, possessed the male determinant gene OGI, implying the presence of genetic factors controlling gynoecia development in genetically male D. oleifera. Importantly, discrepancies in the OGI/MeGI module were found in diploid monoecious D. oleifera compared with polyploid monoecious D. kaki, including no Kali insertion on the promoter of OGI, no different abundance of smRNAs targeting MeGI (a counterpart of OGI), and no different expression of MeGI between female and male floral buds. On the contrary, in both single- and co-sexual plants, female function was expressed in the presence of a genome-wide decrease in methylation levels, along with sexually distinct regulatory networks of smRNAs and their targets. Furthermore, a genome-wide association study (GWAS) identified a genomic region and a DUF247 gene cluster strongly associated with the monoecious phenotype and several regions that may contribute to andromonoecy.

CONCLUSIONS

Collectively, our findings demonstrate stable breakdown of the dioecious system in D. oleifera, presumably also a result of genomic features of the Y-linked region.

Origin of purple asparagus cultivar 'Pacific Purple' based on the sequence of sex determination gene

Garden asparagus is one of the most important crops worldwide. Since this crop is dioecious and male plants generally have higher yields compared to female plants, several DNA markers for sex identification have been developed for acceleration of asparagus breeding. Among these markers, Asp1-T7sp and MSSTS710 were found to be effective in sex determination for many asparagus cultivars. However, we previously found that these markers were not completely suitable for sex identification in the purple asparagus cultivar 'Pacific Purple'. There are two types of male individuals in this cultivar: One type is PP-m, which is identified the sex type by Asp1-T7sp and MSSTS710 markers, while the other type is PP-m* whose sex type is not identified by these markers. Since the sex identification markers are located on the non-recombining Y region, it was expected that the sequence around this region might be different between PP-m and PP-m*. In this study, the sequence of one of the sex-determining genes, MSE1/AoMYB35/AspTDF1, was analyzed, and a comparative analysis was conducted among PP-m and PP-m* of 'Pacific Purple', A. officinalis and related species A. maritimus. The results revealed that PP-m and PP-m* has the similar sequence of MSE1/AoMYB35/AspTDF1 gene from A. officinalis and A. maritimus, respectively. 'Pacific Purple' is a cultivar developed through polycross hybrid from Italian landrace 'Violetto d'Albenga' (VA), suggesting that VA originated from an interspecific crossing between A. officinalis and A. maritimus and that the pollen parent used in 'Pacific Purple' breeding contained two types of male individuals with different MSE1/AoMYB35/AspTDF1 sequence. As a result, PP-m and PP-m* of 'Pacific Purple' harbors the similar sequences of the MSE1/AoMYB35/AspTDF1 gene from A. officinalis and A. maritimus, respectively.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

Integrated evolutionary pattern analyses reveal multiple origins of steroidal saponins in plants

Steroidal saponins are a class of specialized metabolites essential for plant's response to biotic and abiotic stresses. They are also important raw materials for the industrial production of steroid drugs. Steroidal saponins are present in some monocots, such as Dioscorea and Paris, but their distribution, origin, and evolution in plants remain poorly understood. By reconstructing the evolutionary history of the steroidal saponin-associated module (SSAM) in plants, we reveal that the steroidal saponin pathway has its origin in Asparagus and Dioscorea. Through evaluating the distribution and evolutionary pattern of steroidal saponins in angiosperms, we further show that steroidal saponins originated multiple times in angiosperms, and exist in early diverged lineages of certain monocot lineages including Asparagales, Dioscoreales, and Liliales. In these lineages, steroidal saponins are synthesized through the high copy and/or high expression mechanisms of key genes in SSAM. Together with shifts in gene evolutionary rates and amino acid usage, these molecular mechanisms shape the current distribution and diversity of steroidal saponins in plants. Consequently, our results provide new insights into the distribution, diversity and evolutionary history of steroidal saponins in plants, and enhance our understanding of plants' resistance to abiotic and biotic stresses. Additionally, fundamental understanding of the steroidal saponin biosynthesis will facilitate their industrial production and pharmacological applications.

Genome assembly of the bearded iris, Iris pallida Lam

Irises are perennial plants, representing a large genus with hundreds of species. While cultivated extensively for their ornamental value, commercial interest in irises lies in the secondary metabolites present in their rhizomes. The Dalmatian Iris (Iris pallida Lam.) is an ornamental plant that also produces secondary metabolites with potential value to the fragrance and pharmaceutical industries. In addition to providing base notes for the fragrance industry, iris tissues and extracts possess antioxidant, anti-inflammatory and immunomodulatory effects. However, study of these secondary metabolites has been hampered by a lack of genomic information, requiring difficult extraction and analysis techniques. Here, we report the genome sequence of Iris pallida Lam., generated with Pacific Bioscience long-read sequencing, resulting in a 10.04-Gbp assembly with a scaffold N50 of 14.34 Mbp and 91.8% complete BUSCOs. This reference genome will allow researchers to study the biosynthesis of these secondary metabolites in much greater detail, opening new avenues of investigation for drug discovery and fragrance formulations.

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.

High-quality Cymbidium mannii genome and multifaceted regulation of crassulacean acid metabolism in epiphytes

Epiphytes with crassulacean acid metabolism (CAM) photosynthesis are widespread among vascular plants, and repeated evolution of CAM photosynthesis is a key innovation for micro-ecosystem adaptation. However, we lack a complete understanding of the molecular regulation of CAM photosynthesis in epiphytes. Here, we report a high-quality chromosome-level genome assembly of a CAM epiphyte, Cymbidium mannii (Orchidaceae). The 2.88-Gb orchid genome with a contig N50 of 22.7 Mb and 27 192 annotated genes was organized into 20 pseudochromosomes, 82.8% of which consisted of repetitive elements. Recent expansions of long terminal repeat retrotransposon families have made a major contribution to the evolution of genome size in Cymbidium orchids. We reveal a holistic scenario of molecular regulation of metabolic physiology using high-resolution transcriptomics, proteomics, and metabolomics data collected across a CAM diel cycle. Patterns of rhythmically oscillating metabolites, especially CAM-related products, reveal circadian rhythmicity in metabolite accumulation in epiphytes. Genome-wide analysis of transcript and protein level regulation revealed phase shifts during the multifaceted regulation of circadian metabolism. Notably, we observed diurnal expression of several core CAM genes (especially βCA and PPC) that may be involved in temporal fixation of carbon sources. Our study provides a valuable resource for investigating post-transcription and translation scenarios in C. mannii, an Orchidaceae model for understanding the evolution of innovative traits in epiphytes.

Floral Aroma and Pollinator Relationships in Two Sympatric Late-Summer-Flowering Mediterranean Asparagus Species

This research delves into plant-pollinator relationships within the Mediterranean region, focusing on two synchronous and sympatric asparagus species: A. acutifolius and A. albus. For the first time, the floral scents of the genus Asparagus are reported. We investigate the volatile organic compounds (VOCs) present in their floral scents and their impact on pollinator attraction. Captured flower-emitted VOCs underwent solid-phase microextraction of headspace (SPME-HS) and gas chromatography and mass spectrometry (GC-MS) analysis. The investigation confirms distinctive aroma profiles for each species. A. albus predominantly emits benzene derivatives and sesquiterpenes, while A. acutifolius is characterized by carotenoid derivatives, monoterpenes, and sesquiterpenes. The only shared compounds between the two species are the sesquiterpenes (Z,E)-α-farnesene and (E,E)-α-farnesene. A positive correlation links peak floral aroma intensity (benzenoids in A. albus and ionones in A. acutifolius) with a higher pollinator visit frequency, emphasizing the critical role of intense floral scents in pollinator attraction. The study of reproductive aspects reveals almost complete gynodioecy in A. acutifolius, influencing unique dynamics for the two species. These adaptations hold significant importance within the Mediterranean ecosystem, particularly during the late dry summer period, when a limited number of plant species vie for a shared primary pollinator.

Genome-wide analysis of genetic diversity in a germplasm collection including wild relatives and interspecific clones of garden asparagus

The Asparagus genus includes approximately 240 species, the most important of which is garden asparagus (Asparagus officinalis L.), as this is a vegetable crop cultivated worldwide for its edible spear. Along with garden asparagus, other species are also cultivated (e.g., Asparagus maritimus L.) or have been proposed as untapped sources of variability in breeding programs (e.g., Asparagus acutifolius L.). In the present work, we applied reduced-representation sequencing to examine a panel of 378 diverse asparagus genotypes, including commercial hybrids, interspecific lines, wild relatives of garden asparagus, and doubled haploids currently used in breeding programs, which enabled the identification of more than 200K single-nucleotide polymorphisms (SNPs). These SNPs were used to assess the extent of linkage disequilibrium in the diploid gene pool of asparagus and combined with preliminary phenotypic information to conduct genome-wide association studies for sex and traits tied to spear quality and production. Moreover, using the same phenotypic and genotypic information, we fitted and cross-validated genome-enabled prediction models for the same set of traits. Overall, our analyses demonstrated that, unlike the diversity detected in wild species related to garden asparagus and in interspecific crosses, cultivated and wild genotypes of A. officinalis L. show a narrow genetic basis, which is a contributing factor hampering the genetic improvement of this crop. Estimating the extent of linkage disequilibrium and providing the first example of genome-wide association study and genome-enabled prediction in this species, we concluded that the asparagus panel examined in the present study can lay the foundation for determination of the genetic bases of agronomically important traits and for the implementation of predictive breeding tools to sustain breeding.

Typification and Nomenclature Notes on Twenty-Nine Names in Asparagus (Asparagaceae)

Nomenclatural types for twenty-nine names belonging to the genus Asparagus are typified and discussed. The following names are lectotypified: A. altiscandens Engl. & Gilg, A. altissimus Munby, A. baumii Engl. & Gilg, A. benguellensis Baker, A. burchellii Baker, A. curillus Buch.-Ham. ex Roxb., A. deflexus Baker, A. duchesnei L.Linden, A. equisetoides Welw. ex Baker, A. fasciculatus Thunb., A. griffithii, Baker, A. homblei De Wild., A. kaessneri De Wild., A. lecardii De Wild., A. longicladus N.E.Br., A. longiflorus Franch., A. monophyllus Baker, A. palaestinus Baker, A. pastorianus Webb & Berthel., A. persicus Baker, A. poissonii H.Perrier, A. psilurus Welw. ex Baker, A. ritschardii De Wild., A. sapinii De Wild., A. scandens Thunb., A. schumanianus Schltr. ex H.Perrier, A. stellatus Baker, A. subfalcatus De Wild., and A. undulatus (L.f.) Thunb. (synonym of Dracaena undulata L.f.). A new name, Asparagus neofasciculatus, is proposed as a replacement name for A. fasciculatus Thunb., which is an illegitimate later homonym of A. fasciculatus R.Br. The original protologue of these names and the original materials are evaluated. Nomenclature remarks discussing the selection of type specimens are given for each name, and known isotypes or isolectotypes are also cited. This information could be utilized as a reference for future taxonomic and systematic studies on Asparagus around the world.

The genome of Acorus deciphers insights into early monocot evolution

Acorales is the sister lineage to all the other extant monocot plants. Genomic resource enhancement of this genus can help to reveal early monocot genomic architecture and evolution. Here, we assemble the genome of Acorus gramineus and reveal that it has ~45% fewer genes than the majority of monocots, although they have similar genome size. Phylogenetic analyses based on both chloroplast and nuclear genes consistently support that A. gramineus is the sister to the remaining monocots. In addition, we assemble a 2.2 Mb mitochondrial genome and observe many genes exhibit higher mutation rates than that of most angiosperms, which could be the reason leading to the controversies of nuclear genes- and mitochondrial genes-based phylogenetic trees existing in the literature. Further, Acorales did not experience tau (τ) whole-genome duplication, unlike majority of monocot clades, and no large-scale gene expansion is observed. Moreover, we identify gene contractions and expansions likely linking to plant architecture, stress resistance, light harvesting, and essential oil metabolism. These findings shed light on the evolution of early monocots and genomic footprints of wetland plant adaptations.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

Identification of novel sex determination loci in Japanese weedy melon

KEY MESSAGE

Japanese weedy melon exhibits unique sex expression with interactions between previously reported sex determination genes and two novel loci. Sex expression contributes to fruit quality and yield in the Cucurbitaceae. In melon, orchestrated regulation by sex determination genes explains the mechanism of sex expression, resulting in a great variety of sexual morphologies. In this study, we examined the Japanese weedy melon UT1, which does not follow the reported model of sex expression. We conducted QTL analysis using F₂ plants for flower sex on the main stem and the lateral branch and mapped "occurrence of pistil-bearing flower on the main stem" locus on Chr. 3 (Opbf3.1) and "type of pistil-bearing flower" (female or bisexual) loci on Chr. 2 (tpbf2.1) and Chr. 8 (tpbf8.1). The Opbf3.1 included the known sex determination gene CmACS11. Sequence comparison of CmACS11 between parental lines revealed three nonsynonymous SNPs. A CAPS marker developed from one of the SNPs was closely linked to the occurrence of pistil-bearing flowers on the main stem in two F₂ populations with different genetic backgrounds. The UT1 allele on Opbf3.1 was dominant in F₁ lines from crosses between UT1 and diverse cultivars and breeding lines. This study suggests that Opbf3.1 and tpbf8.1 may promote the development of pistil and stamen primordia by inhibiting CmWIP1 and CmACS-7 functions, respectively, making the UT1 plants hermaphrodite. The results of this study provide new insights into the molecular mechanisms of sex determination in melons and considerations for the application of femaleness in melon breeding.

Transcriptome Sequencing of Agave amaniensis Reveals Shoot-Related Expression Patterns of Expansin A Genes in Agave

Agave species are widely planted for fiber production. However, the molecular basis of agave fiber development has not been well understood. In this study, we performed a transcriptomic analysis in A. amaniensi, a well-known variety with high-quality fiber production. Approximately 43.87 million clean reads were obtained using Illumina sequencing. The de novo assembly produced 66,746 unigrams, 54% of which were annotated in a public database. In the Nr database, 21,490 unigenes of A. amaniensis were shown to be most closely related to Asparagus officinalis. Nine expansin A orthologs with full coding regions were obtained, which were named EXP1a, EXP1b, EXP2, EXP3, EXP4a, EXP4b, EXP11, EXP12, and EXP13. The maximum likelihood phylogenetic tree revealed the species-specific expansion of expansin genes in Arabidopsis, rice and agave. The expression analysis suggested the negative correlation between the expression of expansin genes and the leaf growth rate, except AhEXP11. Moreover, expansin genes were differentially affected by abiotic and biotic stresses. Notably, AhEXP2 expression level was highly upgraded after the infection of Phytophthora nicotiana. Nutrient deficiency also influent expansin genes expression. Together, our research will benefit future studies related to fiber development, disease resistance and nutrient usage in agave.

Disruption of a DUF247 Containing Protein Alters Cell Wall Polysaccharides and Reduces Growth in Arabidopsis

Plant cell wall biosynthesis is a complex process that requires proteins and enzymes from glycan synthesis to wall assembly. We show that disruption of At3g50120 (DUF247-1), a member of the DUF247 multigene family containing 28 genes in Arabidopsis, results in alterations to the structure and composition of cell wall polysaccharides and reduced growth and plant size. An ELISA using cell wall antibodies shows that the mutants also exhibit ~50% reductions in xyloglucan (XyG), glucuronoxylan (GX) and heteromannan (HM) epitopes in the NaOH fraction and ~50% increases in homogalacturonan (HG) epitopes in the CDTA fraction. Furthermore, the polymer sizes of XyGs and GXs are reduced with concomitant increases in short-chain polymers, while those of HGs and mHGs are slightly increased. Complementation using 35S:DUF247-1 partially recovers the XyG and HG content, but not those of GX and HM, suggesting that DUF247-1 is more closely associated with XyGs and HGs. DUF247-1 is expressed throughout Arabidopsis, particularly in vascular and developing tissues, and its disruption affects the expression of other gene members, indicating a regulatory control role within the gene family. Our results demonstrate that DUF247-1 is required for normal cell wall composition and structure and Arabidopsis growth.

The proposed role of MSL-lncRNAs in causing sex lability of female poplars

Labile sex expression is frequently observed in dioecious plants, but the underlying genetic mechanism remains largely unknown. Sex plasticity is also observed in many Populus species. Here we carried out a systematic study on a maleness-promoting gene, MSL, detected in the Populus deltoides genome. Our results showed that both strands of MSL contained multiple cis-activating elements, which generated long non-coding RNAs (lncRNAs) promoting maleness. Although female P. deltoides did not have the male-specific MSL gene, a large number of partial sequences with high sequence similarity to this gene were detected in the female poplar genome. Based on sequence alignment, the MSL sequence could be divided into three partial sequences, and heterologous expression of these partial sequences in Arabidopsis confirmed that they could promote maleness. Since activation of the MSL sequences can only result in female sex lability, we propose that MSL-lncRNAs might play a role in causing sex lability of female poplars.

AoSAP8-P encoding A20 and/or AN1 type zinc finger protein in asparagus officinalis L. Improving stress tolerance in transgenic Nicotiana sylvestris

The full length CDS of an A20 and AN1 type zinc finger gene (named AoSAP8-P), located nearby the male specific Y chromosome (MSY) region of Asparagus officinalis (garden asparagus) was amplified by RT-PCR from purple passion. This gene, predicted as the stress associated protein (SAPs) gene families, encodes 172 amino acids with multiple cis elements including light, stress response box, MYB and ERF binding sites on its promoter. To analyze its function, the gene expression of different organs in different asparagus gender were analyzed and the overexpressed transgenic Nicotiana sylvestris lines were generated. The results showed the gene was highly expressed in both flower and root of male garden asparagus; the germination rate of seeds of the T2 transgenic lines (T2-5-4 and T2-7-1) under the stress conditions of 125 mM NaCl and 150 mM mannitol were significantly higher than the wild type (WT) respectively. When the potted T2-5-4, T2-7-1 lines and WT were subjected to drought stress for 24 days and the leaf discs immerged into 20 % PEG6000 and 300 mM NaCl solution for 48 h respectively, the T2-5-4 and T2-7-1 with AoSAP8-P expression showed stronger drought, salt and osmotic stress tolerance. When compared, the effects of AoSAP8-P overexpression on productive development showed that the flowering time of transgenic lines, were ∼ 9 day earlier with larger but fewer pollens than its WT counterparts. However, there were no significant differences in anthers, stigmas and pollen viability between the transgenic lines and WT. Our results suggested that, the AoSAP8-P gene plays a role in improving the stress resistance and shortening seeds generation time for perianal survival during the growth and development of garden asparagus.

Apostasia Mitochondrial Genome Analysis and Monocot Mitochondria Phylogenomics

Apostasia shenzhenica belongs to the subfamily Apostasioideae and is a primitive group located at the base of the Orchidaceae phylogenetic tree. However, the A. shenzhenica mitochondrial genome (mitogenome) is still unexplored, and the phylogenetic relationships between monocots mitogenomes remain unexplored. In this study, we discussed the genetic diversity of A. shenzhenica and the phylogenetic relationships within its monocotyledon mitogenome. We sequenced and assembled the complete mitogenome of A. shenzhenica, resulting in a circular mitochondrial draft of 672,872 bp, with an average read coverage of 122× and a GC content of 44.4%. A. shenzhenica mitogenome contained 36 protein-coding genes, 16 tRNAs, two rRNAs, and two copies of nad4L. Repeat sequence analysis revealed a large number of medium and small repeats, accounting for 1.28% of the mitogenome sequence. Selection pressure analysis indicated high mitogenome conservation in related species. RNA editing identified 416 sites in the protein-coding region. Furthermore, we found 44 chloroplast genomic DNA fragments that were transferred from the chloroplast to the mitogenome of A. shenzhenica, with five plastid-derived genes remaining intact in the mitogenome. Finally, the phylogenetic analysis of the mitogenomes from A. shenzhenica and 28 other monocots showed that the evolution and classification of most monocots were well determined. These findings enrich the genetic resources of orchids and provide valuable information on the taxonomic classification and molecular evolution of monocots.

Steroidal saponin profiles and their key genes for synthesis and regulation in Asparagus officinalis L. by joint analysis of metabolomics and transcriptomics

BACKGROUND

Asparagus officinalis L. is a worldwide cultivated vegetable enrichened in both nutrient and steroidal saponins with multiple pharmacological activities. The upstream biosynthetic pathway of steroidal saponins (USSP) for cholesterol (CHOL) synthesis has been studied, while the downstream pathway of steroidal saponins (DSSP) starting from cholesterol and its regulation in asparagus remains unknown.

RESULTS

Metabolomics, Illumina RNAseq, and PacBio IsoSeq strategies were applied to different organs of both cultivated green and purple asparagus to detect the steroidal metabolite profiles & contents and to screen their key genes for biosynthesis and regulation. The results showed that there is a total of 427 compounds, among which 18 steroids were detected with fluctuated concentrations in roots, spears and flowering twigs of two garden asparagus cultivars. The key genes of DSSP include; steroid-16-hydroxylase (S16H), steroid-22-hydroxylase (S22H) and steroid-22-oxidase-16-hydroxylase (S22O-16H), steroid-26-hydroxylase (S26H), steroid-3-β-glycosyltransferase (S3βGT) and furostanol glycoside 26-O-beta-glucosidases (F26GHs) which were correlated with the contents of major steroidal saponins were screened, and the transcriptional factors (TFs) co-expressing with the resulted from synthetic key genes, including zinc fingers (ZFs), MYBs and WRKYs family genes were also screened.

CONCLUSIONS

Based on the detected steroidal chemical structures, profiles and contents which correlated to the expressions of screened synthetic and TFs genes, the full steroidal saponin synthetic pathway (SSP) of asparagus, including its key regulation networks was proposed for the first time.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Evolution of sex in crops: recurrent scrap and rebuild

Sexuality is the main strategy for maintaining genetic diversity within a species. In flowering plants (angiosperms), sexuality is derived from ancestral hermaphroditism and multiple sexualities can be expressed in an individual. The mechanisms conferring chromosomal sex determination in plants (or dioecy) have been studied for over a century by both biologists and agricultural scientists, given the importance of this field for crop cultivation and breeding. Despite extensive research, the sex determining gene(s) in plants had not been identified until recently. In this review, we dissect plant sex evolution and determining systems, with a focus on crop species. We introduced classic studies with theoretical, genetic, and cytogenic approaches, as well as more recent research using advanced molecular and genomic techniques. Plants have undergone very frequent transitions into, and out of, dioecy. Although only a few sex determinants have been identified in plants, an integrative viewpoint on their evolutionary trends suggests that recurrent neofunctionalization events are potentially common, in a "scrap and (re)build" cycle. We also discuss the potential association between crop domestication and transitions in sexual systems. We focus on the contribution of duplication events, which are particularly frequent in plant taxa, as a trigger for the creation of new sexual systems.

The garden asparagus (Asparagus officinalis L.) mitochondrial genome revealed rich sequence variation throughout whole sequencing data

Garden asparagus (Asparagus officinalis L.) is a horticultural crop with high nutritional and medical value, considered an ideal plant for sex determination research among many dioecious plants, whose genomic information can support genetic analysis and breeding programs. In this research, the entire mitochondrial genome of A. officinalis was sequenced, annotated and assembled using a mixed Illumina and PacBio data. The garden asparagus circular mitochondrial genome measures 492,062 bp with a GC value of 45.9%. Thirty-six protein-coding genes, 17 tRNA and 6 rRNA genes were annotated, among which 8 protein-coding genes contained 16 introns. In addition, 254 SSRs with 10 complete tandem repeats and 293 non-tandem repeats were identified. It was found that the codons of edited sites located in the amino acids showed a leucine-formation trend, and RNA editing sites mainly caused the mutual transformation of amino acids with the same properties. Furthermore, 72 sequence fragments accounting for 20,240 bp, presentating 4.11% of the whole mitochondrial genome, were observed to migrate from chloroplast to mitochondrial genome of A. officinalis. The phylogenetic analysis showed that the closest genetic relationship between A. officinalis with onion (Allium cepa) inside the Liliaceae family. Our results demonstrated that high percentage of protein-coding genes had evolutionary conservative properties, with Ka/Ks values less than 1. Therefore, this study provides a high-quality garden asparagus mitochondrial genome, useful to promote better understanding of gene exchange between organelle genomes.

Male-linked gene TsRPL10a' in androdioecious tree Tapiscia sinensis: implications for sex differentiation by influencing gynoecium development

The mechanism of sex differentiation in androdioecy is of great significance for illuminating the origin and evolution of dioecy. Tapiscia sinensis Oliv. is a functionally androdioecious species with both male and hermaphroditic individuals. Male flowers of T. sinensis lack the ovules of gynoecia compared with hermaphrodites. To identify sex simply and accurately, and further find the potential determinants of sex differentiation in T. sinensis, we found that TsRPL10a', a duplicate of TsRPL10a, was a male-linked gene. The promoter (5' untranslated region and the first intron) of TsRPL10a' can be used to accurately identify sex in T. sinensis. TsRPL10a is a ribosomal protein that is involved in gynoecium development, and sufficient ribosomal levels are necessary for female gametogenesis. The expression level of TsRPL10a was significantly downregulated in male flower primordia compared with hermaphrodites. The RNA fluorescence in situ hybridization (FISH) assay demonstrated that TsRPL10a was almost undetectable in male gynoecia at the gynoecial ridge stage, which was a key period of ovule formation by scanning electron microscope observation. In male flowers, although the promoter activity of TsRPL10a was significantly higher than TsRPL10a' verified by transgenic Arabidopsis thaliana, the transcriptional expression ratio of TsRPL10a was obviously lower than TsRPL10a' and reached its lowest at the gynoecial ridge stage, indicating the existence of a female suppressor. The promoter similarity of TsRPL10a and TsRPL10a' was only 45.29%; the genomic sequence similarity was 89.8%; four amino acids were altered in TsRPL10a'. The secondary structure of TsRPL10a' was different from TsRPL10a, and TsRPL10a' did not exhibit FISH and GUS expression in the gynoecium the way TsRPL10a did. From the perspective of RT-qPCR, its high expression level, followed by the low expression level of TsRPL10a in male flowers, indicates its antagonism function with TsRPL10a. The evolutionary analysis, subcellular localization and flower expression pattern suggested that TsRPL10a might be functionally conserved with AtRPL10aA, AtRPL10aB and AtRPL10aC in A. thaliana. Overall, we speculated that TsRPL10a and its duplicate TsRPL10a' might be involved in sex differentiation by influencing gynoecium development in T. sinensis.

Identifying Genes Associated with Female Flower Development of Phellodendron amurense Rupr. Using a Transcriptomics Approach

Phellodendron amurense Rupr., a species of Rutaceae, is a nationally protected and valuable medicinal plant. It is generally considered to be dioecious. With the discovery of monoecious P. amurense, the phenomenon that its sex development is regulated by epigenetics has been revealed, but the way epigenetics affects the sex differentiation of P. amurense is still unclear. In this study, we investigated the effect of DNA methylation on the sexual development of P. amurense. The young inflorescences of male plants were treated with the demethylation agent 5-azaC, and the induced female flowers were obtained. The induced female flowers’ morphological functions and transcriptome levels were close to those of normally developed plants. Genes associated with the development of female flowers were studied by comparing the differences in transcriptome levels between the male and female flowers. Referring to sex-related genes reported in other plants, 188 candidate genes related to the development of female flowers were obtained, including sex-regulating genes, genes related to the formation and development of sexual organs, genes related to biochemical pathways, and hormone-related genes. RPP0W, PAL3, MCM2, MCM6, SUP, PIN1, AINTEGUMENTA, AINTEGUMENTA-LIKE6, AGL11, SEUSS, SHI-RELATED SEQUENCE 5, and ESR2 were preliminarily considered the key genes for female flower development. This study has demonstrated that epigenetics was involved in the sex regulation of P. amurense, with DNA methylation as one of its regulatory modes. Moreover, some candidate genes related to the sexual differentiation of P. amurense were obtained with analysis. These results are of great significance for further exploring the mechanism of sex differentiation of P. amurense and studying of sex differentiation of plants.

A divergent RWP-RK transcription factor determines mating type in heterothallic Closterium

The Closterium peracerosum-strigosum-littorale complex (Closterium, Zygnematophyceae) has an isogamous mating system. Members of the Zygnematophyceae are the closest relatives to extant land plants and are distantly related to chlorophytic models, for which a genetic basis of mating type (MT) determination has been reported. We thus investigated MT determination in Closterium. We sequenced genomes representing the two MTs, mt+ and mt-, in Closterium and identified CpMinus1, a gene linked to the mt- phenotype. We analyzed its function using reverse genetics methods. CpMinus1 encodes a divergent RWP-RK domain-containing-like transcription factor and is specifically expressed during gamete differentiation. Introduction of CpMinus1 into an mt+ strain was sufficient to convert it to a phenotypically mt- strain, while CpMinus1-knockout mt- strains were phenotypically mt+. We propose that CpMinus1 is the major MT determinant that acts by evoking the mt- phenotype and suppressing the mt+ phenotype in heterothallic Closterium. CpMinus1 likely evolved independently in the Zygnematophyceae lineage, which lost an egg-sperm anisogamous mating system. mt- specific regions possibly constitute an MT locus flanked by common sequences that undergo some recombination.

Regulatory mechanism of MeGI on sexuality in Diospyros oleifera

Dioecy system is an important strategy for maintaining genetic diversity. The transcription factor MeGI, contributes to dioecy by promoting gynoecium development in Diospyros lotus and D. kaki. However, the function of MeGI in D. oleifera has not been identified. In this study, we confirmed that MeGI, cloned from D. oleifera, repressed the androecium development in Arabidopsis thaliana. Subsequently, chromatin immunoprecipitation-sequencing (ChIP-seq), DNA affinity purification-sequencing (DAP-seq), and RNA-seq were used to uncover the gene expression response to MeGI. The results showed that the genes upregulated and downregulated in response to MeGI were mainly enriched in the circadian rhythm-related and flavonoid biosynthetic pathways, respectively. Additionally, the WRKY DNA-binding protein 28 (WRKY28) gene, which was detected by ChIP-seq, DAP-seq, and RNA-seq, was emphasized. WRKY28 has been reported to inhibit salicylic acid (SA) biosynthesis and was upregulated in MeGI-overexpressing A. thaliana flowers, suggesting that MeGI represses the SA level by increasing the expression level of WRKY28. This was confirmed that SA level was lower in D. oleifera female floral buds than male. Overall, our findings indicate that the MeGI mediates its sex control function in D. oleifera mainly by regulating genes in the circadian rhythm, SA biosynthetic, and flavonoid biosynthetic pathways.

Genomic profiling of dioecious Amaranthus species provides novel insights into species relatedness and sex genes

BACKGROUND

Amaranthus L. is a diverse genus consisting of domesticated, weedy, and non-invasive species distributed around the world. Nine species are dioecious, of which Amaranthus palmeri S. Watson and Amaranthus tuberculatus (Moq.) J.D. Sauer are troublesome weeds of agronomic crops in the USA and elsewhere. Shallow relationships among the dioecious Amaranthus species and the conservation of candidate genes within previously identified A. palmeri and A. tuberculatus male-specific regions of the Y (MSYs) in other dioecious species are poorly understood. In this study, seven genomes of dioecious amaranths were obtained by paired-end short-read sequencing and combined with short reads of seventeen species in the family Amaranthaceae from NCBI database. The species were phylogenomically analyzed to understand their relatedness. Genome characteristics for the dioecious species were evaluated and coverage analysis was used to investigate the conservation of sequences within the MSY regions.

RESULTS

We provide genome size, heterozygosity, and ploidy level inference for seven newly sequenced dioecious Amaranthus species and two additional dioecious species from the NCBI database. We report a pattern of transposable element proliferation in the species, in which seven species had more Ty3 elements than copia elements while A. palmeri and A. watsonii had more copia elements than Ty3 elements, similar to the TE pattern in some monoecious amaranths. Using a Mash-based phylogenomic analysis, we accurately recovered taxonomic relationships among the dioecious Amaranthus species that were previously identified based on comparative morphology. Coverage analysis revealed eleven candidate gene models within the A. palmeri MSY region with male-enriched coverages, as well as regions on scaffold 19 with female-enriched coverage, based on A. watsonii read alignments. A previously reported FLOWERING LOCUS T (FT) within A. tuberculatus MSY contig was also found to exhibit male-enriched coverages for three species closely related to A. tuberculatus but not for A. watsonii reads. Additional characterization of the A. palmeri MSY region revealed that 78% of the region is made of repetitive elements, typical of a sex determination region with reduced recombination.

CONCLUSIONS

The results of this study further increase our understanding of the relationships among the dioecious species of the Amaranthus genus as well as revealed genes with potential roles in sex function in the species.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Dirigent protein subfamily function and structure in terrestrial plant phenol metabolism

Aquatic plant transition to land, and subsequent terrestrial plant species diversification, was accompanied by the emergence and massive elaboration of plant phenol chemo-diversity. Concomitantly, dirigent protein (DP) and dirigent-like protein subfamilies, derived from large multigene families, emerged and became extensively diversified. DP biochemical functions as gateway entry points into new and diverse plant phenol skeletal types then markedly expanded. DPs have at least eight non-uniformly distributed subfamilies, with different DP subfamily members of known biochemical/physiological function now implicated as gateway entries to lignan, lignin, aromatic diterpenoid, pterocarpan and isoflavene pathways. While some other DP subfamily members have jacalin domains, both these and indeed the majority of DPs throughout the plant kingdom await discovery of their biochemical roles. Methods and approaches were developed to discover DP biochemical function as gateway entry points to distinct plant phenol skeletal types in land plants. Various DP 3D X-ray structural determinations enabled structure-based comparative sequence analysis and modeling to understand similarities and differences among the different DP subfamilies. We consider that the core DP β-barrel fold and associated characteristics are likely common to all DPs, with several residues conserved and nearly invariant. There is also considerable variation in residue composition and topography of the putative substrate binding pockets, as well as substantial differences in several loops, such as the β1-β2 loop. All DPs likely bind and stabilize quinone methide intermediates, while guiding distinctive regio- and/or stereo-chemical entry into Nature's chemo-diverse land plant phenol metabolic classes.

A chromosome-level genome assembly for Dracaena cochinchinensis reveals the molecular basis of its longevity and formation of dragon's blood

Dracaena, a remarkably long-lived and slowly maturing species of plant, is world famous for its ability to produce dragon's blood, a precious traditional medicine used by different cultures since ancient times. However, there is no detailed and high-quality genome available for this species at present; thus, the molecular mechanisms that underlie its important traits are largely unknown. These factors seriously limit the protection and regeneration of this rare and endangered plant resource. Here, we sequenced and assembled the genome of Dracaena cochinchinensis at the chromosome level. The D. cochinchinensis genome covers 1.21 Gb with a scaffold N50 of 50.06 Mb and encodes 31 619 predicted protein-coding genes. Analysis showed that D. cochinchinensis has undergone two whole-genome duplications and two bursts of long terminal repeat insertions. The expansion of two gene classes, cis-zeatin O-glucosyltransferase and small auxin upregulated RNA, were found to account for its longevity and slow growth. Two transcription factors (bHLH and MYB) were found to be core regulators of the flavonoid biosynthesis pathway, and reactive oxygen species were identified as the specific signaling molecules responsible for the injury-induced formation of dragon's blood. Our study provides high-quality genomic information relating to D. cochinchinensis and significant insight into the molecular mechanisms responsible for its longevity and formation of dragon's blood. These findings will facilitate resource protection and sustainable utilization of Dracaena.

Chromosome-level genome assembly of bunching onion illuminates genome evolution and flavor formation in Allium crops

The Allium genus is cultivated globally as vegetables, condiments, or medicinal plants and is characterized by large genomes and strong pungency. However, the genome evolution and genomic basis underlying their unique flavor formation remain poorly understood. Herein, we report an 11.27-Gb chromosome-scale genome assembly for bunching onion (A. fistulosum). The uneven bursts of long-terminal repeats contribute to diversity in genome constituents, and dispersed duplication events largely account for gene expansion in Allium genomes. The extensive duplication and differentiation of alliinase and lachrymatory factor synthase manifest as important evolutionary events during flavor formation in Allium crops. Furthermore, differential selective preference for flavor-related genes likely lead to the variations in isoalliin content in bunching onions. Moreover, we reveal that China is the origin and domestication center for bunching onions. Our findings provide insights into Allium genome evolution, flavor formation and domestication history and enable future genome-assisted breeding of important traits in these crops.

A tandem CCCH type zinc finger protein gene CpC3H3 from Chimonanthus praecox promotes flowering and enhances drought tolerance in Arabidopsis

BACKGROUND

CCCH-type zinc finger proteins play important roles in plant development and biotic/abiotic stress responses. Wintersweet (Chimonanthus praecox) is a popular ornamental plant with strong resistance to various stresses, which is a good material for exploring gene resource for stress response. In this study, we isolated a CCCH type zinc finger protein gene CpC3H3 (MZ964860) from flower of wintersweet and performed functional analysis with a purpose of identifying gene resource for floral transition and stress tolerance.

RESULTS

CpC3H3 was predicted a CCCH type zinc finger protein gene encoding a protein containing 446 amino acids with five conserved C-X₈-C-X₅-C-X₃-H motifs. CpC3H3 was localized in the cell membrane but with a nuclear export signal at the N-terminal. Transcripts of CpC3H3 were significantly accumulated in flower buds at floral meristem formation stage, and were induced by polyethylene glycol. Overexpression of CpC3H3 promoted flowering, and enhanced drought tolerance in transgenic A. thaliana. CpC3H3 overexpression affects the expression level of genes involved in flower inducement and stress responses. Further comparative studies on physiological indices showed the contents of proline and soluble sugar, activity of peroxidase and the rates of electrolyte leakage were significantly increased and the content of malondialdehyde and osmotic potential was significantly reduced in transgenic A. thaliana under PEG stress.

CONCLUSION

Overall, CpC3H3 plays a role in flowering inducement and drought tolerance in transgenic A. thaliana. The CpC3H3 gene has the potential to be used to promote flowering and enhance drought tolerance in plants.

A CLAVATA3-like Gene Acts as a Gynoecium Suppression Function in White Campion

How do separate sexes originate and evolve? Plants provide many opportunities to address this question as they have diverse mating systems and separate sexes (dioecy) that evolved many times independently. The classic "two-factor" model for evolution of separate sexes proposes that males and females can evolve from hermaphrodites via the spread of male and female sterility mutations that turn hermaphrodites into females and males, respectively. This widely accepted model was inspired by early genetic work in dioecious white campion (Silene latifolia) that revealed the presence of two sex-determining factors on the Y-chromosome, though the actual genes remained unknown. Here, we report identification and functional analysis of the putative sex-determining gene in S. latifolia, corresponding to the gynoecium suppression factor (GSF). We demonstrate that GSF likely corresponds to a Y-linked CLV3-like gene that is specifically expressed in early male flower buds and encodes the protein that suppresses gynoecium development in S. latifolia. Interestingly, GSFY has a dysfunctional X-linked homolog (GSFX) and their synonymous divergence (dS = 17.9%) is consistent with the age of sex chromosomes in this species. We propose that female development in S. latifolia is controlled via the WUSCHEL-CLAVATA feedback loop, with the X-linked WUSCHEL-like and Y-linked CLV3-like genes, respectively. Evolution of dioecy in the S. latifolia ancestor likely involved inclusion of ancestral GSFY into the nonrecombining region on the nascent Y-chromosome and GSFX loss of function, which resulted in disbalance of the WUSCHEL-CLAVATA feedback loop between the sexes and ensured gynoecium suppression in males.

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.

Repeated turnovers keep sex chromosomes young in willows

Background

Salicaceae species have diverse sex determination systems and frequent sex chromosome turnovers. However, compared with poplars, the diversity of sex determination in willows is poorly understood, and little is known about the evolutionary forces driving their turnover. Here, we characterized the sex determination in two Salix species, S. chaenomeloides and S. arbutifolia, which have an XY system on chromosome 7 and 15, respectively.

Results

Based on the assemblies of their sex determination regions, we found that the sex determination mechanism of willows may have underlying similarities with poplars, both involving intact and/or partial homologs of a type A cytokinin response regulator (RR) gene. Comparative analyses suggested that at least two sex turnover events have occurred in Salix, one preserving the ancestral pattern of male heterogamety, and the other changing heterogametic sex from XY to ZW, which could be partly explained by the “deleterious mutation load” and “sexually antagonistic selection” theoretical models. We hypothesize that these repeated turnovers keep sex chromosomes of willow species in a perpetually young state, leading to limited degeneration.

Conclusions

Our findings further improve the evolutionary trajectory of sex chromosomes in Salicaceae species, explore the evolutionary forces driving the repeated turnovers of their sex chromosomes, and provide a valuable reference for the study of sex chromosomes in other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02769-w.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

Differential timing of gene expression and recruitment in independent origins of CAM in the Agavoideae (Asparagaceae)

Crassulacean acid metabolism (CAM) photosynthesis has evolved repeatedly across the plant tree of life, however our understanding of the genetic convergence across independent origins remains hampered by the lack of comparative studies. Here, we explore gene expression profiles in eight species from the Agavoideae (Asparagaceae) encompassing three independent origins of CAM. Using comparative physiology and transcriptomics, we examined the variable modes of CAM in this subfamily and the changes in gene expression across time of day and between well watered and drought-stressed treatments. We further assessed gene expression and the molecular evolution of genes encoding phosphoenolpyruvate carboxylase (PPC), an enzyme required for primary carbon fixation in CAM. Most time-of-day expression profiles are largely conserved across all eight species and suggest that large perturbations to the central clock are not required for CAM evolution. By contrast, transcriptional response to drought is highly lineage specific. Yucca and Beschorneria have CAM-like expression of PPC2, a copy of PPC that has never been shown to be recruited for CAM in angiosperms. Together the physiological and transcriptomic comparison of closely related C₃ and CAM species reveals similar gene expression profiles, with the notable exception of differential recruitment of carboxylase enzymes for CAM function.

A convergent mechanism of sex determination in dioecious plants: Distinct sex-determining genes display converged regulation on floral B-class genes

Sex determination in dioecious plants has been broadly and progressively studied with the blooming of genome sequencing and editing techniques. This provides us with a great opportunity to explore the evolution and genetic mechanisms underlining the sex-determining system in dioecious plants. In this study, comprehensively reviewing advances in sex-chromosomes, sex-determining genes, and floral MADS-box genes in dioecious plants, we proposed a convergent model that governs plant dioecy across divergent species using a cascade regulation pathway connecting sex-determining genes and MADS-box genes e.g., B-class genes. We believe that this convergent mechanism of sex determination in dioecious plants will shed light on our understanding of gene regulation and evolution of plant dioecy. Perspectives concerning the evolutionary pathway of plant dioecy are also suggested.

Chromosome-level Genomes Reveal the Genetic Basis of Descending Dysploidy and Sex Determination in Morus Plants

Multiple plant lineages have independently evolved sex chromosomes and variable karyotypes to maintain their sessile lifestyles through constant biological innovation. Morus notabilis, a dioecious mulberry species, has the fewest chromosomes among Morus spp., but the genetic basis of sex determination and karyotype evolution in this species has not been identified. In this study, three high-quality genome assemblies were generated for Morus spp. [including dioecious M. notabilis (male and female) and Morus yunnanensis (female)] with genome sizes of 301-329 Mb and were grouped into six pseudochromosomes. Using a combination of genomic approaches, we found that the putative ancestral karyotype of Morus species was close to 14 protochromosomes, and that several chromosome fusion events resulted in descending dysploidy (2n = 2x = 12). We also characterized a ∼ 6.2-Mb sex-determining region on chromosome 3. Four potential male-specific genes, a partially duplicatedDNA helicase gene (named MSDH) and three Ty3_Gypsy long terminal repeat retrotransposons (named MSTG1/2/3), were identified in the Y-linked area and considered to be strong candidate genes for sex determination or differentiation. Population genomic analysis showed that Guangdong accessions in China were genetically similar to Japanese accessions of mulberry. In addition, genomic areas containing selective sweeps that distinguish domesticated mulberry from wild populations in terms of flowering and disease resistance were identified. Our findings provide an important genetic resource for sex identification research and molecular breeding in mulberry.

Remarkable Divergence of the Sex-Linked Region between Two Wild Spinach Progenitors, Spinacia turkestanica and Spinacia tetrandra

The sex-linked region (SLR) plays an important role in determining the sex of a plant. The SLR of the Y chromosome, composed of a 14.1-Mb inversion and a 10-Mb Y-duplication region (YDR), was deciphered in Spinacia oleracea previously. However, our understanding of the SLR in its wild relatives, S. turkestanica and S. tetrandra, remains limited. In this study, we used 63 resequencing data from the three Spinacia species to infer the evolution of the SLR among the Spinacia species. In the SLR, all the cultivated spinach and S. turkestanica accessions were clustered into two distinct categories with both sexes, while the S. tetrandra accessions of both sexes were grouped. This suggests that S. oleracea shared a similar SLR with S. turkestanica, but not with S. tetrandra, which was further confirmed based on the population structure and principal component analysis. Furthermore, we identified 3910 fully sex-linked SNPs in S. oleracea and 92.82% of them were available in S. turkestanica, while none of the SNPs were adopted in S. tetrandra. Genome coverage in males and females supported the hypothesis that the YDR increasingly expanded during its evolution. Otherwise, we identified 13 sex-linked transposable element insertion polymorphisms within the inversion in both S. oleracea and S. turkestanica, demonstrating that the transposable element insertions might have occurred before the recombination suppression event of the inversion. The SLR was conserved compared with the pseudoautosomal region given that the genetic hitchhiking process occurred in the SLR during its evolution. Our findings will significantly advance our understanding of the characteristics and evolution of the SLR in Spinacia species.

Recent expansion of the non-recombining sex-linked region on Silene latifolia sex chromosomes

Evolution of a non-recombining sex-specific region on the Y (or W) chromosome (NRY) is a key step in sex chromosome evolution, but how recombination suppression evolves is not well understood. Studies in many different organisms indicated that NRY evolution often involves several expansion steps. Why such NRY expansions occur remains unclear, although it is though that they are likely driven by sexually antagonistic selection. This paper describes a recent NRY expansion due to shift of the pseudoautosomal boundary on the sex chromosomes of a dioecious plant Silene latifolia. The shift resulted in inclusion of at least 16 pseudoautosomal genes into the NRY. This region is pseudoautosomal in closely related Silene dioica and Silene diclinis, indicating that the NRY expansion occurred in S. latifolia after it speciated from the other species ~120 thousand years ago. As S. latifolia and S. dioica actively hybridise across Europe, interspecific gene flow could blur the PAR boundary in these species. The pseudoautosomal genes have significantly elevated genetic diversity (π ~ 3% at synonymous sites), which is consistent with balancing selection maintaining diversity in this region. The recent shift of the PAR boundary in S. latifolia offers an opportunity to study the process of on-going NRY expansion.