Rapid evolution of a geographic cline in size in an introduced fly

Clinal patterns of desiccation and starvation resistance in ancestral and invading populations of Drosophila subobscura

As invading species expand, they eventually encounter physical and biotic stressors that limit their spread. We examine latitudinal and climatic variation in physiological tolerance in one native and two invading populations of Drosophila subobscura. These flies are native to the Palearctic region, but invaded both South and North America around 1980 and spread rapidly across 15° of latitude on each continent. Invading flies rapidly evolved latitudinal clines in chromosome inversion frequencies and in wing size that parallel those of native populations in the Old World. Here we investigate whether flies on all three continents have evolved parallel clines in desiccation and starvation tolerance, such that flies in low-latitude regions (hot, dry) might have increased stress resistance. Starvation tolerance does not vary with latitude or climate on any continent. In contrast, desiccation tolerance varies clinally with latitude on all three continents, although not in parallel. In North American and Europe, desiccation tolerance is inversely related to latitude, as expected. But in South America, desiccation tolerance increases with latitude and is greatest in relatively cool and wet areas. Differences among continents in latitudinal patterns of interspecific-competition potentially influence clinal selection for physiological resistance, but no simple pattern is evident on these continents.

The global expansion of a single ant supercolony

Ants are among the most damaging invasive species, and their success frequently arises from the widespread cooperation displayed by introduced populations, often across hundreds of kilometers. Previous studies of the invasive Argentine ant (Linepithema humile) have shown that introduced populations on different continents each contain a single, vast supercolony and, occasionally, smaller secondary colonies. Here, we perform inter-continental behavioral analyses among supercolonies in North America, Europe, Asia, Hawaii, New Zealand and Australia and show that these far-flung supercolonies also recognize and accept each other as if members of a single, globally distributed supercolony. Furthermore, populations also possess similar genetic and chemical profiles. However, these ants do show aggression toward ants from South Africa and the smaller secondary colonies that occur in Hawaii and California. Thus, the largest and most dominant introduced populations are likely descended from the same ancestral colony and, despite having been established more than 100 years ago, have diverged very little. This apparent evolutionary stasis is surprising because, in other species, some of the most rapid rates of evolutionary change have occurred in introduced populations. Given the spatial extent of the Argentine ant society we report here, there can be little doubt that this intercontinental supercolony represents the most populous known animal society.

Within-population genetic effects of mtDNA on metabolic rate in Drosophila subobscura

A growing body of research supports the view that within-species sequence variation in the mitochondrial genome (mtDNA) is functional, in the sense that it has important phenotypic effects. However, most of this empirical foundation is based on comparisons across populations, and few studies have addressed the functional significance of mtDNA polymorphism within populations. Here, using mitonuclear introgression lines, we assess differences in whole-organism metabolic rate of adult Drosophila subobscura fruit flies carrying either of three different sympatric mtDNA haplotypes. We document sizeable, up to 20%, differences in metabolic rate across these mtDNA haplotypes. Further, these mtDNA effects are to some extent sex specific. We found no significant nuclear or mitonuclear genetic effects on metabolic rate, consistent with a low degree of linkage disequilibrium between mitochondrial and nuclear genes within populations. The fact that mtDNA haplotype variation within a natural population affects metabolic rate, which is a key physiological trait with important effects on life-history traits, adds weight to the emergent view that mtDNA haplotype variation is under natural selection and it revitalizes the question as to what processes act to maintain functional mtDNA polymorphism within populations.

Biological invasions, climate change and genomics

The rate of biological invasions is expected to increase as the effects of climate change on biological communities become widespread. Climate change enhances habitat disturbance which facilitates the establishment of invasive species, which in turn provides opportunities for hybridization and introgression. These effects influence local biodiversity that can be tracked through genetic and genomic approaches. Metabarcoding and metagenomic approaches provide a way of monitoring some types of communities under climate change for the appearance of invasives. Introgression and hybridization can be followed by the analysis of entire genomes so that rapidly changing areas of the genome are identified and instances of genetic pollution monitored. Genomic markers enable accurate tracking of invasive species' geographic origin well beyond what was previously possible. New genomic tools are promoting fresh insights into classic questions about invading organisms under climate change, such as the role of genetic variation, local adaptation and climate pre-adaptation in successful invasions. These tools are providing managers with often more effective means to identify potential threats, improve surveillance and assess impacts on communities. We provide a framework for the application of genomic techniques within a management context and also indicate some important limitations in what can be achieved.

Overcoming Allee effects through evolutionary, genetic, and demographic rescue

Despite the amplified threats of extinction facing small founder populations, successful colonization sometimes occurs, bringing devastating ecological and economic consequences. One explanation may be rapid evolution, which can increase mean fitness in populations declining towards extinction, permitting persistence and subsequent expansion. Such evolutionary rescue may be particularly important, given Allee effects. When a population is introduced at low density, individuals often experience a reduction in one or more components of fitness due to novel selection pressures that arise from diminished intraspecific interactions and positive density dependence (i.e. component Allee effects). A population can avoid extinction if it can adapt and recover on its own (i.e. evolutionary rescue), or if additional immigration sustains the population (i.e. demographic rescue) or boosts its genetic variation that facilitates adaptation (i.e. genetic rescue). These various forms of rescue have often been invoked as possible mechanisms for specific invasions, but their relative importance to invasion is not generally understood. Within a spatially explicit modelling framework, we consider the relative impact of each type of rescue on the probability of successful colonization, when there is evolution of a multi-locus quantitative trait that influences the strength of component Allee effects. We demonstrate that when Allee effects are important, the effect of demographic rescue via recurrent immigration overall provides the greatest opportunity for success. While highlighting the role of evolution in the invasion process, we underscore the importance of the ecological context influencing the persistence of small founder populations.

