Molecular signatures of alternative reproductive strategies in a facultatively social hover wasp.

Social insect reproductives and non-reproductives represent ideal models with which to understand the expression and regulation of alternative phenotypes. Most research in this area has focused on the developmental regulation of reproductive phenotypes in obligately social taxa such as honey bees, while relatively few studies have addressed the molecular correlates of reproductive differentiation in species in which the division of reproductive labour is established only in plastic dominance hierarchies. To address this knowledge gap, we generate the first genome for any stenogastrine wasp and analyse brain transcriptomic data for non-reproductives and reproductives of the facultatively social species Liostenogaster flavolineata, a representative of one of the simplest forms of social living. By experimentally manipulating the reproductive 'queues' exhibited by social colonies of this species, we show that reproductive division of labour in this species is associated with transcriptomic signatures that are more subtle and variable than those observed in social taxa in which colony living has become obligate; that variation in gene expression among non-reproductives reflects their investment into foraging effort more than their social rank; and that genes associated with reproductive division of labour overlap to some extent with those underlying division of labour in the separate polistine origin of wasp sociality but only explain a small portion of overall variation in this trait. These results indicate that broad patterns of within-colony transcriptomic differentiation in this species are similar to those in Polistinae but offer little support for the existence of a strongly conserved 'toolkit' for sociality.

A complete absence of indirect genetic effects on brain gene expression in a highly social context.

Genes not only control traits of their carrier organism (known as direct genetic effects or DGEs) but also shape their carrier's physical environment and the phenotypes of their carrier's social partners (known as indirect genetic effects or IGEs). Theoretical research has shown that the effects that genes exert on social partners can have profound consequences, potentially altering heritability and the direction of trait evolution. Complementary empirical research has shown that in various contexts (particularly in animal agriculture) IGEs can explain a large proportion of variation in specific traits. However, little is known about the general prevalence of IGEs. We conducted a reciprocal cross-fostering experiment with two genetic lineages of the clonal raider ant Ooceraea biroi to quantify the relative contribution of DGEs and IGEs to variation in brain gene expression (which underlies behavioural variation). We found that thousands of genes are differentially expressed by DGEs but not a single gene is differentially expressed by IGEs. This is surprising given the highly social context of ant colonies and given that individual behaviour varies according to the genotypic composition of the social environment in O. biroi. Overall, these findings indicate that we have a lot to learn about how the magnitude of IGEs varies across species and contexts.

Genome-wide sexually antagonistic variants reveal long-standing constraints on sexual dimorphism in fruit flies.

The evolution of sexual dimorphism is constrained by a shared genome, leading to 'sexual antagonism', in which different alleles at given loci are favoured by selection in males and females. Despite its wide taxonomic incidence, we know little about the identity, genomic location, and evolutionary dynamics of antagonistic genetic variants. To address these deficits, we use sex-specific fitness data from 202 fully sequenced hemiclonal Drosophila melanogaster fly lines to perform a genome-wide association study (GWAS) of sexual antagonism. We identify approximately 230 chromosomal clusters of candidate antagonistic single nucleotide polymorphisms (SNPs). In contradiction to classic theory, we find no clear evidence that the X chromosome is a hot spot for sexually antagonistic variation. Characterising antagonistic SNPs functionally, we find a large excess of missense variants but little enrichment in terms of gene function. We also assess the evolutionary persistence of antagonistic variants by examining extant polymorphism in wild D. melanogaster populations and closely related species. Remarkably, antagonistic variants are associated with multiple signatures of balancing selection across the D. melanogaster distribution range and in their sister species D. simulans, indicating widespread and evolutionarily persistent (about 1 million years) genomic constraints on the evolution of sexual dimorphism. Based on our results, we propose that antagonistic variation accumulates because of constraints on the resolution of sexual conflict over protein coding sequences, thus contributing to the long-term maintenance of heritable fitness variation.

A non-coding indel polymorphism in the fruitless gene of Drosophila melanogaster exhibits antagonistically pleiotropic fitness effects.

