The Origin and Evolution of Sex Peptide and Sex Peptide Receptor Interactions

Post-mating responses play a vital role in successful reproduction across diverse species. In fruit flies, sex peptide binds to the sex peptide receptor, triggering a series of post-mating responses. However, the origin of sex peptide receptor predates the emergence of sex peptide. The evolutionary origins of the interactions between sex peptide and sex peptide receptor and the mechanisms by which they interact remain enigmatic. In this study, we used ancestral sequence reconstruction, AlphaFold2 predictions, and molecular dynamics simulations to study sex peptide-sex peptide receptor interactions and their origination. Using AlphaFold2 and long-time molecular dynamics simulations, we predicted the structure and dynamics of sex peptide-sex peptide receptor interactions. We show that sex peptide potentially binds to the ancestral states of Diptera sex peptide receptor. Notably, we found that only a few amino acid changes in sex peptide receptor are sufficient for the formation of sex peptide-sex peptide receptor interactions. Ancestral sequence reconstruction and molecular dynamics simulations further reveal that sex peptide receptor interacts with sex peptide through residues that are mostly involved in the interaction interface of an ancestral ligand, myoinhibitory peptides. We propose a potential mechanism whereby sex peptide-sex peptide receptor interactions arise from the preexisting myoinhibitory peptides-sex peptide receptor interface as well as early chance events both inside and outside the preexisting interface that created novel sex peptide-specific sex peptide-sex peptide receptor interactions. Our findings provide new insights into the origin and evolution of sex peptide-sex peptide receptor interactions and their relationship with myoinhibitory peptides-sex peptide receptor interactions.

The Origin and Evolution of Sex Peptide and Sex Peptide Receptor Interactions

Post-mating responses play a vital role in successful reproduction across diverse species. In fruit flies, sex peptide (SP) binds to the sex peptide receptor (SPR), triggering a series of post-mating responses. However, the origin of SPR predates the emergence of SP. The evolutionary origins of the interactions between SP and SPR and the mechanisms by which they interact remain enigmatic. In this study, we used ancestral sequence reconstruction, AlphaFold2 predictions, and molecular dynamics simulations to study SP-SPR interactions and their origination. Using AlphaFold2 and long-time molecular dynamics (MD) simulations, we predicted the structure and dynamics of SP-SPR interactions. We show that SP potentially binds to the ancestral states of Diptera SPR. Notably, we found that only a few amino acid changes in SPR are sufficient for the formation of SP-SPR interactions. Ancestral sequence reconstruction and MD simulations further reveal that SPR interacts with SP through residues that are mostly involved in the interaction interface of an ancestral ligand, myoinhibitory peptides (MIPs). We propose a potential mechanism whereby SP-SPR interactions arise from the pre-existing MIP-SPR interface as well as early chance events both inside and outside the pre-existing interface that created novel SP-specific SP-SPR interactions. Our findings provide new insights into the origin and evolution of SP-SPR interactions and their relationship with MIP-SPR interactions.

The evolution and mutational robustness of chromatin accessibility in Drosophila

BACKGROUND

The evolution of genomic regulatory regions plays a critical role in shaping the diversity of life. While this process is primarily sequence-dependent, the enormous complexity of biological systems complicates the understanding of the factors underlying regulation and its evolution. Here, we apply deep neural networks as a tool to investigate the sequence determinants underlying chromatin accessibility in different species and tissues of Drosophila.

RESULTS

We train hybrid convolution-attention neural networks to accurately predict ATAC-seq peaks using only local DNA sequences as input. We show that our models generalize well across substantially evolutionarily diverged species of insects, implying that the sequence determinants of accessibility are highly conserved. Using our model to examine species-specific gains in accessibility, we find evidence suggesting that these regions may be ancestrally poised for evolution. Using in silico mutagenesis, we show that accessibility can be accurately predicted from short subsequences in each example. However, in silico knock-out of these sequences does not qualitatively impair classification, implying that accessibility is mutationally robust. Subsequently, we show that accessibility is predicted to be robust to large-scale random mutation even in the absence of selection. Conversely, simulations under strong selection demonstrate that accessibility can be extremely malleable despite its robustness. Finally, we identify motifs predictive of accessibility, recovering both novel and previously known motifs.