A highly pleiotropic amino acid polymorphism in the Drosophila insulin receptor contributes to life-history adaptation

Finding the specific nucleotides that underlie adaptive variation is a major goal in evolutionary biology, but polygenic traits pose a challenge because the complex genotype-phenotype relationship can obscure the effects of individual alleles. However, natural selection working in large wild populations can shift allele frequencies and indicate functional regions of the genome. Previously, we showed that the two most common alleles of a complex amino acid insertion-deletion polymorphism in the Drosophila insulin receptor show independent, parallel clines in frequency across the North American and Australian continents. Here, we report that the cline is stable over at least a five-year period and that the polymorphism also demonstrates temporal shifts in allele frequency concurrent with seasonal change. We tested the alleles for effects on levels of insulin signaling, fecundity, development time, body size, stress tolerance, and life span. We find that the alleles are associated with predictable differences in these traits, consistent with patterns of Drosophila life-history variation across geography that likely reflect adaptation to the heterogeneous climatic environment. These results implicate insulin signaling as a major mediator of life-history adaptation in Drosophila, and suggest that life-history trade-offs can be explained by extensive pleiotropy at a single locus.

Geometric morphometric analysis of Colombian Anopheles albimanus (Diptera: Culicidae) reveals significant effect of environmental factors on wing traits and presence of a metapopulation

Anopheles albimanus is a major malaria mosquito vector in Colombia. In the present study, wing variability (size and shape) in An. albimanus populations from Colombian Maracaibo and Chocó bio-geographical eco-regions and the relationship of these phenotypic traits with environmental factors were evaluated. Microsatellite and morphometric data facilitated a comparison of the genetic and phenetic structure of this species. Wing size was influenced by elevation and relative humidity, whereas wing shape was affected by these two variables and also by rainfall, latitude, temperature and eco-region. Significant differences in mean shape between populations and eco-regions were detected, but they were smaller than those at the intra-population level. Correct assignment based on wing shape was low at the population level (<58%) and only slightly higher (>70%) at the eco-regional level, supporting the low population structure inferred from microsatellite data. Wing size was similar among populations with no significant differences between eco-regions. Population relationships in the genetic tree did not agree with those from the morphometric data; however, both datasets consistently reinforced a panmictic population of An. albimanus. Overall, site-specific population differentiation is not strongly supported by wing traits or genotypic data. We hypothesize that the metapopulation structure of An. albimanus throughout these Colombian eco-regions is favoring plasticity in wing traits, a relevant characteristic of species living under variable environmental conditions and colonizing new habitats.

Laboratory selection quickly erases historical differentiation

The roles of history, chance and selection have long been debated in evolutionary biology. Though uniform selection is expected to lead to convergent evolution between populations, contrasting histories and chance events might prevent them from attaining the same adaptive state, rendering evolution somewhat unpredictable. The predictability of evolution has been supported by several studies documenting repeatable adaptive radiations and convergence in both nature and laboratory. However, other studies suggest divergence among populations adapting to the same environment. Despite the relevance of this issue, empirical data is lacking for real-time adaptation of sexual populations with deeply divergent histories and ample standing genetic variation across fitness-related traits. Here we analyse the real-time evolutionary dynamics of Drosophila subobscura populations, previously differentiated along the European cline, when colonizing a new common environment. By analysing several life-history, physiological and morphological traits, we show that populations quickly converge to the same adaptive state through different evolutionary paths. In contrast with other studies, all analysed traits fully converged regardless of their association with fitness. Selection was able to erase the signature of history in highly differentiated populations after just a short number of generations, leading to consistent patterns of convergent evolution.

Aridity shapes cyanogenesis cline evolution in white clover (Trifolium repens L.)

Adaptive differentiation between populations is often proposed to be the product of multiple interacting selective pressures, although empirical support for this is scarce. In white clover, populations show adaptive differentiation in frequencies of cyanogenesis, the ability to produce hydrogen cyanide after tissue damage. This polymorphism arises through independently segregating polymorphisms for the presence/absence of two required cyanogenic components, cyanogenic glucosides and their hydrolysing enzyme. White clover populations worldwide have evolved a series of recurrent, climate-associated clines, with higher frequencies of cyanogenic plants in warmer locations. These clines have traditionally been hypothesized to reflect a fitness trade-off between chemical defence in herbivore-rich areas (warmer climates) and energetic costs of producing cyanogenic components in areas of low herbivore pressure (cooler climates). Recent observational studies suggest that cyanogenic components may also be beneficial in water-stressed environments. We investigated fitness trade-offs associated with temperature-induced water stress in the cyanogenesis system using manipulative experiments in growth chambers and population surveys across a longitudinal precipitation gradient in the central United States. We find that plants producing cyanogenic glucosides have higher relative fitness in treatments simulating a moderate, persistent drought stress. In water-neutral treatments, there are energetic costs to producing cyanogenic components, but only in treatments with nutrient stress. These fitness trade-offs are consistent with cyanogenesis frequencies in natural populations, where we find clinal variation in the proportion of plants producing cyanogenic glucosides along the precipitation gradient. These results suggest that multiple selective pressures interact to maintain this adaptive polymorphism and that modelling adaptation will require knowledge of environment-specific fitness effects.

