Latitudinal and cold-tolerance variation associate with DNA repeat-number variation in the hsr-omega RNA gene of Drosophila melanogaster

Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen

Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.

From Heterochromatin to Long Noncoding RNAs in Drosophila: Expanding the Arena of Gene Function and Regulation

Recent years have witnessed a remarkable interest in exploring the significance of pervasive noncoding transcripts in diverse eukaryotes. Classical cytogenetic studies using the Drosophila model system unraveled the perplexing attributes and "functions" of the "gene"-poor heterochromatin. Recent molecular studies in the fly model are likewise revealing the very diverse and significant roles played by long noncoding RNAs (lncRNAs) in development, gene regulation, chromatin organization, cell and nuclear architecture, etc. There has been a rapid increase in the number of identified lncRNAs, although a much larger number still remains unknown. The diversity of modes of actions and functions of the limited number of Drosophila lncRNAs, which have been examined, already reflects the profound roles of such RNAs in generating and sustaining the biological complexities of eukaryotes. Several of the known Drosophila lncRNAs originate as independent sense or antisense transcripts from promoter or intergenic, intronic, or 5'/3'-UTR regions, while many of them are independent genes that produce only lncRNAs or coding as well as noncoding RNAs. The different lncRNAs affect chromatin organization (local or large-scale pan-chromosomal), transcription, RNA processing/stability, or translation either directly through interaction with their target DNA sequences or indirectly by acting as intermediary molecules for specific regulatory proteins or may act as decoys/sinks, or storage sites for specific proteins or groups of proteins, or may provide a structural framework for the assembly of substructures in nucleus/cytoplasm. It is interesting that many of the "functions" alluded to heterochromatin in earlier cytogenetic studies appear to find correlates with the known subtle as well as far-reaching actions of the different small and long noncoding RNAs. Further studies exploiting the very rich and powerful genetic and molecular resources available for the Drosophila model are expected to unravel the mystery underlying the long reach of ncRNAs.

CRISPR-induced null alleles show that Frost protects Drosophila melanogaster reproduction after cold exposure

The ability to survive and reproduce after cold exposure is important in all kingdoms of life. However, even in a sophisticated genetic model system like Drosophila melanogaster, few genes have been identified as functioning in cold tolerance. The accumulation of the Frost (Fst) gene transcript increases after cold exposure, making it a good candidate for a gene that has a role in cold tolerance. However, despite extensive RNAi knockdown analysis, no role in cold tolerance has been assigned to Fst. CRISPR is an effective technique for completely knocking down genes, and less likely to produce off-target effects than GAL4-UAS RNAi systems. We have used CRISPR-mediated homologous recombination to generate Fst null alleles, and these Fst alleles uncovered a requirement for FST protein in maintaining female fecundity following cold exposure. However, FST does not have a direct role in survival following cold exposure. FST mRNA accumulates in the Malpighian tubules, and the FST protein is a highly disordered protein with a putative signal peptide for export from the cell. Future work is needed to determine whether FST is exported from the Malpighian tubules and directly interacts with female reproductive tissues post-cold exposure, or if it is required for other repair/recovery functions that indirectly alter energy allocation to reproduction.

Cold acclimation wholly reorganizes the Drosophila melanogaster transcriptome and metabolome

Cold tolerance is a key determinant of insect distribution and abundance, and thermal acclimation can strongly influence organismal stress tolerance phenotypes, particularly in small ectotherms like Drosophila. However, there is limited understanding of the molecular and biochemical mechanisms that confer such impressive plasticity. Here, we use high-throughput mRNA sequencing (RNA-seq) and liquid chromatography – mass spectrometry (LC-MS) to compare the transcriptomes and metabolomes of D. melanogaster acclimated as adults to warm (rearing) (21.5 °C) or cold conditions (6 °C). Cold acclimation improved cold tolerance and led to extensive biological reorganization: almost one third of the transcriptome and nearly half of the metabolome were differentially regulated. There was overlap in the metabolic pathways identified via transcriptomics and metabolomics, with proline and glutathione metabolism being the most strongly-supported metabolic pathways associated with increased cold tolerance. We discuss several new targets in the study of insect cold tolerance (e.g. dopamine signaling and Na⁺-driven transport), but many previously identified candidate genes and pathways (e.g. heat shock proteins, Ca²⁺ signaling, and ROS detoxification) were also identified in the present study, and our results are thus consistent with and extend the current understanding of the mechanisms of insect chilling tolerance.