CONCLUSIONS

These results demonstrate the conservation of the sequence determinants of accessibility and the general robustness of chromatin accessibility, as well as the power of deep neural networks to explore fundamental questions in regulatory genomics and evolution.

The evolution and mutational robustness of chromatin accessibility in Drosophila

The evolution of regulatory regions in the genome plays a critical role in shaping the diversity of life. While this process is primarily sequence-dependent, the enormous complexity of biological systems has made it difficult to understand the factors underlying regulation and its evolution. Here, we apply deep neural networks as a tool to investigate the sequence determinants underlying chromatin accessibility in different tissues of Drosophila. We train hybrid convolution-attention neural networks to accurately predict ATAC-seq peaks using only local DNA sequences as input. We show that a model trained in one species has nearly identical performance when tested in another species, implying that the sequence determinants of accessibility are highly conserved. Indeed, model performance remains excellent even in distantly-related species. By using our model to examine species-specific gains in chromatin accessibility, we find that their orthologous inaccessible regions in other species have surprisingly similar model outputs, suggesting that these regions may be ancestrally poised for evolution. We then use in silico saturation mutagenesis to reveal evidence of selective constraint acting specifically on inaccessible chromatin regions. We further show that chromatin accessibility can be accurately predicted from short subsequences in each example. However, in silico knock-out of these sequences does not qualitatively impair classification, implying that chromatin accessibility is mutationally robust. Subsequently, we demonstrate that chromatin accessibility is predicted to be robust to large-scale random mutation even in the absence of selection. We also perform in silico evolution experiments under the regime of strong selection and weak mutation (SSWM) and show that chromatin accessibility can be extremely malleable despite its mutational robustness. However, selection acting in different directions in a tissue-specific manner can substantially slow adaptation. Finally, we identify motifs predictive of chromatin accessibility and recover motifs corresponding to known chromatin accessibility activators and repressors. These results demonstrate the conservation of the sequence determinants of accessibility and the general robustness of chromatin accessibility, as well as the power of deep neural networks as tools to answer fundamental questions in regulatory genomics and evolution.

Intermolecular interactions drive protein adaptive and co-adaptive evolution at both species and population levels

Proteins are the building blocks for almost all the functions in cells. Understanding the molecular evolution of proteins and the forces that shape protein evolution is essential in understanding the basis of function and evolution. Previous studies have shown that adaptation frequently occurs at the protein surface, such as in genes involved in host–pathogen interactions. However, it remains unclear whether adaptive sites are distributed randomly or at regions associated with particular structural or functional characteristics across the genome, since many proteins lack structural or functional annotations. Here, we seek to tackle this question by combining large-scale bioinformatic prediction, structural analysis, phylogenetic inference, and population genomic analysis of Drosophila protein-coding genes. We found that protein sequence adaptation is more relevant to function-related rather than structure-related properties. Interestingly, intermolecular interactions contribute significantly to protein adaptation. We further showed that intermolecular interactions, such as physical interactions, may play a role in the coadaptation of fast-adaptive proteins. We found that strongly differentiated amino acids across geographic regions in protein-coding genes are mostly adaptive, which may contribute to the long-term adaptive evolution. This strongly indicates that a number of adaptive sites tend to be repeatedly mutated and selected throughout evolution in the past, present, and maybe future. Our results highlight the important roles of intermolecular interactions and coadaptation in the adaptive evolution of proteins both at the species and population levels.

Behavioral and Genomic Sensory Adaptations Underlying the Pest Activity of Drosophila suzukii

