Plant contributions to our understanding of sex chromosome evolution

Evolution of a plant sex chromosome driven by expanding pericentromeric recombination suppression

Recombination suppression around sex-determining gene(s) is a key step in evolution of sex chromosomes, but it is not well understood how it evolves. Recently evolved sex-linked regions offer an opportunity to understand the mechanisms of recombination cessation. This paper analyses such a region on Silene latifolia (Caryophyllaceae) sex chromosomes, where recombination was suppressed in the last 120 thousand years ("stratum 3"). Locating the boundaries of the stratum 3 in S. latifolia genome sequence revealed that this region is far larger than assumed previously-it is about 14 Mb long and includes 202 annotated genes. A gradient of X:Y divergence detected in the stratum 3, with divergence increasing proximally, indicates gradual recombination cessation, possibly caused by expansion of pericentromeric recombination suppression (PRS) into the pseudoautosomal region. Expansion of PRS was also the likely cause for the formation of the older stratum 2 on S. latifolia sex chromosomes. The role of PRS in sex chromosome evolution has been underappreciated, but it may be a significant factor, especially in the species with large chromosomes where PRS is often extensive.

Differentially methylated genes involved in reproduction and ploidy levels in recent diploidized and tetraploidized Eragrostis curvula genotypes

Epigenetics studies changes in gene activity without changes in the DNA sequence. Methylation is an epigenetic mechanism important in many pathways, such as biotic and abiotic stresses, cell division, and reproduction. Eragrostis curvula is a grass species reproducing by apomixis, a clonal reproduction by seeds. This work employed the MCSeEd technique to identify deferentially methylated positions, regions, and genes in the CG, CHG, and CHH contexts in E. curvula genotypes with similar genomic backgrounds but with different reproductive modes and ploidy levels. In this way, we focused the analysis on the cvs. Tanganyika INTA (4x, apomictic), Victoria (2x, sexual), and Bahiense (4x, apomictic). Victoria was obtained from the diploidization of Tanganyika INTA, while Bahiense was produced from the tetraploidization of Victoria. This study showed that polyploid/apomictic genotypes had more differentially methylated positions and regions than the diploid sexual ones. Interestingly, it was possible to observe fewer differentially methylated positions and regions in CG than in the other contexts, meaning CG methylation is conserved across the genotypes regardless of the ploidy level and reproductive mode. In the comparisons between sexual and apomictic genotypes, we identified differentially methylated genes involved in the reproductive pathways, specifically in meiosis, cell division, and fertilization. Another interesting observation was that several differentially methylated genes between the diploid and the original tetraploid genotype recovered their methylation status after tetraploidization, suggesting that methylation is an important mechanism involved in reproduction and ploidy changes.

Evolution of a ZW sex chromosome system in willows

Transitions in the heterogamety of sex chromosomes (e.g., XY to ZW or vice versa) fundamentally alter the genetic basis of sex determination, however the details of these changes have been studied in only a few cases. In an XY to ZW transition, the X is likely to give rise to the W because they both carry feminizing genes and the X is expected to harbour less genetic load than the Y. Here, using a new reference genome for Salix exigua, we trace the X, Y, Z, and W sex determination regions during the homologous transition from an XY system to a ZW system in willow (Salix). We show that both the W and the Z arose from the Y chromosome. We find that the new Z chromosome shares multiple homologous putative masculinizing factors with the ancestral Y, whereas the new W lost these masculinizing factors and gained feminizing factors. The origination of both the W and Z from the Y was permitted by an unexpectedly low genetic load on the Y and this indicates that the origins of sex chromosomes during homologous transitions may be more flexible than previously considered.

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

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

A rapid method for assembly of single chromosome and identification of sex determination region based on single-chromosome sequencing

The sex-determining-region (SDR) may offer the best prospects for studying sex-determining gene, recombination suppression, and chromosome heteromorphism. However, current progress of SDR identification and cloning showed following shortcomings: large near-isogenic lines need to be constructed, and a relatively large population is needed; the cost of whole-genome sequencing and assembly is high. Herein, the X/Y chromosomes of Spinacia oleracea L. subsp. turkestanica were successfully microdissected and assembled using single-chromosome sequencing. The assembly length of X and Y chromosome is c. 192.1 and 195.2 Mb, respectively. Three large inversions existed between X and Y chromosome. The SDR size of X and Y chromosome is c. 13.2 and 24.1 Mb, respectively. MSY region and six male-biased genes were identified. A Y-chromosome-specific marker in SDR was constructed and used to verify the chromosome assembly quality at cytological level via fluorescence in situ hybridization. Meanwhile, it was observed that the SDR located on long arm of Y chromosome and near the centromere. Overall, a technical system was successfully established for rapid cloning the SDR and it is also applicable to rapid assembly of specific chromosome in other plants. Furthermore, this study laid a foundation for studying the molecular mechanism of sex chromosome evolution in spinach.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S. latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S. latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.

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.

The male-heterogametic sex determination system on chromosome 15 of Salix triandra and Salix arbutifolia reveals ancestral male heterogamety and subsequent turnover events in the genus Salix

Dioecious Salix evolved more than 45 million years ago, but have homomorphic sex chromosomes, suggesting that turnover event(s) prevented major differentiation. Sex chromosome turnover events have been inferred in the sister genus Populus. The genus Salix includes two main clades, Salix and Vetrix, with several previously studied Vetrix clade species having female-heterogametic (ZW) or male-heterogametic (XY) sex-determining systems (SDSs) on chromosome 15, while three Salix clade species have XY SDSs on chromosome 7. We here studied two basal taxa of the Vetrix clade, S. arbutifolia and S. triandra using S. purpurea as the reference genome. Analyses of whole genome resequencing data for genome-wide associations (GWAS) with the sexes and genetic differentiation between the sexes (FST values) showed that both species have male heterogamety with a sex-determining locus on chromosome 15, suggesting an early turnover event within the Vetrix clade, perhaps promoted by sexually antagonistic or (and) sex-ratio selection. Changepoint analysis based on FST values identified small sex-linked regions of ~3.33 Mb and ~2.80 Mb in S. arbutifolia and S. triandra, respectively. The SDS of S. arbutifolia was consistent with recent results that used its own genome as reference. Ancestral state reconstruction of SDS suggests that at least two turnover events occurred in Salix.