Latitudinal clines in heat tolerance, protein synthesis rate and transcript level of a candidate gene in Drosophila melanogaster

The occurrence of climatic adaptation in Drosophila melanogaster is highlighted by the presence of latitudinal clines in several quantitative traits, particularly clines in adult heat knockdown tolerance that is higher in tropical populations. However the presence of latitudinal patterns in physiological characteristics that may underlie these traits have rarely been assessed. Protein synthesis has been implicated as an important physiological process that influences thermal tolerance, and this has not been examined in a clinal context. Here, we characterise latitudinal variation in D. melanogaster from eastern Australia in both adult heat knockdown tolerance and rates of protein synthesis following rearing at both 25 °C, approximating summer conditions, and 18 °C approximating winter development. We also examined clinal variation in the predominant nuclear transcript of the heat-inducible RNA gene hsr-omega, which has been implicated in regulating protein synthesis. We find significant clines in heat-hardened tolerance when cultured at both 18 and 25 °C - tolerance increased towards the low latitude tropics. Rates of protein synthesis measured in ovarian tissue also associated negatively with latitude, however the presence of the clines depended on rearing temperature and heat stress conditions. Finally, omega-n levels measured without heat stress showed a positive linear cline. When measured after a mild heat stress higher levels of omega-n were detected and the clinal pattern became parabolic - mid-latitude populations had lower levels of the transcript. While congruent latitudinal trends were detected for these three traits, only a low level of positive association was detected between protein synthesis and thermal tolerance providing little evidence that these traits are related at the level of cellular physiology. However the new clinal patterns of protein synthesis and hsr-omega variation suggest that these variables exert important influences on traits involved with latitudinal climatic adaptation.

Dietary cholesterol modulates pathogen blocking by Wolbachia

The bacterial endosymbiont Wolbachia pipientis protects its hosts from a range of pathogens by limiting their ability to form infections inside the insect. This “pathogen blocking” could be explained by innate immune priming by the symbiont, competition for host-derived resources between pathogens and Wolbachia, or the direct modification of the cell or cellular environment by Wolbachia. Recent comparative work in Drosophila and the mosquito Aedes aegypti has shown that an immune response is not required for pathogen blocking, implying that there must be an additional component to the mechanism. Here we have examined the involvement of cholesterol in pathogen blocking using a system of dietary manipulation in Drosophila melanogaster in combination with challenge by Drosophila C virus (DCV), a common fly pathogen. We observed that flies reared on cholesterol-enriched diets infected with the Wolbachia strains wMelPop and wMelCS exhibited reduced pathogen blocking, with viral-induced mortality occurring 2–5 days earlier than flies reared on Standard diet. This shift toward greater virulence in the presence of cholesterol also corresponded to higher viral copy numbers in the host. Interestingly, an increase in dietary cholesterol did not have an effect on Wolbachia density except in one case, but this did not directly affect the strength of pathogen blocking. Our results indicate that host cholesterol levels are involved with the ability of Wolbachia-infected flies to resist DCV infections, suggesting that cholesterol contributes to the underlying mechanism of pathogen blocking.

Wolbachia pipientis is an intracellular bacterium that naturally infects many insect species. These bacteria can block the replication and dissemination of a variety of pathogens that coinfect the insect. In mosquitoes this effect applies to viruses including dengue and the parasite that causes malaria. This makes Wolbachia a promising method of controlling a number of insect-transmitted diseases of humans. Some studies suggest that pathogen blocking is due to Wolbachia stimulating the host innate immune system, however in the fly Drosophila melanogaster, blocking occurs in the absence of such an immune effect. This suggests that there is an additional component to the mechanism. Host nutrients such as cholesterol are a limiting factor in virus replication and may serve as a source of competition that underlies pathogen blocking. Here we show that competition over cholesterol contributes to the pathogen blocking effect, with higher available cholesterol levels causing less effective blocking, and increased viral titres in Wolbachia-infected D. melanogaster challenged with the pathogenic Drosophila C virus. These results suggest that competition over cholesterol is an important part of the pathogen blocking mechanism.

Long non-coding RNAs coordinate cellular responses to stress

