Sexual dimorphism in the dioecious willow Salix purpurea

Floral Morphometry and Sexual System Determination in Pink Pepper (Schinus terebinthifolia - Anacardiaceae)

The increasing global importance of pink peppertree (Schinus terebinthifolia, Anacardiaceae) as a high-value commercial crop and its potential for expansion in production demand appropriate management due to uncertainties regarding its sexual system. This study focused on evaluating the morphology of sterile and fertile floral whorls, as well as analyzing the sexual system of pink pepper in two populations in northeastern Brazil. The results revealed no significant differences in the morphological characteristics of the flowers between the studied areas, suggesting that the species possesses notable adaptability to environmental conditions. However, a significant difference in the proportion of staminate individuals was observed in both areas, representing over 88% and 72%, respectively. A correlation was observed between the size of the stamens and the presence of apparently atrophied pistils (r=0.275; df=178; p<0.001), along with the occurrence of fruits in these hermaphroditic plants. In this context, the species should be considered gynodioecious due to the presence of plants with hermaphroditic flowers and plants with pistillate flowers. However, further research is essential to elucidate the role of pollinators, especially bees and wasps, and to better understand the fruiting process in hermaphroditic flowers. These insights have the potential to significantly enhance management aiming for efficient fruit production, promoting its economic and ecological relevance.

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.

Structural variation of a sex-linked region confers monoecy and implicates GATA15 as a master regulator of sex in Salix purpurea

The Salicaceae, including Populus and Salix, are dioecious perennials that utilize different sex determination systems. This family provides a useful system to better understand the evolution of dioecy and sex chromosomes. Here, a rare monoecious genotype of Salix purpurea, 94003, was self- and cross-pollinated and progeny sex ratios were used to test hypotheses on possible mechanisms of sex determination. To delimit genomic regions associated with monoecious expression, the 94003 genome sequence was assembled and DNA- and RNA-Seq of progeny inflorescences was performed. Based on alignments of progeny shotgun DNA sequences to the haplotype-resolved monoecious 94003 genome assembly and reference male and female genomes, a 1.15 Mb sex-linked region on Chr15W was confirmed to be absent in monecious plants. Inheritance of this structural variation is responsible for the loss of a male-suppressing function in what would otherwise be genetic females (ZW), resulting in monoecy (ZW^H or WW^H ), or lethality, if homozygous (W^H W^H ). We present a refined, two-gene sex determination model for Salix purpurea, mediated by ARR17 and GATA15 that is different from the single-gene ARR17-mediated system in the related genus Populus.

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.

Genetic mapping of sexually dimorphic volatile and non-volatile floral secondary chemistry of a dioecious willow

Secondary chemistry often differs between sexes in dioecious plant species, a pattern attributed to its possible role in the evolution and/or maintenance of dioecy. We used GC-MS to measure floral volatiles emitted from, and LC-MS to quantitate non-volatile secondary compounds contained in, female and male Salix purpurea willow catkins from an F2 family. Using the abundance of these chemicals, we then performed quantitative trait locus (QTL) mapping to locate them on the genome, identified biosynthetic candidate genes in the QTL intervals, and examined expression patterns of candidate genes using RNA-seq. Male flowers emitted more total terpenoids than females, but females produced more benzenoids. Male tissue contained greater amounts of phenolic glycosides, but females had more chalcones and flavonoids. A flavonoid pigment and a spermidine derivative were found only in males. Male catkins were almost twice the mass of females. Forty-two QTL were mapped for 25 chemical traits and catkin mass across 16 of the 19 S. purpurea chromosomes. Several candidate genes were identified, including a chalcone isomerase associated with seven compounds. A better understanding of the genetic basis of the sexually dimorphic chemistry of a dioecious species may shed light on how chemically mediated ecological interactions may have helped in the evolution and maintenance of dioecy.

Sex-based metabolic and microbiota differences in roots and rhizosphere soils of dioecious papaya (Carica papaya L.)