The amount of genetic variation for fitness within populations tends to exceed that expected under mutation-selection-drift balance. Several mechanisms have been proposed to actively maintain polymorphism and account for this discrepancy, including antagonistic pleiotropy (AP), where allelic variants have opposing effects on different components of fitness. Here, we identify a non-coding indel polymorphism in the fruitless gene of Drosophila melanogaster and measure survival and reproductive components of fitness in males and females of replicate lines carrying each respective allele. Expressing the fruitless region in a hemizygous state reveals a pattern of AP, with one allele generating greater reproductive fitness and the other conferring greater survival to adulthood. Different fitness effects were observed in an alternative genetic background, which may reflect dominance reversal and/or epistasis. Our findings link sequence-level variation at a single locus with complex effects on a range of fitness components, thus helping to explain the maintenance of genetic variation for fitness. Transcription factors, such as fruitless, may be prime candidates for targets of balancing selection since they interact with multiple target loci and their associated phenotypic effects.

Evolving Plastic Responses to External and Genetic Environments.

Phenotypic plasticity can mitigate adaptive trade-offs in fluctuating environments but how plasticity arises is little known. New research documents this process in a bacterial system. We highlight remarkable parallels to the evolution of sexual dimorphism and argue that their approach can aid our understanding of adaptive conflicts between the sexes.

Sexual antagonism drives the displacement of polymorphism across gene regulatory cascades.

Males and females have different reproductive roles and are often subject to contrasting selection pressures. This sexual antagonism can lead, at a given locus, to different alleles being favoured in each sex and, consequently, to genetic variation being maintained in a population. Although the presence of sexually antagonistic (SA) polymorphisms has been documented across a range of species, their evolutionary dynamics remain poorly understood. Here, we study SA selection on gene expression, which is fundamental to sexual dimorphism, via the evolution of regulatory binding sites. We show that for sites longer than 1 nucleotide, expression polymorphism is maintained only when intermediate expression levels are deleterious to both sexes. We then show that, in a regulatory cascade, expression polymorphism tends to become displaced over evolutionary time from the target of SA selection to upstream regulators. Our results have consequences for understanding the evolution of sexual dimorphism, and provide specific empirical predictions for the regulatory architecture of genes under SA selection.

Sex and genotype effects on nutrient-dependent fitness landscapes in Drosophila melanogaster.

The sexes perform different reproductive roles and have evolved sometimes strikingly different phenotypes. One focal point of adaptive divergence occurs in the context of diet and metabolism, and males and females of a range of species have been shown to require different nutrients to maximize their fitness. Biochemical analyses in Drosophila melanogaster have confirmed that dimorphism in dietary requirements is associated with molecular sex differences in metabolite titres. In addition, they also showed significant within-sex genetic variation in the metabolome. To date however, it is unknown whether this metabolic variation translates into differences in reproductive fitness. The answer to this question is crucial to establish whether genetic variation is selectively neutral or indicative of constraints on sex-specific physiological adaptation and optimization. Here we assay genetic variation in consumption and metabolic fitness effects by screening male and female fitness of thirty D. melanogaster genotypes across four protein-to-carbohydrate ratios. In addition to confirming sexual dimorphism in consumption and fitness, we find significant genetic variation in male and female dietary requirements. Importantly, these differences are not explained by feeding responses and probably reflect metabolic variation that, in turn, suggests the presence of genetic constraints on metabolic dimorphism.

Under-dominance constrains the evolution of negative autoregulation in diploids.

Regulatory networks have evolved to allow gene expression to rapidly track changes in the environment as well as to buffer perturbations and maintain cellular homeostasis in the absence of change. Theoretical work and empirical investigation in Escherichia coli have shown that negative autoregulation confers both rapid response times and reduced intrinsic noise, which is reflected in the fact that almost half of Escherichia coli transcription factors are negatively autoregulated. However, negative autoregulation is rare amongst the transcription factors of Saccharomyces cerevisiae. This difference is surprising because E. coli and S. cerevisiae otherwise have similar profiles of network motifs. In this study we investigate regulatory interactions amongst the transcription factors of Drosophila melanogaster and humans, and show that they have a similar dearth of negative autoregulation to that seen in S. cerevisiae. We then present a model demonstrating that this striking difference in the noise reduction strategies used amongst species can be explained by constraints on the evolution of negative autoregulation in diploids. We show that regulatory interactions between pairs of homologous genes within the same cell can lead to under-dominance--mutations which result in stronger autoregulation, and decrease noise in homozygotes, paradoxically can cause increased noise in heterozygotes. This severely limits a diploid's ability to evolve negative autoregulation as a noise reduction mechanism. Our work offers a simple and general explanation for a previously unexplained difference between the regulatory architectures of E. coli and yeast, Drosophila and humans. It also demonstrates that the effects of diploidy in gene networks can have counter-intuitive consequences that may profoundly influence the course of evolution.