On the stability of the Bayenv method in assessing human SNP-environment associations

Background

Phenotypic variation along environmental gradients has been documented among and within many species, and in some cases, genetic variation has been shown to be associated with these gradients. Bayenv is a relatively new method developed to detect patterns of polymorphisms associated with environmental gradients. Using a Bayesian Markov Chain Monte Carlo (MCMC) approach, Bayenv evaluates whether a linear model relating population allele frequencies to environmental variables is more probable than a null model based on observed frequencies of neutral markers. Although this method has been used to detect environmental adaptation in a number of species, including humans, plants, fish, and mosquitoes, stability between independent runs of this MCMC algorithm has not been characterized. In this paper, we explore the variability of results between runs and the factors contributing to it.

Results

Independent runs of the Bayenv program were carried out using genome-wide single-nucleotide polymorphism (SNP) data from samples from 60 worldwide human populations following previous applications of the Bayenv method. To assess factors contributing to the method's stability, we used varying numbers of MCMC iterations and also analyzed a second modified data set that excluded two Siberian populations with extreme climate variables. Between any two runs, correlations between Bayes factors and the overlap of SNPs in the empirical p value tails were surprisingly low. Enrichments of genic versus non-genic SNPs in the empirical tails were more robust than the empirical p values; however, the significance of the enrichments for some environmental variables still varied among runs, contradicting previously published conclusions. Runs with a greater number of MCMC iterations slightly reduced run-to-run variability, and excluding the Siberian populations did not have a large effect on the stability of the runs.

Conclusions

Because of high run-to-run variability, we advise against making conclusions about genome-wide patterns of adaptation based on only one run of the Bayenv algorithm and recommend caution in interpreting previous studies that have used only one run. Moving forward, we suggest carrying out multiple independent runs of Bayenv and averaging Bayes factors between runs to produce more stable and reliable results. With these modifications, future discoveries of environmental adaptation within species using the Bayenv method will be more accurate, interpretable, and easily compared between studies.

Thermodynamics constrains allometric scaling of optimal development time in insects

Development time is a critical life-history trait that has profound effects on organism fitness and on population growth rates. For ectotherms, development time is strongly influenced by temperature and is predicted to scale with body mass to the quarter power based on 1) the ontogenetic growth model of the metabolic theory of ecology which describes a bioenergetic balance between tissue maintenance and growth given the scaling relationship between metabolism and body size, and 2) numerous studies, primarily of vertebrate endotherms, that largely support this prediction. However, few studies have investigated the allometry of development time among invertebrates, including insects. Abundant data on development of diverse insects provides an ideal opportunity to better understand the scaling of development time in this ecologically and economically important group. Insects develop more quickly at warmer temperatures until reaching a minimum development time at some optimal temperature, after which development slows. We evaluated the allometry of insect development time by compiling estimates of minimum development time and optimal developmental temperature for 361 insect species from 16 orders with body mass varying over nearly 6 orders of magnitude. Allometric scaling exponents varied with the statistical approach: standardized major axis regression supported the predicted quarter-power scaling relationship, but ordinary and phylogenetic generalized least squares did not. Regardless of the statistical approach, body size alone explained less than 28% of the variation in development time. Models that also included optimal temperature explained over 50% of the variation in development time. Warm-adapted insects developed more quickly, regardless of body size, supporting the "hotter is better" hypothesis that posits that ectotherms have a limited ability to evolutionarily compensate for the depressing effects of low temperatures on rates of biological processes. The remaining unexplained variation in development time likely reflects additional ecological and evolutionary differences among insect species.

Ecological genomics of local adaptation

It is increasingly important to improve our understanding of the genetic basis of local adaptation because of its relevance to climate change, crop and animal production, and conservation of genetic resources. Phenotypic patterns that are generated by spatially varying selection have long been observed, and both genetic mapping and field experiments provided initial insights into the genetic architecture of adaptive traits. Genomic tools are now allowing genome-wide studies, and recent theoretical advances can help to design research strategies that combine genomics and field experiments to examine the genetics of local adaptation. These advances are also allowing research in non-model species, the adaptation patterns of which may differ from those of traditional model species.

Searching for the bull's eye: agents and targets of selection vary among geographically disparate cyanogenesis clines in white clover (Trifolium repens L.)

