Integrative genomics reveals paths to sex dimorphism in Salix purpurea L

Abstract

Sex dimorphism and gene expression were studied in developing catkins in 159 F₂ 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.

References

[1] Harkess, A. et al. Sex determination by two Y-Linked genes in garden asparagus. Plant Cell 32, 1790-1796 (2020).

[2] Akagi, T., Henry, I. M., Tao, R. & Comai, L. A Y-chromosome-encoded small RNA acts as a sex determinant in persimmons. Science 346, 646-650 (2014).

[3] Akagi, T. et al. Two Y-chromosome-encoded genes determine sex in kiwifruit. Nat. Plants 5, 801-809 (2019).

[4] Muller, N. A. et al. A single gene underlies the dynamic evolution of poplar sex determination. Nat. Plants 6, 630-637 (2020).

[5] Yang, H. W., Akagi, T., Kawakatsu, T. & Tao, R. Gene networks orchestrated by MeGI: a single-factor mechanism underlying sex determination in persimmon. Plant J. 98, 97-111 (2019).

[6] Charlesworth, D., Charlesworth, B. & Marais, G. Steps in the evolution of heteromorphic sex chromosomes. Heredity 95, 118-128 (2005).

[7] Akagi, T. & Charlesworth, D. Pleiotropic effects of sex-determining genes in the evolution of dioecy in two plant species. Proc. R. Soc. B Biol. Sci. 286, 20191805 (2019).

[8] Barrett, S. C. & Hough, J. Sexual dimorphism in flowering plants. J. Exp. Bot. 64, 67-82 (2013).

[9] Albert, A. Y. K. & Otto, S. P. Sexual selection can resolve sex-linked sexual antagonism. Science 310, 119-121 (2005).

[10] Montalvão, A. P. L., Kersten, B., Fladung, M. & Müller, N. A. The diversity and dynamics of sex determination in dioecious plants. Front. Plant Sci. 11, 580488 (2020).

[11] Feng, G. et al. Pathways to sex determination in plants: how many roads lead to Rome? Curr. Opin. Plant Biol. 54, 61-68 (2020).

[12] Sanderson, B. J. et al. Sex determination through X-Y heterogamety in Salix nigra. Heredity 126, 630-639 (2021).

[13] Yang, G. et al. Sex-related differences in growth, herbivory, and defense of two Salix species. Forests, https://doi.org/10.3390/f11040450 (2020).

[14] Zhou, R. et al. A willow sex chromosome reveals convergent evolution of complex palindromic repeats. Genome Biol. 21, 38 (2020).

[15] Pucholt, P., Ronnberg-Wastljung, A. C. & Berlin, S. Single locus sex determination and female heterogamety in the basket willow (Salix viminalis L.). Heredity 114, 575-583 (2015).

[16] Liu, J. et al. Transcriptome analysis of the differentially expressed genes in the male and female shrub willows (Salix suchowensis). PLoS ONE 8, e60181 (2013).

[17] Tuskan, G. A. et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science 313, 1596-1604 (2006).

[18] Paolucci, I. et al. Genetic linkage maps of Populus alba L. and comparative mapping analysis of sex determination across Populus species. Tree Genet. Genomes 6, 863-875 (2010).

[19] Yang, W. et al. A general model to explain repeated turnovers of sex determination in the Salicaceae. Molec. Biol. Evol. 38, 968-980 (2020).

[20] Gammerdinger, W. J. & Kocher, T. D. Unusual diversity of sex chromosomes in african cichlid fishes. Genes 9, 480 (2018).

[21] Jeffries, D. L. et al. A rapid rate of sex-chromosome turnover and non-random transitions in true frogs. Nat. Commun. 9, 4088 (2018).

[22] Dickmann, D. I. & Kuzovkina, J. Poplars and willows of the world, with emphasis on silviculturally important species. In Poplars and willows: Trees for society and the environment Vol 22 CABI (2014).

[23] Gouker, F. E. et al. Sexual dimorphism in the dioecious willow Salix purpurea. Amer. J. Bot. https://doi.org/10.1002/ajb2.1704 (2021).

[24] Fussel, U., Dotterl, S., Jurgens, A. & Aas, G. Inter- and intraspecific variation in floral scent in the genus Salix and its implication for pollination. J. Chem. Ecol. 33, 749-765 (2007).

[25] Ashman, T. L. Sniffing out patterns of sexual dimorphism in floral scent. Funct. Ecol. 23, 852-862 (2009).