Heterochiasmy and Sex Chromosome Evolution in Silene

The evolution of a non-recombining sex-specific region is a key step in sex chromosome evolution. Suppression of recombination between the (proto-) X- and Y-chromosomes in male meiosis creates a non-recombining Y-linked region (NRY), while the X-chromosome continues to recombine in females. Lack of recombination in the NRY defines its main properties—genetic degeneration and accumulation of repetitive DNA, making X and Y chromosomes very different from each other. How and why recombination suppression on sex chromosomes evolves remains controversial. A strong difference in recombination rates between the sexes (heterochiasmy) can facilitate or even cause recombination suppression. In the extreme case—complete lack of recombination in the heterogametic sex (achiasmy)—the entire sex-specific chromosome is automatically non-recombining. In this study, I analyse sex-specific recombination rates in a dioecious plant Silene latifolia (Caryophyllaceae), which evolved separate sexes and sex chromosomes ~11 million years ago. I reconstruct high-density RNAseq-based genetic maps including over five thousand genic markers for the two sexes separately. The comparison of the male and female maps reveals only modest heterochiasmy across the genome, with the exception of the sex chromosomes, where recombination is suppressed in males. This indicates that heterochiasmy likely played only a minor, if any, role in NRY evolution in S. latifolia, as recombination suppression is specific to NRY rather than to the entire genome in males. Other mechanisms such as structural rearrangements and/or epigenetic modifications were likely involved, and comparative genome analysis and genetic mapping in multiple Silene species will help to shed light on the mechanism(s) of recombination suppression that led to the evolution of sex chromosomes.

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.

Does polyploidy inhibit sex chromosome evolution in angiosperms?

Dioecy is rare in flowering plants (5-6% of species), but is often controlled genetically by sex-linked regions (SLRs). It has so far been unclear whether, polyploidy affects sex chromosome evolution, as it does in animals, though polyploidy is quite common in angiosperms, including in dioecious species. Plants could be different, as, unlike many animal systems, degenerated sex chromosomes, are uncommon in plants. Here we consider sex determination in plants and plant-specific factors, and propose that constraints created at the origin of polyploids limit successful polyploidization of species with SLRs. We consider the most likely case of a polyploid of a dioecious diploid with an established SLR, and discuss the outcome in autopolyploids and allopolyploids. The most stable system possibly has an SLR on just one chromosome, with a strongly dominant genetic factor in the heterogametic sex (e.g., xxxY male in a tetraploid). If recombination occurs with its homolog, this will prevent Y chromosome degeneration. Polyploidy may also allow for reversibility of multiplied Z or X chromosomes into autosomes. Otherwise, low dosage of Y-linked SLRs compared to their multiple homologous x copies may cause loss of reliable sex-determination at higher ploidy levels. We discuss some questions that can be studied using genome sequencing, chromosome level-assemblies, gene expression studies and analysis of loci under selection.

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.

Sheltering of deleterious mutations explains the stepwise extension of recombination suppression on sex chromosomes and other supergenes

Many organisms have sex chromosomes with large nonrecombining regions that have expanded stepwise, generating "evolutionary strata" of differentiation. The reasons for this remain poorly understood, but the principal hypotheses proposed to date are based on antagonistic selection due to differences between sexes. However, it has proved difficult to obtain empirical evidence of a role for sexually antagonistic selection in extending recombination suppression, and antagonistic selection has been shown to be unlikely to account for the evolutionary strata observed on fungal mating-type chromosomes. We show here, by mathematical modeling and stochastic simulation, that recombination suppression on sex chromosomes and around supergenes can expand under a wide range of parameter values simply because it shelters recessive deleterious mutations, which are ubiquitous in genomes. Permanently heterozygous alleles, such as the male-determining allele in XY systems, protect linked chromosomal inversions against the expression of their recessive mutation load, leading to the successive accumulation of inversions around these alleles without antagonistic selection. Similar results were obtained with models assuming recombination-suppressing mechanisms other than chromosomal inversions and for supergenes other than sex chromosomes, including those without XY-like asymmetry, such as fungal mating-type chromosomes. However, inversions capturing a permanently heterozygous allele were found to be less likely to spread when the mutation load segregating in populations was lower (e.g., under large effective population sizes or low mutation rates). This may explain why sex chromosomes remain homomorphic in some organisms but are highly divergent in others. Here, we model a simple and testable hypothesis explaining the stepwise extensions of recombination suppression on sex chromosomes, mating-type chromosomes, and supergenes in general.

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'.

Identification of Sex Differentiation-Related microRNAs in Spinach Female and Male Flower

Sex determination and differentiation is an important biological process for unisexual flower development. Spinach is a model plant to study the mechanism of sex determination and differentiation of dioecious plant. Till now, little is known about spinach sex determination and differentiation mechanism. MicroRNAs are key factors in flower development. Herein, small RNA sequencing was performed to explore the roles of microRNAs in spinach sex determination and differentiation. As a result, 92 known and 3402 novel microRNAs were identified in 18 spinach female and male flower samples. 74 differentially expressed microRNAs were identified between female and male flowers, including 20 female-biased and 48 male-biased expression microRNAs. Target prediction identified 22 sex-biased microRNA-target pairs, which may be involved in spinach sex determination or differentiation. Among the differentially expressed microRNAs between FNS and M03, 55 microRNAs were found to reside in sex chromosome; one of them, sol-miR2550n, was functionally studied via genetic transformation. Silencing of sol-miR2550n resulted in abnormal anther while overexpression of sol-miR2550n induced early flowering, indicating sol-miR2550n was a male-promoting factor and validating the reliability of our small RNA sequencing data. Conclusively, this work can supply valuable information for exploring spinach sex determination and differentiation and provide a new insight in studying unisexual flower development.

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.