Studying how novel phenotypes originate and evolve is fundamental to the field of evolutionary biology as it allows us to understand how organismal diversity is generated and maintained. However, determining the basis of novel phenotypes is challenging as it involves orchestrated changes at multiple biological levels. Here, we aim to overcome this challenge by using a comparative species framework combining behavioral, gene expression, and genomic analyses to understand the evolutionary novel egg-laying substrate-choice behavior of the invasive pest species Drosophila suzukii. First, we used egg-laying behavioral assays to understand the evolution of ripe fruit oviposition preference in D. suzukii compared with closely related species D. subpulchrella and D. biarmipes as well as D. melanogaster. We show that D. subpulchrella and D. biarmipes lay eggs on both ripe and rotten fruits, suggesting that the transition to ripe fruit preference was gradual. Second, using two-choice oviposition assays, we studied how D. suzukii, D. subpulchrella, D. biarmipes, and D. melanogaster differentially process key sensory cues distinguishing ripe from rotten fruit during egg-laying. We found that D. suzukii's preference for ripe fruit is in part mediated through a species-specific preference for stiff substrates. Last, we sequenced and annotated a high-quality genome for D. subpulchrella. Using comparative genomic approaches, we identified candidate genes involved in D. suzukii's ability to seek out and target ripe fruits. Our results provide detail to the stepwise evolution of pest activity in D. suzukii, indicating important cues used by this species when finding a host, and the molecular mechanisms potentially underlying their adaptation to a new ecological niche.

Transcription Factors Drive Opposite Relationships between Gene Age and Tissue Specificity in Male and Female Drosophila Gonads

Evolutionarily young genes are usually preferentially expressed in the testis across species. Although it is known that older genes are generally more broadly expressed than younger genes, the properties that shaped this pattern are unknown. Older genes may gain expression across other tissues uniformly, or faster in certain tissues than others. Using Drosophila gene expression data, we confirmed previous findings that younger genes are disproportionately testis biased and older genes are disproportionately ovary biased. We found that the relationship between gene age and expression is stronger in the ovary than any other tissue and weakest in testis. We performed ATAC-seq on Drosophila testis and found that although genes of all ages are more likely to have open promoter chromatin in testis than in ovary, promoter chromatin alone does not explain the ovary bias of older genes. Instead, we found that upstream transcription factor (TF) expression is highly predictive of gene expression in ovary but not in testis. In the ovary, TF expression is more predictive of gene expression than open promoter chromatin, whereas testis gene expression is similarly influenced by both TF expression and open promoter chromatin. We propose that the testis is uniquely able to express younger genes controlled by relatively few TFs, whereas older genes with more TF partners are broadly expressed with peak expression most likely in the ovary. The testis allows widespread baseline expression that is relatively unresponsive to regulatory changes, whereas the ovary transcriptome is more responsive to trans-regulation and has a higher ceiling for gene expression.

The evolution of sex-biased gene expression in the Drosophila brain

Genes with sex-biased expression in Drosophila are thought to underlie sexually dimorphic phenotypes and have been shown to possess unique evolutionary properties. However, the forces and constraints governing the evolution of sex-biased genes in the somatic tissues of Drosophila are largely unknown. By using population-scale RNA sequencing data, we show that sex-biased genes in the Drosophila brain are highly enriched on the X Chromosome and that most are biased in a species-specific manner. We show that X-linked male-biased genes, and to a lesser extent female-biased genes, are enriched for signatures of directional selection at the gene expression level. By examining the evolutionary properties of gene-flanking regions on the X Chromosome, we find evidence that adaptive cis-regulatory changes are more likely to drive the expression evolution of X-linked male-biased genes than other X-linked genes. Finally, we examine whether constraint owing to broad expression across multiple tissues and genetic constraint owing to the largely shared male and female genomes could be responsible for the observed patterns of gene expression evolution. We find that expression breadth does not constrain the directional evolution of gene expression in the brain. Additionally, we find that the shared genome between males and females imposes a substantial constraint on the expression evolution of sex-biased genes. Overall, these results significantly advance our understanding of the patterns and forces shaping the evolution of sexual dimorphism in the Drosophila brain.

Testis single-cell RNA-seq reveals the dynamics of de novo gene transcription and germline mutational bias in Drosophila

