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

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.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

Genome-Wide Identification and Analysis of the NAC Transcription Factor Gene Family in Garden Asparagus (Asparagus officinalis)

As a large plant-specific gene family, the NAC (NAM, ATAF1/2, and CUC2) transcription factor is related to plant growth, development, and response to abiotic stresses. Although the draft genome of garden asparagus (Asparagus officinalis) has been released, the genome-wide investigation of the NAC gene family is still unavailable. In this study, a total of 85 A. officinalis NAC genes were identified, and a comprehensive analysis of the gene family was performed, including physicochemical properties, phylogenetic relationship, chromosome localization, gene structure, conserved motifs, intron/exon, cis-acting elements, gene duplication, syntenic analysis, and differential gene expression analysis. The phylogenetic analysis demonstrated that there were 14 subgroups in both A. officinalis and Arabidopsis thaliana, and the genes with a similar gene structure and motif distribution were clustered in the same group. The cis-acting regulatory analysis of AoNAC genes indicated four types of cis-acting elements were present in the promoter regions, including light-responsive, hormone-responsive, plant-growth-and-development-related, and stress-responsive elements. The chromosomal localization analysis found that 81 NAC genes in A. officinalis were unevenly distributed on nine chromosomes, and the gene duplication analysis showed three pairs of tandem duplicated genes and five pairs of segmental duplications, suggesting that gene duplication is possibly associated with the amplification of the A. officinalis NAC gene family. The differential gene expression analysis revealed one and three AoNAC genes that were upregulated and downregulated under different types of salinity stress, respectively. This study provides insight into the evolution, diversity, and characterization of NAC genes in garden asparagus and will be helpful for future understanding of their biological roles and molecular mechanisms in plants.

Phylogeny and Expression Atlas of the NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY in Agave

Agave species are widely cultivated crassulacean acid metabolism (CAM) plants for alcoholic beverages, food and fiber production. Among these, the Agave hybrid H11648 ((A. amaniensis × A. angustifolia) × A. amaniensis) is the main cultivar for sisal fiber in the tropical areas of Brazil, China, and African countries. The plants of Agave hybrid H11648 have a long life cycle and large leaves, which require a huge amount of nitrogen nutrient. However, the molecular basis of nitrogen transport and allocation has not been well understood in agave. In this study, we identified 19 NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY(NPF) genes (called AhNPFs) with full-length coding sequences in Agave hybrid H11648. Our analysis of gene expression in various types of tissues revealed the tissue-specific expression pattern of AhNPFs. We further examined their expression patterns at different leaf developmental stages, under abiotic/biotic stresses and nutrient deficiency. The results reveal several candidate regulators in the agave NPF family, including AhNPF4.3/5.2/7.1. We first characterized the NPF genes in agave based on published leaf transcriptome datasets and emphasized their potential functions. The study will benefit future studies related to nitrogen nutrient in agave.

Aspects of In Vitro Plant Tissue Culture and Breeding of Asparagus: A Review

In vitro plant tissue culture and biotechnology used to assist and support the development of plant breeding when classical methods of propagation must be accelerated or it was necessary to overcome barriers inaccessible by classical approaches. In asparagus, to improve multiple breeding tasks, a high number of in vitro methods have been used, such as plant regeneration methods through organogenesis, embryogenesis, manipulation of ploidy, protoplast isolation, genetic manipulation (protoplast fusion, genetic transformation), embryo rescue and germplasm preservation (in vitro, in vitro slow growth, cryopreservation). Plant tissue culture methods can overcome multiple problems in asparagus breeding such as, barriers of self and cross-incompatibility between asparagus species through embryo rescue of interspecific hybrids and protoplast fusion or genetic transformation, introgression of new genes, clonal propagation of elite genotypes of asparagus, mass screening, and the generation of haploid and polyploid genotypes, among others, becoming the tool of choice for asparagus breeding programs. Some of these in vitro methods are still under development.

Representing sex chromosomes in genome assemblies

Sex chromosomes have evolved hundreds of independent times across eukaryotes. As genome sequencing, assembly, and scaffolding techniques rapidly improve, it is now feasible to build fully phased sex chromosome assemblies. Despite technological advances enabling phased assembly of whole chromosomes, there are currently no standards for representing sex chromosomes when publicly releasing a genome. Furthermore, most computational analysis tools are unable to efficiently investigate their unique biology relative to autosomes. We discuss a diversity of sex chromosome systems and consider the challenges of representing sex chromosome pairs in genome assemblies. By addressing these issues now as technologies for full phasing of chromosomal assemblies are maturing, we can collectively ensure that future genome analysis toolkits can be broadly applied to all eukaryotes with diverse types of sex chromosome systems. Here we provide best practice guidelines for presenting a genome assembly that contains sex chromosomes. These guidelines can also be applied to other non-recombining genomic regions, such as S-loci in plants and mating-type loci in fungi and algae.

Are plant and animal sex chromosomes really all that different?

Sex chromosomes in plants have often been contrasted with those in animals with the goal of identifying key differences that can be used to elucidate fundamental evolutionary properties. For example, the often homomorphic sex chromosomes in plants have been compared to the highly divergent systems in some animal model systems, such as birds, Drosophila and therian mammals, with many hypotheses offered to explain the apparent dissimilarities, including the younger age of plant sex chromosomes, the lesser prevalence of sexual dimorphism, or the greater extent of haploid selection. Furthermore, many plant sex chromosomes lack complete sex chromosome dosage compensation observed in some animals, including therian mammals, Drosophila, some poeciliids, and Anolis, and plant dosage compensation, where it exists, appears to be incomplete. Even the canonical theoretical models of sex chromosome formation differ somewhat between plants and animals. However, the highly divergent sex chromosomes observed in some animal groups are actually the exception, not the norm, and many animal clades are far more similar to plants in their sex chromosome patterns. This begs the question of how different are plant and animal sex chromosomes, and which of the many unique properties of plants would be expected to affect sex chromosome evolution differently than animals? In fact, plant and animal sex chromosomes exhibit more similarities than differences, and it is not at all clear that they differ in terms of sexual conflict, dosage compensation, or even degree of divergence. Overall, the largest difference between these two groups is the greater potential for haploid selection in plants compared to animals. This may act to accelerate the expansion of the non-recombining region at the same time that it maintains gene function within it. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Recurrent allopolyploidization, Y-chromosome introgression and the evolution of sexual systems in the plant genus Mercurialis