Sexual dimorphism in the dioecious willow Salix purpurea

PREMISE

The evolution of sex chromosomes is driven by sexual dimorphism, yet it can be challenging to document sexually dimorphic traits in dioecious plant species. At the genetic level, sexual dimorphism can be identified through sequence variation between females and males associated with sexually antagonistic traits and different fitness optima. This study aims to examine sexual dimorphism for 26 traits in three populations of Salix purpurea (a diversity panel and F₁ and F₂ populations) and determine the effect of the traits on biomass yield, a key trait in Salix bioenergy crops across multiple years, locations, and under manipulated growth conditions.

METHODS

Sexual dimorphism was evaluated for morphological, phenological, physiological, and wood composition traits in a diversity panel of unrelated S. purpurea accessions and in full-sib F₁ and F₂ families produced through controlled cross pollinations and grown in replicated field trials.

RESULTS

We observed sexual dimorphism in the timing of development for several traits that were highly predictive of biomass yield across three populations of S. purpurea. Across all populations and years surveyed, males had significantly shallower branching angle. Male plants highly predictive of biomass yield across three populations of S. purpurea also accumulated more nitrogen under fertilizer amendment as measured by SPAD in the diversity panel and had greater susceptibility to the rust fungus Melampsora americana in the F₂ family. Allometric modelling of biomass yield showed an effect of sex and of location on the interaction between yield and stem height.

CONCLUSIONS

These results provide evidence of sexual dimorphism for certain traits in S. purpurea that may be involved in sex chromosome 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.

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.

A General Model to Explain Repeated Turnovers of Sex Determination in the Salicaceae

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms by which sex systems evolve and their commonalities among plant species remain poorly understood. With both XY and ZW sex systems, the family Salicaceae provides a system to uncover the evolutionary forces driving sex chromosome turnovers. In this study, we performed a genome-wide association study to characterize sex determination in two Populus species, P. euphratica and P. alba. Our results reveal an XY system of sex determination on chromosome 14 of P. euphratica, and a ZW system on chromosome 19 of P. alba. We further assembled the corresponding sex-determination regions, and found that their sex chromosome turnovers may be driven by the repeated translocations of a Helitron-like transposon. During the translocation, this factor may have captured partial or intact sequences that are orthologous to a type-A cytokinin response regulator gene. Based on results from this and other recently published studies, we hypothesize that this gene may act as a master regulator of sex determination for the entire family. We propose a general model to explain how the XY and ZW sex systems in this family can be determined by the same RR gene. Our study provides new insights into the diversification of incipient sex chromosomes in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems.

Study of Sexual-Linked Genes (OGI and MeGI) on the Performance of Androecious Persimmons (Diospyros kaki Thunb.)

It is reported that the production of floral sexual phenotype in hexaploid monoecious persimmon (Diospyros kaki) is closely related to a pseudogene called OGI, and a short interspersed nuclear element (SINE)-like insertion (named Kali) in the OGI promoter leads to the gene silence. As a result, DNA methylation level of MeGI promoter determines the development of male or female flowers. However, the molecular mechanism in androecious D. kaki, which only bear male flowers, remains elusive. Here, real-time quantitative polymerase chain reaction (RT-qPCR), molecular cloning, and bisulfite PCR sequencing technique were carried out using 87 materials, including 56 androecious resources, 15 monoecious, and 16 gynoecious cultivars, to investigate the performance of OGI and MeGI on the specific androecious type of D. kaki in China. In conclusion, the Kali insertion was exactly located in the OGI promoter region, and the OGI gene and the Kali sequence were existing and conserved in androecious D. kaki. Meanwhile, we also demonstrated that the MeGI gene was widespread in our investigated samples. Ultimately, our result convincingly provided evidence that the low expression of OGI is probably ascribed to the presence of Kali displaying strong methylation in the OGI promoter, and low expression of MeGI, as well as high DNA methylation level, in the promoter was closely connected with the production of male flowers; this result was consistent with the monoecious persimmon model. Our findings provide predominant genetic aspects for investigation into androecious D. kaki, and future perfecting the sex-determining mechanisms in persimmon.

Comparative transcriptome analysis of male and female flowers in Spinacia oleracea L

Background

Dioecious spinach (Spinacia oleracea L.), a commercial and nutritional vegetable crop, serves as a model for studying the mechanisms of sex determination and differentiation in plants. However, this mechanism is still unclear. Herein, based on PacBio Iso-seq and Illumina RNA-seq data, comparative transcriptome analysis of male and female flowers were performed to explore the sex differentiation mechanism in spinach.

Results

Compared with published genome of spinach, 10,800 transcripts were newly annotated; alternative splicing, alternative polyadenylation and lncRNA were analyzed for the first time, increasing the diversity of spinach transcriptome. A total of 2965 differentially expressed genes were identified between female and male flowers at three early development stages. The differential expression of RNA splicing-related genes, polyadenylation-related genes and lncRNAs suggested the involvement of alternative splicing, alternative polyadenylation and lncRNA in sex differentiation. Moreover, 1946 male-biased genes and 961 female-biased genes were found and several candidate genes related to gender development were identified, providing new clues to reveal the mechanism of sex differentiation. In addition, weighted gene co-expression network analysis showed that auxin and gibberellin were the common crucial factors in regulating female or male flower development; however, the closely co-expressed genes of these two factors were different between male and female flower, which may result in spinach sex differentiation.

Conclusions

In this study, 10,800 transcripts were newly annotated, and the alternative splicing, alternative polyadenylation and long-noncoding RNA were comprehensively analyzed for the first time in spinach, providing valuable information for functional genome study. Moreover, candidate genes related to gender development were identified, shedding new insight on studying the mechanism of sex determination and differentiation in plant.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07277-4.

How to Become an Apomixis Model: The Multifaceted Case of Paspalum

