Relaxed Selection Research Articles

Genomic stasis over millions of years in subseafloor sediment.

One of the significant challenges in microbiology is to understand the extent and mechanisms of evolution within life beneath the surface of the Earth. The population bottleneck that microbes in deep marine sediment experience implies that mutational and population genetic forces could lead to higher levels of relaxed selection and an increase in pseudogenes. To investigate this hypothesis, a group of Thalassospira strains were isolated from subseafloor sediment that is 3 to 6 million years old, as reported by Orsi and colleagues in 2021. These isolates, representing lineages that have been buried for millions of years, offer an excellent opportunity to study the evolution of life beneath the seafloor over a long period. The existence of closely related strains from environments on the surface of the Earth enabled us to examine the impact of selection within each group. We discovered that isolates from beneath the seafloor show lineage-specific similarities to Thalassospira from the surface world, both in the overall intensity of selection on the genome and in the specific genes affected by mutation. We found no signs of increased relaxed selection or other notable genomic changesin the genomes of the Thalassospira isolates from beneath the seafloor, suggesting that these subseafloor isolates were awakened from a million-year near-stasis. The unique genomic characteristics of each Thalassospira lineage from beneath the seafloor must then reflect genetic changes that surface-inhabiting decendants acquired in the past 3-6 million years. Remarkably, Thalassospira lineages beneath the surface appear to have stably maintained their genomes in the midst of metabolic dormancy and extremely long generation times.

Plastid genome evolution in leafless members of the orchid subfamily Orchidoideae, with a focus on Degranvillea dermaptera.

Leafless, heterotrophic plants are prime examples of organismal modification, the genomic consequences of which have received considerable interest. In particular, plastid genomes (plastomes) are being sequenced at a high rate, allowing continual refinement of conceptual models of reductive evolution in heterotrophs. However, numerous sampling gaps exist, hindering the ability to conduct comprehensive phylogenomic analyses in these plants. Using floral tissue from an herbarium specimen, we sequenced and analyzed the plastome of Degranvillea dermaptera, a rarely collected, leafless orchid species from South America about which little is known, including its phylogenetic affinities. The plastome is the most reduced of those sequenced among the orchid subfamily Orchidoideae. In Degranvillea, it has lost the majority of genes found in leafy autotrophic species, is structurally rearranged, and has similar gene content to the most reduced plastomes among the orchids. We found strong evidence for the placement of Degranvillea within the subtribe Spiranthinae using models that explicitly account for heterotachy, or lineage-specific evolutionary rate variation over time. We further found evidence of relaxed selection on several genes and of correlations among substitution rates and several other "traits" of the plastome among leafless members of orchid subfamily Orchidoideae. Our findings advance knowledge on the phylogenetic relationships and paths of plastid genome evolution among the orchids, which have experienced more independent transitions to heterotrophy than any other plant family. This study demonstrates the importance of herbarium collections in comparative genomics of poorly known species of conservation concern.

Comparative Genomic Analysis of the Pattern of Evolution of Male and Female Reproductive Proteins in Seed Beetles.

Male seminal fluid proteins often show signs of positive selection and divergent evolution, believed to reflect male-female coevolution. Yet, our understanding of the predicted concerted evolution of seminal fluid proteins and female reproductive proteins is limited. We sequenced, assembled, and annotated the genome of two species of seed beetles allowing a comparative analysis of four closely related species of these herbivorous insects. We compare the general pattern of evolution in genes encoding seminal fluid proteins and female reproductive proteins with those in digestive protein genes and well-conserved reference genes. We found that female reproductive proteins showed an overall ratio of nonsynonymous to synonymous substitutions (ω) similar to that of conserved genes, while seminal fluid proteins and digestive proteins exhibited higher overall ω values. Further, seminal fluid proteins and digestive proteins showed a higher proportion of sites putatively under positive selection, and explicit tests showed no difference in relaxed selection between protein types. Evolutionary rate covariation analyses showed that evolutionary rates among seminal fluid proteins were on average more closely correlated with those in female reproductive proteins than with either digestive or conserved genes. Gene expression showed the expected negative covariation with ω values, except for male-biased genes where this negative relationship was reversed. In conclusion, seminal fluid proteins showed relatively rapid evolution and signs of positive selection. In contrast, female reproductive proteins evolved at a lower rate under selective constraints, on par with genes known to be well conserved. Although our findings provide support for concerted evolution of seminal fluid proteins and female reproductive proteins, they also suggest that these two classes of proteins evolve under partly distinct selective regimes.

