Single locus sex determination and female heterogamety in the basket willow (Salix viminalis L.)

A haplotype-resolved genome assembly and gene expression map of Cushion willow

Salix brachista, commonly known as Cushion willow, is a common component of subnival alpine assemblages and a dioecious or monoecious plant with a creeping stem and numerous lateral branches. Cushion willow takes cuttings more easier and has a specific sex system, making it a suitable system for studying the evolution of plant sex determination, adaptive evolution of alpine plants, and mining stress resistance gene resource that cope with the hostile alpine environment. Therefore, Cushion willow has potential value in genetic improvements for willows used as bioenergy crops, in gardening, and as ornamental plants. However, the genome of Cushion willow still contains some un-assembled repetitive sequences, and there is limited availability of a gene expression atlas, which hinders its potential use for the aforementioned purposes. Here, we updated the genome of Cushion willow to be haplotype-resolved and near telomere-to-telomere, and obtained a high-quality transcriptomic map. Our research provides a potential model species for alpine adaptive research, sex determination evolution studies, and improving willow crops.

Gender specific SNP markers in Coscinium fenestratum (Gaertn.) Colebr. for resource augmentation

BACKGROUND

Unregulated extraction of highly traded medicinal plant species results in drastic decline of the natural resources and alters viable sex ratio of populations. Conservation and long-term survival of such species, require gender specific restoration programs to ensure reproductive success. However, it is often difficult to differentiate sex of individuals before reaching reproductive maturity. C. fenestratum is one of the medicinally important and overexploited dioecious woody liana, with a reproductive maturity of 15 years. Currently, no information is available to identify sex of C. fenestratum in seedling stage while augmenting the resources. Thus, the current study envisages to utilize transcriptomics approach for gender differentiation which is imperative for undertaking viable resource augmentation programmes.

METHODS AND RESULTS

Gender specific SNPs with probable role in sexual reproduction/sex determination was located using comparative transcriptomics approach (sampling male and female individuals), alongside gene ontology and annotation. Nine sets of primers were synthesized from 7 transcripts (involved in sexual reproduction/other biological process) containing multiple SNP variants. Out of the nine primer pairs, only one SNP locus with no available information of its role in reproduction, showed consistent and accurate results (males-heterozygous and females-homozygous), in the analyzed 40 matured individuals of known sexes. Thus validated the efficiency of this SNP marker in differentiating male and female individuals.

CONCLUSIONS

The study could identify SNPs linked to the loci with apparent role in gender differentiation. This SNP marker can be used for early sexing of seedlings for in-situ conservation and resource augmentation of C. fenestratum in Kerala, India.

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.

Functional analysis of Salix purpurea genes support roles for ARR17 and GATA15 as master regulators of sex determination

The Salicaceae family is of growing interest in the study of dioecy in plants because the sex determination region (SDR) has been shown to be highly dynamic, with differing locations and heterogametic systems between species. Without the ability to transform and regenerate Salix in tissue culture, previous studies investigating the mechanisms regulating sex in the genus Salix have been limited to genome resequencing and differential gene expression, which are mostly descriptive in nature, and functional validation of candidate sex determination genes has not yet been conducted. Here, we used Arabidopsis to functionally characterize a suite of previously identified candidate genes involved in sex determination and sex dimorphism in the bioenergy shrub willow Salix purpurea. Six candidate master regulator genes for sex determination were heterologously expressed in Arabidopsis, followed by floral proteome analysis. In addition, 11 transcription factors with predicted roles in mediating sex dimorphism downstream of the SDR were tested using DAP-Seq in both male and female S. purpurea DNA. The results of this study provide further evidence to support models for the roles of ARR17 and GATA15 as master regulator genes of sex determination in S. purpurea, contributing to a regulatory system that is notably different from that of its sister genus Populus. Evidence was also obtained for the roles of two transcription factors, an AP2/ERF family gene and a homeodomain-like transcription factor, in downstream regulation of sex dimorphism.

Salix gordejevii females exhibit more resistance against wind erosion than males under aeolian environment

Effects of wind erosion on growth and adaptability have been widely reported in many plants, but little attention has been paid to dioecious plants. Recent studies have shown that sex-specific responses to environmental changes in many plants exist. To explore sexual differences in response to wind erosion, female and male Salix gordejevii saplings growing on inter-dune land (no erosion) and on the windward slope of the dune (20cm wind erosion) in Hunshandake Sandy Land were chosen and their morphology, biomass and physiological traits were investigated, respectively. Wind erosion significantly reduced plant growth, biomass accumulation, gas exchange and chlorophyll fluorescence, and obviously disrupted osmotic regulation function and antioxidant enzyme system in both sexes, especially in males. Under wind erosion condition, females exhibited higher sapling height (SH), basal diameter (BD), leaf dry mass (LDM), root dry mass (RDM), total dry mass (TDM), root percentage in total dry mass, net photosynthesis rate (P _n), maximum efficiency of photosystem II (F _v/F _m), effective quantum yield of PSII (Φ _PSII), relative water content (RWC) of leaves, superoxide dismutase (SOD) and peroxidase (POD) activities, but lower malondialdehyde (MDA), proline as well as soluble sugar content than did males. However, no significant sexual differences in most of these traits were observed under no erosion condition. Our results demonstrated that females possess a greater resistance to wind erosion than do males, with females having a better photosynthetic capacity, stronger water retention capacity and more efficient antioxidant system to alleviate negative effects caused by aeolian environment.

Mapping the sex determination region in the Salix F1 hybrid common parent population confirms a ZW system in six diverse species

Within the genus Salix, there are approximately 350 species native primarily to the northern hemisphere and adapted to a wide range of habitats. This diversity can be exploited to mine novel alleles conferring variation important for production as a bioenergy crop, but also to identify evolutionarily important genes, such as those involved in sex determination. To leverage this diversity, we created a mapping population by crossing 6 Salix species (Salix viminalis, Salix suchowensis, Salix integra, Salix koriyanagi, Salix udensis, and Salix alberti) to common male and female Salix purpurea parents. Each family was genotyped via genotyping-by-sequencing and assessed for kinship and population structure as well as the construction of 16 backcross linkage maps to be used as a genetic resource for breeding and selection. Analyses of population structure resolved both the parents and F₁ progeny to their respective phylogenetic section and indicated that the S. alberti parent was misidentified and was most likely S.suchowensis. Sex determining regions were identified on Salix chromosome 15 in the female-informative maps for seven of the eight families indicating that these species share a common female heterogametic ZW sex system. The eighth family, S. integra × S. purpurea, was entirely female and had a truncated chromosome 15. Beyond sex determination, the Salix F₁ hybrid common parent population (Salix F₁ HCP) introduced here will be useful in characterizing genetic factors underlying complex traits, aid in marker-assisted selection, and support genome assemblies for this promising bioenergy crop.

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

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.