In the past decades, the grasses of the Paspalum genus have emerged as a versatile model allowing evolutionary, genetic, molecular, and developmental studies on apomixis as well as successful breeding applications. The rise of such an archetypal system progressed through integrative phases, which were essential to draw conclusions based on solid standards. Here, we review the steps adopted in Paspalum to establish the current body of knowledge on apomixis and provide model breeding programs for other agronomically important apomictic crops. In particular, we discuss the need for previous detailed cytoembryological and cytogenetic germplasm characterization; the establishment of sexual and apomictic materials of identical ploidy level; the development of segregating populations useful for inheritance analysis, positional mapping, and epigenetic control studies; the development of omics data resources; the identification of key molecular pathways via comparative gene expression studies; the accurate molecular characterization of genomic loci governing apomixis; the in-depth functional analysis of selected candidate genes in apomictic and model species; the successful building of a sexual/apomictic combined breeding scheme.

Identification of the First Oomycete Mating-type Locus Sequence in the Grapevine Downy Mildew Pathogen, Plasmopara viticola

Mating types are self-incompatibility systems that promote outcrossing in plants, fungi, and oomycetes. Mating-type genes have been widely studied in plants and fungi but have yet to be identified in oomycetes, eukaryotic organisms closely related to brown algae that cause many destructive animal and plant diseases. We identified the mating-type locus of Plasmopara viticola, the oomycete responsible for grapevine downy mildew, one of the most damaging grapevine diseases worldwide. Using a genome-wide association approach, we identified a 570-kb repeat-rich non-recombining region controlling mating types, with two highly divergent alleles. We showed that one mating type was homozygous, whereas the other was heterozygous at this locus. The mating-type locus encompassed 40 genes, including one encoding a putative hormone receptor. Functional studies will, however, be required to validate the function of these genes and find the actual determinants of mating type. Our findings have fundamental implications for our understanding of the evolution of mating types, as they reveal a unique determinism involving an asymmetry of heterozygosity, as in sex chromosomes and unlike other mating-type systems. This identification of the mating-type locus in such an economically important crop pathogen also has applied implications, as outcrossing facilitates rapid evolution and resistance to harsh environmental conditions.

Sequencing and Analysis of the Sex Determination Region of Populus trichocarpa

The ages and sizes of a sex-determination region (SDR) are difficult to determine in non-model species. Due to the lack of recombination and enrichment of repetitive elements in SDRs, the quality of assembly with short sequencing reads is universally low. Unique features present in the SDRs help provide clues about how SDRs are established and how they evolve in the absence of recombination. Several Populus species have been reported with a male heterogametic configuration of sex (XX/XY system) mapped on chromosome 19, but the exact location of the SDR has been inconsistent among species, and thus far, none of these SDRs has been fully assembled in a genomic context. Here we identify the Y-SDR from a Y-linked contig directly from a long-read PacBio assembly of a Populus trichocarpa male individual. We also identified homologous gene sequences in the SDR of P. trichocarpa and the SDR of the W chromosome in Salix purpurea. We show that inverted repeats (IRs) found in the Y-SDR and the W-SDR are lineage-specific. We hypothesize that, although the two IRs are derived from the same orthologous gene within each species, they likely have independent evolutionary histories. Furthermore, the truncated inverted repeats in P. trichocarpa may code for small RNAs that target the homologous gene for RNA-directed DNA methylation. These findings support the hypothesis that diverse sex-determining systems may be achieved through similar evolutionary pathways, thereby providing a possible mechanism to explain the lability of sex-determination systems in plants in general.

Molecular Sex Identification in the Hardy Rubber Tree (Eucommia ulmoides Oliver) via ddRAD Markers

Eucommia ulmoides, also known as the industrially and medicinally important hardy rubber tree, is the sole species of Eucommiaceae. Nevertheless, its dioecious property hinders sex recognition by traditional morphological observation at very early developmental stages, thus inhibiting breeding and economic cropping. In this study, double-digest restriction site-associated DNA sequencing (ddRAD-seq) was applied to screen sex-linked molecular markers for sex identification and investigation of the sex determination system in 20 male and female E. ulmoides individual plants, respectively. In consequence, five candidate male-specific loci but no female-specific loci were predicated among the 183,752 male and 147,122 female catalogue loci by bioinformatics analysis. Subsequent PCR (polymerase chain reaction) amplification and Sanger sequencing examinations were performed on another 24 individuals, 12 for each sex, from a separate population. One ideal sex-linked locus, MSL4, was identified among the five putative male-specific loci that were found using ddRAD data. MSL4 is 479 bp in length and highly conserved in all the male individuals, suggesting its feature of being stable and repeatable. Our results also indicated that the sex of E. ulmoides is likely determined genetically. In short, this study provides a consistent and reproducible ddRAD marker (MSL4) that is able to discriminate male from female seedlings in E. ulmoides, which will be valuable for rapid breeding practice and better commercial production of this economically important tree.

The Location of the Pseudoautosomal Boundary in Silene latifolia

Y-chromosomes contain a non-recombining region (NRY), and in many organisms it was shown that the NRY expanded over time. How and why the NRY expands remains unclear. Young sex chromosomes, where NRY expansion occurred recently or is on-going, offer an opportunity to study the causes of this process. Here, we used the plant Silene latifolia, where sex chromosomes evolved ~11 million years ago, to study the location of the boundary between the NRY and the recombining pseudoautosomal region (PAR). The previous work devoted to the NRY/PAR boundary in S. latifolia was based on a handful of genes with locations approximately known from the genetic map. Here, we report the analysis of 86 pseudoautosomal and sex-linked genes adjacent to the S. latifolia NRY/PAR boundary to establish the location of the boundary more precisely. We take advantage of the dense genetic map and polymorphism data from wild populations to identify 20 partially sex-linked genes located in the “fuzzy boundary”, that rarely recombines in male meiosis. Genes proximal to this fuzzy boundary show no evidence of recombination in males, while the genes distal to this partially-sex-linked region are actively recombining in males. Our results provide a more accurate location for the PAR boundary in S. latifolia, which will help to elucidate the causes of PAR boundary shifts leading to NRY expansion over time.