The plant genus Mercurialis includes dioecious, monoecious and androdioecious species (where males coexist with hermaphrodites). Its diversification involved reticulate evolution via hybridization and polyploidization. The Y chromosome of the diploid species Mercurialis annua shows only mild signs of degeneration. We used sequence variation at a Y-linked locus in several species and at multiple autosomal and pseudoautosomal loci to investigate the origin and evolution of the Y chromosome across the genus. Our study provides evidence for further cases of allopolyploid speciation. It also reveals that all lineages with separate sexes (with one possible exception) share the same ancestral Y chromosome. Surprisingly, males in androdioecious populations of hexaploid M. annua carry a Y chromosome that is not derived from either of its two putative progenitor lineages but from a more distantly related perennial dioecious lineage via introgression. These results throw new light on the evolution of sexual systems and polyploidy in Mercurialis and secure it as a promising model for further study of plant sex chromosomes. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Evolution of the sex-determining region in Ginkgo biloba

Sex chromosomes or sex-determining regions (SDR) have been discovered in many dioecious plant species, including the iconic 'living fossil' Ginkgo biloba, though the location and size of the SDR in G. biloba remain contradictory. Here we resolve these controversies and analyse the evolution of the SDR in this species. Based on transcriptome sequencing data from four genetic crosses we reconstruct male- and female-specific genetic maps and locate the SDR to the middle of chromosome 2. Integration of the genetic maps with the genome sequence reveals that recombination in and around the SDR is suppressed in a region of about 50 Mb in both males and females. However, occasional recombination does occur except a small, less than 5 Mb long region that does not recombine in males. Based on synonymous divergence between homologous X- and Y-linked genes in this region, we infer that the Ginkgo SDR is fairly old-at least of Cretaceous origin. The analysis of substitution rates and gene expression reveals only slight Y-degeneration. These results are consistent with findings in other dioecious plants with homomorphic sex chromosomes, where the SDR is typically small and evolves in a region with pre-existing reduced recombination, surrounded by long actively recombining pseudoautosomal regions. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

Chemical genetics in Silene latifolia elucidate regulatory pathways involved in gynoecium development

Dioecious plants possess diverse sex determination systems and unique mechanisms of reproductive organ development; however, little is known about how sex-linked genes shape the expression of regulatory cascades that lead to developmental differences between sexes. In Silene latifolia, a dioecious plant with stable dimorphism in floral traits, early experiments suggested that female-regulator genes act on the factors that determine the boundaries of the flower whorls. To identify these regulators, we sequenced the transcriptome of male flowers with fully developed gynoecia, induced by rapid demethylation in the parental generation. Eight candidates were found to have a positive role in gynoecium promotion, floral organ size, and whorl boundary, and affect the expression of class B MADS-box flower genes. To complement our transcriptome analysis, we closely examined the floral organs in their native state using field emission environmental scanning electron microscopy, and examined the differences between females and androhermaphrodites in their placenta and ovule organization. Our results reveal the regulatory pathways potentially involved in sex-specific flower development in the classical model of dioecy, S. latifolia. These pathways include previously hypothesized and unknown female-regulator genes that act on the factors that determine the flower boundaries, and a negative regulator of anther development, SUPERMAN-like (SlSUP).

Chromosome evolution and the genetic basis of agronomically important traits in greater yam

The nutrient-rich tubers of the greater yam, Dioscorea alata L., provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an orphan crop. Here, we address this resource gap by presenting a highly contiguous chromosome-scale genome assembly of D. alata combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient allotetraploidization in the Dioscorea lineage, followed by extensive genome-wide reorganization. Using the genomic tools, we find quantitative trait loci for resistance to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

The Gastrodia menghaiensis (Orchidaceae) genome provides new insights of orchid mycorrhizal interactions

BACKGROUND

To illustrate the molecular mechanism of mycoheterotrophic interactions between orchids and fungi, we assembled chromosome-level reference genome of Gastrodia menghaiensis (Orchidaceae) and analyzed the genomes of two species of Gastrodia.

RESULTS

Our analyses indicated that the genomes of Gastrodia are globally diminished in comparison to autotrophic orchids, even compared to Cuscuta (a plant parasite). Genes involved in arbuscular mycorrhizae colonization were found in genomes of Gastrodia, and many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids. The highly expressed genes for fatty acid and ammonium root transporters suggest that fungi receive material from orchids, although most raw materials flow from the fungi. Many nuclear genes (e.g. biosynthesis of aromatic amino acid L-tryptophan) supporting plastid functions are expanded compared to photosynthetic orchids, an indication of the importance of plastids even in totally mycoheterotrophic species.

CONCLUSION

Gastrodia menghaiensis has the smallest proteome thus far among angiosperms. Many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids.

Genomes of leafy and leafless Platanthera orchids illuminate the evolution of mycoheterotrophy

To improve our understanding of the origin and evolution of mycoheterotrophic plants, we here present the chromosome-scale genome assemblies of two sibling orchid species: partially mycoheterotrophic Platanthera zijinensis and holomycoheterotrophic Platanthera guangdongensis. Comparative analysis shows that mycoheterotrophy is associated with increased substitution rates and gene loss, and the deletion of most photoreceptor genes and auxin transporter genes might be linked to the unique phenotypes of fully mycoheterotrophic orchids. Conversely, trehalase genes that catalyse the conversion of trehalose into glucose have expanded in most sequenced orchids, in line with the fact that the germination of orchid non-endosperm seeds needs carbohydrates from fungi during the protocorm stage. We further show that the mature plant of P. guangdongensis, different from photosynthetic orchids, keeps expressing trehalase genes to hijack trehalose from fungi. Therefore, we propose that mycoheterotrophy in mature orchids is a continuation of the protocorm stage by sustaining the expression of trehalase genes. Our results shed light on the molecular mechanism underlying initial, partial and full mycoheterotrophy.

What Was Old Is New Again: The Pennate Diatom Haslea ostrearia (Gaillon) Simonsen in the Multi-Omic Age