Genome-wide association mapping uncovers sex-associated copy number variation markers and female hemizygous regions on the W chromosome in Salix viminalis

Background

Sex chromosomes are in some species largely undifferentiated (homomorphic) with restricted sex determination regions. Homomorphic but different sex chromosomes are found in the closely related genera Populus and Salix indicating flexible sex determination systems, ideal for studies of processes involved in sex chromosome evolution. We have performed genome-wide association studies of sex and analysed sex chromosomes in a population of 265 wild collected Salix viminalis accessions and studied the sex determining locus.

Results

A total of 19,592 markers were used in association analyses using both Fisher’s exact tests and a single-marker mixed linear model, which resulted in 48 and 41 sex-associated (SA) markers respectively. Across all 48 SA markers, females were much more often heterozygous than males, which is expected if females were the heterogametic sex. The majority of the SA markers were, based on positions in the S. purpurea genome, located on chromosome 15, previously demonstrated to be the sex chromosome. Interestingly, when mapping the genotyping-by-sequencing sequence tag harbouring the two SA markers with the highest significance to the S. viminalis genomic scaffolds, five regions of very high similarity were found: three on a scaffold that represents a part of chromosome 15, one on a scaffold that represents a part of chromosome 9 and one on a scaffold not anchored to the genome. Based on segregation differences of the alleles at the two marker positions and on differences in PCR amplification between females and males we conclude that females had multiple copies of this DNA fragment (chromosome 9 and 15), whereas males only had one (chromosome 9). We therefore postulate that the female specific sequences have been copied from chromosome 9 and inserted on chromosome 15, subsequently developing into a hemizygous W chromosome linked region.

Conclusions

Our results support that sex determination in S. viminalis is controlled by one locus on chromosome 15. The segregation patterns observed at the SA markers furthermore confirm that S. viminalis females are the heterogametic sex. We also identified a translocation from chromosome 9 to the W chromosome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08021-2.

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.

Integrative genomics reveals paths to sex dimorphism in Salix purpurea L

Sex dimorphism and gene expression were studied in developing catkins in 159 F2 individuals from the bioenergy crop Salix purpurea, and potential mechanisms and pathways for regulating sex development were explored. Differential expression, eQTL, bisulfite sequencing, and network analysis were used to characterize sex dimorphism, detect candidate master regulator genes, and identify pathways through which the sex determination region (SDR) may mediate sex dimorphism. Eleven genes are presented as candidates for master regulators of sex, supported by gene expression and network analyses. These include genes putatively involved in hormone signaling, epigenetic modification, and regulation of transcription. eQTL analysis revealed a suite of transcription factors and genes involved in secondary metabolism and floral development that were predicted to be under direct control of the sex determination region. Furthermore, data from bisulfite sequencing and small RNA sequencing revealed strong differences in expression between males and females that would implicate both of these processes in sex dimorphism pathways. These data indicate that the mechanism of sex determination in Salix purpurea is likely different from that observed in the related genus Populus. This further demonstrates the dynamic nature of SDRs in plants, which involves a multitude of mechanisms of sex determination and a high rate of turnover.

Sex determination through X-Y heterogamety in Salix nigra

The development of non-recombining sex chromosomes has radical effects on the evolution of discrete sexes and sexual dimorphism. Although dioecy is rare in plants, sex chromosomes have evolved repeatedly throughout the diversification of angiosperms, and many of these sex chromosomes are relatively young compared to those found in vertebrates. In this study, we designed and used a sequence capture array to identify a novel sex-linked region (SLR) in Salix nigra, a basal species in the willow clade, and demonstrated that this species has XY heterogamety. We did not detect any genetic overlap with the previously characterized ZW SLRs in willows, which map to a different chromosome. The S. nigra SLR is characterized by strong recombination suppression across a 2 MB region and an excess of low-frequency alleles, resulting in a low Tajima's D compared to the remainder of the genome. We speculate that either a recent bottleneck in population size or factors related to positive or background selection generated this differential pattern of Tajima's D on the X and autosomes. This discovery provides insights into factors that may influence the evolution of sex chromosomes in plants and contributes to a large number of recent observations that underscore their dynamic nature.

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.

The Diversity and Dynamics of Sex Determination in Dioecious Plants

The diversity of inflorescences among flowering plants is captivating. Such charm is not only due to the variety of sizes, shapes, colors, and flowers displayed, but also to the range of reproductive systems. For instance, hermaphrodites occur abundantly throughout the plant kingdom with both stamens and carpels within the same flower. Nevertheless, 10% of flowering plants have separate unisexual flowers, either in different locations of the same individual (monoecy) or on different individuals (dioecy). Despite their rarity, dioecious plants provide an excellent opportunity to investigate the mechanisms involved in sex expression and the evolution of sex-determining regions (SDRs) and sex chromosomes. The SDRs and the evolution of dioecy have been studied in many species ranging from Ginkgo to important fruit crops. Some of these studies, for example in asparagus or kiwifruit, identified two sex-determining genes within the non-recombining SDR and may thus be consistent with the classical model for the evolution of dioecy from hermaphroditism via gynodioecy, that predicts two successive mutations, the first one affecting male and the second one female function, becoming linked in a region of suppressed recombination. On the other hand, aided by genome sequencing and gene editing, single factor sex determination has emerged in other species, such as persimmon or poplar. Despite the diversity of sex-determining mechanisms, a tentative comparative analysis of the known sex-determining genes and candidates in different species suggests that similar genes and pathways may be employed repeatedly for the evolution of dioecy. The cytokinin signaling pathway appears important for sex determination in several species regardless of the underlying genetic system. Additionally, tapetum-related genes often seem to act as male-promoting factors when sex is determined via two genes. We present a unified model that synthesizes the genetic networks of sex determination in monoecious and dioecious plants and will support the generation of hypothesis regarding candidate sex determinants in future studies.

Male Linked Genomic Region Determines Sex in Dioecious Amaranthus palmeri

