Genetic variation within a stick insect species associated with community-level traits.

Phenotypic variation within species can affect the ecological dynamics of populations and communities. Characterizing the genetic variation underlying such effects can help parse the roles of genetic evolution and plasticity in 'eco-evolutionary dynamics' and inform how genetic variation may shape patterns of evolution. Here we employ genome-wide association (GWA) methods in Timema cristinae stick insects and their co-occurring arthropod communities to identify genetic variation associated with community-level traits. Previous studies have shown that maladaptation (i.e., imperfect crypsis) of T. cristinae can reduce the abundance and species richness of other arthropods due to an increase in bird predation. Whether genetic variation that is independent from crypsis has similar effects is unknown and was tested here using genome-wide genotyping-by-sequencing data of stick insects, arthropod community information, and GWA mapping with Bayesian sparse linear mixed models. We find associations between genetic variation in stick insects and arthropod community traits. However, these associations disappeared when host-plant traits are accounted for. We thus use path analysis to disentangle interrelationships among stick-insect genetic variation, host-plant traits and community traits. This revealed that host-plant size has large effects on arthropod communities, while genetic variation in stick insects has a smaller, but still significant effect. Our findings demonstrate that: (1) genetic variation in a species can be associated with community-level traits, but that (2) interrelationships among multiple factors may need to be analyzed to disentangle whether such associations represent causal relationships. This work helps to build a framework for genomic studies of eco-evolutionary dynamics.

Divergence and introgression among the virilis group of Drosophila.

Speciation with gene flow is now widely regarded as common. However, the frequency of introgression between recently diverged species and the evolutionary consequences of gene flow are still poorly understood. The virilis group of Drosophila contains 12 species that are geographically widespread and show varying levels of prezygotic and postzygotic isolation. Here, we use de novo genome assemblies and whole-genome sequencing data to resolve phylogenetic relationships and describe patterns of introgression and divergence across the group. We suggest that the virilis group consists of three, rather than the traditional two, subgroups. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and desiccation tolerance, and may be related to the evolution of sexual isolation and adaptation. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and more recent gene flow between closely related species. When assessing patterns of genome-wide divergence in species pairs across the group, we found no consistent genomic evidence of a disproportionate role for the X chromosome as has been found in other systems. Our results show how ancient and recent introgressions confuse phylogenetic reconstruction, but may play an important role during early radiation of a group.

DrosoPhyla: Resources for Drosophilid Phylogeny and Systematics.

The vinegar fly Drosophila melanogaster is a pivotal model for invertebrate development, genetics, physiology, neuroscience, and disease. The whole family Drosophilidae, which contains over 4,400 species, offers a plethora of cases for comparative and evolutionary studies. Despite a long history of phylogenetic inference, many relationships remain unresolved among the genera, subgenera, and species groups in the Drosophilidae. To clarify these relationships, we first developed a set of new genomic markers and assembled a multilocus data set of 17 genes from 704 species of Drosophilidae. We then inferred a species tree with highly supported groups for this family. Additionally, we were able to determine the phylogenetic position of some previously unplaced species. These results establish a new framework for investigating the evolution of traits in fruit flies, as well as valuable resources for systematics.

Multiple paths to cold tolerance: the role of environmental cues, morphological traits and the circadian clock gene vrille.

BACKGROUND: Tracing the association between insect cold tolerance and latitudinally and locally varying environmental conditions, as well as key morphological traits and molecular mechanisms, is essential for understanding the processes involved in adaptation. We explored these issues in two closely-related species, Drosophila montana and Drosophila flavomontana, originating from diverse climatic locations across several latitudes on the coastal and mountainous regions of North America. We also investigated the association between sequence variation in one of the key circadian clock genes, vrille, and cold tolerance in both species. Finally, we studied the impact of vrille on fly cold tolerance and cold acclimation ability by silencing it with RNA interference in D. montana. RESULTS: We performed a principal component analysis (PCA) on variables representing bioclimatic conditions on the study sites and used latitude as a proxy of photoperiod. PC1 separated the mountainous continental sites from the coastal ones based on temperature variability and precipitation, while PC2 arranged the sites based on summer and annual mean temperatures. Cold tolerance tests showed D. montana to be more cold-tolerant than D. flavomontana and chill coma resistance (CTmin) of this species showed an association with PC2. Chill coma recovery time (CCRT) of both species improved towards northern latitudes, and in D. flavomontana this trait was also associated with PC1. D. flavomontana flies were darkest in the coast and in the northern mountainous populations, but coloration showed no linkage with cold tolerance. Body size decreased towards cold environments in both species, but only within D. montana populations largest flies showed fastest recovery from cold. Finally, both the sequence analysis and RNAi study on vrille suggested this gene to play an essential role in D. montana cold resistance and acclimation, but not in recovery time. CONCLUSIONS: Our study demonstrates the complexity of insect cold tolerance and emphasizes the need to trace its association with multiple environmental variables and morphological traits to identify potential agents of natural selection. It also shows that a circadian clock gene vrille is essential both for short- and long-term cold acclimation, potentially elucidating the connection between circadian clock system and cold tolerance.