An efficient RNA-seq-based segregation analysis identifies the sex chromosomes of Cannabis sativa

Cannabissativa–derived tetrahydrocannabinol (THC) production is increasing very fast worldwide. C. sativa is a dioecious plant with XY Chromosomes, and only females (XX) are useful for THC production. Identifying the sex chromosome sequence would improve early sexing and better management of this crop; however, the C. sativa genome projects have failed to do so. Moreover, as dioecy in the Cannabaceae family is ancestral, C. sativa sex chromosomes are potentially old and thus very interesting to study, as little is known about old plant sex chromosomes. Here, we RNA-sequenced a C. sativa family (two parents and 10 male and female offspring, 576 million reads) and performed a segregation analysis for all C. sativa genes using the probabilistic method SEX-DETector. We identified >500 sex-linked genes. Mapping of these sex-linked genes to a C. sativa genome assembly identified the largest chromosome pair being the sex chromosomes. We found that the X-specific region (not recombining between X and Y) is large compared to other plant systems. Further analysis of the sex-linked genes revealed that C. sativa has a strongly degenerated Y Chromosome and may represent the oldest plant sex chromosome system documented so far. Our study revealed that old plant sex chromosomes can have large, highly divergent nonrecombining regions, yet still be roughly homomorphic.

An efficient RNA-seq-based segregation analysis identifies the sex chromosomes of Cannabis sativa

Cannabis sativa-derived tetrahydrocannabinol (THC) production is increasing very fast worldwide. C. sativa is a dioecious plant with XY Chromosomes, and only females (XX) are useful for THC production. Identifying the sex chromosome sequence would improve early sexing and better management of this crop; however, the C. sativa genome projects have failed to do so. Moreover, as dioecy in the Cannabaceae family is ancestral, C. sativa sex chromosomes are potentially old and thus very interesting to study, as little is known about old plant sex chromosomes. Here, we RNA-sequenced a C. sativa family (two parents and 10 male and female offspring, 576 million reads) and performed a segregation analysis for all C. sativa genes using the probabilistic method SEX-DETector. We identified >500 sex-linked genes. Mapping of these sex-linked genes to a C. sativa genome assembly identified the largest chromosome pair being the sex chromosomes. We found that the X-specific region (not recombining between X and Y) is large compared to other plant systems. Further analysis of the sex-linked genes revealed that C. sativa has a strongly degenerated Y Chromosome and may represent the oldest plant sex chromosome system documented so far. Our study revealed that old plant sex chromosomes can have large, highly divergent nonrecombining regions, yet still be roughly homomorphic.

Sex is determined by XY chromosomes across the radiation of dioecious Nepenthes pitcher plants

Species with separate sexes (dioecy) are a minority among flowering plants, but dioecy has evolved multiple times independently in their history. The sex-determination system and sex-linked genomic regions are currently identified in a limited number of dioecious plants only. Here, we study the sex-determination system in a genus of dioecious plants that lack heteromorphic sex chromosomes and are not amenable to controlled breeding: Nepenthes pitcher plants. We genotyped wild populations of flowering males and females of three Nepenthes taxa using ddRAD-seq and sequenced a male inflorescence transcriptome. We developed a statistical tool (privacy rarefaction) to distinguish true sex specificity from stochastic noise in read coverage of sequencing data from wild populations and identified male-specific loci and XY-patterned single nucleotide polymorphsims (SNPs) in all three Nepenthes taxa, suggesting the presence of homomorphic XY sex chromosomes. The male-specific region of the Y chromosome showed little conservation among the three taxa, except for the essential pollen development gene DYT1 that was confirmed as male specific by PCR in additional Nepenthes taxa. Hence, dioecy and part of the male-specific region of the Nepenthes Y-chromosomes likely have a single evolutionary origin.

Somatic Mutation and Evolution in Plants

Somatic mutations are common in plants, and they may accumulate and be passed on to gametes. The determinants of somatic mutation accumulation include the intraorganismal selective effect of mutations, the number of cell divisions that separate the zygote from the formation of gametes, and shoot apical meristem structure and branching. Somatic mutations can promote the evolution of diploidy, polyploidy, sexual recombination, outcrossing, clonality, and separate sexes, and they may contribute genetic variability in many other traits. The amplification of beneficial mutations via intraorganismal selection may relax selection to reduce the genomic mutation rate or to protect the germline in plants. The total rate of somatic mutation, the distribution of selective effects and fates in the plant body, and the degree to which the germline is sheltered from somatic mutations are still poorly understood. Our knowledge can be improved through empirical estimates of mutation rates and effects on cell lineages and whole organisms, such as estimates of the reduction in fitness of progeny produced by within- versus between-flower crosses on the same plant, mutation coalescent studies within the canopy, and incorporation of somatic mutation into theoretical models of plant evolutionary genetics.

How to identify sex chromosomes and their turnover

Although sex is a fundamental component of eukaryotic reproduction, the genetic systems that control sex determination are highly variable. In many organisms the presence of sex chromosomes is associated with female or male development. Although certain groups possess stable and conserved sex chromosomes, others exhibit rapid sex chromosome evolution, including transitions between male and female heterogamety, and turnover in the chromosome pair recruited to determine sex. These turnover events have important consequences for multiple facets of evolution, as sex chromosomes are predicted to play a central role in adaptation, sexual dimorphism, and speciation. However, our understanding of the processes driving the formation and turnover of sex chromosome systems is limited, in part because we lack a complete understanding of interspecific variation in the mechanisms by which sex is determined. New bioinformatic methods are making it possible to identify and characterize sex chromosomes in a diverse array of non-model species, rapidly filling in the numerous gaps in our knowledge of sex chromosome systems across the tree of life. In turn, this growing data set is facilitating and fueling efforts to address many of the unanswered questions in sex chromosome evolution. Here, we synthesize the available bioinformatic approaches to produce a guide for characterizing sex chromosome system and identity simultaneously across clades of organisms. Furthermore, we survey our current understanding of the processes driving sex chromosome turnover, and highlight important avenues for future research.

