Cold adaptation drives population genomic divergence in the ecological specialist, Drosophila montana.

Detecting signatures of ecological adaptation in comparative genomics is challenging, but analysing population samples with characterised geographic distributions, such as clinal variation, can help identify genes showing covariation with important ecological variation. Here, we analysed patterns of geographic variation in the cold-adapted species Drosophila montana across phenotypes, genotypes and environmental conditions and tested for signatures of cold adaptation in population genomic divergence. We first derived the climatic variables associated with the geographic distribution of 24 populations across two continents to trace the scale of environmental variation experienced by the species, and measured variation in the cold tolerance of the flies of six populations from different geographic contexts. We then performed pooled whole genome sequencing of these six populations, and used Bayesian methods to identify SNPs where genetic differentiation is associated with both climatic variables and the population phenotypic measurements, while controlling for effects of demography and population structure. The top candidate SNPs were enriched on the X and fourth chromosomes, and they also lay near genes implicated in other studies of cold tolerance and population divergence in this species and its close relatives. We conclude that ecological adaptation has contributed to the divergence of D. montana populations throughout the genome and in particular on the X and fourth chromosomes, which also showed highest interpopulation FST . This study demonstrates that ecological selection can drive genomic divergence at different scales, from candidate genes to chromosome-wide effects.

The genetic architecture of sexually selected traits in two natural populations of Drosophila montana.

We investigated the genetic architecture of courtship song and cuticular hydrocarbon traits in two phygenetically distinct populations of Drosophila montana. To study natural variation in these two important traits, we analysed within-population crosses among individuals sampled from the wild. Hence, the genetic variation analysed should represent that available for natural and sexual selection to act upon. In contrast to previous between-population crosses in this species, no major quantitative trait loci (QTLs) were detected, perhaps because the between-population QTLs were due to fixed differences between the populations. Partitioning the trait variation to chromosomes suggested a broadly polygenic genetic architecture of within-population variation, although some chromosomes explained more variation in one population compared with the other. Studies of natural variation provide an important contrast to crosses between species or divergent lines, but our analysis highlights recent concerns that segregating variation within populations for important quantitative ecological traits may largely consist of small effect alleles, difficult to detect with studies of moderate power.

How consistent are the transcriptome changes associated with cold acclimation in two species of the Drosophila virilis group?

For many organisms the ability to cold acclimate with the onset of seasonal cold has major implications for their fitness. In insects, where this ability is widespread, the physiological changes associated with increased cold tolerance have been well studied. Despite this, little work has been done to trace changes in gene expression during cold acclimation that lead to an increase in cold tolerance. We used an RNA-Seq approach to investigate this in two species of the Drosophila virilis group. We found that the majority of genes that are differentially expressed during cold acclimation differ between the two species. Despite this, the biological processes associated with the differentially expressed genes were broadly similar in the two species. These included: metabolism, cell membrane composition, and circadian rhythms, which are largely consistent with previous work on cold acclimation/cold tolerance. In addition, we also found evidence of the involvement of the rhodopsin pathway in cold acclimation, a pathway that has been recently linked to thermotaxis. Interestingly, we found no evidence of differential expression of stress genes implying that long-term cold acclimation and short-term stress response may have a different physiological basis.

Photoperiodic regulation of life-history traits before and after eclosion: egg-to-adult development time, juvenile body mass and reproductive diapause in Drosophila montana.

Photoperiod is the main environmental cue used by northern insects to predict the forthcoming seasonal changes and to adjust their life-history traits to fit these changes. We studied the effects of photoperiod on egg-to-adult development time, juvenile body mass and female reproductive diapause in two northern Drosophila montana populations with different patterns of voltinism. The most interesting findings were consistent between the populations: (1) when maintained before eclosion in short day conditions, representing early autumn, the flies developed faster and were lighter than when maintained in long day conditions, representing early summer, (2) photoperiodic time measurement is apparently reset after eclosion, adjusting the flies' development according to post eclosion conditions, (3) the sensitive period for diapause induction took place after eclosion and (4) there was no direct connection between females' egg-to-adult development time and their reproductive state at adulthood, which suggests that these traits can be determined by photoperiodic cues through different time measurement systems. Independence of photoperiodic regulation of life-history traits before and after eclosion enable D. montana flies to respond to changing photoperiods on a short time scale and match their life-history traits according to seasons.