Functional Genomics Offers New Tests of Speciation Hypotheses.

Speciation is a fundamental process shaping biodiversity. However, existing empirical methods often cannot provide key genetic and functional details required to validate speciation theory. New gene modification technologies can verify the causal functionality of genes with astonishing accuracy, helping resolve questions about how reproductive isolation evolves during speciation.

Adaptation of gene loci to heterochromatin in the course of Drosophila evolution is associated with insulator proteins.

Pericentromeric heterochromatin is generally composed of repetitive DNA forming a transcriptionally repressive environment. Dozens of genes were embedded into pericentromeric heterochromatin during evolution of Drosophilidae lineage while retaining activity. However, factors that contribute to insusceptibility of gene loci to transcriptional silencing remain unknown. Here, we find that the promoter region of genes that can be embedded in both euchromatin and heterochromatin exhibits a conserved structure throughout the Drosophila phylogeny and carries motifs for binding of certain chromatin remodeling factors, including insulator proteins. Using ChIP-seq data, we demonstrate that evolutionary gene relocation between euchromatin and pericentric heterochromatin occurred with preservation of sites of insulation of BEAF-32 in evolutionarily distant species, i.e. D. melanogaster and D. virilis. Moreover, promoters of virtually all protein-coding genes located in heterochromatin in D. melanogaster are enriched with insulator proteins BEAF-32, GAF and dCTCF. Applying RNA-seq of a BEAF-32 mutant, we show that the impairment of BEAF-32 function has a complex effect on gene expression in D. melanogaster, affecting even those genes that lack BEAF-32 association in their promoters. We propose that conserved intrinsic properties of genes, such as sites of insulation near the promoter regions, may contribute to adaptation of genes to the heterochromatic environment and, hence, facilitate the evolutionary relocation of genes loci between euchromatin and heterochromatin.

Diapause affects cuticular hydrocarbon composition and mating behavior of both sexes in Drosophila montana.

Environmental cues, mainly photoperiod and temperature, are known to control female adult reproductive diapause in several insect species. Diapause enhances female survival during adverse conditions and postpones progeny production to the favorable season. Male diapause (a reversible inability to inseminate receptive females) has been studied much less than female diapause. However, if the males maximized their chances to fertilize females while minimizing their energy expenditure, they would be expected to be in diapause at the same time as females. We investigated Drosophila montana male mating behavior under short-day conditions that induce diapause in females and found the males to be reproductively inactive. We also found that males reared under long-day conditions (reproducing individuals) court reproducing postdiapause females, but not diapausing ones. The diapausing flies of both sexes had more long-chain and less short-chain hydrocarbons on their cuticle than the reproducing ones, which presumably increase their survival under stressful conditions, but at the same time decrease their attractiveness. Our study shows that the mating behavior of females and males is well coordinated during and after overwintering and it also gives support to the dual role of insect cuticular hydrocarbons in adaptation and mate choice.

Inter and Intraspecific Genomic Divergence in Drosophila montana Shows Evidence for Cold Adaptation.

The genomes of species that are ecological specialists will likely contain signatures of genomic adaptation to their niche. However, distinguishing genes related to ecological specialism from other sources of selection and more random changes is a challenge. Here, we describe the genome of Drosophila montana, which is the most extremely cold-adapted Drosophila species known. We use branch tests to identify genes showing accelerated divergence in contrasts between cold- and warm-adapted species and identify about 250 genes that show differences, possibly driven by a lower synonymous substitution rate in cold-adapted species. We also look for evidence of accelerated divergence between D. montana and D. virilis, a previously sequenced relative, but do not find strong evidence for divergent selection on coding sequence variation. Divergent genes are involved in a variety of functions, including cuticular and olfactory processes. Finally, we also resequenced three populations of D. montana from across its ecological and geographic range. Outlier loci were more likely to be found on the X chromosome and there was a greater than expected overlap between population outliers and those genes implicated in cold adaptation between Drosophila species, implying some continuity of selective process at these different evolutionary scales.

Localization of quantitative trait loci for diapause and other photoperiodically regulated life history traits important in adaptation to seasonally varying environments.

Seasonally changing environments at high latitudes present great challenges for the reproduction and survival of insects, and photoperiodic cues play an important role in helping them to synchronize their life cycle with prevalent and forthcoming conditions. We have mapped quantitative trait loci (QTL) responsible for the photoperiodic regulation of four life history traits, female reproductive diapause, cold tolerance, egg-to-eclosion development time and juvenile body weight in Drosophila montana strains from different latitudes in Canada and Finland. The F2 progeny of the cross was reared under a single photoperiod (LD cycle 16:8), which the flies from the Canadian population interpret as early summer and the flies from the Finnish population as late summer. The analysis revealed a unique QTL for diapause induction on the X chromosome and several QTL for this and the other measured traits on the 4th chromosome. Flies' cold tolerance, egg-to-eclosion development time and juvenile body weight had several QTL also on the 2nd, 3rd and 5th chromosome, some of the peaks overlapping with each other. These results suggest that while the downstream output of females' photoperiodic diapause response is partly under a different genetic control from that of the other traits in the given day length, all traits also share some QTL, possibly involving genes with pleiotropic effects and/or multiple tightly linked genes. Nonoverlapping QTL detected for some of the traits also suggest that the traits are potentially capable of independent evolution, even though this may be restricted by epistatic interactions and/or correlations and trade-offs between the traits.