YY males of the dioecious plant Mercurialis annua are fully viable but produce largely infertile pollen

The suppression of recombination during sex-chromosome evolution is thought to be favoured by linkage between the sex-determining locus and sexually antagonistic loci, and leads to the degeneration of the chromosome restricted to the heterogametic sex. Despite substantial evidence for genetic degeneration at the sequence level, the phenotypic effects of the earliest stages of sex-chromosome evolution are poorly known. Here, we compare the morphology, viability and fertility between XY and YY individuals produced by crossing seed-producing males in the dioecious plant Mercurialis annua, which has young sex chromosomes with limited X-Y sequence divergence. We found no significant difference in viability or vegetative morphology between XY and YY males. However, electron microscopy revealed clear differences in pollen anatomy, and YY males were significantly poorer sires in competition with their XY counterparts. Our study suggests either that the X chromosome is required for full male fertility in M. annua, or that male fertility is sensitive to the dosage of relevant Y-linked genes. We discuss the possibility that the maintenance of male-fertility genes on the X chromosome might have been favoured in recent population expansions that selected for the ability of females to produce pollen in the absence of males.

Insights into the genetic architecture of sexual dimorphism from an interspecific cross between two diverging Silene (Caryophyllaceae) species

The evolution of sexual dimorphism in species with separate sexes is influenced by the resolution of sexual conflicts creating sex differences through genetic linkage or sex-biased expression. Plants with different degrees of sexual dimorphism are thus ideal to study the genetic basis of sexual dimorphism. In this study we explore the genetic architecture of sexual dimorphism between Silene latifolia and Silene dioica. These species have chromosomal sex determination and differ in the extent of sexual dimorphism. To test whether QTL for sexually dimorphic traits have accumulated on the sex chromosomes and to quantify their contribution to species differences, we create a linkage map and performed QTL analysis of life history, flower and vegetative traits using an unidirectional interspecific F2 hybrid cross. We found support for an accumulation of QTL on the sex chromosomes and that sex differences explained a large proportion of the variance between species, suggesting that both natural and sexual selection contributed to species divergence. Sexually dimorphic traits that also differed between species displayed transgressive segregation. We observed a reversal in sexual dimorphism in the F2 population, where males tended to be larger than females, indicating that sexual dimorphism is constrained within populations but not in recombinant hybrids. This study contributes to the understanding of the genetic basis of sexual dimorphism and its evolution in Silene.

Two Y-chromosome-encoded genes determine sex in kiwifruit

Dioecy, the presence of male and female individuals, has evolved independently in multiple flowering plant lineages^1-3. Although theoretical models for the evolution of dioecy, such as the 'two-mutations' model, are well established^4,5, little is known about the specific genes determining sex and their evolutionary history³. Kiwifruit, a major tree crop consumed worldwide, is a dioecious species. In kiwifruit we previously identified a Y-encoded sex-determinant candidate gene acting as the suppressor of feminization (SuF), named Shy Girl (SyGI)⁶. Here, we identify a second Y-encoded sex-determinant that we named Friendly Boy (FrBy), which exhibits strong expression in tapetal cells. Gene-editing and complementation analyses in Arabidopsis thaliana and Nicotiana tabacum indicated that FrBy acts for the maintenance of male (M) functions, independently of SyGI, and that these functions are conserved across angiosperm species. We further characterized the genomic architecture of the small (<1 megabase pairs (Mb)) male-specific region of the Y chromosome (MSY), which harbours only two genes expressed extensively in developing gynoecia and androecia, respectively: SyGI and FrBy. Re-sequencing of the genome of a natural hermaphrodite kiwifruit revealed that this individual is genetically male but carries deletion(s) of parts of the Y chromosome, including SyGI. Additionally, expression of FrBy in female kiwifruit resulted in hermaphrodite plants. These results clearly indicate that Y-encoded SyGI and FrBy act independently as the SuF and M factors in kiwifruit, respectively, and provide insight into not only the evolutionary path leading to a two-factor sex-determination system, but also a new breeding approach for dioecious species.

Towards sex identification of Asian Palmyra palm (Borassus flabellifer L.) by DNA fingerprinting, suppression subtractive hybridization and de novo transcriptome sequencing

BACKGROUND

Asian Palmyra palm, the source of palm-sugar, is dioecious with a long juvenile period requiring at least 12 years to reach its maturity. To date, there is no reliable molecular marker for identifying sexes before the first bloom, limiting crop designs and utilization. We aimed to identify sex-linked markers for this palm using PCR-based DNA fingerprinting, suppression subtractive hybridization (SSH) and transcriptome sequencing.

METHODS

DNA fingerprints were generated between males and females based on RAPD, AFLP, SCoT, modified SCoT, ILP, and SSR techniques. Large-scale cloning and screening of SSH libraries and de novo transcriptome sequencing of male and female cDNA from inflorescences were performed to identify sex-specific genes for developing sex-linked markers.

RESULTS

Through extensive screening and re-testing of the DNA fingerprints (up to 1,204 primer pairs) and transcripts from SSH (>10,000 clones) and transcriptome data, however, no sex-linked marker was identified. Although de novo transcriptome sequencing of male and female inflorescences provided ∼32 million reads and 187,083 assembled transcripts, PCR analysis of selected sex-highly represented transcripts did not yield any sex-linked marker. This result may suggest the complexity and small sex-determining region of the Asian Palmyra palm. To this end, we provide the first global transcripts of male and female inflorescences of Asian Palmyra palm. Interestingly, sequence annotation revealed a large proportion of transcripts related to sucrose metabolism, which corresponds to the sucrose-rich sap produced in the inflorescences, and these transcripts will be useful for further understanding of sucrose production in sugar crop plants. Provided lists of sex-specific and differential-expressed transcripts would be beneficial to the further study of sexual development and sex-linked markers in palms and related species.