The marine pennate diatom Haslea ostrearia has long been known for its characteristic blue pigment marennine, which is responsible for the greening of invertebrate gills, a natural phenomenon of great importance for the oyster industry. For two centuries, this taxon was considered unique; however, the recent description of a new blue Haslea species revealed unsuspected biodiversity. Marennine-like pigments are natural blue dyes that display various biological activities—e.g., antibacterial, antioxidant and antiproliferative—with a great potential for applications in the food, feed, cosmetic and health industries. Regarding fundamental prospects, researchers use model organisms as standards to study cellular and physiological processes in other organisms, and there is a growing and crucial need for more, new and unconventional model organisms to better correspond to the diversity of the tree of life. The present work, thus, advocates for establishing H. ostrearia as a new model organism by presenting its pros and cons—i.e., the interesting aspects of this peculiar diatom (representative of benthic-epiphytic phytoplankton, with original behavior and chemodiversity, controlled sexual reproduction, fundamental and applied-oriented importance, reference genome, and transcriptome will soon be available); it will also present the difficulties encountered before this becomes a reality as it is for other diatom models (the genetics of the species in its infancy, the transformation feasibility to be explored, the routine methods needed to cryopreserve strains of interest).

The spinach YY genome reveals sex chromosome evolution, domestication, and introgression history of the species

Background

Spinach (Spinacia oleracea L.) is a dioecious species with an XY sex chromosome system, but its Y chromosome has not been fully characterized. Our knowledge about the history of its domestication and improvement remains limited.

Results

A high-quality YY genome of spinach is assembled into 952 Mb in six pseudo-chromosomes. By a combination of genetic mapping, Genome-Wide Association Studies, and genomic analysis, we characterize a 17.42-Mb sex determination region (SDR) on chromosome 1. The sex chromosomes of spinach evolved when an insertion containing sex determination genes occurred, followed by a large genomic inversion about 1.98 Mya. A subsequent burst of SDR-specific repeats (0.1–0.15 Mya) explains the large size of this SDR. We identify a Y-specific gene, NRT1/PTR 6.4 which resides in this insertion, as a strong candidate for the sex determination or differentiation factor. Resequencing of 112 spinach genomes reveals a severe domestication bottleneck approximately 10.87 Kya, which dates the domestication of spinach 7000 years earlier than the archeological record. We demonstrate that a strong selection signal associated with internode elongation and leaf area expansion is associated with domestication of edibility traits in spinach. We find that several strong genomic introgressions from the wild species Spinacia turkestanica and Spinacia tetrandra harbor desirable alleles of genes related to downy mildew resistance, frost resistance, leaf morphology, and flowering-time shift, which likely contribute to spinach improvement.

Conclusions

Analysis of the YY genome uncovers evolutionary forces shaping nascent sex chromosome evolution in spinach. Our findings provide novel insights about the domestication and improvement of spinach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02633-x.

Reinvention of hermaphroditism via activation of a RADIALIS-like gene in hexaploid persimmon

In flowering plants, different lineages have independently transitioned from the ancestral hermaphroditic state into and out of various sexual systems¹. Polyploidizations are often associated with this plasticity in sexual systems^2,3. Persimmons (the genus Diospyros) have evolved dioecy via lineage-specific palaeoploidizations. More recently, hexaploid D. kaki has established monoecy and also exhibits reversions from male to hermaphrodite flowers in response to natural environmental signals (natural hermaphroditism, NH), or to artificial cytokinin treatment (artificial hermaphroditism, AH). We sought to identify the molecular pathways underlying these polyploid-specific reversions to hermaphroditism. Co-expression network analyses identified regulatory pathways specific to NH or AH transitions. Surprisingly, the two pathways appeared to be antagonistic, with abscisic acid and cytokinin signalling for NH and AH, respectively. Among the genes common to both pathways leading to hermaphroditic flowers, we identified a small-Myb RADIALIS-like gene, named DkRAD, which is specifically activated in hexaploid D. kaki. Consistently, ectopic overexpression of DkRAD in two model plants resulted in hypergrowth of the gynoecium. These results suggest that production of hermaphrodite flowers via polyploidization depends on DkRAD activation, which is not associated with a loss-of-function within the existing sex determination pathway, but rather represents a new path to (or reinvention of) hermaphroditism.

Identification of structural variation and polymorphisms of a sex co-segregating scaffold in spinach

Spinach is a common vegetable, and dioecy is maintained by a pair of XY sex chromosomes. Due to limited genomic resources and its highly repetitive genome, limited studies were conducted to investigate the genomic landscape of the region near sex-determining loci. In this study, we screened the structure variations (SVs) between Y-linked contigs and a 1.78-Mb X scaffold (Super_scaffold 66), which enabled the development of 12 sex co-segregating DNA markers. These markers were tested in one F₁ mapping population and 40 spinach accessions, which comprised 692 individual plants with the strong sex linkage pattern. In addition, we found that Super_scaffold 66 was highly repetitive along with the enriched LTR-RTs insertions and decreased microsatellite distribution compared with the rest genome, which matches extremely low gene density featured by only nine annotated genes. Synteny analysis between Y contigs and Superscaffold_66 revealed a 340-Kb accumulative Y contig (non-continuous) and a 500-Kb X counterpart along with SVs and wide-spread tandem duplications. Among the nine genes, one ABC transporter gene revealed noticeable SVs between Y contig and X counterpart, as an approximate 5-Kb recent Gypsy LTR-RT insertion in the Y-linked allele, but not the X allele. The gene paucity, SVs, and sex-linked polymorphisms attributed to the recombination suppression. We proposed that Super_scaffold 66 is part of the non-recombining region containing the sex determination genes. The spread of 12 sex co-segregating markers from this 1.78 Mb genomic region indicated the existence and expansion of sex determination region during progression of the Y chromosome.

Comparative genomics reveal the convergent evolution of CYP82D and CYP706X members related to flavone biosynthesis in Lamiaceae and Asteraceae