Locally adaptive inversions in structured populations.

Inversions have been proposed to facilitate local adaptation, by linking together locally coadapted alleles at different loci. Prior work addressing this question theoretically has considered the spread of inversions in "continent-island" scenarios in which there is a unidirectional flow of maladapted migrants into the island population. In this setting, inversions capturing locally adaptive haplotypes are most likely to invade when selection is weak, because stronger local selection (i) more effectively purges maladaptive alleles and (ii) generates linkage disequilibrium between adaptive alleles, thus lessening the advantage of inversions. We show this finding only holds under limited conditions by studying the establishment of inversions in a more general two-deme model, which explicitly considers the dynamics of allele frequencies in both populations linked by bidirectional migration. In this model, the level of symmetry between demes can be varied from complete asymmetry (continent-island) to complete symmetry. For symmetric selection and migration, strong selection increases the allele frequency divergence between demes thereby increasing the frequency of maladaptive alleles in migrants, favoring inversions-the opposite of the pattern seen in the asymmetric continent-island scenario. We also account for the likelihood that a new inversion captures an adaptive haplotype in the first instance. When considering the combined process of capture and invasion in "continent island" and symmetric scenarios, relatively strong selection increases inversion establishment probability. Migration must also be low enough that the inversion is likely to capture an adaptive allele combination, but not so low as to eliminate the inversion's advantage. Overall, our analysis suggests that inversions are likely to harbor larger effect alleles that experience relatively strong selection.

On the genetic architecture of cytoplasmic incompatibility: inference from phenotypic data.

Numerous insects carry intracellular bacteria that manipulate the insects' reproduction and thus facilitate their own spread. Cytoplasmic incompatibility (CI) is a common form of such manipulation, where a (currently uncharacterized) bacterial modification of male sperm induces the early death of embryos unless the fertilized eggs carry the same bacteria, inherited from the mother. The death of uninfected embryos provides an indirect selective advantage to infected ones, thus enabling the spread of the bacteria. Here we use and expand recently developed algorithms to infer the genetic architecture underlying the complex incompatibility data from the mosquito Culex pipiens. We show that CI requires more genetic determinants than previously believed and that quantitative variation in gene products potentially contributes to the observed CI patterns. In line with population genetic theory of CI, our analysis suggests that toxin factors (those inducing embryo death) are present in fewer copies in the bacterial genomes than antitoxin factors (those ensuring that infected embryos survive). In combination with comparative genomics, our approach will provide helpful guidance to identify the genetic basis of CI and more generally of other toxin/antitoxin systems that can be conceptualized under the same framework.

Molecular evolution of Drosophila Sex-lethal and related sex determining genes.

