Stochastic Modeling of Gene Expression Evolution Uncovers Tissue- and Sex-Specific Properties of Expression Evolution in the Drosophila Genus

Stochastic modeling of single-cell gene expression adaptation reveals non-genomic contribution to evolution of tumor subclones

Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present a formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Sublines previously observed to be resistant to anti-CTLA4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression. A record of this paper's transparent peer review process is included in the supplemental information.

cis- and trans-regulatory contributions to a hierarchy of factors influencing gene expression variation

Gene expression variation results from numerous sources including genetic, environmental, life stage, and even the environment experienced by previous generations. While the importance of each has been demonstrated in diverse organisms, their relative contributions remain understudied because few investigations have simultaneously determined each within a single experiment. Here we quantified genome-wide gene expression traits in Drosophila, quantified the contribution of multiple different sources of trait variation and determined the molecular mechanisms underlying observed variation. Our results show that there is a clear hierarchy in our data with genome and developmental stage contributing on average considerably more than current and finally previous generation environmental effects. We also determined the role of cis and trans-regulatory changes across different sources of trait variation, highlighting their importance in adaptation and environmental responses and showing unexpectedly that transgenerational effects herein were predominantly associated with changes in trans-regulation.

Stochastic modelling of single-cell gene expression adaptation reveals non-genomic contribution to evolution of tumor subclones

Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present the first formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Interestingly, sublines previously observed to be resistant to anti-CTLA-4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression.

Evolution of chemosensory tissues and cells across ecologically diverse Drosophilids

Chemosensory tissues exhibit significant between-species variability, yet the evolution of gene expression and cell types underlying this diversity remain poorly understood. To address these questions, we conducted transcriptomic analyses of five chemosensory tissues from six Drosophila species and integrated the findings with single-cell datasets. While stabilizing selection predominantly shapes chemosensory transcriptomes, thousands of genes in each tissue have evolved expression differences. Genes that have changed expression in one tissue have often changed in multiple other tissues but at different past epochs and are more likely to be cell type-specific than unchanged genes. Notably, chemosensory-related genes have undergone widespread expression changes, with numerous species-specific gains/losses including novel chemoreceptors expression patterns. Sex differences are also pervasive, including a D. melanogaster-specific excess of male-biased expression in sensory and muscle cells in its forelegs. Together, our analyses provide new insights for understanding evolutionary changes in chemosensory tissues at both global and individual gene levels.

The seminal odorant binding protein Obp56g is required for mating plug formation and male fertility in Drosophila melanogaster

In Drosophila melanogaster and other insects, the seminal fluid proteins (SFPs) and male sex pheromones that enter the female with sperm during mating are essential for fertility and induce profound post-mating effects on female physiology. The SFPs in D. melanogaster and other taxa include several members of the large gene family known as odorant binding proteins (Obps). Work in Drosophila has shown that some Obp genes are highly expressed in the antennae and can mediate behavioral responses to odorants, potentially by binding and carrying these molecules to odorant receptors. These observations have led to the hypothesis that the seminal Obps might act as molecular carriers for pheromones or other compounds important for male fertility, though functional evidence in any species is lacking. Here, we used functional genetics to test the role of the seven seminal Obps in D. melanogaster fertility and the post-mating response (PMR). We found that Obp56g is required for male fertility and the induction of the PMR, whereas the other six genes are dispensable. We found males lacking Obp56g fail to form a mating plug in the mated female's reproductive tract, leading to ejaculate loss and reduced sperm storage, likely due to its expression in the male ejaculatory bulb. We also examined the evolutionary history of these seminal Obp genes, as several studies have documented rapid evolution and turnover of SFP genes across taxa. We found extensive lability in gene copy number and evidence of positive selection acting on two genes, Obp22a and Obp51a. Comparative RNAseq data from the male reproductive tract of multiple Drosophila species revealed that Obp56g shows high male reproductive tract expression in a subset of taxa, though conserved head expression across the phylogeny. Together, these functional and expression data suggest that Obp56g may have been co-opted for a reproductive function over evolutionary time.

The seminal odorant binding protein Obp56g is required for mating plug formation and male fertility in Drosophila melanogaster

In Drosophila melanogaster and other insects, the seminal fluid proteins (SFPs) and male sex pheromones that enter the female with sperm during mating are essential for fertility and induce profound post-mating effects on female physiology and behavior. The SFPs in D. melanogaster and other taxa include several members of the large gene family known as odorant binding proteins (Obps). Previous work in Drosophila has shown that some Obp genes are highly expressed in the antennae and can mediate behavioral responses to odorants, potentially by binding and carrying these molecules to odorant receptors. These observations have led to the hypothesis that the seminal Obps might act as molecular carriers for pheromones or other compounds important for male fertility in the ejaculate, though functional evidence in any species is lacking. Here, we used RNAi and CRISPR/Cas9 generated mutants to test the role of the seven seminal Obps in D. melanogaster fertility and the post-mating response (PMR). We found that Obp56g is required for male fertility and the induction of the PMR, whereas the other six genes had no effect on fertility when mutated individually. Obp56g is expressed in the male's ejaculatory bulb, an important tissue in the reproductive tract that synthesizes components of the mating plug. We found males lacking Obp56g fail to form a mating plug in the mated female's reproductive tract, leading to ejaculate loss and reduced sperm storage. We also examined the evolutionary history of these seminal Obp genes, as several studies have documented rapid evolution and turnover of SFP genes across taxa. We found extensive lability in gene copy number and evidence of positive selection acting on two genes, Obp22a and Obp51a. Comparative RNAseq data from the male reproductive tract of multiple Drosophila species revealed that Obp56g shows high male reproductive tract expression only in species of the melanogaster and obscura groups, though conserved head expression in all species tested. Together, these functional and expression data suggest that Obp56g may have been co-opted for a reproductive function over evolutionary time.