Distant species producing the same secondary metabolites is an interesting and common phenomenon in nature. A classic example of this is scutellarein whose derivatives have been used clinically for more than 30 years. Scutellarein occurs in significant amounts in species of two different orders, Scutellaria baicalensis and Erigeron breviscapus, which diverged more than 100 million years ago. Here, according to the genome-wide selection and functional identification of 39 CYP450 genes from various angiosperms, we confirmed that only seven Scutellaria-specific CYP82D genes and one Erigeron CYP706X gene could perform the catalytic activity of flavone 6-hydroxylase (F6H), suggesting that the convergent evolution of scutellarein production in these two distant species was caused by two independently evolved CYP450 families. We also identified seven Scutellaria-specific CYP82D genes encoding flavone 8-hydroxylase (F8H). The evolutionary patterns of CYP82 and CYP706 families via kingdom-wide comparative genomics highlighted the evolutionary diversity of CYP82D and the specificity of CYP706X in angiosperms. Multi-collinearity and phylogenetic analysis of CYP82D in Scutellaria confirmed that the function of F6H evolved from F8H. Furthermore, the SbaiCYP82D1^A319D , EbreCYP706X^R130A , EbreCYP706X^F312D and EbreCYP706X^A318D mutants can significantly decrease the catalytic activity of F6H, revealing the contribution of crucial F6H amino acids to the scutellarein biosynthesis of distant species. This study provides important insights into the multi-origin evolution of the same secondary metabolite biosynthesis in the plant kingdom.

Comparative Transcriptome Analysis Reveals Sex-Biased Expression of Hormone-Related Genes at an Early Stage of Sex Differentiation in Red Bayberry (Morella rubra)

The molecular mechanism of sex development and differentiation in the economically important dioecious fruit tree, red bayberry (Morella rubra), was revealed using next-generation transcriptome sequencing (NGS), and comparative analyses were used to identify differentially expressed genes (DEGs) in female and male flower buds. A total of 7,029 of these DEGs were identified at two early development stages. KEGG pathway enrichment analysis revealed that plant hormone signal transduction was significantly overrepresented, and 91 genes related to hormones were identified. An analysis of 7,029 DEGs revealed 161 hormone-related genes, with the 42 related to auxin and 26 related to ethylene being the most highly represented. A total of 62 genes were significantly up-regulated in females and 29 were in males, with 18 of them specifically expressed in females and 10 in males. A total of 415 transcription factors were identified, with 129 genes up-regulated in females and 53 in males. Moreover, 38 had female-specific expression and 18 had male-specific expression. Using weighted gene co-expression network analysis (WGCNA), two modules were found to be associated with sexual type. In the module coded light-green, there were five genes related to hormones, one to flower development and ten transcription factors with four genes specifically expressed in the males and four in females. The hub gene in the light-green module is MR0TCONS_00017483.1 (ACO), which is involved in ethylene biosynthesis and had male-specific expression. Among the transcription factors, three of the four male-specific expressed genes involved in flavonoid biosynthesis, including the MYB gene MR1TCONS_00020658.1 and two BHLH genes, MR6G001563.1 and MR8G020751.1, played important roles in male floral differentiation. In the dark-cyan module, six hormone-related genes, five transcription factors and three flower development genes were identified with the hub gene MR1G019545.1 (ETR1), which participates in the ethylene signaling pathway, and MR4G023618.1, which encodes the C3H zinc finger transcription factor. These results indicate that ethylene is the key hormone that interacts with other hormones and transcription factors to regulate sex differentiation in the red bayberry, which also provides new insights into the mechanism of sex determination and differentiation in the red bayberry.

A derived ZW chromosome system in Amborella trichopoda, representing the sister lineage to all other extant flowering plants

The genetic basis and evolution of sex determination in dioecious plants is emerging as an active area of research with exciting advances in genome sequencing and analysis technologies. As the sole species within the sister lineage to all other extant flowering plants, Amborella trichopoda is an important model for understanding the evolution and development of flowers. Plants typically produce only male or female flowers, but sex determination mechanisms are unknown for the species. Sequence data derived from plants of natural origin and an F1 mapping population were used to identify sex-linked genes and the nonrecombining region. Amborella trichopoda has a ZW sex determination system. Analysis of genes in a 4 Mb nonrecombining sex-determination region reveals recent divergence of Z and W gametologs, and few Z- and W-specific genes. The sex chromosomes of A. trichopoda evolved less than 16.5 Myr ago, long after the divergence of the extant angiosperms.

Insights from the analysis of draft genome sequence of Crocus sativus L

Saffron (Crocus sativus L.) is the low yielding plant of medicinal and economic importance. Therefore, it is of interest to report the draft genome sequence of C. sativus. The draft genome of C. sativus has been assembled using Illumina sequencing and is 3.01 Gb long covering 84.24% of genome. C. sativus genome annotation identified 53,546 functional genes (including 5726 transcription factors), 862,275 repeats and 964,231 SSR markers. The genes involved in the apocarotenoids biosynthesis pathway (crocin, crocetin, picrocrocin, and safranal) were found in the draft genome analysis.

Chromosome-scale haplotype-phased genome assemblies of the male and female lines of wild asparagus (Asparagus kiusianus), a dioecious plant species

Asparagus kiusianus is a disease-resistant dioecious plant species and a wild relative of garden asparagus (Asparagus officinalis). To enhance A. kiusianus genomic resources, advance plant science, and facilitate asparagus breeding, we determined the genome sequences of the male and female lines of A. kiusianus. Genome sequence reads obtained with a linked-read technology were assembled into four haplotype-phased contig sequences (∼1.6 Gb each) for the male and female lines. The contig sequences were aligned onto the chromosome sequences of garden asparagus to construct pseudomolecule sequences. Approximately 55,000 potential protein-encoding genes were predicted in each genome assembly, and ∼70% of the genome sequence was annotated as repetitive. Comparative analysis of the genomes of the two species revealed structural and sequence variants between the two species as well as between the male and female lines of each species. Genes with high sequence similarity with the male-specific sex determinant gene in A. officinalis, MSE1/AoMYB35/AspTDF1, were presented in the genomes of the male line but absent from the female genome assemblies. Overall, the genome sequence assemblies, gene sequences, and structural and sequence variants determined in this study will reveal the genetic mechanisms underlying sexual differentiation in plants, and will accelerate disease-resistance breeding in garden asparagus.