Dioecy, the separation of reproductive organs on different individuals, has evolved repeatedly in different plant families. Several evolutionary paths to dioecy have been suggested, but the mechanisms behind sex determination is not well understood. The diploid dioecious Amaranthus palmeri represents a well-suited model system to study sex determination in plants. Despite the agricultural importance of the species, the genetic control and evolutionary state of dioecy in A. palmeri is currently unknown. Early cytogenetic experiments did not identify heteromorphic chromosomes. Here, we used whole-genome sequencing of male and female pools from 2 independent populations to elucidate the genetic control of dioecy in A. palmeri. Read alignment to a close monoecious relative and allele frequency comparisons between male and female pools did not reveal significant sex-linked genes. Consequently, we employed an alignment-free k-mer comparison which enabled us to identify a large number of male-specific k-mers. We assembled male-specific contigs comprising a total of almost 2 Mb sequence, proposing a XY sex-determination system in the species. We were able to identify the potential Y chromosome in the A. palmeri draft genome sequence as 90% of our male-specific sequence aligned to a single scaffold. Based on our findings, we suggest an intermediate evolutionary state of dioecy with a young Y chromosome in A. palmeri. Our findings give insight into the evolution of sex chromosomes in plants and may help to develop sustainable strategies for weed management.

Evidences for a role of two Y-specific genes in sex determination in Populus deltoides

Almost all plants in the genus Populus are dioecious (i.e. trees are either male or female), but it is unknown whether dioecy evolved in a common ancestor or independently in different subgenera. Here, we sequence the small peritelomeric X- and Y-linked regions of P. deltoides chromosome XIX. Two genes are present only in the Y-linked region. One is a duplication of a non-Y-linked, female-specifically expressed response regulator, which produces siRNAs that block this gene's expression, repressing femaleness. The other is an LTR/Gypsy transposable element family member, which generates long non-coding RNAs. Overexpression of this gene in A. thaliana promotes androecium development. We also find both genes in the sex-determining region of P. simonii, a different poplar subgenus, which suggests that they are both stable components of poplar sex-determining systems. By contrast, only the duplicated response regulator gene is present in the sex-linked regions of P. davidiana and P. tremula. Therefore, findings in our study suggest dioecy may have evolved independently in different poplar subgenera.

Transcriptome analysis of contrasting resistance to herbivory by Empoasca fabae in two shrub willow species and their hybrid progeny

Short rotation woody biomass cultivars developed from fast-growing shrub species of willow (Salix spp.) have superior properties as perennial energy crops for the Northeast and Midwest US. However, the insect pest potato leafhopper (PLH) Empoasca fabae (Harris) can cause serious damage and reduce yield of susceptible genotypes. Currently, the willow cultivars in use display varying levels of susceptibility under PLH infestation. However, genes and markers for resistance to PLH are not yet available for marker-assisted selection in breeding. In this study, transcriptome differences between a resistant genotype 94006 (S. purpurea) and a susceptible cultivar ‘Jorr’ (S. viminalis), and their hybrid progeny were determined. Over 600 million RNA-Seq reads were generated and mapped to the Salix purpurea reference transcriptome. Gene expression analyses revealed the unique defense mechanism in resistant genotype 94006 that involves PLH-induced secondary cell wall modification. In the susceptible genotypes, genes involved in programed cell death were highly expressed, explaining the necrosis symptoms after PLH feeding. Overall, the discovery of resistance genes and defense mechanisms provides new resources for shrub willow breeding and research in the future.

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.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Fine mapping of the sex locus in Salix triandra confirms a consistent sex determination mechanism in genus Salix

Salix triandra belongs to section Amygdalinae in genus Salix, which is in a different section from the willow species in which sex determination has been well studied. Studying sex determination in distantly related willow species will help to clarify whether the sexes of different willows arise through a common sex determination system. For this purpose, we generated an intraspecific full-sib F1 population for S. triandra and constructed high-density genetic linkage maps for the crossing parents using restriction site-associated DNA sequencing and following a two-way pseudo-testcross strategy. With the established maps, the sex locus was positioned in linkage group XV only in the maternal map, and no sex linkage was detected in the paternal map. Consistent with previous findings in other willow species, our study showed that chromosome XV was the incipient sex chromosome and that females were the heterogametic sex in S. triandra. Therefore, sex in this willow species is also determined through a ZW sex determination system. We further performed fine mapping in the vicinity of the sex locus with SSR markers. By comparing the physical and genetic distances for the target interval encompassing the sex determination gene confined by SSRs, severe recombination repression was revealed in the sex determination region in the female map. The recombination rate in the confined interval encompassing the sex locus was approximately eight-fold lower than the genome-wide average. This study provides critical information relevant to sex determination in S. triandra.

Sex-Related Differences in Growth, Herbivory, and Defense of Two Salix Species

Sex-related differences in sex ratio, growth, and herbivory are widely documented in many dioecious plants. The common pattern is for males to grow faster than females and to be less well-defended against herbivores, but Salix is an exception. To study sex-related differences in the patterns of resource allocation for growth and defense in willows, we conducted a large-scale field experiment to investigate the flowering sex ratio, mortality, growth traits, insect herbivory and content of defensive substances in three Salix populations comprising two species. Results demonstrate that the two Salix suchowensis Cheng populations have a female bias in the sex ratio, whereas no bias is found in the S. triandra L. population. Male individuals in the S. suchowensis populations have significantly higher mortality rates than females. However, the mortality rate of S. triandra population has no gender difference. This finding may be one of the explanations for the difference in sex ratio between the two species. The females are larger in height, ground diameter, and biomass, and have a higher nutritional quality (N concentration) than males in both species. Nevertheless, slow-growing males have a higher concentration of the defense chemical (total phenol) and lower degrees of insect herbivory than females. Additionally, biomass is positively correlated with herbivory and negatively correlated with defense in the two willow species. It is concluded that the degrees of herbivory would have a great influence on resource allocation for growth and defense. Meanwhile, it also provides important implications for understanding the evolution of dioecy.

Construction of an anchoring SSR marker genetic linkage map and detection of a sex-linked region in two dioecious populations of red bayberry