The recurrent evolution of adaptive clines within a species can be used to elucidate the selective factors and genetic responses that underlie adaptation. White clover is polymorphic for cyanogenesis (HCN release with tissue damage), and climate-associated cyanogenesis clines have evolved throughout the native and introduced species range. This polymorphism arises through two independently segregating Mendelian polymorphisms for the presence/absence of two required components: cyanogenic glucosides and their hydrolyzing enzyme linamarase. Cyanogenesis is commonly thought to function in herbivore defense; however, the individual cyanogenic components may also serve other physiological functions. To test whether cyanogenesis clines have evolved in response to the same selective pressures acting on the same genetic targets, we examined cyanogenesis cline shape and its environmental correlates in three world regions: southern New Zealand, the central United States and the US Pacific Northwest. For some regional comparisons, cline shapes are remarkably similar despite large differences in the spatial scales over which clines occur (40-1600 km). However, we also find evidence for major differences in both the agents and targets of selection among the sampled clines. Variation in cyanogenesis frequency is best predicted using a combination of minimum winter temperature and aridity variables. Together, our results provide evidence that recurrent adaptive clines do not necessarily reflect shared adaptive processes.

Antagonistic selection factors induce a continuous population divergence in a polymorphism

Understanding the relative importance of selection and stochastic factors in population divergence of adaptive traits is a classical topic in evolutionary biology. However, it is difficult to separate these factors and detect the effects of selection when two or more contrasting selective factors are simultaneously acting on a single locus. In the damselfly Ischnura senegalensis, females exhibit color dimorphism and morph frequencies change geographically. We here evaluated the role of selection and stochastic factors in population divergence of morph frequencies by comparing the divergences in color locus and neutral loci. Comparisons between population pairwise FST for neutral loci and for the color locus did not detect any stochastic factors affecting color locus. Although comparison between population divergence in color and neutral loci using all populations detected only divergent selection, we detected two antagonistic selective factors acting on the color locus, that is, balancing and divergent selection, when considering geographical distance between populations. Our results suggest that a combination of two antagonistic selective factors, rather than stochastic factors, establishes the geographic cline in morph frequency in this system.

Rapid adaptation to climate facilitates range expansion of an invasive plant

Adaptation to climate, evolving over contemporary time scales, could facilitate rapid range expansion across environmental gradients. Here, we examine local adaptation along a climatic gradient in the North American invasive plant Lythrum salicaria. We show that the evolution of earlier flowering is adaptive at the northern invasion front where it increases fitness as much as, or more than, the effects of enemy release and the evolution of increased competitive ability. However, early flowering decreases investment in vegetative growth, which reduces fitness by a factor of 3 in southern environments where the North American invasion commenced. Our results demonstrate that local adaptation can evolve quickly during range expansion, overcoming environmental constraints on propagule production.

High genetic diversity is not essential for successful introduction

Some introduced populations thrive and evolve despite the presumed loss of diversity at introduction. We aimed to quantify the amount of genetic diversity retained at introduction in species that have shown evidence of adaptation to their introduced environments. Samples were taken from native and introduced ranges of Arctotheca populifolia and Petrorhagia nanteuilii. Using microsatellite data, we identified the source for each introduction, estimated genetic diversity in native and introduced populations, and calculated the amount of diversity retained in introduced populations. These values were compared to those from a literature review of diversity in native, confamilial populations and to estimates of genetic diversity retained at introduction. Gene diversity in the native range of both species was significantly lower than for confamilials. We found that, on average, introduced populations showing evidence of adaptation to their new environments retained 81% of the genetic diversity from the native range. Introduced populations of P. nanteuilii had higher genetic diversity than found in the native source populations, whereas introduced populations of A. populifolia retained only 14% of its native diversity in one introduction and 1% in another. Our literature review has shown that most introductions demonstrating adaptive ability have lost diversity upon introduction. The two species studied here had exceptionally low native range genetic diversity. Further, the two introductions of A. populifolia represent the largest percentage loss of genetic diversity in a species showing evidence of substantial morphological change in the introduced range. While high genetic diversity may increase the likelihood of invasion success, the species examined here adapted to their new environments with very little neutral genetic diversity. This finding suggests that even introductions founded by small numbers of individuals have the potential to become invasive.

Genetic architecture of parallel pelvic reduction in ninespine sticklebacks

Teleost fish genomes are known to be evolving faster than those of other vertebrate taxa. Thus, fish are suited to address the extent to which the same vs. different genes are responsible for similar phenotypic changes in rapidly evolving genomes of evolutionary independent lineages. To gain insights into the genetic basis and evolutionary processes behind parallel phenotypic changes within and between species, we identified the genomic regions involved in pelvic reduction in Northern European ninespine sticklebacks (Pungitius pungitius) and compared them to those of North American ninespine and threespine sticklebacks (Gasterosteus aculeatus). To this end, we conducted quantitative trait locus (QTL) mapping using 283 F₂ progeny from an interpopulation cross. Phenotypic analyses indicated that pelvic reduction is a recessive trait and is inherited in a simple Mendelian fashion. Significant QTL for pelvic spine and girdle lengths were identified in the region of the Pituitary homeobox transcription factor 1 (Pitx1) gene, also responsible for pelvic reduction in threespine sticklebacks. The fact that no QTL was observed in the region identified in the mapping study of North American ninespine sticklebacks suggests that an alternative QTL for pelvic reduction has emerged in this species within the past 1.6 million years after the split between Northern European and North American populations. In general, our study provides empirical support for the view that alternative genetic mechanisms that lead to similar phenotypes can evolve over short evolutionary time scales.

Two invaders instead of one: the true identity of species under the name Ceracis cucullatus (Coleoptera: Ciidae)