Transcriptomic Analysis of Sex-Associated DEGs in Female and Male Flowers of Kiwifruit (Actinidia deliciosa [A. Chev] C. F. Liang & A. R. Ferguson)

Kiwifruit (Actinidia deliciosa [A. Chev.], C.V. Liang & A. R. Ferguson, 1984) is a perennial plant, with morphologically hermaphroditic and functionally dioecious flowers. Fruits of this species are berries of great commercial and nutritional importance. Nevertheless, few studies have analyzed the molecular mechanisms involved in sexual differentiation in this species. To determine these mechanisms, we performed RNA-seq in floral tissue at stage 60 on the BBCH scale in cultivar ‘Hayward’ (H, female) and a seedling from ‘Green Light’ × ‘Tomuri’ (G × T, male). From these analyses, we obtained expression profiles of 24,888 (H) and 27,027 (G × T) genes, of which 6413 showed differential transcript abundance. Genetic ontology (GO) and KEGG analysis revealed activation of pathways associated with the translation of hormonal signals, plant-pathogen interaction, metabolism of hormones, sugars, and nucleotides. The analysis of the protein-protein interaction network showed that the genes ERL1, AG, AGL8, LFY, WUS, AP2, WRKY, and CO, are crucial elements in the regulation of the hormonal response for the formation and development of anatomical reproductive structures and gametophytes. On the other hand, genes encoding four Putative S-adenosyl-L-methionine-dependent methyltransferases (Achn201401, Achn281971, Achn047771 and Achn231981) were identified, which were up-regulated mainly in the male flowers. Moreover, the expression profiles of 15 selected genes through RT-qPCR were consistent with the results of RNA-seq. Finally, this work provides gene expression-based interactions between transcription factors and effector genes from hormonal signaling pathways, development of floral organs, biological and metabolic processes or even epigenetic mechanisms which could be involved in the kiwi sex-determination. Thus, in order to decode the nature of these interactions, it could be helpful to propose new models of flower development and sex determination in the Actinidia genus.

VviPLATZ1 is a major factor that controls female flower morphology determination in grapevine

Plant genetic sex determinants that mediate the transition to dioecy are predicted to be diverse, as this type of mating system independently evolved multiple times in angiosperms. Wild Vitis species are dioecious with individuals producing morphologically distinct female or male flowers; whereas, modern domesticated Vitis vinifera cultivars form hermaphrodite flowers capable of self-pollination. Here, we identify the VviPLATZ1 transcription factor as a key candidate female flower morphology factor that localizes to the Vitis SEX-DETERMINING REGION. The expression pattern of this gene correlates with the formation reflex stamens, a prominent morphological phenotype of female flowers. After generating CRISPR/Cas9 gene-edited alleles in a hermaphrodite genotype, phenotype analysis shows that individual homozygous lines produce flowers with reflex stamens. Taken together, our results demonstrate that loss of VviPLATZ1 function is a major factor that controls female flower morphology in Vitis.

Unlike wild Vitis species, which produce either female or male flowers, modern grapevine cultivars form hermaphrodite flowers for self-pollination. Here, the authors report that the VviPLATZ1 (plant AT-rich sequence-and zinc-binding protein1) transcription factor functions in controlling female flower morphology determination.

Evolution and expression of the MADS-box flowering transition genes AGAMOUS-like 24/SHORT VEGETATIVE PHASE with emphasis in selected Neotropical orchids

In angiosperms the reproductive transition results in the transformation of a vegetative apical meristem (SAM) into an inflorescence meristem (IM), capable of forming floral meristems (FM). Two key players in the flowering transition are AGAMOUS-like 24 (AGL24) and SHORT VEGETATIVE PHASE (SVP). They are eudicot MADS-box paralogs performing opposite roles, as AGL24 positively regulates flowering while SVP represses the reproductive transition in Arabidopsis. We confirm that the Arabidopsis functional reference cannot be readily extrapolated to all eudicots as there are additional duplications of AGL24 in early divergent eudicots and core eudicots with significant sequence variation. In addition, we found that in monocots, two additional independent duplication events have resulted in at least three clades of AGL24/SVP homologs, some only found in Orchidaceae. Protein sequence analyses and comparative evolutionary rates point to higher rates of relaxed negative selection in the Core Eudicot AGL24 B and the Orch SVP-like B clades, in eudicots and monocots respectively. On the other hand, expression data points to plesiomorphic pleiotropic roles of AGL24/SVP genes likely similar to SVP core eudicot genes, and the acquisition of new roles as flowering positive regulators in Core Eudicot AGL24 A genes. Our research presents evidence on the diversification and recruitment of AGL24/SVP homologs in flowering transition in orchids. Although, broad expression of most copies does not allow to determine if they act as flowering repressors or promoters, the restricted expression of some homologs in the SAM suggests putative roles in maintaining the vegetative phase. If so studying in detail the function of AGL24/SVP homologs in orchids is critical to identify putative flowering repressors in a lineage where other canonical repressors remain elusive.

High-quality reference genome sequences of two coconut cultivars provide insights into evolution of monocot chromosomes and differentiation of fiber content and plant height

BACKGROUND

Coconut is an important tropical oil and fruit crop whose evolutionary position renders it a fantastic species for the investigation of the evolution of monocot chromosomes and the subsequent differentiation of ancient plants.

RESULTS

Here, we report the assembly and annotation of reference-grade genomes of Cn. tall and Cn. dwarf, whose genome sizes are 2.40 Gb and 2.39 Gb, respectively. The comparative analysis reveals that the two coconut subspecies diverge about 2-8 Mya while the conserved Arecaceae-specific whole-genome duplication (ω WGD) occurs approximately 47-53 Mya. It additionally allows us to reconstruct the ancestral karyotypes of the ten ancient monocot chromosomes and the evolutionary trajectories of the 16 modern coconut chromosomes. Fiber synthesis genes in Cn. tall, related to lignin and cellulose synthesis, are found at a higher copy number and expression level than dwarf coconuts. Integrated multi-omics analysis reveals that the difference in coconut plant height is the result of altered gibberellin metabolism, with both the GA20ox copy number and a single-nucleotide change in the promoter together leading to the difference in plant height between Cn. tall and Cn. dwarf.