The testis is a peculiar tissue in many respects. It shows patterns of rapid gene evolution and provides a hotspot for the origination of genetic novelties such as de novo genes, duplications and mutations. To investigate the expression patterns of genetic novelties across cell types, we performed single-cell RNA-sequencing of adult Drosophila testis. We found that new genes were expressed in various cell types, the patterns of which may be influenced by their mode of origination. In particular, lineage-specific de novo genes are commonly expressed in early spermatocytes, while young duplicated genes are often bimodally expressed. Analysis of germline substitutions suggests that spermatogenesis is a highly reparative process, with the mutational load of germ cells decreasing as spermatogenesis progresses. By elucidating the distribution of genetic novelties across spermatogenesis, this study provides a deeper understanding of how the testis maintains its core reproductive function while being a hotbed of evolutionary innovation.

Seasonal cues induce phenotypic plasticity of Drosophila suzukii to enhance winter survival

BACKGROUND

As global climate change and exponential human population growth intensifies pressure on agricultural systems, the need to effectively manage invasive insect pests is becoming increasingly important to global food security. Drosophila suzukii is an invasive pest that drastically expanded its global range in a very short time since 2008, spreading to most areas in North America and many countries in Europe and South America. Preliminary ecological modeling predicted a more restricted distribution and, for this reason, the invasion of D. suzukii to northern temperate regions is especially unexpected. Investigating D. suzukii phenology and seasonal adaptations can lead to a better understanding of the mechanisms through which insects express phenotypic plasticity, which likely enables invasive species to successfully colonize a wide range of environments.

RESULTS

We describe seasonal phenotypic plasticity in field populations of D. suzukii. Specifically, we observed a trend of higher proportions of flies with the winter morph phenotype, characterized by darker pigmentation and longer wing length, as summer progresses to winter. A laboratory-simulated winter photoperiod and temperature (12:12 L:D and 10 °C) were sufficient to induce the winter morph phenotype in D. suzukii. This winter morph is associated with increased survival at 1 °C when compared to the summer morph, thus explaining the ability of D. suzukii to survive cold winters. We then used RNA sequencing to identify gene expression differences underlying seasonal differences in D. suzukii physiology. Winter morph gene expression is consistent with known mechanisms of cold-hardening such as adjustments to ion transport and up-regulation of carbohydrate metabolism. In addition, transcripts involved in oogenesis and DNA replication were down-regulated in the winter morph, providing the first molecular evidence of a reproductive diapause in D. suzukii.

CONCLUSIONS

To date, D. suzukii cold resistance studies suggest that this species cannot overwinter in northern locations, e.g. Canada, even though they are established pests in these regions. Combining physiological investigations with RNA sequencing, we present potential mechanisms by which D. suzukii can overwinter in these regions. This work may contribute to more accurate population models that incorporate seasonal variation in physiological parameters, leading to development of better management strategies.

The Adaptive Significance of Natural Genetic Variation in the DNA Damage Response of Drosophila melanogaster

Despite decades of work, our understanding of the distribution of fitness effects of segregating genetic variants in natural populations remains largely incomplete. One form of selection that can maintain genetic variation is spatially varying selection, such as that leading to latitudinal clines. While the introduction of population genomic approaches to understanding spatially varying selection has generated much excitement, little successful effort has been devoted to moving beyond genome scans for selection to experimental analysis of the relevant biology and the development of experimentally motivated hypotheses regarding the agents of selection; it remains an interesting question as to whether the vast majority of population genomic work will lead to satisfying biological insights. Here, motivated by population genomic results, we investigate how spatially varying selection in the genetic model system, Drosophila melanogaster, has led to genetic differences between populations in several components of the DNA damage response. UVB incidence, which is negatively correlated with latitude, is an important agent of DNA damage. We show that sensitivity of early embryos to UVB exposure is strongly correlated with latitude such that low latitude populations show much lower sensitivity to UVB. We then show that lines with lower embryo UVB sensitivity also exhibit increased capacity for repair of damaged sperm DNA by the oocyte. A comparison of the early embryo transcriptome in high and low latitude embryos provides evidence that one mechanism of adaptive DNA repair differences between populations is the greater abundance of DNA repair transcripts in the eggs of low latitude females. Finally, we use population genomic comparisons of high and low latitude samples to reveal evidence that multiple components of the DNA damage response and both coding and non-coding variation likely contribute to adaptive differences in DNA repair between populations.

Parallel Gene Expression Differences between Low and High Latitude Populations of Drosophila melanogaster and D. simulans