Photoperiodic regulation of cold tolerance and expression levels of regucalcin gene in Drosophila montana.

Temperature-induced plasticity of cold tolerance has been reported in many insect species, but cold tolerance can also be affected by changes in day (or night) length. In the present study we elucidate the direct and indirect effects of photoperiod on the cold tolerance of females of two Drosophila montana strains--one which possesses a robust photoperiodic diapause and another which does not. In the diapause-strain the time needed for recovery from chill coma showed a positive correlation with day length, but diapause itself played only a minor role in photoperiodic acclimation. The strain that was not able to enter to diapause as a response to day length also lacked photoperiodic cold acclimation ability indicating that this strain has deficiencies in its photoperiodic time measurement system. In the diapause-strain, the expression level of regucalcin gene was more than two times higher in diapausing than in non-diapausing females maintained in a single photoperiod, but day length per se did not cause significant changes in expression levels of this gene in either of the strains. In the non-diapausing strain this gene showed no expression changes in any comparison. Overall, the study shows that a decrease in day length can induce cold acclimation in D. montana, while changes in regucalcin expression are linked with photoperiodic diapause.

Two distinct genomic regions, harbouring the period and fruitless genes, affect male courtship song in Drosophila montana.

Acoustic signals often have a significant role in pair formation and in species recognition. Determining the genetic basis of signal divergence will help to understand signal evolution by sexual selection and its role in the speciation process. An earlier study investigated quantitative trait locus for male courtship song carrier frequency (FRE) in Drosophila montana using microsatellite markers. We refined this study by adding to the linkage map markers for 10 candidate genes known to affect song production in Drosophila melanogaster. We also extended the analyses to additional song characters (pulse train length (PTL), pulse number (PN), interpulse interval, pulse length (PL) and cycle number (CN)). Our results indicate that loci in two different regions of the genome control distinct features of the courtship song. Pulse train traits (PTL and PN) mapped to the X chromosome, showing significant linkage with the period gene. In contrast, characters related to song pulse properties (PL, CN and carrier FRE) mapped to the region of chromosome 2 near the candidate gene fruitless, identifying these genes as suitable loci for further investigations. In previous studies, the pulse train traits have been found to vary substantially between Drosophila species, and so are potential species recognition signals, while the pulse traits may be more important in intra-specific mate choice.

Cold tolerance and cold-induced modulation of gene expression in two Drosophila virilis group species with different distributions.

The importance of high and low temperature tolerance in adaptation to changing environmental conditions has evoked new interest in modulations in gene expression and metabolism linked with stress tolerance. We investigated the effects of rapid cold hardening and cold acclimatization on the chill coma recovery times of two Drosophila virilis group species, Drosophila montana and D. virilis, with different distributions and utilized a candidate gene approach to trace changes in their gene expression during and after the cold treatments. The study showed that cold acclimatization clearly decreases chill coma recovery times in both species, whereas rapid cold hardening did not have a significant effect. Microarray analysis revealed several genes showing expression changes during different stages of cold response. Amongst the 219 genes studied, two genes showed rather consistent expression changes: hsr-omega, which was up-regulated in both study species during cold acclimatization, and Eip71CD, which was down-regulated in nearly all of the cold treatments. In addition, 29 genes showed expression changes that were more treatment- and/or species specific. Overall, different stages of cold response elicited changes mainly in genes involved in heat shock response, circadian rhythm and metabolism.

Multiple quantitative trait loci influence intra-specific variation in genital morphology between phylogenetically distinct lines of Drosophila montana.