BACKGROUND: Sex determining mechanisms are evolutionarily labile and related species often use different primary signals and gene regulatory networks. This is well illustrated by the sex determining cascade of Drosophila fruitflies, which have recruited Sex-lethal as the master switch and cellular memory of sexual identity, a role performed in other insects by the gene transformer. Here we investigate the evolutionary change in the coding sequences of sex determining genes associated with the recruitment of Sex-lethal. We analyze sequences of Sex-lethal itself, its Drosophila paralogue sister-or-Sex-lethal and downstream targets transformer and doublesex. RESULTS: We find that the recruitment of sister-or-Sex-lethal was associated with a number of adaptive amino acid substitutions, followed by a tightening of purifying selection within the Drosophila clade. Sequences of the paralogue sister-or-Sex-lethal, in contrast, show a signature of rampant positive selection and relaxation of purifying selection. The recruitment of Sex-lethal as top regulator and memory gene is associated with a significant release from purifying selection in transformer throughout the Drosophila clade. In addition, doublesex shows a signature of positive selection and relaxation of purifying selection in the Drosophila clade. A similar pattern is seen in sequences from the sister Tephritidae clade. CONCLUSIONS: The pattern of molecular evolution we observe for Sex-lethal and its paralogue sister-or-Sex-lethal is not characteristic of a duplication followed by neo-functionalization. Rather, evidence suggests a sub-functionalization scenario achieved through the evolution of sophisticated splicing. As expected, we find that transformer evolves under relaxed purifying selection after the recruitment of Sex-lethal in Drosophila. Finally, the observation of doublesex adaptation in both Drosophila and Tephritidae suggests that these changes are due to ongoing adaptation of downstream sex-specific regulation, rather than being associated the recruitment of Sex-lethal and the resulting change in the topology of the sex determining cascade.

Evolutionary genetics: Dissecting a sexually antagonistic polymorphism.

Males and females experience divergent selection on many shared traits, which can lead to 'sexual antagonism' - opposing fitness effects of genetic variants in each sex. A new study in the fly Drosophila serrata links sexually antagonistic selection on cuticular hydrocarbons to a single major-effect gene.

Sex differences in deleterious mutational effects in Drosophila melanogaster: combining quantitative and population genetic insights.

Fitness effects of deleterious mutations can differ between females and males due to: (i) sex differences in the strength of purifying selection; and (ii) sex differences in ploidy. Although sex differences in fitness effects have important broader implications (e.g., for the evolution of sex and lifespan), few studies have quantified their scope. Those that have belong to one of two distinct empirical traditions: (i) quantitative genetics, which focusses on multi-locus genetic variances in each sex, but is largely agnostic about their genetic basis; and (ii) molecular population genetics, which focusses on comparing autosomal and X-linked polymorphism, but is poorly suited for inferring contemporary sex differences. Here, we combine both traditions to present a comprehensive analysis of female and male adult reproductive fitness among 202 outbred, laboratory-adapted, hemiclonal genomes of Drosophila melanogaster. While we find no clear evidence for sex differences in the strength of purifying selection, sex differences in ploidy generate multiple signals of enhanced purifying selection for X-linked loci. These signals are present in quantitative genetic metrics-i.e., a disproportionate contribution of the X to male (but not female) fitness variation-and population genetic metrics-i.e., steeper regressions of an allele's average fitness effect on its frequency, and proportionally less nonsynonymous polymorphism on the X than autosomes. Fitting our data to models for both sets of metrics, we infer that deleterious alleles are partially recessive. Given the often-large gap between quantitative and population genetic estimates of evolutionary parameters, our study showcases the benefits of combining genomic and fitness data when estimating such parameters.

Mother's curse is pervasive across a large mitonuclear Drosophila panel.

The maternal inheritance of mitochondrial genomes entails a sex-specific selective sieve, whereby mutations in mitochondrial DNA can only respond to selection acting on females. In theory, this enables male-harming mutations to accumulate in mitochondrial genomes as long as they are neutral, beneficial, or only slightly deleterious to females. Ultimately, this bias could drive the evolution of male-specific mitochondrial mutation loads, an idea known as mother's curse. Earlier work on this hypothesis has mainly used small Drosophila panels, in which naturally sourced mitochondrial genomes were coupled to an isogenic nuclear background. The lack of nuclear genetic variation in these designs has precluded robust generalization. Here, we test the predictions of mother's curse using a large Drosophila mitonuclear genetic panel, comprising nine isogenic nuclear genomes coupled to nine mitochondrial haplotypes, giving a total of 81 different mitonuclear genotypes. Following a predictive framework, we tested the mother's curse hypothesis by screening our panel for wing size. This trait is tightly correlated with overall body size and is sexually dimorphic in Drosophila. Moreover, growth is heavily reliant on metabolism and mitochondrial function, making wing size an ideal trait for the study of the impact of mitochondrial variation. We detect high levels of mitonuclear epistasis, and more importantly, we report that mitochondrial genetic variance is larger in male than female Drosophila for eight out of the nine nuclear genetic backgrounds used. These results demonstrate that the maternal inheritance of mitochondrial DNA does indeed modulate male life history traits in a more generalisable way than previously demonstrated.