Red bayberry (Morella rubra) is an evergreen fruit tree found in southern China whose whole-genome sequence has recently been published. We updated the linkage map of the species by adding 118 SSR markers and the female-specific marker MrFT2_BD-SEX. The integrated map included eight linkage groups and spanned 491 cM. Eleven sex-associated markers were identified, six of which were located in linkage group 8, in agreement with the previously reported location of the sex-determining region. The MrFT2_BD-SEX marker was genotyped in 203 cultivated accessions. Among the females of the accessions, we found two female-specific alleles, designated W-b (151 bp) and W-d (129 bp). We previously found that 'Dongkui', a female cultivar, could produce viable pollen (we refer to such plants 'Dongkui-male') and serve as the paternal parent in crosses. The genotypes of the MrFT2_BD-SEX marker were W-b/Z in 'Biqi' and W-d/Z in 'Dongkui-male'. The progeny of a cross between these parents produced a 3:1 female (W-) to male (ZZ) ratio and the expected 1:1:1:1 ratio of W-b/W-d: W-b/Z: W-d/Z: Z/Z. In addition, the flowering and fruiting phenotypes of all the F1 progeny fit their genotypes. Our results confirm the existence of ZW sex determination and show that the female phenotype is controlled by a single dominant locus (W) in a small genomic region (59 kb and less than 3.3 cM). Furthermore, we have produced a homozygous "super female" (WW) that should produce all-female offspring in the F2 generation, providing a foundation for commercial use and presenting great potential for use in modern breeding programs.

A willow sex chromosome reveals convergent evolution of complex palindromic repeats

Background

Sex chromosomes have arisen independently in a wide variety of species, yet they share common characteristics, including the presence of suppressed recombination surrounding sex determination loci. Mammalian sex chromosomes contain multiple palindromic repeats across the non-recombining region that show sequence conservation through gene conversion and contain genes that are crucial for sexual reproduction. In plants, it is not clear if palindromic repeats play a role in maintaining sequence conservation in the absence of homologous recombination.

Results

Here we present the first evidence of large palindromic structures in a plant sex chromosome, based on a highly contiguous assembly of the W chromosome of the dioecious shrub Salix purpurea. The W chromosome has an expanded number of genes due to transpositions from autosomes. It also contains two consecutive palindromes that span a region of 200 kb, with conspicuous 20-kb stretches of highly conserved sequences among the four arms that show evidence of gene conversion. Four genes in the palindrome are homologous to genes in the sex determination regions of the closely related genus Populus, which is located on a different chromosome. These genes show distinct, floral-biased expression patterns compared to paralogous copies on autosomes.

Conclusion

The presence of palindromes in sex chromosomes of mammals and plants highlights the intrinsic importance of these features in adaptive evolution in the absence of recombination. Convergent evolution is driving both the independent establishment of sex chromosomes as well as their fine-scale sequence structure.

Joint linkage and association mapping of complex traits in shrub willow (Salix purpurea L.)

BACKGROUND AND AIMS

Increasing energy demands and the necessity to reduce greenhouse gas emissions are key motivating factors driving the development of lignocellulosic crops as an alternative to non-renewable energy sources. The effects of global climate change will require a better understanding of the genetic basis of complex adaptive traits to breed more resilient bioenergy feedstocks, like willow (Salix spp.). Shrub willow is a sustainable and dedicated bioenergy crop, bred to be fast-growing and high-yielding on marginal land without competing with food crops. In a rapidly changing climate, genomic advances will be vital for the sustained improvement of willow and other non-model bioenergy crops. Here, joint genetic mapping was used to exploit genetic variation garnered from both recent and historical recombination events in S. purpurea.

METHODS

A panel of North American naturalized S. purpurea accessions and full-sib F2S. purpurea population were genotyped and phenotyped for a suite of morphological, physiological, pest and disease resistance, and wood chemical composition traits, collected from multi-environment and multi-year replicated field trials. Controlling for population stratification and kinship in the association panel and spatial variation in the F2, a comprehensive mixed model analysis was used to dissect the complex genetic architecture and plasticity of these important traits.

KEY RESULTS

Individually, genome-wide association (GWAS) models differed in terms of power, but the combined approach, which corrects for yearly and environmental co-factors across datasets, improved the overall detection and resolution of associated loci. Although there were few significant GWAS hits located within support intervals of QTL for corresponding traits in the F2, many large-effect QTL were identified, as well as QTL hotspots.

CONCLUSIONS

This study provides the first comparison of linkage analysis and linkage disequilibrium mapping approaches in Salix, and highlights the complementarity and limits of these two methods for elucidating the genetic architecture of complex bioenergy-related traits of a woody perennial breeding programme.

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.

Sex-specific polymorphism of MET1 and ARR17 genes in Populus × sibirica

The genus Populus is an effective model in tree genetics. This genus includes dioecious species and, recently, whole genome resequencing of P. trichocarpa and P. balsamifera enabled the identification of sex-linked regions and sex-associated single nucleotide polymorphisms (SNPs). These results created new opportunities to study sex determination in poplars. In the present work, we performed deep sequencing of genes encoding METHYLTRANSFERASE1 (MET1) and homolog of ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17), which are localized in a sex-linked region of Populus genome and contain a number of sex-associated SNPs. Amplicon libraries for 38 samples of P. × sibirica (19 males and 19 females) were sequenced on MiSeq Illumina (300 nt paired-end reads) and approximately 4000× coverage was obtained for each sample. In total, from 80 to 179 SNPs were detected in poplar individuals for MET1, and from 16 to 49 SNPs were detected for ARR17. We identified 17 sex-specific SNPs (11 in MET1 and 6 in ARR17) - they were present in all males but absent in all females. For identified sex-specific SNP sites, females were homozygous, while males were heterozygous. Moreover, colocation of sex-specific SNPs confirming the XY sex-determination system of poplars was revealed: in one allelic variant, males had the same nucleotides as females, while in the other, sex-specific SNPs were present. Based on the data obtained, we developed and successfully applied a high-resolution melting-based approach for sex identification in poplars. The developed molecular markers are useful for distinguishing between male and female poplars in scientific research and can also be applied to select male-only genotypes for use in city landscaping and production of paper, pulp, and biofuel.

The red bayberry genome and genetic basis of sex determination

Morella rubra, red bayberry, is an economically important fruit tree in south China. Here, we assembled the first high-quality genome for both a female and a male individual of red bayberry. The genome size was 313-Mb, and 90% sequences were assembled into eight pseudo chromosome molecules, with 32 493 predicted genes. By whole-genome comparison between the female and male and association analysis with sequences of bulked and individual DNA samples from female and male, a 59-Kb region determining female was identified and located on distal end of pseudochromosome 8, which contains abundant transposable element and seven putative genes, four of them are related to sex floral development. This 59-Kb female-specific region was likely to be derived from duplication and rearrangement of paralogous genes and retained non-recombinant in the female-specific region. Sex-specific molecular markers developed from candidate genes co-segregated with sex in a genetically diverse female and male germplasm. We propose sex determination follow the ZW model of female heterogamety. The genome sequence of red bayberry provides a valuable resource for plant sex chromosome evolution and also provides important insights for molecular biology, genetics and modern breeding in Myricaceae family.

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.