CONCLUSION

We provide high-quality coconut genomes and reveal the genetic basis of trait differences between two coconuts through multi-omics analysis. We also reveal that the selection of plant height has been targeted for the same gene for millions of years, not only in natural selection of ancient plant as illustrated in coconut, but also for artificial selection in cultivated crops such as rice and maize.

Transcriptome analysis and genetic diversity of Allium victorialis germplasms from the Changbai Mountains

The Changbai Mountains comprise one of the main distribution areas of A. victorialis in China, and this species is endangered owing to habitat changes and overexploitation. However, A. victorialis germplasms have not been systematically collected and studied. The aims of this study were to obtain some detailed genetic information, analyze the genetic diversity, and further promote the protection of A. victorialis germplasms from the Changbai Mountains. Transcriptomic analysis was performed with six A. victorialis samples collected from the Changbai Mountains. At least 146,759 genes for each sample were obtained after performing de novo assembly of the RNA-seq data, and at least 92% of these genes were found to have only one mRNA isoform. These sequences and their functional annotations provided a large-scale genetic resource of this species. Phylogenetic analysis showed that A. victorialis was genetically distant from some related species, e.g. Allium sativum, Allium fistulosum, and Allium cepa, but genetically close to Allium tuberosum. The two A. victorialis var. listera samples were phylogenetically separated from the other four samples, and these two samples should be regarded as Allium listera. In addition, two KASP markers for discriminating the Dongfeng samples from the other four A. victorialis samples were successfully developed. This study lays the foundation for future studies on the genetic diversity and evolution of Allium species, as well as for the conservation of A. victorialis germplasms from the Changbai Mountains and other populations of this species.

Insights from the first genome assembly of Onion (Allium cepa)

Onion is an important vegetable crop with an estimated genome size of 16 Gb. We describe the de novo assembly and ab initio annotation of the genome of a doubled haploid onion line DHCU066619, which resulted in a final assembly of 14.9 Gb with an N50 of 464 Kb. Of this, 2.4 Gb was ordered into eight pseudomolecules using four genetic linkage maps. The remainder of the genome is available in 89.6 K scaffolds. Only 72.4% of the genome could be identified as repetitive sequences and consist, to a large extent, of (retro) transposons. In addition, an estimated 20% of the putative (retro) transposons had accumulated a large number of mutations, hampering their identification, but facilitating their assembly. These elements are probably already quite old. The ab initio gene prediction indicated 540,925 putative gene models, which is far more than expected, possibly due to the presence of pseudogenes. Of these models, 47,066 showed RNASeq support. No gene rich regions were found, genes are uniformly distributed over the genome. Analysis of synteny with Allium sativum (garlic) showed collinearity but also major rearrangements between both species. This assembly is the first high-quality genome sequence available for the study of onion and will be a valuable resource for further research.

Insights from the first genome assembly of Onion (Allium cepa)

Onion is an important vegetable crop with an estimated genome size of 16 Gb. We describe the de novo assembly and ab initio annotation of the genome of a doubled haploid onion line DHCU066619, which resulted in a final assembly of 14.9 Gb with an N50 of 464 Kb. Of this, 2.4 Gb was ordered into eight pseudomolecules using four genetic linkage maps. The remainder of the genome is available in 89.6 K scaffolds. Only 72.4% of the genome could be identified as repetitive sequences and consist, to a large extent, of (retro) transposons. In addition, an estimated 20% of the putative (retro) transposons had accumulated a large number of mutations, hampering their identification, but facilitating their assembly. These elements are probably already quite old. The ab initio gene prediction indicated 540,925 putative gene models, which is far more than expected, possibly due to the presence of pseudogenes. Of these models, 47,066 showed RNASeq support. No gene rich regions were found, genes are uniformly distributed over the genome. Analysis of synteny with Allium sativum (garlic) showed collinearity but also major rearrangements between both species. This assembly is the first high-quality genome sequence available for the study of onion and will be a valuable resource for further research.

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.

Plant genera Cannabis and Humulus share the same pair of well-differentiated sex chromosomes

We recently described, in Cannabis sativa, the oldest sex chromosome system documented so far in plants (12-28 Myr old). Based on the estimated age, we predicted that it should be shared by its sister genus Humulus, which is known also to possess XY chromosomes. Here, we used transcriptome sequencing of an F₁ family of H. lupulus to identify and study the sex chromosomes in this species using the probabilistic method SEX-DETector. We identified 265 sex-linked genes in H. lupulus, which preferentially mapped to the C. sativa X chromosome. Using phylogenies of sex-linked genes, we showed that a region of the sex chromosomes had already stopped recombining in an ancestor of both species. Furthermore, as in C. sativa, Y-linked gene expression reduction is correlated to the position on the X chromosome, and highly Y degenerated genes showed dosage compensation. We report, for the first time in Angiosperms, a sex chromosome system that is shared by two different genera. Thus, recombination suppression started at least 21-25 Myr ago, and then (either gradually or step-wise) spread to a large part of the sex chromosomes (c. 70%), leading to a degenerated Y chromosome.

The timing of genetic degeneration of sex chromosomes

Genetic degeneration is an extraordinary feature of sex chromosomes, with the loss of functions of Y-linked genes in species with XY systems, and W-linked genes in ZW systems, eventually affecting almost all genes. Although degeneration is familiar to most biologists, important aspects are not yet well understood, including how quickly a Y or W chromosome can become completely degenerated. I review the current understanding of the time-course of degeneration. Degeneration starts after crossing over between the sex chromosome pair stops, and theoretical models predict an initially fast degeneration rate and a later much slower one. It has become possible to estimate the two quantities that the models suggest are the most important in determining degeneration rates-the size of the sex-linked region, and the time when recombination became suppressed (which can be estimated using Y-X or W-Z sequence divergence). However, quantifying degeneration is still difficult. I review evidence on gene losses (based on coverage analysis) or loss of function (by classifying coding sequences into functional alleles and pseudogenes). I also review evidence about whether small genome regions degenerate, or only large ones, whether selective constraints on the genes in a sex-linked region also strongly affect degeneration rates, and about how long it takes before all (or almost all) genes are lost. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Sex-linked gene expression and the emergence of hermaphrodites in Carica papaya