Involvement of circadian oscillation(s) in the photoperiodic time measurement and the induction of reproductive diapause in a northern Drosophila species.

An ability to predict forthcoming changes in environmental conditions and get prepared for them in advance is crucial for the survival and reproduction of organisms living in a seasonally changing environment. We have studied the possible involvement of circadian oscillator(s) in the photoperiodic timer controlling seasonal responses by tracing Drosophila montana females' diapause induction in constant darkness and in a classical Nanda-Hamner experiment. Nearly all females developed ovaries in continuous darkness, which shows the direct development to be their default developmental pathway in the absence of photoperiods. In Nanda-Hamner experiment the females' diapause incidence was close to zero in light:dark cycle 12:4 (photoperiod 16 h) and increased to nearly 100% in 12:8 and 12:12 (photoperiods 20 and 24 h). In longer photoperiods (28-72 h) the females' diapause percentages decreased gradually along with an increase in the length of the dark period, showing no peaks of high diapause incidence in the multiples of 24h. These findings suggest that the photoperiodic timer of D. montana is based on heavily damping circadian oscillator(s) or that it lacks strong oscillators. Damping of the photoperiodic timer under prolonged nights and constant darkness fits well with our earlier finding that these flies lose their locomotor activity rhythm in constant darkness, and suggests that the mechanisms underlying females' photoperiodic diapause response and their free-running locomotor activity rhythm may be partly based on same oscillators.

Northern Drosophila montana flies show variation both within and between cline populations in the critical day length evoking reproductive diapause.

Reproductive diapause, and its correct timing, plays an important role in the life cycle of many insect species living in a seasonally varying environment at high latitudes. In the present paper we have documented variation in the critical day length (CDL) for adult reproductive diapause and the steepness of photoperiodic response curves (PPRCs) in seven clinal populations of Drosophila montana in Finland between the latitudes 61 and 67°N, paying special attention to variation in these traits within and between cline populations. The isofemale lines representing these populations showed a sharp transition from 0% to 100% in females' diapause incidence in the shortening day lengths, indicated by steep PPRCs. The mean CDL showed a clear latitudinal cline decreasing by 1.6h from North to South regardless of the age of the lines, variation within the populations (i.e. among lines) in this trait being up to 3h. The steepness of the PPRCs correlated with the age of the line and this trait showed no clear latitudinal cline. Further studies on a large number of lines from one D. montana population confirmed that while maintaining the flies in diapause preventing conditions in the laboratory has no effect on CDL, older lines had steeper PPRCs. High variation in CDL within and between D. montana cline populations is likely to be heritable and provide a good potential for the evolution of photoperiodic responses. Information on genetic variation in life-history traits, such as diapause, is of utmost importance for predicting the ability of insects to survive in seasonally changing environmental conditions and to respond to long term changes in the length of the growing period e.g. by postponing the timing of diapause towards shorter day length and later calendar date.

What do we need to know about speciation?

Speciation has been a major focus of evolutionary biology research in recent years, with many important advances. However, some of the traditional organising principles of the subject area no longer provide a satisfactory framework, such as the classification of speciation mechanisms by geographical context into allopatric, parapatric and sympatry classes. Therefore, we have asked where speciation research should be directed in the coming years. Here, we present a distillation of questions about the mechanisms of speciation, the genetic basis of speciation and the relationship between speciation and diversity. Our list of topics is not exhaustive; rather we aim to promote discussion on research priorities and on the common themes that underlie disparate speciation processes.

Adaptation to a seasonally varying environment: a strong latitudinal cline in reproductive diapause combined with high gene flow in Drosophila montana.

Adaptation to seasonal changes in the northern hemisphere includes an ability to predict the forthcoming cold season from gradual changes in environmental cues early enough to prepare for the harsh winter conditions. The magnitude and speed of changes in these cues vary between the latitudes, which induces strong selection pressures for local adaptation.We studied adaptation to seasonal changes in Drosophila montana, a northern maltfly, by defining the photoperiodic conditions leading to adult reproductive diapause along a latitudinal cline in Finland and by measuring genetic differentiation and the amount of gene flow between the sampling sites with microsatellites. Our data revealed a clear correlation between the latitude and the critical day length (CDL), in which half of the females of different cline populations enter photoperiodic reproductive diapause. There was no sign of limited gene flow between the cline populations, even though these populations showed isolation by distance. Our results show that local adaptation may occur even in the presence of high gene flow, when selection for locally adaptive life-history traits is strong. A wide range of variation in the CDLs of the fly strains within and between the cline populations may be partly due to gene flow and partly due to the opposing selection pressures for fly reproduction and overwinter survival. This variation in the timing of diapause will enhance populations' survival over the years that differ in the severity of the winter and in the length of the warm period and may also help them respond to long-term changes in environmental conditions.