The evolution of animal genitalia has gained renewed interest because of their potential roles during sexual selection and early stages of species formation. Although central to understanding the evolutionary process, knowledge of the genetic basis of natural variation in genital morphology is limited to a very few species. Using an outbred cross between phylogenetically distinct lines of Drosophila montana, we characterized quantitative trait loci (QTLs) affecting the size and shape of the distiphallus, a prominent part of the male intromittent organ. Our microsatellite-based linkage analysis shows that intra-specific variation in the distiphallus involves several QTLs of largely additive effect and that a highly significant QTL co-localizes with the same inversion where we have earlier localized a large QTL for a sexually selected courtship song trait. The latter indicates that inversions can play an important role in shaping the evolution of rapidly evolving traits with a potential influence on speciation.

A microsatellite linkage map for Drosophila montana shows large variation in recombination rates, and a courtship song trait maps to an area of low recombination.

Current advances in genetic analysis are opening up our knowledge of the genetics of species differences, but challenges remain, particularly for out-bred natural populations. We constructed a microsatellite-based linkage map for two out-bred lines of Drosophila montana derived from divergent populations by taking advantage of the Drosophila virilis genome and available cytological maps of both species. Although the placement of markers was quite consistent with cytological predictions, the map indicated large heterogeneity in recombination rates along chromosomes. We also performed a quantitative trait locus (QTL) analysis on a courtship song character (carrier frequency), which differs between populations and is subject to strong sexual selection. Linkage mapping yielded two significant QTLs, which explained 3% and 14% of the variation in carrier frequency, respectively. Interestingly, as in other recent studies of traits which can influence speciation, the strongest QTL mapped to a genomic region partly covered by an inversion polymorphism.

The extent of variation in male song, wing and genital characters among allopatric Drosophila montana populations.

Drosophila montana, a species of the Drosophila virilis group, has distributed around the northern hemisphere. Phylogeographic analyses of two North American and one Eurasian population of this species offer a good background for the studies on the extent of variation in phenotypic traits between populations as well as for tracing the selection pressures likely to play a role in character divergence. In the present paper, we studied variation in the male courtship song, wing and genital characters among flies from Colorado (USA), Vancouver (Canada) and Oulanka (Finland) populations. The phenotypic divergence among populations did not coincide with the extent of their genetic divergence, suggesting that the characters are not evolving neutrally. Divergence in phenotypic traits was especially high between the Colorado and Vancouver populations, which are closer to each other in terms of their mtDNA genotypes than they are to the Oulanka population. The males of the Colorado population showed high divergence especially in song traits and the males of the Vancouver population in wing characters. Among the male song traits, two characters known to be under sexual selection and a trait important in species recognition differed clearly between populations, implying a history of directional and/or diversifying rather than balancing selection. The population divergence in wing characters is likely to have been enhanced by natural selection associated with environmental factors, whereas the male genitalia traits may have been influenced by sexual selection and/or sexual conflict.

Phylogeographic patterns in Drosophila montana.

The Drosophila virilis species group offers valuable opportunities for studying the roles of chromosomal re-arrangements and mating signals in speciation. The 13 species are divided into two subgroups, the montana and virilis 'phylads'. There is greater differentiation among species within the montana phylad in both karyotype and acoustic signals than exists among members of the virilis phylad. Drosophila montana is a divergent species which is included in the montana phylad. Here, we analyse the phylogeography of D. montana to provide a framework for understanding divergence of acoustic signals among populations. We analysed mitochondrial sequences corresponding to the cytochrome oxidase I and cytochrome oxidase II genes, as well as 16 microsatellite loci, from 108 lines of D. montana covering most of the species' range. The species shows a clear genetic differentiation between North American and Scandinavian populations. Microsatellite allele frequencies and mitochondrial DNA haplotypes gave significant FST values between populations from Canada, USA and Finland. A Bayesian analysis of population structure based on the microsatellite frequencies showed four genetically distinct groups, corresponding to these three populations plus a small sample from Japan. A network based on mitochondrial haplotypes showed two Finnish clades of very different shape and variability, and another clade with all sequences from North America and Japan. All D. montana populations showed evidence of demographic expansion but the patterns inferred by coalescent analysis differed between populations. The divergence times between Scandinavian and North American clades were estimated to range from 450,000 to 900,000 years with populations in Canada and the USA possibly representing descendants of different refugial populations. Long-term separation of D. montana populations could have provided the opportunity for differentiation observed in male signal traits, especially carrier frequency of the song, but relaxation of sexual selection during population expansion may have been necessary.