Uneven selection pressure accelerating divergence of Populus and Salix

Populus (poplars) and Salix (willows) are sister genera in the Salicaceae family that arise from a common tetraploid ancestor. The karyotypes of these two lineages are distinguished by two major interchromosomal and some minor intrachromosomal rearrangements, but which one is evolutionarily more primitive remains debatable. In this study, we compare the selection pressure acting on the paralogous genes resulting from salicoid duplication (PGRS) within and between the genomes of the two lineages. Purifying selection was determined to act more strongly on the PGRS in willow than on those in poplar, which would cause a faster loss of paralogous duplicates in willow. Therefore, Salix species are supposed to evolve faster than Populus species, which is consistent with the observation that the former are taxonomically and morphologically more diverse than the latter. In these two lineages, different autosomes were found to have been evolving into sex chromosomes. Examining the ω ratio and the PGRS in the sex determination regions in willow and poplar revealed higher convergent selection pressure and a faster loss of PGRS in the sex determination regions of both lineages. At the chromosome level, the sex chromosome in poplar is characterized by the lowest gene density among all chromosome members, while this feature is not observed on the sex chromosome in willow, suggesting that Populus species may inherit the more incipient sex chromosome from their progenitor. Taken together, Salix is supposed to be the nascent lineage arising from the additional round of genome reorganization that distinguishes the karyotypes of the two sister genera. In this study, assessment of ω ratios also detected a list of paralogous genes under unusual selection pressure, which could have special consequences for the adaptive evolution of Salicaceae species. In conclusion, the results of this study provide unique information for better understanding the genetic mechanism accelerating the divergence of these two closely related lineages.

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.

Hardwood Tree Genomics: Unlocking Woody Plant Biology

Woody perennial angiosperms (i.e., hardwood trees) are polyphyletic in origin and occur in most angiosperm orders. Despite their independent origins, hardwoods have shared physiological, anatomical, and life history traits distinct from their herbaceous relatives. New high-throughput DNA sequencing platforms have provided access to numerous woody plant genomes beyond the early reference genomes of Populus and Eucalyptus, references that now include willow and oak, with pecan and chestnut soon to follow. Genomic studies within these diverse and undomesticated species have successfully linked genes to ecological, physiological, and developmental traits directly. Moreover, comparative genomic approaches are providing insights into speciation events while large-scale DNA resequencing of native collections is identifying population-level genetic diversity responsible for variation in key woody plant biology across and within species. Current research is focused on developing genomic prediction models for breeding, defining speciation and local adaptation, detecting and characterizing somatic mutations, revealing the mechanisms of gender determination and flowering, and application of systems biology approaches to model complex regulatory networks underlying quantitative traits. Emerging technologies such as single-molecule, long-read sequencing is being employed as additional woody plant species, and genotypes within species, are sequenced, thus enabling a comparative ("evo-devo") approach to understanding the unique biology of large woody plants. Resource availability, current genomic and genetic applications, new discoveries and predicted future developments are illustrated and discussed for poplar, eucalyptus, willow, oak, chestnut, and pecan.

RAD-seq reveals genetic structure of the F2-generation of natural willow hybrids (Salix L.) and a great potential for interspecific introgression

BACKGROUND

Hybridization of species with porous genomes can eventually lead to introgression via repeated backcrossing. The potential for introgression between species is reflected by the extent of segregation distortion in later generation hybrids. Here we studied a population of hybrids between Salix purpurea and S. helvetica that has emerged within the last 30 years on a glacier forefield in the European Alps due to secondary contact of the parental species. We used 5758 biallelic SNPs produced by RAD sequencing with the aim to ascertain the predominance of backcrosses (F1 hybrid x parent) or F2 hybrids (F1 hybrid x F1 hybrid) among hybrid offspring. Further, the SNPs were used to study segregation distortion in the second hybrid generation.

RESULTS

The analyses in STRUCTURE and NewHybrids revealed that the population consisted of parents and F1 hybrids, whereas hybrid offspring consisted mainly of backcrosses to either parental species, but also some F2 hybrids. Although there was a clear genetic differentiation between S. purpurea and S. helvetica (FST = 0.24), there was no significant segregation distortion in the backcrosses or the F2 hybrids. Plant height of the backcrosses resembled the respective parental species, whereas F2 hybrids were more similar to the subalpine S. helvetica.

CONCLUSIONS

The co-occurrence of the parental species and the hybrids on the glacier forefield, the high frequency of backcrossing, and the low resistance to gene flow via backcrossing make a scenario of introgression in this young hybrid population highly likely, potentially leading to the transfer of adaptive traits. We further suggest that this willow hybrid population may serve as a model for the evolutionary processes initiated by recent global warming.

Repeated translocation of a gene cassette drives sex-chromosome turnover in strawberries

Turnovers of sex-determining systems represent important diversifying forces across eukaryotes. Shifts in sex chromosomes—but conservation of the master sex-determining genes—characterize distantly related animal lineages. Yet in plants, in which separate sexes have evolved repeatedly and sex chromosomes are typically homomorphic, we do not know whether such translocations drive sex-chromosome turnovers within closely related taxonomic groups. This phenomenon can only be demonstrated by identifying sex-associated nucleotide sequences, still largely unknown in plants. The wild North American octoploid strawberries (Fragaria) exhibit separate sexes (dioecy) with homomorphic, female heterogametic (ZW) inheritance, yet sex maps to three different chromosomes in different taxa. To characterize these turnovers, we identified sequences unique to females and assembled their reads into contigs. For most octoploid Fragaria taxa, a short (13 kb) sequence was observed in all females and never in males, implicating it as the sex-determining region (SDR). This female-specific “SDR cassette” contains both a gene with a known role in fruit and pollen production and a novel retrogene absent on Z and autosomal chromosomes. Phylogenetic comparison of SDR cassettes revealed three clades and a history of repeated translocation. Remarkably, the translocations can be ordered temporally due to the capture of adjacent sequence with each successive move. The accumulation of the “souvenir” sequence—and the resultant expansion of the hemizygous SDR over time—could have been adaptive by locking genes into linkage with sex. Terminal inverted repeats at the insertion borders suggest a means of movement. To our knowledge, this is the first plant SDR shown to be translocated, and it suggests a new mechanism (“move-lock-grow”) for expansion and diversification of incipient sex chromosomes.