Following the initial discovery of the heat shock RNA omega (hsrω) gene of Drosophila melanogaster to be non-coding (nc) and also inducible by cell stress, other stress-inducible long non-coding RNAs (lncRNA) have been described in diverse organisms. In view of the rapid sequence divergence of lncRNAs, present knowledge of stress trasncriptome is limited and fragmented. Several known stress-related lncRNAs, associated with specific nuclear speckled domains or nucleolus, provide structural base for sequestering diverse RNA-processing/regulatory proteins. Others have roles in transcriptional or translational inhibition during stress or in signaling pathways; functions of several other lncRNAs are not yet known. Most stress-related lncRNAs act primarily by modulating activity of the proteins to which they bind or by sequestering specific sets of proteins away from the active pool. A common emerging theme is that a given lncRNA targets one or more protein/s with key role/s in the cascade of events triggered by the stress and therefore has a widespread integrative effect. Since proteins associate with RNA through short sequence motifs, the overall base sequence of functionally similar ncRNAs is often not conserved except for specific motifs. The rapid evolvability of ncRNA sequences provides elegant modules for adaptability to changing environment as binding of one or the other protein to ncRNA can alter its structure and functions in distinct ways. Thus the stress-related lncRNAs act as hubs in the cellular networks to coordinate activities of the members within and between different networks to maintain cellular homeostasis for survival or to trigger cell death.

Genetics of climate change adaptation

The rapid rate of current global climate change is having strong effects on many species and, at least in some cases, is driving evolution, particularly when changes in conditions alter patterns of selection. Climate change thus provides an opportunity for the study of the genetic basis of adaptation. Such studies include a variety of observational and experimental approaches, such as sampling across clines, artificial evolution experiments, and resurrection studies. These approaches can be combined with a number of techniques in genetics and genomics, including association and mapping analyses, genome scans, and transcription profiling. Recent research has revealed a number of candidate genes potentially involved in climate change adaptation and has also illustrated that genetic regulatory networks and epigenetic effects may be particularly relevant for evolution driven by climate change. Although genetic and genomic data are rapidly accumulating, we still have much to learn about the genetic architecture of climate change adaptation.

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.

Biogeography of species richness gradients: linking adaptive traits, demography and diversification

Here we review how adaptive traits contribute to the emergence and maintenance of species richness gradients through their influence on demographic and diversification processes. We start by reviewing how demographic dynamics change along species richness gradients. Empirical studies show that geographical clines in population parameters and measures of demographic variability are frequent along latitudinal and altitudinal gradients. Demographic variability often increases at the extremes of regional species richness gradients and contributes to shape these gradients. Available studies suggest that adaptive traits significantly influence demographic dynamics, and set the limits of species distributions. Traits related to thermal tolerance, resource use, phenology and dispersal seem to play a significant role. For many traits affecting demography and/or diversification processes, complex mechanistic approaches linking genotype, phenotype and fitness are becoming progressively available. In several taxa, species can be distributed along adaptive trait continuums, i.e. a main axis accounting for the bulk of inter-specific variation in some correlated adaptive traits. It is shown that adaptive trait continuums can provide useful mechanistic frameworks to explain demographic dynamics and diversification in species richness gradients. Finally, we review the existence of sequences of adaptive traits in phylogenies, the interactions of adaptive traits and community context, the clinal variation of traits across geographical gradients, and the role of adaptive traits in determining the history of dispersal and diversification of clades. Overall, we show that the study of demographic and evolutionary mechanisms that shape species richness gradients clearly requires the explicit consideration of adaptive traits. To conclude, future research lines and trends in the field are briefly outlined.

Genetic variation in senescence marker protein-30 is associated with natural variation in cold tolerance in Drosophila

A comprehensive understanding of the genetic basis of phenotypic adaptation in nature requires the identification of the functional allelic variation underlying adaptive phenotypes. The manner in which organisms respond to temperature extremes is an adaptation in many species. In the current study, we investigate the role of molecular variation in senescence marker protein-30 (Smp-30) on natural phenotypic variation in cold tolerance in Drosophila melanogaster. Smp-30 encodes a product that is thought to be involved in the regulation of Ca2+ ion homeostasis and has been shown previously to be differentially expressed in response to cold stress. Thus, we sought to assess whether molecular variation in Smp-30 was associated with natural phenotypic variation in cold tolerance in a panel of naturally derived inbred lines from a population in Raleigh, North Carolina. We identified four non-coding polymorphisms that were strongly associated with natural phenotypic variation in cold tolerance. Interestingly, two polymorphisms that were in close proximity to one another (2 bp apart) exhibited opposite phenotypic effects. Consistent with the maintenance of a pair of antagonistically acting polymorphisms, tests of molecular evolution identified a significant excess of maintained variation in this region, suggesting balancing selection is acting to maintain this variation. These results suggest that multiple mutations in non-coding regions can have significant effects on phenotypic variation in adaptive traits within natural populations, and that balancing selection can maintain polymorphisms with opposite effects on phenotypic variation.

