Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant

  1. Lei Zhao
  2. Wei Zhou
  3. Jun He
  4. De-Zhu Li  Is a corresponding author
  5. Hong-Tao Li  Is a corresponding author
  1. Germplasm Bank of Wild Species & Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, China
  2. Kunming College of Life Science, University of Chinese Academy of Sciences, China
5 figures and 16 additional files

Figures

Floral buds and flowers at anthesis of females (A, B) and males (C) in Trichosanthes pilosa.
Figure 2 with 1 supplement
Sex-biased gene expression for floral buds and flowers at anthesis in males and females of Trichosanthes pilosa.

Volcano plots of average expression between female-biased, male-biased and unbiased genes in floral buds (A) and flowers at anthesis (B). M1 and F1 indicate male and female floral buds; M2 and F2 indicate male and female flowers at anthesis. The value of y coordinate represents -log10(FDR), and the value of x coordinate represents log2(Fold Change) identified by DESeq2. Heatmap of sex-biased gene expression (C) using hierarchical clustering analysis. Hierarchical gene clustering is based on Euclidean distance with an average of log2(FPKM) for differentially expressed genes. The color gradient represents from high to low (from red to green) gene expression.

Figure 2—figure supplement 1
The length distribution of unigenes (A) and the e-value distribution (B) and the species distribution (C) of BLAST hits for each unigene.
Figure 3 with 1 supplement
The overlap between sex-biased and tissue-biased genes in two types of sexes and tissues.

Male-biased genes in floral buds (M1BGs) (A) or flowers at anthesis (M2BGs) (B) overlapped with tissue-biased genes in floral buds (M1TGs) and flowers at anthesis (M2TGs). Female-biased genes in floral buds (F1BGs) (C) or flowers at anthesis (F2BGs) (D) overlapped with tissue-biased genes in floral buds (F1TGs) and flowers at anthesis (F2TGs).

Figure 3—figure supplement 1
The overlap between male-biased genes with faster evolutionary rates and tissue-biased genes in floral buds.

M1BGs indicate male-biased genes with faster evolutionary rates in male floral buds (M1-biased genes). M1TGs indicate tissue-biased genes in male floral buds (M1-tissue-biased genes). 343 M1BGs overlapped with 1755 M1TGs (A). 27 out of 343 M1BGs & M1TGs overlapped with M1BGs_98 (B). 27 out of 343 M1BGs & M1TGs overlapped with M1BGs_18 (C).

Figure 4 with 3 supplements
Violin plots of dN/dS values (0<ω<2) of female-biased, male-biased, and unbiased genes in floral buds and flowers at anthesis of Trichosanthes pilosa.

White dot indicates the median of dN/dS values for sex-biased and unbiased genes. Wilcoxon rank sum tests are used to test for significant differences (***p<0.0005, **p<0.005, and *p<0.05). The distributions of dN/dS values for female-biased, male-biased and unbiased genes in floral buds (A) and flowers at anthesis (B) using ‘two-ratio’ branch model. The distributions of dN/dS values for female-biased, male-biased and unbiased genes in floral buds (C) and flowers at anthesis (D) using ‘free-ratio’ branch model.

Figure 4—figure supplement 1
Boxplot of dN/dS values (including all ω values) of female-biased, male-biased, and unbiased genes in floral buds and flowers at anthesis of Trichosanthes pilosa.

White dot indicates the median of dN/dS values for sex-biased and unbiased genes. Wilcoxon rank sum tests are used to test for significant differences (***p<0.0005, **p<0.005, and *p<0.05). The distributions of dN/dS values for female-biased, male-biased, and unbiased genes in floral buds (A) and flowers at anthesis (B) using ‘free-ratio’ branch model. The distributions of dN/dS values for female-biased, male-biased, and unbiased genes in floral buds (C) and flowers at anthesis (D) using ‘two-ratio’ branch model.

Figure 4—figure supplement 2
Violin plots of effective number of codons (ENCs) values of female-biased, male-biased and unbiased genes in floral buds.

Significant differences using Wilcoxon rank sum tests are represented by * (***p<0.0005 and *p<0.05).

Figure 4—figure supplement 3
Boxplot of dS values of female-biased, male-biased, and unbiased genes in floral buds of dioecious Trichosanthes pilosa using ‘free-ratio’ (A) and ‘two-ratio’ (B) branch model.

Significant differences are represented by * in Wilcoxon rank sum tests (***p<0.0005).

Figure 5 with 3 supplements
Venn diagrams of male-biased genes detected to be under positive selection using aBSREL, BUSTED, CodeML, and RELAX in floral buds.
Figure 5—figure supplement 1
Venn diagrams of female-biased genes under positive selection in floral buds.