PREMISE

One evolutionary path from hermaphroditism to dioecy is via a gynodioecious intermediate. The evolution of dioecy may also coincide with the formation of sex chromosomes that possess sex-determining loci that are physically linked in a region of suppressed recombination. Dioecious papaya (Carica papaya) has an XY chromosome system, where the presence of a Y chromosome determines maleness. However, in cultivation, papaya is gynodioecious, due to the conversion of the male Y chromosome to a hermaphroditic Y^h chromosome during its domestication.

METHODS

We investigated gene expression linked to the X, Y, and Y^h chromosomes at different floral developmental stages to identify differentially expressed genes that may be involved in the sexual transition of males to hermaphrodites.

RESULTS

We identified 309 sex-biased genes found on the sex chromosomes, most of which are found in the pseudoautosomal regions. Female (XX) expression in the sex-determining region was almost double that of X-linked expression in males (XY) and hermaphrodites (XY^h ), which rules out dosage compensation for most sex-linked genes; although, an analysis of hemizygous X-linked loci found evidence of partial dosage compensation. Furthermore, we identified a candidate gene associated with sex determination and the transition to hermaphroditism, a homolog of the MADS-box protein SHORT VEGETATIVE PHASE.

CONCLUSIONS

We identified a pattern of partial dosage compensation for hemizygous genes located in the papaya sex-determining region. Furthermore, we propose that loss-of-expression of the Y-linked SHORT VEGETATIVE PHASE homolog facilitated the transition from males to hermaphrodites in papaya.

TGIF-DB: terse genomics interface for developing botany

Objectives

Pearl millet (Pennisetum glaucum) is a staple cereal crop for semi-arid regions. Its whole genome sequence and deduced putative gene sequences are available. However, the functions of many pearl millet genes are unknown. Situations are similar for other crop species such as garden asparagus (Asparagus officinalis), chickpea (Cicer arietinum) and Tartary buckwheat (Fagopyrum tataricum). The objective of the data presented here was to improve functional annotations of genes of pearl millet, garden asparagus, chickpea and Tartary buckwheat with gene annotations of model plants, to systematically provide such annotations as well as their sequences on a website, and thereby to promote genomics for those crops.

Data description