Early Sex-Chromosome Evolution in the Diploid Dioecious Plant Mercurialis annua

Suppressed recombination allows divergence between homologous sex chromosomes and the functionality of their genes. Here, we reveal patterns of the earliest stages of sex-chromosome evolution in the diploid dioecious herb Mercurialis annua on the basis of cytological analysis, de novo genome assembly and annotation, genetic mapping, exome resequencing of natural populations, and transcriptome analysis. The genome assembly contained 34,105 expressed genes, of which 10,076 were assigned to linkage groups. Genetic mapping and exome resequencing of individuals across the species range both identified the largest linkage group, LG1, as the sex chromosome. Although the sex chromosomes of M. annua are karyotypically homomorphic, we estimate that about one-third of the Y chromosome, containing 568 transcripts and spanning 22.3 cM in the corresponding female map, has ceased recombining. Nevertheless, we found limited evidence for Y-chromosome degeneration in terms of gene loss and pseudogenization, and most X- and Y-linked genes appear to have diverged in the period subsequent to speciation between M. annua and its sister species M. huetii, which shares the same sex-determining region. Taken together, our results suggest that the M. annua Y chromosome has at least two evolutionary strata: a small old stratum shared with M. huetii, and a more recent larger stratum that is probably unique to M. annua and that stopped recombining ∼1 MYA. Patterns of gene expression within the nonrecombining region are consistent with the idea that sexually antagonistic selection may have played a role in favoring suppressed recombination.

Origin and evolution of sex-determination systems in the brown algae

Sexual reproduction is a nearly universal feature of eukaryotic organisms. Meiosis appears to have had a single ancient origin, but the mechanisms underlying male or female sex determination are diverse and have emerged repeatedly and independently in the different eukaryotic groups. The brown algae are a group of multicellular photosynthetic eukaryotes that have a distinct evolutionary history compared with animals and plants, as they have been evolving independently for over 1 billion yr. Here, we review recent work using the brown alga Ectocarpus as a model organism to study haploid sex chromosomes, and highlight how the diversity of reproductive and life cycle features of the brown algae offer unique opportunities to characterize the evolutionary forces and the mechanisms underlying the evolution of sex determination.

A reference high-density genetic map of greater yam (Dioscorea alata L.)

This study generated the first high-density genetic map for D. alata based on genotyping-by-sequencing and provides new insight on sex determination in yam. Greater yam (Dioscorea alata L.) is a major staple food in tropical and subtropical areas. This study aimed to produce the first reference genetic map of this dioecious species using genotyping-by-sequencing. In this high-density map combining information of two F1 outcrossed populations, 20 linkage groups were resolved as expected and 1579 polymorphic markers were ordered. The consensus map length was 2613.5 cM with an average SNP interval of 1.68 cM. An XX/XY sex determination system was identified on LG6 via the study of sex ratio, homology of parental linkage groups and the identification of a major QTL for sex determination. Homology with the sequenced D. rotundata is described, and the median physical distance between SNPs was estimated at 139.1 kb. The effects of segregation distortion and the presence of heteromorphic sex chromosomes are discussed. This D. alata linkage map associated with the available genomic resources will facilitate quantitative trait mapping, marker-assisted selection and evolutionary studies in the important yet scarcely studied yam species.

Genetic and cytological analyses reveal the recombination landscape of a partially differentiated plant sex chromosome in kiwifruit

BACKGROUND

Angiosperm sex chromosomes, where present, are generally recently evolved. The key step in initiating the development of sex chromosomes from autosomes is the establishment of a sex-determining locus within a region of non-recombination. To better understand early sex chromosome evolution, it is important to determine the process by which recombination is suppressed around the sex determining genes. We have used the dioecious angiosperm kiwifruit Actinidia chinensis var. chinensis, which has an active-Y sex chromosome system, to study recombination rates around the sex locus, to better understand key events in the development of sex chromosomes.

RESULTS

We have confirmed the sex-determining region (SDR) in A. chinensis var. chinensis, using a combination of high density genetic mapping and fluorescent in situ hybridisation (FISH) of Bacterial Artificial Chromosomes (BACs) linked to the sex markers onto pachytene chromosomes. The SDR is a subtelomeric non-recombining region adjacent to the nucleolar organiser region (NOR). A region of restricted recombination of around 6 Mbp in size in both male and female maps spans the SDR and covers around a third of chromosome 25.

CONCLUSIONS

As recombination is suppressed over a similar region between X chromosomes and between and X and Y chromosomes, we propose that recombination is suppressed in this region because of the proximity of the NOR and the centromere, with both the NOR and centromere suppressing recombination, and this predates suppressed recombination due to differences between X and Y chromosomes. Such regions of suppressed recombination in the genome provide an opportunity for the evolution of sex chromosomes, if a sex-determining locus develops there or translocates into this region.

Sex Chromosome Turnover in Moths of the Diverse Superfamily Gelechioidea

Sex chromosomes play a central role in genetics of speciation and their turnover was suggested to promote divergence. In vertebrates, sex chromosome-autosome fusions resulting in neo-sex chromosomes occur frequently in male heterogametic taxa (XX/XY), but are rare in groups with female heterogamety (WZ/ZZ). We examined sex chromosomes of seven pests of the diverse lepidopteran superfamily Gelechioidea and confirmed the presence of neo-sex chromosomes in their karyotypes. Two synteny blocks, which correspond to autosomes 7 (LG7) and 27 (LG27) in the ancestral lepidopteran karyotype exemplified by the linkage map of Biston betularia (Geometridae), were identified as sex-linked in the tomato leafminer, Tuta absoluta (Gelechiidae). Testing for sex-linkage performed in other species revealed that while LG7 fused to sex chromosomes in a common ancestor of all Gelechioidea, the second fusion between the resulting neo-sex chromosome and the other autosome is confined to the tribe Gnoreschemini (Gelechiinae). Our data accentuate an emerging pattern of high incidence of neo-sex chromosomes in Lepidoptera, the largest clade with WZ/ZZ sex chromosome system, which suggest that the paucity of neo-sex chromosomes is not an intrinsic feature of female heterogamety. Furthermore, LG7 contains one of the major clusters of UDP-glucosyltransferases, which are involved in the detoxification of plant secondary metabolites. Sex chromosome evolution in Gelechioidea thus supports an earlier hypothesis postulating that lepidopteran sex chromosome-autosome fusions can be driven by selection for association of Z-linked preference or host-independent isolation genes with larval performance and thus can contribute to ecological specialization and speciation of moths.