[26] Melnikova, N. V. et al. Sex-specific polymorphism of MET1 and ARR17 genes in Populus x sibirica. Biochimie 162, 26-32 (2019).

[27] Cronk, Q. & Müller, N. A. Default sex and single gene sex determination in dioecious plants. Front. Plant Sci., https://doi.org/10.3389/fpls.2020.01162 (2020).

[28] Robinson, K. M. et al. Populus tremula (European aspen) shows no evidence of sexual dimorphism. BMC Plant Biol. 14, 276 (2014).

[29] Carlson, C. H. et al. Dominance and sexual dimorphism pervade the Salix purpurea L. transcriptome. Genome Biol. Evol. 9, 2377-2394 (2017).

[30] Wu, G. et al. The Sequential action of miR156 and miR172 regulates developmental timing in Arabidopsis. Cell 138, 750-759 (2009).

[31] Akagi, T. et al. A Y-Encoded suppressor of feminization arose via lineage-specific duplication of a cytokinin response regulator in kiwifruit. Plant Cell 30, 780-795 (2018).

[32] Fedoroff, N. V. Transposable elements, epigenetics, and genome evolution. Science 338, 758-767 (2012).

[33] Janoušek, B., Široký, J. & Vyskot, B. Epigenetic control of sexual phenotype in a dioecious plant Melandrium album. Mol. Gen. Genet. MGG 250, 483-490 (1996).

[34] Eamens, A. L., Smith, N. A., Curtin, S. J., Wang, M.-B. & Waterhouse, P. M. The Arabidopsis thaliana double-stranded RNA binding protein DRB1 directs guide strand selection from microRNA duplexes. RNA 15, 2219-2235 (2009).

[35] Ranftl, Q. L., Bastakis, E., Klermund, C. & Schwechheimer, C. LLM-domain containing B-GATA factors control different aspects of cytokinin-regulated development in Arabidopsis thaliana. Plant Physiol. 170, 2295-2311 (2016).

[36] Ocarez, N. & Mejía, N. Suppression of the D-class MADS-box AGL11 gene triggers seedlessness in fleshy fruits. Plant Cell Rep. 35, 239-254 (2016).

[37] Verelst, W., Saedler, H. & Münster, T. MIKC* MADS-protein complexes bind motifs enriched in the proximal region of late pollen-specific Arabidopsis Promoters. Plant Physiol. 143, 447-460 (2007).

[38] Dreni, L. & Zhang, D. Flower development: the evolutionary history and functions of the AGL6 subfamily MADS-box genes. J. Exp. Bot. 67, 1625-1638 (2016).

[39] The Arabidopsis Information Resource (TAIR), https://www.arabidopsis.org/tools/bulk/go/index.jsp, on www.arabidopsis.org, last accessed on 23 May 2021.

[40] Lin, J., Gibbs, J. P. & Smart, L. B. Population genetic structure of native versus naturalized sympatric shrub willows (Salix; Salicaceae). Am. J. Bot. 96, 771-785 (2009).

[41] Gouker, F. E., DiFazio, S. P., Bubner, B., Zander, M. & Smart, L. B. Genetic diversity and population structure of native, naturalized, and cultivated Salix purpurea. Tree Genetics Genomes, https://doi.org/10.1007/s11295-019-1359-0 (2019).

[42] Carlson, C. H. et al. Joint linkage and association mapping of complex traits in shrub willow (Salix purpurea L.). Ann. Bot. 124, 701-716 (2019).

[43] Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15-21 (2013).

[44] Liao, Y., Smyth, G. K. & Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923-930 (2014).

[45] Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 1-21 (2014).

[46] Axtell, M. J. Classification and comparison of small RNAs from plants. Annu Rev. Plant Biol. 64, 137-159 (2013).

[47] Krueger, F. & Andrews, S. R. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics 27, 1571-1572 (2011).

[48] Gaspar, J. M. & Hart, R. P. DMRfinder: efficiently identifying differentially methylated regions from MethylC-seq data. BMC Bioinforma. 18, 1-8 (2017).

[49] Bradbury, P. J. et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633-2635 (2007).

[50] Shabalin, A. A. Matrix eQTL: ultra fast eQTL analysis via large matrix operations. Bioinformatics 28, 1353-1358 (2012).

[51] Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139-140 (2010).

[52] Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma. 9, 1-13 (2008).

[53] Maere, S., Heymans, K. & Kuiper, M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21, 3448-3449 (2005).

[54] Blázquez, M. A. Flower development pathways. J. Cell Sci. 113, 3547-3548, https://doi.org/10.1242/jcs.113.20.3547 (2000).