Sex chromosomes frequently restructure themselves during organismal evolution, often becoming highly differentiated. This dynamic process is poorly understood for most taxa, especially during the early stages typical of many dioecious flowering plants. We show that in wild strawberries, a female-specific region of DNA is associated with sex and has repeatedly changed its genomic location, each time increasing the size of the hemizygous female-specific sequence on the W sex chromosome. This observation shows, for the first time to our knowledge, that plant sex regions can “jump” and suggests that this phenomenon may be adaptive by gathering and locking new genes into linkage with sex. This conserved and presumed causal sex-determining sequence, which varies in both genomic location and degree of differentiation, will facilitate future studies to understand how sex chromosomes first begin to differentiate.

Characterization of a large sex determination region in Salix purpurea L. (Salicaceae)

Dioecy has evolved numerous times in plants, but heteromorphic sex chromosomes are apparently rare. Sex determination has been studied in multiple Salix and Populus (Salicaceae) species, and P. trichocarpa has an XY sex determination system on chromosome 19, while S. suchowensis and S. viminalis have a ZW system on chromosome 15. Here we use whole genome sequencing coupled with quantitative trait locus mapping and a genome-wide association study to characterize the genomic composition of the non-recombining portion of the sex determination region. We demonstrate that Salix purpurea also has a ZW system on chromosome 15. The sex determination region has reduced recombination, high structural polymorphism, an abundance of transposable elements, and contains genes that are involved in sex expression in other plants. We also show that chromosome 19 contains sex-associated markers in this S. purpurea assembly, along with other autosomes. This raises the intriguing possibility of a translocation of the sex determination region within the Salicaceae lineage, suggesting a common evolutionary origin of the Populus and Salix sex determination loci.

Haploid selection, sex ratio bias, and transitions between sex-determining systems

Sex determination is remarkably dynamic; many taxa display shifts in the location of sex-determining loci or the evolution of entirely new sex-determining systems. Predominant theories for why we observe such transitions generally conclude that novel sex-determining systems are favoured by selection if they equalise the sex ratio or increase linkage with a locus that experiences different selection in males versus females. We use population genetic models to extend these theories in two ways: (1) We consider the dynamics of loci very tightly linked to the ancestral sex-determining loci, e.g., within the nonrecombining region of the ancestral sex chromosomes. Variation at such loci can favour the spread of new sex-determining systems in which the heterogametic sex changes (XY to ZW or ZW to XY) and the new sex-determining region is less closely linked (or even unlinked) to the locus under selection. (2) We consider selection upon haploid genotypes either during gametic competition (e.g., pollen competition) or meiosis (i.e., nonmendelian segregation), which can cause the zygotic sex ratio to become biased. Haploid selection can drive transitions between sex-determining systems without requiring selection to act differently in diploid males versus females. With haploid selection, we find that transitions between male and female heterogamety can evolve so that linkage with the sex-determining locus is either strengthened or weakened. Furthermore, we find that sex ratio biases may increase or decrease with the spread of new sex chromosomes, which implies that transitions between sex-determining systems cannot be simply predicted by selection to equalise the sex ratio. In fact, under many conditions, we find that transitions in sex determination are favoured equally strongly in cases in which the sex ratio bias increases or decreases. Overall, our models predict that sex determination systems should be highly dynamic, particularly when haploid selection is present, consistent with the evolutionary lability of this trait in many taxa.

Size and Content of the Sex-Determining Region of the Y Chromosome in Dioecious Mercurialis annua, a Plant with Homomorphic Sex Chromosomes

Dioecious plants vary in whether their sex chromosomes are heteromorphic or homomorphic, but even homomorphic sex chromosomes may show divergence between homologues in the non-recombining, sex-determining region (SDR). Very little is known about the SDR of these species, which might represent particularly early stages of sex-chromosome evolution. Here, we assess the size and content of the SDR of the diploid dioecious herb Mercurialis annua, a species with homomorphic sex chromosomes and mild Y-chromosome degeneration. We used RNA sequencing (RNAseq) to identify new Y-linked markers for M. annua. Twelve of 24 transcripts showing male-specific expression in a previous experiment could be amplified by polymerase chain reaction (PCR) only from males, and are thus likely to be Y-linked. Analysis of genome-capture data from multiple populations of M. annua pointed to an additional six male-limited (and thus Y-linked) sequences. We used these markers to identify and sequence 17 sex-linked bacterial artificial chromosomes (BACs), which form 11 groups of non-overlapping sequences, covering a total sequence length of about 1.5 Mb. Content analysis of this region suggests that it is enriched for repeats, has low gene density, and contains few candidate sex-determining genes. The BACs map to a subset of the sex-linked region of the genetic map, which we estimate to be at least 14.5 Mb. This is substantially larger than estimates for other dioecious plants with homomorphic sex chromosomes, both in absolute terms and relative to their genome sizes. Our data provide a rare, high-resolution view of the homomorphic Y chromosome of a dioecious plant.

Recent Sex Chromosome Divergence despite Ancient Dioecy in the Willow Salix viminalis

Sex chromosomes can evolve when recombination is halted between a pair of chromosomes, and this can lead to degeneration of the sex-limited chromosome. In the early stages of differentiation sex chromosomes are homomorphic, and even though homomorphic sex chromosomes are very common throughout animals and plants, we know little about the evolutionary forces shaping these types of sex chromosomes. We used DNA- and RNA-Seq data from females and males to explore the sex chromosomes in the female heterogametic willow, Salix viminalis, a species with ancient dioecy but with homomorphic sex chromosomes. We detected no major sex differences in read coverage in the sex determination (SD) region, indicating that the W region has not significantly degenerated. However, single nucleotide polymorphism densities in the SD region are higher in females compared with males, indicating very recent recombination suppression, followed by the accumulation of sex-specific single nucleotide polymorphisms. Interestingly, we identified two female-specific scaffolds that likely represent W-chromosome-specific sequence. We show that genes located in the SD region display a mild excess of male-biased expression in sex-specific tissue, and we use allele-specific gene expression analysis to show that this is the result of masculinization of expression on the Z chromosome rather than degeneration of female-expression on the W chromosome. Together, our results demonstrate that insertion of small DNA fragments and accumulation of sex-biased gene expression can occur before the detectable decay of the sex-limited chromosome.

The Evolution of Sex Chromosomes and Dosage Compensation in Plants