Dioecious plant species have a high genetic variation that is important for coping with or adapting to environmental stress through natural selection. Intensive studies have reported dimorphism morphism in morphology, physiology, as well as biotic and abiotic stress responses in dioecious plants. Here, we demonstrated the dimorphism of metabolic profile and the preference of some microorganisms in the roots and rhizosphere soils of male and female papaya. The metabolic composition of roots were significantly different between the males and females. Some sex hormones occurred in the differential metabolites in roots and rhizosphere soils. For example, testosterone was up-regulated in male papaya roots and rhizosphere soils, whereas norgestrel was up-regulated in the female papaya roots, indicating a possible balance in papaya roots to control the sexual differentiation. Plant hormones such as BRs, JAs, SA and GAs were also detected among the differential metabolites in the roots and rhizosphere soils of dioecious papaya. In addition, some metabolites that have medicinal values, such as ecliptasaponin A, crocin, berberine and sapindoside A were also expressed differentially between the two sexes. Numerous differential metabolites from the papaya roots were secreted in the soil, resulting in the differences in microbial community structure in the roots and rhizosphere soils. Some nitrogen-fixing bacteria such as Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium, Brevundimonas and Microvirga were enriched in the male papaya roots or rhizosphere soils. While Candidatus Solibacter and Tumebacillus, which utilize organic matters, were enriched in the roots or rhizosphere soils of the female papaya. Some differences in the fungi abundance were also observed in both male and female papaya roots. These findings uncovered the effect of sex types on the metabolic and microbiota differences in roots and rhizosphere soils in papaya and will lead to investigations of underlining genomic and molecular mechanisms.

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

Population Variability of Almond-Leaved Willow (Salix triandra L.) Based on the Leaf Morphometry: Isolation by Distance and Environment Explain Phenotypic Diversity

Almond-leaved willow (Salix triandra L., Salicaceae) is a dioecious shrub, rarely a small tree that grows under various environmental conditions. We examined the population structure of 12 populations of almond-leaved willow using nine leaf morphological traits and specific leaf area. Populations were selected from a range of habitats, from continental to the sub-Mediterranean zone, to examine the influence of environmental conditions (climate and altitude) and geographic distance on leaf variability. Significant differences were confirmed among all populations for all traits, with significant correlations between geographic location of populations and morphological traits, and between environmental conditions and morphological traits. Large-leaved populations were found in continental and sub-Mediterranean climates, while small-leaved populations were found in higher elevations and smaller karstic rivers. In addition, populations from floodplains showed greater variability than populations from the karstic habitats, indicating a positive influence of lowland habitats and possible underlying differences in gene pool size. In conclusion, we found that environmental conditions and geographical distances in addition to genetic drift, are the main influences on the variability in almond-leaved willow, with the species showing a high level of plasticity and adaptation to local environmental conditions.

Nonadditive gene expression is correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow (Salix spp.)

Many studies have highlighted the complex and diverse basis for heterosis in inbred crops. Despite the lack of a consensus model, it is vital that we turn our attention to understanding heterosis in undomesticated, heterozygous, and polyploid species, such as willow (Salix spp.). Shrub willow is a dedicated energy crop bred to be fast-growing and high yielding on marginal land without competing with food crops. A trend in willow breeding is the consistent pattern of heterosis in triploids produced from crosses between diploid and tetraploid species. Here, we test whether differentially expressed genes are associated with heterosis in triploid families derived from diploid Salix purpurea, diploid Salix viminalis, and tetraploid Salix miyabeana parents. Three biological replicates of shoot tips from all family progeny and parents were collected after 12 weeks in the greenhouse and RNA extracted for RNA-Seq analysis. This study provides evidence that nonadditive patterns of gene expression are correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow. Expression-level dominance was most correlated with heterosis for biomass yield traits and was highly enriched for processes involved in starch and sucrose metabolism. In addition, there was a global dosage effect of parent alleles in triploid hybrids, with expression proportional to copy number variation. Importantly, differentially expressed genes between family parents were most predictive of heterosis for both field and greenhouse collected traits. Altogether, these data will be used to progress models of heterosis to complement the growing genomic resources available for the improvement of heterozygous perennial bioenergy crops.