Gene expression variation within species is relatively common, however, the role of natural selection in the maintenance of this variation is poorly understood. Here we investigate low and high latitude populations of Drosophila melanogaster and its sister species, D. simulans, to determine whether the two species show similar patterns of population differentiation, consistent with a role for spatially varying selection in maintaining gene expression variation. We compared at two temperatures the whole male transcriptome of D. melanogaster and D. simulans sampled from Panama City (Panama) and Maine (USA). We observed a significant excess of genes exhibiting differential expression in both species, consistent with parallel adaptation to heterogeneous environments. Moreover, the majority of genes showing parallel expression differentiation showed the same direction of differential expression in the two species and the magnitudes of expression differences between high and low latitude populations were correlated across species, further bolstering the conclusion that parallelism for expression phenotypes results from spatially varying selection. However, the species also exhibited important differences in expression phenotypes. For example, the genomic extent of genotype × environment interaction was much more common in D. melanogaster. Highly differentiated SNPs between low and high latitudes were enriched in the 3’ UTRs and CDS of the geographically differently expressed genes in both species, consistent with an important role for cis-acting variants in driving local adaptation for expression-related phenotypes.

While gene expression variation in natural populations is common, the population genetic processes responsible for the maintenance of this variation remain obscure. Here we study geographic differences in gene expression in recently established low and high latitude populations of two closely related species of Drosophila. We observe substantial parallelism in expression differences and expression plasticity between populations, which supports the idea that spatially varying selection correlated with latitude contributes to the maintenance of gene expression variation in these species. Comparison of inter-population sequence differentiation and expression differentiation suggests that cis-acting variants play a role in geographic expression differentiation.

Evidence that natural selection maintains genetic variation for sleep in Drosophila melanogaster

Background

Drosophila melanogaster often shows correlations between latitude and phenotypic or genetic variation on different continents, which suggests local adaptation with respect to a heterogeneous environment. Previous phenotypic analyses of latitudinal clines have investigated mainly physiological, morphological, or life-history traits. Here, we studied latitudinal variation in sleep in D. melanogaster populations from North and Central America. In parallel, we used RNA-seq to identify interpopulation gene expression differences.

Results

We found that in D. melanogaster the average nighttime sleep bout duration exhibits a latitudinal cline such that sleep bouts of equatorial populations are roughly twice as long as those of temperate populations. Interestingly, this pattern of latitudinal variation is not observed for any daytime measure of activity or sleep. We also found evidence for geographic variation for sunrise anticipation. Our RNA-seq experiment carried out on heads from a low and high latitude population identified a large number of gene expression differences, most of which were time dependent. Differentially expressed genes were enriched in circadian regulated genes and enriched in genes potentially under spatially varying selection.

Conclusion

Our results are consistent with a mechanistic and selective decoupling of nighttime and daytime activity. Furthermore, the present study suggests that natural selection plays a major role in generating transcriptomic variation associated with circadian behaviors. Finally, we identified genomic variants plausibly causally associated with the observed behavioral and transcriptomic variation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0316-2) contains supplementary material, which is available to authorized users.

Gene-environment interplay in Drosophila melanogaster: chronic food deprivation in early life affects adult exploratory and fitness traits

Early life adversity has known impacts on adult health and behavior, yet little is known about the gene-environment interactions (GEIs) that underlie these consequences. We used the fruit fly Drosophila melanogaster to show that chronic early nutritional adversity interacts with rover and sitter allelic variants of foraging (for) to affect adult exploratory behavior, a phenotype that is critical for foraging, and reproductive fitness. Chronic nutritional adversity during adulthood did not affect rover or sitter adult exploratory behavior; however, early nutritional adversity in the larval period increased sitter but not rover adult exploratory behavior. Increasing for gene expression in the mushroom bodies, an important center of integration in the fly brain, changed the amount of exploratory behavior exhibited by sitter adults when they did not experience early nutritional adversity but had no effect in sitters that experienced early nutritional adversity. Manipulation of the larval nutritional environment also affected adult reproductive output of sitters but not rovers, indicating GEIs on fitness itself. The natural for variants are an excellent model to examine how GEIs underlie the biological embedding of early experience.