Plant sex chromosomes can be vastly different from those of the few historical animal model organisms from which most of our understanding of sex chromosome evolution is derived. Recently, we have seen several advancements from studies on green algae, brown algae, and land plants that are providing a broader understanding of the variable ways in which sex chromosomes can evolve in distant eukaryotic groups. Plant sex-determining genes are being identified and, as expected, are completely different from those in animals. Species with varying levels of differentiation between the X and Y have been found in plants, and these are hypothesized to be representing different stages of sex chromosome evolution. However, we are also finding that sex chromosomes can remain morphologically unchanged over extended periods of time. Where degeneration of the Y occurs, it appears to proceed similarly in plants and animals. Dosage compensation (a phenomenon that compensates for the consequent loss of expression from the Y) has now been documented in a plant system, its mechanism, however, remains unknown. Research has also begun on the role of sex chromosomes in sexual conflict resolution, and it appears that sex-biased genes evolve similarly in plants and animals, although the functions of these genes remain poorly studied. Because the difficulty in obtaining sex chromosome sequences is increasingly being overcome by methodological developments, there is great potential for further discovery within the field of plant sex chromosome evolution.

Dominance and Sexual Dimorphism Pervade the Salix purpurea L. Transcriptome

The heritability of gene expression is critical in understanding heterosis and is dependent on allele-specific regulation by local and remote factors in the genome. We used RNA-Seq to test whether variation in gene expression among F1 and F2 intraspecific Salix purpurea progeny is attributable to cis- and trans-regulatory divergence. We assessed the mode of inheritance based on gene expression levels and allele-specific expression for F1 and F2 intraspecific progeny in two distinct tissue types: shoot tip and stem internode. In addition, we explored sexually dimorphic patterns of inheritance and regulatory divergence among F1 progeny individuals. We show that in S. purpurea intraspecific crosses, gene expression inheritance largely exhibits a maternal dominant pattern, regardless of tissue type or pedigree. A significantly greater number of cis- and trans-regulated genes coincided with upregulation of the maternal parent allele in the progeny, irrespective of the magnitude, whereas the paternal allele was higher expressed for genes showing cis × trans or compensatory regulation. Importantly, consistent with previous genetic mapping results for sex in shrub willow, we have delimited sex-biased gene expression to a 2 Mb pericentromeric region on S. purpurea chr15 and further refined the sex determination region. Altogether, our results offer insight into the inheritance of gene expression in S. purpurea as well as evidence of sexually dimorphic expression which may have contributed to the evolution of dioecy in Salix.

Detecting the Candidate Gender Determinants by Bioinformatic Prediction of miRNAs and Their Targets from Transcriptome Sequences of the Male and Female Flowers in Salix suchowensis

MicroRNAs (miRNAs) belong to a class of small, noncoding, and endogenous single-stranded RNAs that negatively regulate gene expression at the posttranscriptional level. Potential miRNAs can be identified based on sequence homology since miRNAs are highly conserved in plants. In this study, we aligned the expressed sequence tags derived from flower buds of male and female S. suchowensis to miRNAs in the miRBase, which enable us to identify 34 potential miRNAs from flower buds of the alternate sexes. Among them, 11 were from the female and 23 were from the male. Analyzing sequence complementarity led to identification of 124 and 55 miRNA targets in the male and female flower buds, respectively. By mapping the target genes of the predicted miRNAs to the sequence assemblies of S. suchowensis, a miR156 mediated gene was detected at the gender locus of willow, which was a transcription factor involved in flower development. It is noteworthy that this target is not expressed in male flower, while it is expressed fairly highly in female flower based on the transcriptome data derived from the alternate sexes of willows. This study provides new bioinformatic clue for further exploring the genetic mechanism underlying gender determination in willows.

Allelic incompatibility can explain female biased sex ratios in dioecious plants

Background

Biased sex ratios are common among dioecious plant species despite the theoretical prediction of selective advantage of even sex ratios. Albeit the high prevalence of deviations from even sex ratios, the genetic causes to sex biases are rarely known outside of a few model species. Here we present a mechanism underlying the female biased sex ratio in the dioecious willow species Salix viminalis.

Results

We compared the segregation pattern of genome-wide single nucleotide polymorphism markers in two contrasting bi-parental pedigree populations, the S3 with even sex ratio and the S5 with a female biased sex ratio. With the segregation analysis and comparison between the two populations, we were able to demonstrate that sex determination and sex ratio distortion are controlled by different genetic mechanisms. We furthermore located the sex ratio distorter locus to a Z/W-gametologous region on chromosome 15, which was in close linkage with the sex determination locus. Interestingly, all males in the population with biased sex ratio have in this sex ratio distorter locus the same genotype, meaning that males with the Z₁/Z₃-genotype were missing from the population, thereby creating the 2:1 female biased sex ratio.

Conclusions

We attribute the absence of Z₁/Z₃ males to an allelic incompatibility between maternally and paternally inherited alleles in this sex ratio distorter locus. Due to the tight linkage with the sex determination locus only male individuals are purged from the population at an early age, presumably before or during seed development. We showed that such allelic incompatibility could be stably maintained over evolutionary times through a system of overdominant or pseudooverdominant alleles. Thus, it is possible that the same mechanism generates the female biased sex ratio in natural willow populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3634-5) contains supplementary material, which is available to authorized users.

Identification of Quantitative Trait Loci Conditioning the Main Biomass Yield Components and Resistance to Melampsora spp. in Salix viminalis × Salix schwerinii Hybrids

The biomass of Salix viminalis is the most highly valued source of green energy, followed by S. schwerinii, S. dasyclados and other species. Significant variability in productivity and leaf rust resistance are noted both within and among willow species, which creates new opportunities for improving willow yield parameters through selection of desirable recombinants supported with molecular markers. The aim of this study was to identify quantitative trait loci (QTLs) linked with biomass yield-related traits and the resistance/susceptibility of Salix mapping population to leaf rust. The experimental material comprised a mapping population developed based on S. viminalis × S. schwerinii hybrids. Phenotyping was performed on plants grown in a field experiment that had a balanced incomplete block design with 10 replications. Based on a genetic map, 11 QTLs were identified for plant height, 9 for shoot diameter, 3 for number of shoots and 11 for resistance/susceptibility to leaf rust. The QTLs identified in our study explained 3%–16% of variability in the analyzed traits. Our findings make significant contributions to the development of willow breeding programs and research into shrubby willow crops grown for energy.

A Century of Sex Determination in Flowering Plants