The Neotropical obligate fungivorous beetle Ceracis cucullatus (Mellié) has attracted attention of coleopterists due to the increasing number of records of populations in Africa. Although its disjunct populations have been interpreted as a cohesive taxonomic unity, previous comparisons between African and Neotropical specimens revealed differences in their external morphology, causing uncertainty about the true unity of the species. Here, we compare the external morphology of specimens named Cer. cucullatus from several localities of the Neotropical, Palearctic, Afrotropical, Afrotemperate and Oriental regions. As results, we reverse three previous junior synonymies of Cer. cucullatus, proposing Cer. lamellatus (Pic) and Cer. tabellifer (Mellié), both reinstated status and new combinations, as separate species. We also propose Enn. bilamellatum Pic as a new synonym of Cer. tabellifer. In face of these taxonomic changes, we identify Cer. tabellifer as the actual invasive species on African lands, instead of Cer. cucullatus as was previously accepted. Then, through historical records gathered from scientific collections and literature, and through examination of recently collected specimens from South Africa and Brazil, we provide data on host fungi and geographic distribution of Cer. tabellifer. Based on these data, we discuss possible explanations to the successful invasion of Cer. tabellifer in Africa and elsewhere and its potential threat to native faunas of ciids. This study helps to fulfil an old gap in the literature on biological invasions, with considerably more studies on predatory species, disease vectors or potential pests of agricultural crops, than on non-pest fungivorous organisms.

The Fukushima nuclear accident and the pale grass blue butterfly: evaluating biological effects of long-term low-dose exposures

Background

On August 9th 2012, we published an original research article in Scientific Reports, concluding that artificial radionuclides released from the Fukushima Dai-ichi Nuclear Power Plant exerted genetically and physiologically adverse effects on the pale grass blue butterfly Zizeeria maha in the Fukushima area. Immediately following publication, many questions and comments were generated from all over the world. Here, we have clarified points made in the original paper and answered questions posed by the readers.

Results

The following points were clarified. (1) There are many advantages to using the pale grass blue butterfly as an indicator species. (2) The forewings of the individuals collected in Fukushima were significantly smaller than in the northern and southern localities. (3) We observed growth retardation in the butterflies from the Fukushima area. (4) The aberrant colour patterns in the butterflies obtained in the Fukushima area were different from the colour patterns induced by temperature and sibling crosses but similar to those induced by external and internal exposures to the artificial radionuclides and by a chemical mutagen, suggesting that genetic mutations caused the aberrations. (5) This species of butterfly has been plentiful in Fukushima area for at least half a century. We here present specimens collected from Fukushima Prefecture before the accident. (6) Mutation accumulation was detected by the increase in the abnormality rates from May 2011 to September 2011. (7) The abnormal traits were heritable. (8) Our sampling localities were not affected by the tsunami. (9) We used a high enough number of samples to obtain statistically significant results. (10) The standard rearing method was followed, producing normal adults in the control groups. (11) The exposure experiments successfully reproduced the results of the field work. This species of butterfly is vulnerable to long-term low-dose internal and external exposures; however, insect cells are known to be resistant to short-term high-dose irradiation. This discrepancy is reconcilable based on the differences in the experimental conditions.

Conclusions

We are just beginning to understand the biological effects of long-term low-dose exposures in animals. Further research is necessary to accurately assess the possible biological effects of the accident.

Gene-by-temperature interactions and candidate plasticity genes for morphological traits in Drosophila melanogaster

Understanding the genetic architecture of any quantitative trait requires identifying the genes involved in its expression in different environmental conditions. This goal can be achieved by mutagenesis screens in genetically tractable model organisms such as Drosophila melanogaster. Temperature during ontogenesis is an important environmental factor affecting development and phenotypic variation in holometabolous insects. In spite of the importance of phenotypic plasticity and genotype by environment interaction (GEI) for fitness related traits, its genetic basis has remained elusive. In this context, we analyzed five different adult morphological traits (face width, head width, thorax length, wing size and wing shape) in 42 co-isogenic single P-element insertional lines of Drosophila melanogaster raised at 17°C and 25°C. Our analyses showed that all lines differed from the control for at least one trait in males or females at either temperature. However, no line showed those differences for all traits in both sexes and temperatures simultaneously. In this sense, the most pleiotropic candidate genes were CG34460, Lsd-2 and Spn. Our analyses also revealed extensive genetic variation for all the characters mostly indicated by strong GEIs. Further, our results indicate that GEIs were predominantly explained by changes in ranking order in all cases suggesting that a moderate number of genes are involved in the expression of each character at both temperatures. Most lines displayed a plastic response for at least one trait in either sex. In this regard, P-element insertions affecting plasticity of a large number of traits were associated to the candidate genes Btk29A, CG43340, Drak and jim. Further studies will help to elucidate the relevance of these genes on the morphogenesis of different body structures in natural populations of D. melanogaster.

True polyandry and pseudopolyandry: why does a monandrous fly remate?

BACKGROUND