Mating system variation and gene expression in the male reproductive tract of Peromyscus mice.

Genes involved in reproduction often evolve rapidly at the sequence level due to postcopulatory sexual selection (PCSS) driven by male-male competition and male-female sexual conflict, but the impact of PCSS on gene expression has been under-explored. Further, though multiple tissues contribute to male reproductive success, most studies have focused on the testes. To explore the influence of mating system variation on reproductive tract gene expression in natural populations, we captured adult males from monogamous Peromyscus californicus and polygynandrous P. boylii and P. maniculatus. We generated RNAseq libraries, quantified gene expression in the testis, seminal vesicle, epididymis, and liver, and identified 3627 mating system-associated differentially expressed genes (MS-DEGs), where expression shifted in the same direction in P. maniculatus and P. boylii relative to P. californicus. Gene expression variation was most strongly associated with mating behaviour in the seminal vesicles, where 89% of differentially expressed genes were MS-DEGs, including the key seminal fluid proteins Svs2 and Pate4. We also used published rodent genomes to test for positive and relaxed selection on Peromyscus-expressed genes. Though we did not observe more overlap than expected by chance between MS-DEGs and positively selected genes, 203 MS-DEGs showed evidence of positive selection. Fourteen reproductive genes were under tree-wide positive selection but convergent relaxed selection in P. californicus and Microtus ochrogaster, a distantly related monogamous species. Changes in transcript abundance and gene sequence evolution in association with mating behaviour suggest that male mice may respond to sexual selection intensity by altering aspects of sperm motility, sperm-egg binding and copulatory plug formation.

Possible future evolutionary consequences to Homo as a result of the implementation of biotechnology