The role of male courtship song in species recognition in Drosophila montana.

In many Drosophila species the male courtship song is an essential part of a successful courtship, playing a role in species-recognition and/or in sexual selection exercised by the females on conspecific males. In Drosophila montana, specific traits of the courtship song (the length and the carrier frequency of sound pulses) have earlier been shown to play an important role in intraspecific mate choice. Here we show that the male song (especially the interpulse interval) also plays a role in species-recognition, maintaining sexual isolation between sympatric species. We succeeded in breaking down sexual isolation between D. montana females and D. lummei males by playing the females simulated courtship song with a species-specific interpulse interval while they were courted by a mute (wingless) male. Involvement of different song traits in intra- and interspecific mate choice suggests that the song may be affected by both directional and stabilizing selection.

QTL analysis of variation in male courtship song characters in Drosophila virilis.

We have used a quantitative trait locus (QTL) mapping approach to study the genetic basis of differences between two Drosophila virilis strains representing extreme phenotypes in two song characters, the number of pulses in a pulse train (PN) and the length of a pulse train (PTL). Variation in these characters among 520 F2 males was studied by single-marker analysis and composite interval mapping (CIM) using a recombination linkage map constructed for 26 microsatellite markers. In single-marker analysis, two adjacent microsatellite markers on the third chromosome, msat19 and vir84 explained 13.8 and 12.4% of the variation in PN and 9.9 and 6.5% of the variation in PTL, respectively. CIM analysis revealed significant QTLs affecting PN, located on the X and the second, third and fourth chromosome of D. virilis, while variation in PTL was attributable to QTLs located only on the third chromosome.

Enterolignans.

A review with 114 references about mammalian lignans (enterolignans). Several aspects have been reviewed: the precursors of mammalian lignans and their biosynthesis, biological activities and health effects, metabolism (in vivo and in vitro) in human and animals, some synthetic strategies to obtain enterolignan skeleton types, including the synthesis of haptens and deuterated lignans, and finally an overview of the analytical methods to detect and quantify lignans in biological matrices and foods.

Nucleotide variation at the no-on-transient A gene in Drosophila littoralis.

The no-on-transient A (nonA) gene encodes a putative RNA-binding protein, and mutations in this gene are known to affect vision, male courtship song and viability in Drosophila melanogaster. Here we have sequenced the coding region of the nonA gene of Drosophila littoralis and compared it with those of Drosophila virilis and D. melanogaster. All portions of nonA appeared to be conserved between D. littoralis and D. virilis, while the 5' region of the gene of these two species showed high divergence from that of a more distantly-related species, D. melanogaster. The same was true for the glycine repeat regions. No significant deviation from neutrality was observed in the analysis of intraspecific nucleotide variation in 5' or 3' region of the nonA gene in D. littoralis population. Also, comparison of D. littoralis sequences with homologous sequence of D. virilis suggests that the gene is evolving neutrally in D. virilis group. Divergence of the 5' regions between D. virilis group species and D. melanogaster could be a result of positive selection, but this finding is obscured by the long divergence time of the species groups.

Variability levels, population size and structure of American and European Drosophila montana populations.

The level and patterns of nucleotide diversity have been characterized for two X-linked loci, fused (fu; a region of 2362 bp) and suppressor of sable (su(s); a region of 413 bp), in one European and one American D. montana population. Sequence variation at these loci shows that the two populations are divergent, although they may not be completely isolated. Data on the level of silent site variability at su(s) (1.1% and 0.5% for the European and American populations, respectively) suggest that the effective population sizes of the two populations may be similar. At the fused locus, one European sequence was highly divergent and may have resulted from gene conversion, and was excluded from the analysis. With this sequence removed, the level of silent site variability was significantly lower in the European population (0.28%) than in the American population (2.3%), which suggests a selective sweep at or near fu in the former population.

Characterization of female preference functions for Drosophila montana courtship song and a test of the temperature coupling hypothesis.