The rate of female remating can have important impacts on a species, from affecting conflict and cooperation within families, to population viability and gene flow. However, determining the level of polyandry in a species can be difficult, with information on the mating system of many species being based on a single experiment, or completely absent. Here we investigate the mating system of the fruit fly Drosophila subobscura. Reports from England, Spain and Canada suggest D. subobscura is entirely monandrous, with no females remating. However, work in Greece suggests that 23% of females remate. We examine the willingness of female D. subobscura to remate in the laboratory in a range of conditions, using flies from both Greece and England. We make a distinction between pseudopolyandry, where a female remates after an ineffective first mating that is incapable of fertilising her eggs, and true polyandry, where a female remates even though she has received suitable sperm from a previous mating.

RESULTS

We find a low rate of true polyandry by females (4%), with no difference between populations. The rate of true polyandry is affected by temperature, but not starvation. Pseudopolyandry is three times as common as true polyandry, and most females showing pseudopolyandry mated at their first opportunity after their first failed mating. However, despite the lack of differences in polyandry between the populations, we do find differences in the way males respond to exposure to other males prior to mating. In line with previous work, English flies responded to one or more rivals by increasing their copulation duration, a response previously thought to be driven by sperm competition. Greek males only show increased copulation duration when exposed to four or more rival males. This suggests that the response to rivals in D. subobscura is not related to sperm competition, because sperm competition is rare, and there is no correlation of response to rivals and mating system across the populations.

CONCLUSIONS

These results illustrate the difficulties in determining the mating system of a species, even one that is well known and an excellent laboratory species, with results being highly dependent on the conditions used to assay the behaviour, and the population used.

Pleiotropy of FRIGIDA enhances the potential for multivariate adaptation

An evolutionary response to selection requires genetic variation; however, even if it exists, then the genetic details of the variation can constrain adaptation. In the simplest case, unlinked loci and uncorrelated phenotypes respond directly to multivariate selection and permit unrestricted paths to adaptive peaks. By contrast, 'antagonistic' pleiotropic loci may constrain adaptation by affecting variation of many traits and limiting the direction of trait correlations to vectors that are not favoured by selection. However, certain pleiotropic configurations may improve the conditions for adaptive evolution. Here, we present evidence that the Arabidopsis thaliana gene FRI (FRIGIDA) exhibits 'adaptive' pleiotropy, producing trait correlations along an axis that results in two adaptive strategies. Derived, low expression FRI alleles confer a 'drought escape' strategy owing to fast growth, low water use efficiency and early flowering. By contrast, a dehydration avoidance strategy is conferred by the ancestral phenotype of late flowering, slow growth and efficient water use during photosynthesis. The dehydration avoidant phenotype was recovered when genotypes with null FRI alleles were transformed with functional alleles. Our findings indicate that the well-documented effects of FRI on phenology result from differences in physiology, not only a simple developmental switch.

Wing shape allometry and aerodynamics in calopterygid damselflies: a comparative approach

BACKGROUND

Wing size and shape have important aerodynamic implications on flight performance. We explored how wing size was related to wing shape in territorial males of 37 taxa of the damselfly family Calopterygidae. Wing coloration was also included in the analyses because it is sexually and naturally selected and has been shown to be related to wing shape. We studied wing shape using both the non-dimensional radius of the second moment of wing area (RSM) and geometric morphometrics. Lower values of the RSM result in less energetically demanding flight and wider ranges of flight speed. We also re-analyzed previously published data on other damselflies and dragonflies.

RESULTS

The RSM showed a hump-shaped relationship with wing size. However, after correcting for phylogeny using independent contrast, this pattern changed to a negative linear relationship. The basal genus of the study family, Hetaerina, was mainly driving that change. The obtained patterns were specific for the study family and differed from other damselflies and dragonflies. The relationship between the RSM and wing shape measured by geometric morphometrics was linear, but relatively small changes along the RSM axis can result in large changes in wing shape. Our results also showed that wing coloration may have some effect on RSM.

CONCLUSIONS

We found that RSM showed a complex relationship with size in calopterygid damselflies, probably as a result of other selection pressures besides wing size per se. Wing coloration and specific behavior (e.g. courtship) are potential candidates for explaining the complexity. Univariate measures of wing shape such as RSM are more intuitive but lack the high resolution of other multivariate techniques such as geometric morphometrics. We suggest that the relationship between wing shape and size are taxa-specific and differ among closely-related insect groups.

Latitudinal clines in the grasshopper Dichroplus elongatus: coevolution of the A genome and B chromosomes?

Argentine populations of Dichroplus elongatus (Orthoptera: Acrididae) are polymorphic for B chromosomes. Previous studies showed that B chromosomes affect body size and some fitness components in Northwestern populations. We studied phenotype and B's variation patterns along a latitudinal cline as well as the relationship between karyotype and body size related traits in 17 populations from East. Body size related traits showed a 'saw tooth' pattern of variation being small at low and high latitudes and large at intermediate latitudes in most of the analysed populations. Analyses of variance and principal components demonstrated that in most analysed populations B carrier males are associated with a decrease in body size related traits with respect to individuals with standard karyotype. Accordingly with the relationship between karyotype and body size, an opposite pattern of latitudinal variation in the frequencies of B's with respect to body size variation was observed in this area. i.e. smaller individuals tend to have a higher frequency of B chromosomes. The comparison of the differentiation of both karyotype and body size traits with molecular neutral markers demonstrated the relative importance of selection moulding chromosome and phenotype variation. The observed pattern of phenotypic variation is likely to be the result of local adaptation to season length along the latitudinal gradient. The observed contrary pattern of B's clinal variation may reflect the population ability to maintain this chromosome in relation to the local adaptation. The available evidence indicates that the distribution of B chromosome frequency was shaped by selective factors.