Connecting genetic variation to phenotypic clines

From early allozyme work to recent genome-wide scans, many studies have reported associations between molecular markers and latitude. These geographic patterns are tantalizing because they hint at the possibility of identifying specific mutations responsible for climatic adaptation. Unfortunately, few studies have done so because these exciting first glances often prove extremely challenging to follow up. Many difficulties can hinder connecting genetic and phenotypic variation in this context, and without such links, distinguishing the action of spatially varying selection from the other evolutionary processes capable of generating these patterns can be quite thorny. Nevertheless, two papers in this issue report excellent progress in overcoming these obstacles and provide persuasive evidence supporting the involvement of specific natural variants in clinal adaptation of Drosophila melanogaster populations (Fig. 1). In the first paper, Paaby et al. (2010) describe replicated allele frequency clines for a coding polymorphism in the Insulin-like Receptor (InR) gene on two continents, findings that strongly point to selection acting at this locus and that likely reflect life history adaptation. McKechnie et al. (2010) report compelling functional evidence that cis-regulatory variation in the Dca (drosophila cold acclimation) gene contributes to an adaptive cline in wing size. Notably, these papers employ largely alternative and complementary approaches, and together they exemplify how diverse strategies may be interwoven to draw convincing connections between genotype, phenotype, and evolutionary process.

Characterization of the shsp genes in Drosophila buzzatii and association between the frequency of Valine mutations in hsp23 and climatic variables along a longitudinal gradient in Australia

The small heat shock gene (shsp) cluster of Drosophila buzzatii was sequenced and the gene order and DNA sequence were compared with those of the shsps in Drosophila melanogaster. The D. buzzatii shsp cluster contains an inversion and a duplication of hsp26. A phylogenetic tree was constructed based on hsp26 genes from several Drosophila species of the Sophophora and Drosophila subgenera. The tree shows first a separation of the Sophophora and the Drosophila subgenera and then the Drosophila subgenus is divided into the Hawaiian Drosophila and the repleta/virilis groups. Only the latter contain a duplicated hsp26. Comparing the gene organisation of the shsp cluster shows that all the Drosophila subgenus species contain the inversion. Putative heat shock elements (HSE) were found in the promoters of all the shsp and putative regulator elements for tissue specific expression were found in the promoter of hsp23, hsp27 and one of the hsp26 genes. hsp23 was found to be polymorphic for four non-synonymous changes that all lead to exchange of a Valine. The duplicated hsp26 gene in D. buzzatii (phsp26) was polymorphic for two non-synonymous changes. The allele frequencies of these variants were determined in nine D. buzzatii populations covering most of its distribution in Australia using high-resolution melting curves. The allele frequencies of one of the hsp23 variants showed a significant linear regression with longitude and the pooled frequency of the four Valine changes of hsp23 in the nine populations showed a significant linear regression with longitude and with a composite measure of climatic variables.

Protein synthesis rates in Drosophila associate with levels of the hsr-omega nuclear transcript

Transcripts of the Drosophila hsr-omega gene are known to interact with RNA processing factors and ribosomes and are postulated to aid in co-ordinating nuclear and cytoplasmic activities particularly in stressed cells. However, the significance of these interactions for physiological processes and in turn for whole-organism fitness remains an open question. Because hsr-omega's cellular expression characteristics suggest it may influence protein synthesis, and because both genotypic and expression variation of hsr-omega have been associated with thermotolerance, we characterised 30 lines for variation in the rates of protein synthesis, measured in ovarian tissues, both before and after a mild heat shock, and for basal levels of the two main hsr-omega transcripts, omega-n and omega-c. As expected, the mild heat shock reduced protein synthesis rates. Large variation occurred among lines in levels of omega-n which was negatively associated with rates of basal protein synthesis--a result that supports the model for the cellular function of omega-n. Furthermore, omega-n levels were associated with hsr-omega genotype of the line parents. Little variation occurred among lines for omega-c levels and no associations were detected with protein synthesis or genotype. Since protein synthesis is a fundamental process for growth and development, we characterised the lines for several life-history traits; however, no associations with protein synthesis, omega-n or omega-c levels were detected. Our results are consistent with the idea that natural variation in hsr-omega expression influence rates of protein synthesis in this species.