Identification of X-linked quantitative trait loci affecting cold tolerance in Drosophila melanogaster and fine mapping by selective sweep analysis

Drosophila melanogaster is a cosmopolitan species that colonizes a great variety of environments. One trait that shows abundant evidence for naturally segregating genetic variance in different populations of D. melanogaster is cold tolerance. Previous work has found quantitative trait loci (QTL) exclusively on the second and the third chromosomes. To gain insight into the genetic architecture of cold tolerance on the X chromosome and to compare the results with our analyses of selective sweeps, a mapping population was derived from a cross between substitution lines that solely differed in the origin of their X chromosome: one originates from a European inbred line and the other one from an African inbred line. We found a total of six QTL for cold tolerance factors on the X chromosome of D. melanogaster. Although the composite interval mapping revealed slightly different QTL profiles between sexes, a coherent model suggests that most QTL overlapped between sexes, and each explained around 5-14% of the genetic variance (which may be slightly overestimated). The allelic effects were largely additive, but we also detected two significant interactions. Taken together, this provides evidence for multiple QTL that are spread along the entire X chromosome and whose effects range from low to intermediate. One detected transgressive QTL influences cold tolerance in different ways for the two sexes. While females benefit from the European allele increasing their cold tolerance, males tend to do better with the African allele. Finally, using selective sweep mapping, the candidate gene CG16700 for cold tolerance colocalizing with a QTL was identified.

Effect of laboratory acclimation on the variation of reproduction-related characters in Drosophila melanogaster

The natural variation of sex-specific characters between populations can favor their behavioral isolation, eventually leading to the formation of new species. Marked variations for male courtship, mating and the production of sex pheromones – three complex characters potentially inducing sexual isolation – were found between Drosophila melanogaster populations of various origins acclimated for many generations in research laboratories. However, the natural variation of these three characters between natural populations and their evolution after long-term acclimation in the laboratory remains unknown. We measured many traits involved in these characters in six stocks initiated with distinct populations sampled in a restricted geographic area. Several sex-specific traits varied between stocks freshly brought back to the laboratory. After 100 generations spent in the laboratory without any experimental selection, traits varied in a strain-dependent manner. This variation was not related to a reduction of their variance except for copulation duration. This indicates that reproduction-related characters can diverge between neighboring D. melanogaster populations, and differently adapt to stable laboratory conditions.

Effect of genes, social experience, and their interaction on the courtship behaviour of transgenic Drosophila males

Behaviour depends (a) on genes that specify the neural and non-neural elements involved in the perception of and responses to sensory stimuli and (b) on experience that can modulate the fine development of these elements. We exposed transgenic and control Drosophila melanogaster males, and their hybrids, to male siblings during adult development and measured the contribution of genes and of experience to their courtship behaviour. The transgene CheB42a specifically targets male gustatory sensillae and alters the perception of male inhibitory pheromones which leads to frequent male-male interactions. The age at which social experience occurred and the genotype of tester males induced a variable effect on the intensity of male homo- and heterosexual courtship. The strong interaction between the effects of genes and of social experience reveals the plasticity of the apparently stereotyped elements involved in male courtship behaviour. Finally, a high intensity of homosexual courtship was found only in males that simultaneously carried a mutation in their white gene and the CheB42a transgene.

Social experience and pheromonal perception can change male-male interactions in Drosophila melanogaster

Social interaction with conspecifics can influence the developing brain and behaviour of the exposed animal. This experience can involve the exchange and retention of visual, chemical, acoustic and tactile signals. When several Drosophila melanogaster male flies are associated with mated females in the presence of food, they show frequent aggressive interactions. To measure the role of social experience on male-male interaction, two tester males - naive or exposed to sibling(s) during a variable period of their adult development - were confronted in the absence of female and food. The two males displayed homosexual courtship and aggressive behaviours, the frequency, intensity and directionality of which varied according to their experience. The effect of social experience was greatly enhanced between transgenic males partially defective for pheromonal perception, indicating that male inhibitory pheromones are normally used to repress male-male interaction.