Overlaps of female-biased genes were detected to be under positive selection using aBSREL, BUSTED, CodeML, and RELAX.

Figure 5—figure supplement 2
Venn diagrams of unbiased genes under positive selection in floral buds.

Overlaps of unbiased genes were identified to be under positive selection using aBSREL, BUSTED, CodeML, and RELAX.

Figure 5—figure supplement 3
Scatterplots of KEGG pathway of sex-biased genes in female (A) and male (B) floral buds of the dioecious Trichosanthes pilosa.

The y-axis represents the pathway name, and the x-axis represents the ratio of genes corresponding to the pathways.

Additional files

Supplementary file 1

Overview of sequencing reads from 12 samples of male and female plants in Trichosanthes pilosa.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp1-v1.xls
Supplementary file 2

Numbers of unigenes annotated in public databases.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp2-v1.xls
Supplementary file 3

The mapping rate of reads for each sample in floral buds and flowers at anthesis of Trichosanthes pilosa.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp3-v1.xls
Supplementary file 4

dN, dS, and ω values of each female-biased, male-biased, unbiased orthologous genes of floral buds and flowers at anthesis for each species using ‘two-ratio’ branch model of CodeML in PAML.

Two dioecious species Trichosanthes pilosa and T. kirilowii represent the foreground and two monoecious species T. anguina and Luffa cylindrica represent the background.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp4-v1.xls
Supplementary file 5

dN, dS, and ω values of each female-biased, male-biased, unbiased orthologous genes of floral buds and flowers at anthesis for each species using ‘free-ratio’ branch model of CodeML in PAML.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp5-v1.xls
Supplementary file 6

Genes under positive selection identified by branch-site model of CodeML in PAML and function in NR, KEGG, Swissport, and GO databases for male-biased orthologous genes in floral buds.

Two dioecious species Trichosanthes pilosa and T. kirilowii represent the foreground and two monoecious species, T. anguina and Luffa cylindrica represent the background.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp6-v1.xls
Supplementary file 7

Genes under episodic positive selection tested by aBSREL model in HyPhy and functions in NR, KEGG, Swissport, and GO databases for male-biased orthologous genes in floral buds.

The dioecious species Trichosanthes pilosa represents the foreground and other species represent the background.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp7-v1.xls
Supplementary file 8

Genes under episodic positive selection found by BUSTED model in HyPhy and functions in NR, KEGG, Swissport, and GO databases for male-biased orthologous genes in floral buds.

The dioecious species Trichosanthes pilosa represents the foreground (unconstrained branch) and other species represent the background (constrained branch).

https://cdn.elifesciences.org/articles/89941/elife-89941-supp8-v1.xls
Supplementary file 9

Genes under relaxed selection detected by RELAX model in HyPhy and functions in NR, KEGG, Swissport, and GO databases for male-biased orthologous genes in floral buds.

The dioecious species Trichosanthes pilosa represents the foreground (test) and other species represent the background (reference).

https://cdn.elifesciences.org/articles/89941/elife-89941-supp9-v1.xls
Supplementary file 10

Genes under intensified positive selection identified by RELAX model in HyPhy and functions in NR, KEGG, Swissport, and GO databases for male-biased orthologous genes in floral buds.

The dioecious species Trichosanthes pilosa represents the foreground (test) and other species represent the background (reference).

https://cdn.elifesciences.org/articles/89941/elife-89941-supp10-v1.xls
Supplementary file 11

KEGG pathway enrichment analysis of female-biased and male-biased genes in floral buds of the dioecious Trichosanthes pilosa.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp11-v1.xls
Supplementary file 12

Functions and references associated with abiotic stress and immune responses, organ developments of male-biased genes under significant positive selection (p<0.05) in floral buds.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp12-v1.docx
Supplementary file 13

Functions and references associated with abiotic stress and immune responses, organ developments of male-biased genes under significant relaxed selection (p<0.05) in floral buds.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp13-v1.docx
Supplementary file 14

The expressions and functions of some male-biased genes associated with senescence, raceme inflorescence development and early flowering in floral buds.

https://cdn.elifesciences.org/articles/89941/elife-89941-supp14-v1.xls
MDAR checklist
https://cdn.elifesciences.org/articles/89941/elife-89941-mdarchecklist1-v1.docx
Source data 1

Reference transcriptome, orthology data, and alignments.

https://cdn.elifesciences.org/articles/89941/elife-89941-data1-v1.zip

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Lei Zhao
  2. Wei Zhou
  3. Jun He
  4. De-Zhu Li
  5. Hong-Tao Li
(2024)
Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant
eLife 12:RP89941.
https://doi.org/10.7554/eLife.89941.5