Polymorphism in the neurofibromin gene, Nf1, is associated with antagonistic selection on wing size and development time in Drosophila melanogaster

In many invertebrates, body size shows genetically based clines, with size increasing in colder climates. Large body size is typically associated with prolonged development times. We consider variation in the CNS-specific gene neurofibromin 1 (Nf1) and its association with body size and development time. We identified two major Nf1 haplotypes in natural populations, Nf1-insertion-A and Nf1-deletion-G. These haplotypes are characterized by a 45-base insertion/deletion (INDEL) in Nf1 intron 2 and an A/G synonymous substitution (locus L17277). Linkage disequilibrium (LD) between the INDEL and adjacent sites is high but appears to be restricted within the Nf1 gene interval. In Australia, the frequency of the Nf1-insertion-A haplotype increases with latitude where wing size is larger, independent of the chromosomal inversion In(3R)Payne. Unexpectedly, the Nf1-insertion-A haplotype is negatively associated with wing size. We found that the Nf1-insertion-A haplotype is enriched in females with shorter development time. This suggests that the Nf1 haplotype cline may be driven by selection for development time rather than size; females from southern (higher latitude) D. melanogaster populations maintain a rapid development time despite being relatively larger, and the higher incidence of Nf1-insertion-A in Southern Australia may contribute to this pattern, whereas the effects of the Nf1 haplotypes on size may be countered by other loci with antagonistic effects on size and development time. Our results point to the potential complexity involved in identifying selection on genetic variants exhibiting pleiotropic effects when studies are based on spatial patterns or association studies.

Northern range expansion of European populations of the wasp spider Argiope bruennichi is associated with global warming-correlated genetic admixture and population-specific temperature adaptations

Poleward range expansions are observed for an increasing number of species, which may be an effect of global warming during the past decades. However, it is still not clear in how far these expansions reflect simple geographical shifts of species ranges, or whether new genetic adaptations play a role as well. Here, we analyse the expansion of the wasp spider Argiope bruennichi into Northern Europe during the last century. We have used a range-wide sampling of contemporary populations and historical specimens from museums to trace the phylogeography and genetic changes associated with the range shift. Based on the analysis of mitochondrial, microsatellite and SNP markers, we observe a higher level of genetic diversity in the expanding populations, apparently due to admixture of formerly isolated lineages. Using reciprocal transplant experiments for testing overwintering tolerance, as well as temperature preference and tolerance tests in the laboratory, we find that the invading spiders have possibly shifted their temperature niche. This may be a key adaptation for survival in Northern latitudes. The museum samples allow a reconstruction of the invasion's genetic history. A first, small-scale range shift started around 1930, in parallel with the onset of global warming. A more massive invasion of Northern Europe associated with genetic admixture and morphological changes occurred in later decades. We suggest that the latter range expansion into far Northern latitudes may be a consequence of the admixture that provided the genetic material for adaptations to new environmental regimes. Hence, global warming could have facilitated the initial admixture of populations and this resulted in genetic lineages with new habitat preferences.

Rates of molecular evolution and diversification in plants: chloroplast substitution rates correlate with species-richness in the Proteaceae

BACKGROUND

Many factors have been identified as correlates of the rate of molecular evolution, such as body size and generation length. Analysis of many molecular phylogenies has also revealed correlations between substitution rates and clade size, suggesting a link between rates of molecular evolution and the process of diversification. However, it is not known whether this relationship applies to all lineages and all sequences. Here, in order to investigate how widespread this phenomenon is, we investigate patterns of substitution in chloroplast genomes of the diverse angiosperm family Proteaceae. We used DNA sequences from six chloroplast genes (6278bp alignment with 62 taxa) to test for a correlation between diversification and the rate of substitutions.

RESULTS

Using phylogenetically-independent sister pairs, we show that species-rich lineages of Proteaceae tend to have significantly higher chloroplast substitution rates, for both synonymous and non-synonymous substitutions.

CONCLUSIONS

We show that the rate of molecular evolution in chloroplast genomes is correlated with net diversification rates in this large plant family. We discuss the possible causes of this relationship, including molecular evolution driving diversification, speciation increasing the rate of substitutions, or a third factor causing an indirect link between molecular and diversification rates. The link between the synonymous substitution rate and clade size is consistent with a role for the mutation rate of chloroplasts driving the speed of reproductive isolation. We find no significant differences in the ratio of non-synonymous to synonymous substitutions between lineages differing in net diversification rate, therefore we detect no signal of population size changes or alteration in selection pressures that might be causing this relationship.

Patterns of nucleotide diversity at the regions encompassing the Drosophila insulin-like peptide (dilp) genes: demography vs. positive selection in Drosophila melanogaster