Evolution of sex determination and heterogamety changes in section Otites of the genus Silene

Switches in heterogamety are known to occur in both animals and plants. Although plant sex determination systems probably often evolved more recently than those in several well-studied animals, including mammals, and have had less time for switches to occur, we previously detected a switch in heterogamety in the plant genus Silene: section Otites has both female and male heterogamety, whereas S. latifolia and its close relatives, in a different section of the genus, Melandrium (subgenus Behenantha), all have male heterogamety. Here we analyse the evolution of sex chromosomes in section Otites, which is estimated to have evolved only about 0.55 MYA. Our study confirms female heterogamety in S. otites and newly reveals female heterogamety in S. borysthenica. Sequence analyses and genetic mapping show that the sex-linked regions of these two species are the same, but the region in S. colpophylla, a close relative with male heterogamety, is different. The sex chromosome pairs of S. colpophylla and S. otites each correspond to an autosome of the other species, and both differ from the XY pair in S. latifolia. Silene section Otites species are suitable for detailed studies of the events involved in such changes, and our phylogenetic analysis suggests a possible change from female to male heterogamety within this section. Our analyses suggest a possibility that has so far not been considered, change in heterogamety through hybridization, in which a male-determining chromosome from one species is introgressed into another one, and over-rides its previous sex-determining system.

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

Background

Repetitive sequences, including transposable elements (TEs) and satellite DNAs, occupy a considerable portion of plant genomes. Analysis of the repeat fraction benefits the understanding of genome structure and evolution. Spinach (Spinacia oleracea L.), an important vegetable crop, is also a model dioecious plant species for studying sex determination and sex chromosome evolution. However, the repetitive sequences of the spinach genome have not been fully investigated.

Results

We extensively analyzed the repetitive components of draft spinach genome, especially TEs and satellites, by different strategies. A total of 16,002 full-length TEs were identified. Among the most abundant long terminal repeat (LTR) retrotransposons (REs), Copia elements were overrepresented compared with Gypsy ones. Angela was the most dominating Copia lineage; Ogre/Tat was the most abundant Gypsy lineage. The mean insertion age of LTR-REs was 1.42 million years; approximately 83.7% of these elements were retrotransposed during the last two million years. RepeatMasker totally masked about 64.05% of the spinach genome, with LTR-REs, non-LTR-REs, and DNA transposons occupying 49.2, 2.4, and 5.6%, respectively. Fluorescence in situ hybridization (FISH) analysis showed that most LTR-REs dispersed all over the chromosomes, by contrast, elements of CRM lineage were distributed at the centromeric region of all chromosomes. In addition, Ogre/Tat lineage mainly accumulated on sex chromosomes, and satellites Spsat2 and Spsat3 were exclusively located at the telomeric region of the short arm of sex chromosomes.

Conclusions

We reliably annotated the TE fraction of the draft genome of spinach. FISH analysis indicates that Ogre/Tat lineage and the sex chromosome-specific satellites DNAs might participate in sex chromosome formation and evolution. Based on FISH signals of microsatellites, together with 45S rDNA, a fine karyotype of spinach was established. This study improves our knowledge of repetitive sequence organization in spinach genome and aids in accurate spinach karyotype construction.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0147-6) contains supplementary material, which is available to authorized users.

Testing the translocation hypothesis and Haldane's rule in Rumex hastatulus

The translocation hypothesis regarding the origin of the XX/XY1Y2 sex chromosome system was tested with reference to the F1 hybrids between two chromosomal races of Rumex hastatulus. The hybrids derived from reciprocal crossing between the Texas (T) race and the North Carolina (NC) race were investigated for the first time with respect to the meiotic chromosome configuration in the pollen mother cells, pollen viability, and sex ratio. A sex chromosome trivalent in the NC × T males and two sex chromosome bivalents in the T × NC males were detected. The observed conjugation patterns confirmed the autosomal origin of the extra chromosome segments occurring in the North Carolina neo-sex chromosomes. Decreased pollen viability was found in the T × NC hybrid in contrast to the NC × T hybrid and the parental forms. Moreover, only in the T × NC hybrid sex ratio was significantly female-biased (1:1.72). Thus, Haldane's rule for both male fertility and male rarity was shown in this hybrid. According to the authors' knowledge, R. hastatulus is just the second plant with sex chromosomes in which Haldane's rule was evidenced.

Sex and the flower - developmental aspects of sex chromosome evolution

Background

The evolution of dioecious plants is occasionally accompanied by the establishment of sex chromosomes: both XY and ZW systems have been found in plants. Structural studies of sex chromosomes are now being followed up by functional studies that are gradually shedding light on the specific genetic and epigenetic processes that shape the development of separate sexes in plants.

Scope

This review describes sex determination diversity in plants and the genetic background of dioecy, summarizes recent progress in the investigation of both classical and emerging model dioecious plants and discusses novel findings. The advantages of interspecies hybrids in studies focused on sex determination and the role of epigenetic processes in sexual development are also overviewed.

Conclusions

We integrate the genic, genomic and epigenetic levels of sex determination and stress the impact of sex chromosome evolution on structural and functional aspects of plant sexual development. We also discuss the impact of dioecy and sex chromosomes on genome structure and expression.