Plants have evolved a diverse array of strategies for sexual reproduction, particularly through the modification of male and female organs at distinct points in development. The immense variation in sexual systems across the land plants provides a unique opportunity to study the genetic, epigenetic, phylogenetic, and ecological underpinnings of sex determination. Here, we reflect on more than a century of research into flowering plant sex determination, placing a particular focus on the foundational genetic and cytogenetic observations, experiments, and hypotheses. Building on the seminal work on the genetics of plant sex, modern comparative genomic analyses now allow us to address longstanding questions about sex determination and the origins of sex chromosomes.

Homomorphic ZW chromosomes in a wild strawberry show distinctive recombination heterogeneity but a small sex-determining region

Recombination in ancient, heteromorphic sex chromosomes is typically suppressed at the sex-determining region (SDR) and proportionally elevated in the pseudoautosomal region (PAR). However, little is known about recombination dynamics of young, homomorphic plant sex chromosomes. We examine male and female function in crosses and unrelated samples of the dioecious octoploid strawberry Fragaria chiloensis in order to map the small and recently evolved SDR controlling both traits and to examine recombination patterns on the incipient ZW chromosome. The SDR of this ZW system is located within a 280 kb window, in which the maternal recombination rate is lower than the paternal one. In contrast to the SDR, the maternal PAR recombination rate is much higher than the rates of the paternal PAR or autosomes, culminating in an elevated chromosome-wide rate. W-specific divergence is elevated within the SDR and a single polymorphism is observed in high species-wide linkage disequilibrium with sex. Selection for recombination suppression within the small SDR may be weak, but fluctuating sex ratios could favor elevated recombination in the PAR to remove deleterious mutations on the W. The recombination dynamics of this nascent sex chromosome with a modestly diverged SDR may be typical of other dioecious plants.

Major Chromosomal Rearrangements Distinguish Willow and Poplar After the Ancestral "Salicoid" Genome Duplication

Populus (poplar) and Salix (willow) are sister genera in the Salicaceae family. In both lineages extant species are predominantly diploid. Genome analysis previously revealed that the two lineages originated from a common tetraploid ancestor. In this study, we conducted a syntenic comparison of the corresponding 19 chromosome members of the poplar and willow genomes. Our observations revealed that almost every chromosomal segment had a parallel paralogous segment elsewhere in the genomes, and the two lineages shared a similar syntenic pinwheel pattern for most of the chromosomes, which indicated that the two lineages diverged after the genome reorganization in the common progenitor. The pinwheel patterns showed distinct differences for two chromosome pairs in each lineage. Further analysis detected two major interchromosomal rearrangements that distinguished the karyotypes of willow and poplar. Chromosome I of willow was a conjunction of poplar chromosome XVI and the lower portion of poplar chromosome I, whereas willow chromosome XVI corresponded to the upper portion of poplar chromosome I. Scientists have suggested that Populus is evolutionarily more primitive than Salix. Therefore, we propose that, after the “salicoid” duplication event, fission and fusion of the ancestral chromosomes first give rise to the diploid progenitor of extant Populus species. During the evolutionary process, fission and fusion of poplar chromosomes I and XVI subsequently give rise to the progenitor of extant Salix species. This study contributes to an improved understanding of genome divergence after ancient genome duplication in closely related lineages of higher plants.

Plant Sex Chromosomes

Although individuals in most flowering plant species, and in many haploid plants, have both sex functions, dioecious species-in which individuals have either male or female functions only-are scattered across many taxonomic groups, and many species have genetic sex determination. Among these, some have visibly heteromorphic sex chromosomes, and molecular genetic studies are starting to uncover sex-linked markers in others, showing that they too have fully sex-linked regions that are either too small or are located in chromosomes that are too small to be cytologically detectable from lack of pairing, lack of visible crossovers, or accumulation of heterochromatin. Detailed study is revealing that, like animal sex chromosomes, plant sex-linked regions show evidence for accumulation of repetitive sequences and genetic degeneration. Estimating when recombination stopped confirms the view that many plants have young sex-linked regions, making plants of great interest for studying the timescale of these changes.

Confirmation of Single-Locus Sex Determination and Female Heterogamety in Willow Based on Linkage Analysis

In this study, we constructed high-density genetic maps of Salix suchowensis and mapped the gender locus with an F₁ pedigree. Genetic maps were separately constructed for the maternal and paternal parents by using amplified fragment length polymorphism (AFLP) markers and the pseudo-testcross strategy. The maternal map consisted of 20 linkage groups that spanned a genetic distance of 2333.3 cM; whereas the paternal map contained 21 linkage groups that covered 2260 cM. Based on the established genetic maps, it was found that the gender of willow was determined by a single locus on linkage group LG_03, and the female was the heterogametic gender. Aligned with mapped SSR markers, linkage group LG_03 was found to be associated with chromosome XV in willow. It is noteworthy that marker density in the vicinity of the gender locus was significantly higher than that expected by chance alone, which indicates severe recombination suppression around the gender locus. In conclusion, this study confirmed the findings on the single-locus sex determination and female heterogamety in willow. It also provided additional evidence that validated the previous studies, which found that different autosomes evolved into sex chromosomes between the sister genera of Salix (willow) and Populus (poplar).

What causes female bias in the secondary sex ratios of the dioecious woody shrub Salix sitchensis colonizing a primary successional landscape?

•

PREMISE OF THE STUDY

Females often outnumber males in Salix populations, although the mechanisms behind female bias are not well understood and could be caused by both genetic and ecological factors. We investigated several ecological factors that could bias secondary sex ratios of Salix sitchensis colonizing Mount St. Helens after the 1980 eruption.•

METHODS

We determined whether S. sitchensis secondary sex ratios varied across disturbance zones created by the eruption and across mesic and hydric habitats within each zone. For one population, we tracked adult mortality, whole-plant reproductive allocation, the number of stems, and plant size for 2 years. In a field experiment, we created artificial streams to test whether vegetative reproduction via stem fragments was sex-biased.•

KEY RESULTS

We found a consistent 2:1 female bias in S. sitchensis secondary sex ratios across all disturbance zones and habitats. Despite female plants sometimes allocating more resources (in terms of carbon, nitrogen, and phosphorus) to reproduction than males, we found no evidence of sex-biased mortality. The establishment rate of S. sitchensis experimental stems did not differ between the sexes, indicating that vegetative reproduction was not distorting secondary sex ratios.•

CONCLUSIONS

We hypothesize that S. sitchensis secondary sex ratios depend on either early-acting genetic factors affecting the seed sex ratio or sex-specific germination or survival rates before maturity, as opposed to factors associated with reproduction in adult plants.