Female mate preferences are a major cause of diversity and elaboration in male sexual traits. Here we characterize the shape of female preference functions for pulse length and carrier frequency of the courtship song of Drosophila montana by fitting both parametric and nonparametric functions to the incidence of female receptive gestures to synthetic song. Preference functions for both traits are strongly directional. That for pulse length is linear and favors short pulses, whereas that for carrier frequency is stabilizing in shape, but would exert directional preferences favoring males with high carrier frequency. The preference for carrier frequency has probably evolved under sexual selection, but reasons for the preference for short pulses are less apparent. We also examine the effect of ambient temperature on the carrier frequency of male song and on the preference function for carrier frequency. For many similar acoustic communication systems, temperature coupling, a compensatory effect of temperature on preference functions, is thought to maintain coordination between preferences and signals. However, although the carrier frequency of D. montana song is highly dependent on environmental temperature, there is no temperature coupling of the female preference function. We suggest that temperature coupling may often arise due to a common effect of temperature on song and preference, rather than be an advantageous characteristic whose function is to maintain coordination in temperature-affected communication systems.

Identification of X chromosomal restriction fragment length polymorphism markers and their use in a gene localization study in Drosophila virilis and D. littoralis.

We have identified six restriction fragment length polymorphism (RFLP) markers based on unique gene sequences on the X chromosome of Drosophila virilis and D. littoralis. The markers were localized by in situ hybridization on larval polytene chromosomes, and the conjugation of the X chromosomes of the two species was studied in salivary glands of interspecific hybrid female larvae. The gene arrangement of D. virilis and D. littoralis appeared to be very different at the proximal end of the X chromosome preventing recombination between RFLP markers located in this area. Simple quantitative trait loci (QTL) analysis showed that five of our marker genes (including nonA and Dmca1A, previously found to affect male courtship song in D. melanogaster) are linked with a gene(s) having a major effect on species differences in the male courtship song between D. virilis and D. littoralis. This shows that the song gene(s) may be located inside a large X-chromosomal inversion in D. littoralis (as previously suggested), but that it may also be located on an area between this inversion and the centromere, close to nonA and Dmca1A. Localization of this gene or gene complex will be continued with the aid of our newly identified RFLP markers by making interspecific crosses between D. virilis group species with more similar X chromosomes.

SEM search for sound production and sound perception organs in a variety of Drosophila species.

The discovery of new song types in Hawaiian Drosophila species has raised a question of alternative sound production and sound perception mechanisms in Drosophila. For example D. disjuncta males, which produce song with a carrier frequency of 5,000-6,000 Hz, do not vibrate their wings but only tremble them in a very small amplitude while singing. In our SEM search we found at the wing base of the flies of this species a special structure which could play a part in song production. Our search for auditory sense organs in Drosophila species revealed that the aristae on fly antennae are structurally not as simple as has been assumed. In Hawaiian D. planitibia subgroup species, the fly aristae are more regular than in other Drosophila species and they are also covered with small hairs. We did not find any sign of tympanic hearing organs in any of the studied species, but we discovered on fly legs and wings sensilla which could play a part in sound/vibration perception. Our SEM survey provides a good starting point in the search for alternative pathways for sound production and perception via transmission electron microscopy combined with behavioral studies and electrophysiological recordings in Hawaiian Drosophila species.

Localization of genes affecting species differences in male courtship song between Drosophila virilis and D. littoralis.

The males of six species of the Drosophila virilis group (including D. virilis) keep their wings extended while producing a train of sound pulses, where the pulses follow each other without any pause. The males of the remaining five species of the group produce only one sound pulse during each wing extension/vibration, which results in species-specific songs with long pauses (in D. littoralis about 300 ms) between successive sound pulses. Genetic analyses of the differences between the songs of D. virilis and D. littoralis showed that species-specific song traits are affected by genes on the X chromosome, and for the length of pause, also by genes on chromosomes 3 and 4. The X chromosomal genes having a major impact on pulse and pause length were tightly linked with white, apricot and notched marker genes located at the proximal third of the chromosome. A large inversion in D. littoralis, marked by notched, prevents more precise localization of these genes by classical crossing methods.