In Drosophila, the insulin-signaling pathway controls some life history traits, such as fertility and lifespan, and it is considered to be the main metabolic pathway involved in establishing adult body size. Several observations concerning variation in body size in the Drosophila genus are suggestive of its adaptive character. Genes encoding proteins in this pathway are, therefore, good candidates to have experienced adaptive changes and to reveal the footprint of positive selection. The Drosophila insulin-like peptides (DILPs) are the ligands that trigger the insulin-signaling cascade. In Drosophila melanogaster, there are several peptides that are structurally similar to the single mammalian insulin peptide. The footprint of recent adaptive changes on nucleotide variation can be unveiled through the analysis of polymorphism and divergence. With this aim, we have surveyed nucleotide sequence variation at the dilp1-7 genes in a natural population of D. melanogaster. The comparison of polymorphism in D. melanogaster and divergence from D. simulans at different functional classes of the dilp genes provided no evidence of adaptive protein evolution after the split of the D. melanogaster and D. simulans lineages. However, our survey of polymorphism at the dilp gene regions of D. melanogaster has provided some evidence for the action of positive selection at or near these genes. The regions encompassing the dilp1-4 genes and the dilp6 gene stand out as likely affected by recent adaptive events.

Divergent strategies for adaptations to stress resistance in two tropical Drosophila species: effects of developmental acclimation for D. bipectinata and the invasive species - D. malerkotliana

Previous studies on two tropical Drosophila species (D. malerkotliana and D. bipectinata) have shown lower resistance to stress related traits but the rapid colonization of D. malerkotliana in the past few decades is not consistent with its sensitivity to desiccation and cold stress. We tested the hypothesis whether developmental acclimation at two growth temperatures (17 & 25 °C) can confer adaptations to desiccation and thermal stresses. We found divergence in developmental plastic effects on cuticular traits i.e. a significant increase of body melanisation (~2 fold) and of cuticular lipid mass (~3 fold) in D. malerkotliana but only 1.5 fold higher cuticular lipid mass in D. bipectinata when grown at 17 °C as compared with 25 °C. A comparison of water budget of these two species showed significantly higher effects of developmental acclimation on body water content, reduced rate of water loss and greater dehydration tolerance to confer higher desiccation resistance in D. malerkotliana as compared with D. bipectinata. Under less warmer growth conditions (17 °C), D. malerkotliana has evidenced greater resistance to cold as well as desiccation stress. In contrast, heat resistance of D. bipectinata is higher than D. malerkotliana when grown at 25 °C. These laboratory observations find support from data on seasonally varying populations. Further, adults acclimated to different stresses showed greater increase in D. malerkotliana as compared with D. bipectinata. Thus, significant increase of stress resistance of D. malerkotliana via developmental acclimation may be responsible for its invasion and ecological success on different continents as compared with D. bipectinata.

Genetic and morphometric divergence of an invasive bird: the introduced house sparrow (Passer domesticus) in Brazil

Introduced species are interesting systems for the study of contemporary evolution in new environments because of their spatial and temporal scales. For this study we had three aims: (i) to determine how genetic diversity and genetic differentiation of introduced populations of the house sparrow (Passer domesticus) in Brazil varies with range expansion, (ii) to determine how genetic diversity and differentiation in Brazil compares to ancestral European populations; and (iii) to determine whether selection or genetic drift has been more influential on phenotypic divergence. We used six microsatellite markers to genotype six populations from Brazil and four populations from Europe. We found slightly reduced levels of genetic diversity in Brazilian compared to native European populations. However, among introduced populations of Brazil, we found no association between genetic diversity and time since introduction. Moreover, overall genetic differentiation among introduced populations was low indicating that the expansion took place from large populations in which genetic drift effects would likely have been weak. We found significant phenotypic divergence among sites in Brazil. Given the absence of a spatial genetic pattern, divergent selection and not genetic drift seems to be the main force behind most of the phenotypic divergence encountered. Unravelling whether microevolution (e.g., allele frequency change), phenotypic plasticity, or both mediated phenotypic divergence is challenging and will require experimental work (e.g., common garden experiments or breeding programs).

Genomic conflict drives patterns of X-linked population structure in Drosophila neotestacea

Intragenomic conflict has the potential to cause widespread changes in patterns of genetic diversity and genome evolution. In this study, we investigate the consequences of sex-ratio (SR) drive on the population genetic patterns of the X-chromosome in Drosophila neotestacea. An SR X-chromosome prevents the maturation of Y-bearing sperm during male spermatogenesis and thus is transmitted to ~100% of the offspring, nearly all of which are daughters. Selection on the rest of the genome to suppress SR can be strong, and the resulting conflict over the offspring sex ratio can result in the accumulation of multiple loci on the X-chromosome that are necessary for the expression of drive. We surveyed variation at 12 random X-linked microsatellites across 16 populations of D. neotestacea that range in SR frequency from 0% to 30%. First, every locus was differentiated between SR and wild-type chromosomes, and this drives genetic structure at the X-chromosome. Once the association with SR is accounted for, the patterns of differentiation among populations are similar to the autosomes. Second, within wild-type chromosomes, the relative heterozygosity is reduced in populations with an increased prevalence of drive, and the heterozygosity of SR chromosomes is higher than expected based on its prevalence. The combination of the relatively high prevalence of SR drive and the structuring of polymorphism between the SR and wild-type chromosomes suggests that genetic conflict because of SR drive has had significant consequences on the patterns of X-linked polymorphism and thus also probably affects the tempo of X-chromosome evolution in D. neotestacea.