Biotechnology has become one of the most powerful forces on the planet, since it is capable of altering life processes at a molecular level. Since human bodies are dynamic biological systems, medicine requires to understand the evolutionary antecedents of Homo, especially in relation to neurohormonal regulation. Furthermore, increasing human dependence on biotechnology has led to relaxed natural selection in Homo, with subsequent increase of genetic load. In this paper, we speculate on the possible consequences of the application of parsimoniously derived biotechnologies onto the biological system of humans, with special attention to three areas: 1. human brain augmentation; 2. biotechnology and public health; 3. relaxed natural selection and genetic load. Human ability to manipulate and alter the structure and function of the body may not only make natural selection redundant but will be guided by a teleology whose purpose will seek to improve upon nature’s design.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments.

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family encodes phospholipases essential for the degradation of organelles in the lens of the eye. These phospholipases translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny and selection of plaat1 across bony fishes and tetrapods. We showed that plaat1 (probably ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual-acuity and blind mammals and fishes. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

Genomic signatures associated with recurrent scale loss in cyprinid fish.

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Genetic Causes and Genomic Consequences of Breakdown of Distyly in Linum trigynum.

Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility. Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male self-incompatibility and/or anther height could underlie homostyly and self-compatibility in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.

Impact of pre-copulatory sexual cannibalism on genetic diversity and efficacy of selection.

Factors that increase reproductive variance among individuals act to reduce effective population size (Ne), which accelerates the loss of genetic diversity and decreases the efficacy of purifying selection. These factors include sexual cannibalism, offspring investment and mating system. Pre-copulatory sexual cannibalism, where the female consumes the male prior to mating, exacerbates this effect. We performed comparative transcriptomics in two spider species, the cannibalistic Trechaleoides biocellata and the non-cannibalistic T. keyserlingi, to generate genomic evidence to support these predictions. First, we estimated heterozygosity and found that genetic diversity is relatively lower in the cannibalistic species. Second, we calculated dN/dS ratios as a measure of purifying selection; a higher dN/dS ratio indicated relaxed purifying selection in the cannibalistic species. These results are consistent with the hypothesis that sexual cannibalism impacts operational sex ratio and demographic processes, which interact with evolutionary forces to shape the genetic structure of populations. However, other factors such as the mating system and life-history traits contribute to shaping Ne. Comparative analyses across multiple contrasting species pairs would be required to disentangle these effects. Our study highlights that extreme behaviours such as pre-copulatory cannibalism may have profound eco-evolutionary effects.

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Drought-induced circular RNAs in maize roots: Separating signal from noise.

Circular RNAs (circRNAs) play an important role in diverse biological processes; however, their origin and functions, especially in plants, remain largely unclear. Here, we used 2 maize (Zea mays) inbred lines, as well as 14 of their derivative recombination inbred lines with different drought sensitivity, to systematically characterize 8,790 circRNAs in maize roots under well-watered (WW) and water-stress (WS) conditions. We found that a diverse set of circRNAs expressed at significantly higher levels under WS. Enhanced expression of circRNAs was associated with longer flanking introns and an enrichment of long interspersed nuclear element retrotransposable elements. The epigenetic marks found at the back-splicing junctions of circRNA-producing genes were markedly different from canonical splicing, characterized by increased levels of H3K36me3/H3K4me1, as well as decreased levels of H3K9Ac/H3K27Ac. We found that genes expressing circRNAs are subject to relaxed selection. The significant enrichment of trait-associated sites along their genic regions suggested that genes giving rise to circRNAs were associated with plant survival rate under drought stress, implying that circRNAs play roles in plant drought responses. Furthermore, we found that overexpression of circMED16, one of the drought-responsive circRNAs, enhances drought tolerance in Arabidopsis (Arabidopsis thaliana). Our results provide a framework for understanding the intricate interplay of epigenetic modifications and how they contribute to the fine-tuning of circRNA expression under drought stress.

DNA Damage, Genome Stability, and Adaptation: A Question of Chance or Necessity?

DNA damage causes the mutations that are the principal source of genetic variation. DNA damage detection and repair mechanisms therefore play a determining role in generating the genetic diversity on which natural selection acts. Speciation, it is commonly assumed, occurs at a rate set by the level of standing allelic diversity in a population. The process of speciation is driven by a combination of two evolutionary forces: genetic drift and ecological selection. Genetic drift takes place under the conditions of relaxed selection, and results in a balance between the rates of mutation and the rates of genetic substitution. These two processes, drift and selection, are necessarily mediated by a variety of mechanisms guaranteeing genome stability in any given species. One of the outstanding questions in evolutionary biology concerns the origin of the widely varying phylogenetic distribution of biodiversity across the Tree of Life and how the forces of drift and selection contribute to shaping that distribution. The following examines some of the molecular mechanisms underlying genome stability and the adaptive radiations that are associated with biodiversity and the widely varying species richness and evenness in the different eukaryotic lineages.

Variations and reduction of plastome are associated with the evolution of parasitism in Convolvulaceae.

Parasitic lifestyle can often relax the constraint on the plastome, leading to gene pseudogenization and loss, and resulting in diverse genomic structures and rampant genome degradation. Although several plastomes of parasitic Cuscuta have beenreported, the evolution of parasitism in the family Convolvulaceae which is linked to structural variations and reduction of plastome has not been well investigated. In this study, we assembled and collected 40 plastid genomes belonging to 23 species representing four subgenera of Cuscuta and ten species of autotrophic Convolvulaceae. Our findings revealed nine types of structural variations and six types of inverted repeat (IR) boundary variations in the plastome of Convolvulaceae spp. These structuralvariations were associated with the shift of parasitic lifestyle, and IR boundary shift, as well as the abundance of long repeats. Overall, the degradation of Cuscuta plastome proceeded gradually, with one clade exhibiting an accelerated degradation rate. We observed five stages of gene loss in Cuscuta, including NAD(P)H complex → PEP complex → Photosynthesis-related → Ribosomal protein subunits → ATP synthase complex. Based on our results, we speculated that the shift of parasitic lifestyle in early divergent time promoted relaxed selection on plastomes, leading to the accumulation of microvariations, which ultimately resulted in the plastome reduction. This study provides new evidence towardsa better understanding of plastomic evolution, variation, and reduction in the genus Cuscuta.

Comparative genomics on chloroplasts of Rubus (Rosaceae)

Rubus, the largest genus in Rosaceae, contains over 1400 species that distributed in multiple habitats across the world, with high species diversity in the temperate regions of Northern Hemisphere. Multiple Rubus species are cultivated for their valuable fruits. However, the intrageneric classification and phylogenetic relationships are still poorly understood. In this study, we sequenced, assembled, and characterized 17 plastomes of Rubus, and conducted comparative genomics integrating with 47 previously issued plastomes of this genus. The 64 plastomes of Rubus exhibited typical quadripartite structure with sizes ranging from 155,144 to 156,700 bp, and contained 132 genes including 87 protein-coding genes, 37 tRNA genes and eight rRNA genes. All plastomes are conservative in the gene order, the frequency of different types of long repeats and simple sequence repeats (SSRs), the codon usage, and the selection pressure of protein-coding genes. However, there are also some differences in the Rubus plastomes, including slight contraction and expansion of the IRs, a variation in the numbers of SSRs and long repeats, and some genes in certain clades undergoing intensified or relaxed purifying selection. Phylogenetic analysis based on whole plastomes showed that the monophyly of Rubus was strongly supported and resolved it into six clades corresponding to six subgenera. Moreover, we identified 12 highly variable regions that could be potential molecular markers for phylogenetic, population genetic, and barcoding studies. Overall, our study provided insight into plastomic structure and sequence diversification of Rubus, which could be beneficial for future studies on identification, evolution, and phylogeny in this genus.

Parthenogenetic Stick Insects Exhibit Signatures of Preservation in the Molecular Architecture of Male Reproduction.

After the loss of a trait, theory predicts that the molecular machinery underlying its phenotypic expression should decay. Yet, empirical evidence is contrasting. Here, we test the hypotheses that (i) the molecular ground plan of a lost trait could persist due to pleiotropic effects on other traits and (ii) that gene co-expression network architecture could constrain individual gene expression. Our testing ground has been the Bacillus stick insect species complex, which contains close relatives that are either bisexual or parthenogenetic. After the identification of genes expressed in male reproductive tissues in a bisexual species, we investigated their gene co-expression network structure in two parthenogenetic species. We found that gene co-expression within the male gonads was partially preserved in parthenogens. Furthermore, parthenogens did not show relaxed selection on genes upregulated in male gonads in the bisexual species. As these genes were mostly expressed in female gonads, this preservation could be driven by pleiotropic interactions and an ongoing role in female reproduction. Connectivity within the network also played a key role, with highly connected-and more pleiotropic-genes within male gonad also having a gonad-biased expression in parthenogens. Our findings provide novel insight into the mechanisms which could underlie the production of rare males in parthenogenetic lineages; more generally, they provide an example of the cryptic persistence of a lost trait molecular architecture, driven by gene pleiotropy on other traits and within their co-expression network.

Paradoxical Exception to Island Tameness: Increased Defensiveness in an Insular Population of Rattlesnakes.

Island tameness results largely from a lack of natural predators. Because some insular rattlesnake populations lack functional rattles, presumably the consequence of relaxed selection from reduced predation, we hypothesized that the Santa Catalina Island, California, USA, population of the southern Pacific rattlesnake (Crotalus helleri, which possesses a functional rattle), would exhibit a decrement in defensive behavior relative to their mainland counterparts. Contrary to our prediction, rattlesnakes from the island not only lacked tameness compared to mainland snakes, but instead exhibited measurably greater levels of defensiveness. Island snakes attempted to bite 4.7 times more frequently as we endeavored to secure them by hand, and required 2.1-fold more time to be pinned and captured. When induced to bite a beaker after being grasped, the island snakes also delivered 2.1-fold greater quantities of venom when controlling for body size. The additional venom resulted from 2.1-fold larger pulses of venom ejected from the fangs. We found no effects of duration in captivity (2-36 months), which suggests an absence of long-term habituation of antipredator behaviors. Breeding bird surveys and Christmas bird counts indicated reduced population densities of avian predators on Catalina compared to the mainland. However, historical estimates confirmed that populations of foxes and introduced mammalian predators (cats and pigs) and antagonists (herbivorous ungulates) substantially exceeded those on the mainland in recent centuries, and therefore best explain the paradoxically exaggerated defensive behaviors exhibited by Catalina's rattlesnakes. These findings augment our understanding of anthropogenic effects on the behaviors of island animals and underscore how these effects can negatively affect human safety.

The Global Decline in Human Fertility: The Post-Transition Trap Hypothesis.

Over the past half a century many countries have witnessed a rapid fall in total fertility rates, particularly in the world's most advanced economies including the industrial powerhouses of Eastern Asia and Europe. Such nations have now passed through the first and second demographic transitions and are currently exhibiting fertility rates well below the replacement threshold of 2.1, with no sign of recovery. This paper examines the factors responsible for driving these demographic transitions and considers their impact on both fertility and fecundity (our fundamental capacity to reproduce). I argue that because the first demographic transition was extremely rapid and largely driven by socioeconomic factors, it has had no lasting impact on the genetic/epigenetic underpinnings of human fecundity. However, the second demographic transition will be different. A series of conditions associated with low fertility societies, including relaxed selection pressure for high-fertility genotypes, the indiscriminate use of assisted reproductive technologies to treat human infertility, and environmental contamination with reproductive toxicants, may impact our genetic constitution in ways that compromise the future fecundity of our species. Since any fundamental change in the genetic foundations of human reproduction will be difficult to reverse, we should actively pursue methods to monitor human fecundity, as sub-replacement fertility levels become established across the globe.

Evolution of intraspecific floral variation in a generalist-specialist pollination system.

Intraspecific processes impact macroevolutionary patterns through individual variation, selection, and ecological specialisation. According to the niche variation hypothesis, the broader ecological niche of gen- eralist species results in an increased morphological variation among individuals, either because they are constituted of diversified specialised individuals each exploiting a fraction of the species' niche, or because they are constituted of true generalist individuals that experience relaxed selection. To test this hypoth- esis, we surveyed the individual floral morphology of species of Antillean Gesneriaceae, a group that has transitioned between specialisation for hummingbird pollination and generalisation multiple times throughout its evolutionary history. We characterised the profiles of corollas using geometric morpho- metrics and compared the intraspecific shape variance of specialists and generalists in a phylogenetic context. We used three approaches that differently accounted for the high dimensionality of morphologi- cal traits, the ancestral reconstruction of pollination syndromes over time, and the error associated with the estimation of the intraspecific variance. Our findings provide partial support for the niche variation hypothesis. If considering the whole shape in the analysis corroborated this idea, decomposing the shape into principal components indicated that not all aspects of the corolla exhibit the same pattern of vari- ation. Specifically, pollination generalists tend to display greater intraspecific variation than specialists in terms of tubularity, but not of curvature. Accounting for the error in the variance estimation also reduced the support for the hypothesis, suggesting that larger sample sizes may be required to reach stronger conclusions. This study emphasises the reciprocal influence between plants and their pollinators on floral morphology at different biodiversity scales, and suggests that ecological strategies of species can affect patterns of morphological variation at macroevolutionary scales.

The evolution of seminal fluid gene expression and postmating reproductive isolation in Drosophila.

Seminal fluid protein (Sfp) genes show, in general, a higher rate of sequence divergence than genes from other categories, which is often attributed to forms of postcopulatory sexual selection or sexual conflict. Recently, the relaxation of selective constraints has been proposed as an alternative explanation for the rapid sequence evolution of Sfps and other genes with sex-limited expression. The expression of Sfp genes is a likely target of selection, but the evolution of differences in their expression levels is less understood. Here, we explore both polymorphism and divergence in Sfp gene expression between Drosophila melanogaster and Drosophila simulans, how selection might have influenced their expression, and whether changes in expression might trigger the evolution of reproductive isolating barriers. In our analysis, Sfp genes showed higher divergence, but not higher polymorphism, in expression than the average male reproductive glands gene. Sfp genes with reproductive-tissue-specific expression were enriched for both directional and stabilizing selection, while relaxed selection was the predominant mode of evolution among Sfp genes with any other nonreproductive tissue-specific or nontissue-specific expression. The knockdown of single genes known to affect intraspecific sperm competition, and with patterns of expression divergence and polymorphism suggestive of directional selection, was not enough to break down postmating reproductive isolation barriers between species. Our results identify the expression of male-specific Sfp genes as an enriched target of selection and suggest a complex molecular relationship between postcopulatory sexual selection on a single gene's expression and its effect on the onset of speciation.

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

Sex-biased genes offer insights into the evolution of sexual dimorphism. Sex-biased genes, especially those with male bias, show elevated evolutionary rates of protein sequences driven by positive selection and relaxed purifying selection in animals. Although rapid sequence evolution of sex-biased genes and evolutionary forces have been investigated in animals and brown algae, less is known about evolutionary forces in dioecious angiosperms. In this study, we separately compared the expression of sex-biased genes between female and male floral buds and between female and male flowers at anthesis in dioecious Trichosanthes pilosa (Cucurbitaceae). In floral buds, sex-biased gene expression was pervasive, and had significantly different roles in sexual dimorphism such as physiology. We observed higher rates of sequence evolution for male-biased genes in floral buds compared to female-biased and unbiased genes. Male-biased genes under positive selection were mainly associated with functions to abiotic stress and immune responses, suggesting that high evolutionary rates are driven by adaptive evolution. Additionally, relaxed purifying selection may contribute to accelerated evolution in male-biased genes generated by gene duplication. Our findings, for the first time in angiosperms, suggest evident rapid evolution of male-biased genes, advance our understanding of the patterns and forces driving the evolution of sexual dimorphism in dioecious plants.