The molecular basis of socially mediated phenotypic plasticity in a eusocial paper wasp.

Phenotypic plasticity, the ability to produce multiple phenotypes from a single genotype, represents an excellent model with which to examine the relationship between gene expression and phenotypes. Analyses of the molecular foundations of phenotypic plasticity are challenging, however, especially in the case of complex social phenotypes. Here we apply a machine learning approach to tackle this challenge by analyzing individual-level gene expression profiles of Polistes dominula paper wasps following the loss of a queen. We find that caste-associated gene expression profiles respond strongly to queen loss, and that this change is partly explained by attributes such as age but occurs even in individuals that appear phenotypically unaffected. These results demonstrate that large changes in gene expression may occur in the absence of outwardly detectable phenotypic changes, resulting here in a socially mediated de-differentiation of individuals at the transcriptomic level but not at the levels of ovarian development or behavior.

Male-killing bacteria trigger a cycle of increasing male fatigue and female promiscuity.

Sex-ratio distorters are found in numerous species and can reach high frequencies within populations. Here, we address the compelling, but poorly tested, hypothesis that the sex ratio bias caused by such elements profoundly alters their host's mating system. We compare aspects of female and male reproductive biology between island populations of the butterfly Hypolimnas bolina that show varying degrees of female bias, because of a male-killing Wolbachia infection. Contrary to expectation, female bias leads to an increase in female mating frequency, up to a point where male mating capacity becomes limiting. We show that increased female mating frequency can be explained as a facultative response to the depleted male mating resources in female biased populations. In other words, this system is one where male-killing bacteria trigger a vicious circle of increasing male fatigue and female promiscuity.

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin.

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that-in addition to recombination rate-the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.

Impact of mitonuclear interactions on life-history responses to diet.

Mitochondria are central to both energy metabolism and biosynthesis. Mitochondrial function could therefore influence resource allocation. Critically, mitochondrial function depends on interactions between proteins encoded by the mitochondrial and nuclear genomes. Severe incompatibilities between these genomes can have pervasive effects on both fitness and longevity. How milder deficits in mitochondrial function affect life-history trade-offs is less well understood. Here, we analyse how mitonuclear interactions affect the trade-off between fecundity and longevity in Drosophila melanogaster. We consider a panel of 10 different mitochondrial DNA haplotypes against two contrasting nuclear backgrounds (w1118 (WE) and Zim53 (ZIM)) in response to high-protein versus standard diet. We report strikingly different responses between the two nuclear backgrounds. WE females have higher fecundity and decreased longevity on high protein. ZIM females have much greater fecundity and shorter lifespan than WE flies on standard diet. High protein doubled their fecundity with no effect on longevity. Mitochondrial haplotype reflected nuclear life-history trade-offs, with a negative correlation between longevity and fecundity in WE flies and no correlation in ZIM flies. Mitonuclear interactions had substantial effects but did not reflect genetic distance between mitochondrial haplotypes. We conclude that mitonuclear interactions can have significant impact on life-history trade-offs, but their effects are not predictable by relatedness. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

A test of the null model for 5' UTR evolution based on GC content.

Eukaryotic mRNAs are headed by a stretch of noncoding sequence, the 5' untranslated region (UTR). It has been proposed that the length of 5' UTRs is selectively neutral and evolves under a process of stochastic destruction and recruitment of core promoter elements, combined with selection against the premature initiation of translation. We test this null model by investigating whether 5' UTR length varies with genomic GC content, an implicit prediction of the model. Using simulations, we show that the null model predicts a positive relationship between GC content and UTR length for genes regulated by a TATA box. Although this prediction is borne out qualitatively in genomic data from yeast, fruit flies, and humans, we find marked quantitative discrepancies. We conclude that UTR length may be shaped to some degree by the forces considered in the null model but that the model fails to provide a complete explanation for UTR length evolution.