Sequences of genomes and transcripts of pearl millet, garden asparagus, chickpea and Tartary buckwheat were downloaded from a public database. These transcripts were associated with functional annotations of their Arabidopsis thaliana and rice (Oryza sativa) counterparts identified by BLASTX. Conserved domains in protein sequences of those species were identified by the HMMER scan with the Pfam database. The resulting data was deposited in the figshare repository and can be browsed on the Terse Genomics Interface for Developing Botany (TGIF-DB) website (http://webpark2116.sakura.ne.jp/rlgpr/).

Transcriptome landscapes of multiple tissues highlight the genes involved in the flavor metabolic pathway in Chinese chive (Allium tuberosum)

The unique flavor of Allium tuberosum is primarily associated with the hydrolysis of a series of organosulfur compounds, S-alk(en)yl cysteine sulphoxides (CSOs), upon tissue bruising or maceration. To obtain the tissue-specific transcriptomes, 18 RNA-Seq libraries representing leaf, root, stem, mature flower, inflorescence, and seed tissues of A. tuberosum were sequenced, finally yielding 133.7 Gb clean reads. The de novo assembled transcriptomes enabled the identification of 223,529 unigenes, which were functionally annotated and analyzed for the gene ontology and metabolic pathways. Furthermore, to reveal the flavor metabolic pathways, a total of 205 unigenes involved in the sulfur assimilation and CSO biosynthesis were identified, and their expression profiles were analyzed by RNA-Seq and qRT-PCR. Collectively, this study provides a valuable resource for in-depth molecular and functional researches especially on flavor formation, as well as for the development of molecular markers, and other genetic studies in A. tuberosum.

Reproductive phasiRNA loci and DICER-LIKE5, but not microRNA loci, diversified in monocotyledonous plants

In monocots other than maize (Zea mays) and rice (Oryza sativa), the repertoire and diversity of microRNAs (miRNAs) and the populations of phased, secondary, small interfering RNAs (phasiRNAs) are poorly characterized. To remedy this, we sequenced small RNAs (sRNA) from vegetative and dissected inflorescence tissue in 28 phylogenetically diverse monocots and from several early-diverging angiosperm lineages, as well as publicly available data from 10 additional monocot species. We annotated miRNAs, small interfering RNAs (siRNAs) and phasiRNAs across the monocot phylogeny, identifying miRNAs apparently lost or gained in the grasses relative to other monocot families, as well as a number of transfer RNA fragments misannotated as miRNAs. Using our miRNA database cleaned of these misannotations, we identified conservation at the 8th, 9th, 19th, and 3'-end positions that we hypothesize are signatures of selection for processing, targeting, or Argonaute sorting. We show that 21-nucleotide (nt) reproductive phasiRNAs are far more numerous in grass genomes than other monocots. Based on sequenced monocot genomes and transcriptomes, DICER-LIKE5, important to 24-nt phasiRNA biogenesis, likely originated via gene duplication before the diversification of the grasses. This curated database of phylogenetically diverse monocot miRNAs, siRNAs, and phasiRNAs represents a large collection of data that should facilitate continued exploration of sRNA diversification in flowering plants.

Cloning and physical localization of male-biased repetitive DNA sequences in Spinacia oleracea (Amaranthaceae)

Spinach (Spinacia oleracea Linnaeus, 1753) is an ideal material for studying molecular mechanisms of early-stage sex chromosome evolution in dioecious plants. Degenerate oligonucleotide-primed polymerase chain reaction (DOP-PCR) technique facilitates the retrotransposon-relevant studies by enriching specific repetitive DNA sequences from a micro-dissected single chromosome. We conducted genomic subtractive hybridization to screen sex-biased DNA sequences by using the DOP-PCR amplification products of micro-dissected spinach Y chromosome. The screening yielded 55 male-biased DNA sequences with 30 576 bp in length, of which, 32 DNA sequences (12 049 bp) contained repeat DNA sequences, including LTR/Copia, LTR/Gypsy, simple repeats, and DNA/CMC-EnSpm. Among these repetitive DNA sequences, four DNA sequences that contained a fragment of Ty3-gypsy retrotransposons (SP73, SP75, SP76, and SP77) were selected as fluorescence probes to hybridization on male and female spinach karyotypes. Fluorescence in situ hybridization (FISH) signals of SP73 and SP75 were captured mostly on the centromeres and their surrounding area for each homolog. Hybridization signals primarily appeared near the putative centromeres for each homologous chromosome pair by using SP76 and SP77 probes for FISH, and sporadic signals existed on the long arms. Results can be served as a basis to study the function of repetitive DNA sequences in sex chromosome evolution in spinach.

Epigenetics drive the evolution of sex chromosomes in animals and plants

We review how epigenetics affect sex chromosome evolution in animals and plants. In a few species, sex is determined epigenetically through the action of Y-encoded small RNAs. Epigenetics is also responsible for changing the sex of individuals through time, even in species that carry sex chromosomes, and could favour species adaptation through breeding system plasticity. The Y chromosome accumulates repeats that become epigenetically silenced which leads to an epigenetic conflict with the expression of Y genes and could accelerate Y degeneration. Y heterochromatin can be lost through ageing, which activates transposable elements and lowers male longevity. Y chromosome degeneration has led to the evolution of meiotic sex chromosome inactivation in eutherians (placentals) and marsupials, and dosage compensation mechanisms in animals and plants. X-inactivation convergently evolved in eutherians and marsupials via two independently evolved non-coding RNAs. In Drosophila, male X upregulation by the male specific lethal (MSL) complex can spread to neo-X chromosomes through the transposition of transposable elements that carry an MSL-binding motif. We discuss similarities and possible differences between plants and animals and suggest future directions for this dynamic field of research. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.

Inferring the Genetic Basis of Sex Determination from the Genome of a Dioecious Nightshade

Dissecting the genetic mechanisms underlying dioecy (i.e., separate female and male individuals) is critical for understanding the evolution of this pervasive reproductive strategy. Nonetheless, the genetic basis of sex determination remains unclear in many cases, especially in systems where dioecy has arisen recently. Within the economically important plant genus Solanum (∼2,000 species), dioecy is thought to have evolved independently at least 4 times across roughly 20 species. Here, we generate the first genome sequence of a dioecious Solanum and use it to ascertain the genetic basis of sex determination in this species. We de novo assembled and annotated the genome of Solanum appendiculatum (assembly size: ∼750 Mb scaffold N50: 0.92 Mb; ∼35,000 genes), identified sex-specific sequences and their locations in the genome, and inferred that males in this species are the heterogametic sex. We also analyzed gene expression patterns in floral tissues of males and females, finding approximately 100 genes that are differentially expressed between the sexes. These analyses, together with observed patterns of gene-family evolution specific to S. appendiculatum, consistently implicate a suite of genes from the regulatory network controlling pectin degradation and modification in the expression of sex. Furthermore, the genome of a species with a relatively young sex-determination system provides the foundational resources for future studies on the independent evolution of dioecy in this clade.

Datura genome reveals duplications of psychoactive alkaloid biosynthetic genes and high mutation rate following tissue culture

BACKGROUND

Datura stramonium (Jimsonweed) is a medicinally and pharmaceutically important plant in the nightshade family (Solanaceae) known for its production of various toxic, hallucinogenic, and therapeutic tropane alkaloids. Recently, we published a tissue-culture based transformation protocol for D. stramonium that enables more thorough functional genomics studies of this plant. However, the tissue culture process can lead to undesirable phenotypic and genomic consequences independent of the transgene used. Here, we have assembled and annotated a draft genome of D. stramonium with a focus on tropane alkaloid biosynthetic genes. We then use mRNA sequencing and genome resequencing of transformants to characterize changes following tissue culture.

RESULTS

Our draft assembly conforms to the expected 2 gigabasepair haploid genome size of this plant and achieved a BUSCO score of 94.7% complete, single-copy genes. The repetitive content of the genome is 61%, with Gypsy-type retrotransposons accounting for half of this. Our gene annotation estimates the number of protein-coding genes at 52,149 and shows evidence of duplications in two key alkaloid biosynthetic genes, tropinone reductase I and hyoscyamine 6 β-hydroxylase. Following tissue culture, we detected only 186 differentially expressed genes, but were unable to correlate these changes in expression with either polymorphisms from resequencing or positional effects of transposons.

CONCLUSIONS

We have assembled, annotated, and characterized the first draft genome for this important model plant species. Using this resource, we show duplications of genes leading to the synthesis of the medicinally important alkaloid, scopolamine. Our results also demonstrate that following tissue culture, mutation rates of transformed plants are quite high (1.16 × 10- 3 mutations per site), but do not have a drastic impact on gene expression.

Chromosome-scale assembly of the genome of Salix dunnii reveals a male-heterogametic sex determination system on chromosome 7

Sex determination systems in plants can involve either female or male heterogamety (ZW or XY, respectively). Here we used Illumina short reads, Oxford Nanopore Technologies (ONT) long reads and Hi-C reads to assemble the first chromosome-scale genome of a female willow tree (Salix dunnii), and to predict genes using transcriptome sequences and available databases. The final genome sequence of 328 Mb in total was assembled in 29 scaffolds, and includes 31,501 predicted genes. Analyses of short-read sequence data that included female and male plants suggested a male heterogametic sex-determining factor on chromosome 7, implying that, unlike the female heterogamety of most species in the genus Salix, male heterogamety evolved in the subgenus Salix. The S. dunnii sex-linked region occupies about 3.21 Mb of chromosome 7 in females (representing its position in the X chromosome), probably within a pericentromeric region. Our data suggest that this region is enriched for transposable element insertions, and about one-third of its 124 protein-coding genes were gained via duplications from other genome regions. We detect purifying selection on the genes that were ancestrally present in the region, though some have been lost. Transcriptome data from female and male individuals show more male- than female-biased genes in catkin and leaf tissues, and indicate enrichment for male-biased genes in the pseudo-autosomal regions. Our study provides valuable genomic resources for further studies of sex-determining regions in the family Salicaceae, and sex chromosome evolution.