Genome differentiation of Drosophila melanogaster from a microclimate contrast in Evolution Canyon, Israel

Scattered differentiation of unlinked loci across the genome underlines ecological divergence of the selfing grass Brachypodium stacei

Ecological divergence without geographic isolation, as an early speciation process that may lead finally to reproductive isolation through natural selection, remains a captivating topic in evolutionary biology. However, the pattern of genetic divergence underlying this process across the genome may vary between species and mating systems. Here, we present evidence that Brachypodium stacei, an annual and highly selfing grass model species, has undergone sympatric ecological divergence without geographic isolation. Genomic, transcriptomic, and metabolomic analyses together with lab experiments mimicking the two opposite environmental conditions suggest that diploid B. stacei populations have diverged sympatrically in two slopes characterized by distinct biomes at Evolution Canyon I (ECI), Mount Carmel, Israel. Despite ongoing gene flow, primarily facilitated by seed dispersal, the level of gene flow has progressively decreased over time. This local adaptation involves the scattered divergence of many unlinked loci across the total genome that include both coding genes and noncoding regions. Additionally, we have identified significant differential expressions of genes related to the ABA signaling pathway and contrasting metabolome composition between the arid- vs. forest-adapted B. stacei populations in ECI. These results suggest that multiple small loci involved in environmental responses act additively to account for ecological adaptations by this selfing species in contrasting environments.

Comparative gene retention analysis in barley, wild emmer, and bread wheat pangenome lines reveals factors affecting gene retention following gene duplication

BACKGROUND

Gene duplication is a prevalent phenomenon and a major driving force underlying genome evolution. The process leading to the fixation of gene duplicates following duplication is critical to understand how genome evolves but remains fragmentally understood. Most previous studies on gene retention are based on gene duplicate analyses in single reference genome. No population-based comparative gene retention analysis has been performed to date.

RESULTS

Taking advantage of recently published genomic data in Triticeae, we dissected a divergent homogentisate phytyltransferase (HPT2) lineage caught in the middle stage of gene fixation following duplication. The presence/absence of HPT2 in barley (diploid), wild emmer (tetraploid), and bread wheat (hexaploid) pangenome lines appears to be associated with gene dosage constraint and environmental adaption. Based on these observations, we adopted a phylogeny-based orthology inference approach and performed comparative gene retention analyses across barley, wild emmer, and bread wheat. This led to the identification of 326 HPT2-pattern-like genes at whole genome scale, representing a pool of gene duplicates in the middle stage of gene fixation. Majority of these HPT2-pattern-like genes were identified as small-scale duplicates, such as dispersed, tandem, and proximal duplications. Natural selection analyses showed that HPT2-pattern-like genes have experienced relaxed selection pressure, which is generally accompanied with partial positive selection and transcriptional divergence. Functional enrichment analyses showed that HPT2-pattern-like genes are over-represented with molecular-binding and defense response functions, supporting the potential role of environmental adaption during gene retention. We also observed that gene duplicates from larger gene family are more likely to be lost, implying a gene dosage constraint effect. Further comparative gene retention analysis in barley and bread wheat pangenome lines revealed combined effects of species-specific selection and gene dosage constraint.

CONCLUSIONS

Comparative gene retention analyses at the population level support gene dosage constraint, environmental adaption, and species-specific selection as three factors that may affect gene retention following gene duplication. Our findings shed light on the evolutionary process leading to the retention of newly formed gene duplicates and will greatly improve our understanding on genome evolution via duplication.

Selection of p53 pathway in adaptive evolution and reproductive isolation in incipient sympatric speciation of Drosophila at Evolution Canyon

Sympatric speciation (SS) refers to the origin of new species within a freely breeding population. The ‘Evolution Canyon’ (EC) in Israel is a natural microsite model for SS of species across phylogenies from viruses and bacteria to mammals, adapting to, and speciating in, interslope microclimates. The cosmopolitan Drosophila melanogaster at EC I, Mount Carmel, is undergoing incipient SS in response to sharply divergent interslope microclimate stresses, including solar radiation, temperature, humidity and pathogenicity. We demonstrated here a selective interslope divergence of single nucleotide polymorphism (SNP) distribution in the Drosophila p53 pathway. This involves a total of 71 genes, which are associated with DNA repair, heat response, and fungal and bacterial resistant pathways. This distribution pattern links the previously observed thermotolerance and ageing divergence of D. melanogaster between the opposite canyon slopes: the south-facing slope (SFS, or African slope: tropical, savannoid and dry) and the abutting north-facing slope (NFS, or European slope; temperate, forested, cool and humid). The genes with interslope-significant differential SNPs link the p53 pathway with pathways related to the responses to microclimates through protein-protein interaction. Moreover, for the first time we provide evidence that the p53 pathway is linked to reproductive isolation, and is thus actively participating in incipient SS of D. melanogaster. This is the first demonstration of a link between the p53 pathway and reproductive isolation, thereby contributing to adaptive incipient sympatric speciation.

Conserved ecophysiology despite disparate microclimatic conditions in a gecko

Microscale differences in the habitats organisms occupy can influence selection regimes and promote intraspecific variation of traits. Temperature-dependent traits can be locally adapted to climatic conditions or be highly conserved and insensitive to directional selection under all but the most extreme regimes, and thus be similar across populations. The opposing slopes of Nahal Oren canyon in the Carmel Mountains, Israel, are strikingly different: the south-facing slope receives intensive solar radiation, is hot and supports mostly annual vegetation, whereas the north-facing slope is ~10°C cooler, more humid, and supports Mediterranean woodland. We examined whether these differences manifest in the thermal physiology of a common gecko species Ptyodactylus guttatus in controlled laboratory conditions. We predicted that geckos from the hotter south-facing slope would prefer higher temperatures, have faster gut passage times, lower metabolic and evaporative water loss rates, and start diel activity earlier compared with north-facing slope conspecifics. Contrary to these predictions, there were no differences between any of the ecophysiological traits in geckos from the opposing slopes. Nevertheless, our data showed that individuals from the north-facing slope were generally more active in earlier hours of the afternoon compared with south-facing individuals. We suggest that P. guttatus individuals disperse between the slopes and either gene-flow or behavioral plasticity deter local adaptation, resulting in similar physiological traits. Perhaps a stronger contrast in climatic conditions and a stronger barrier are needed to result in interpopulation divergence in temperature-dependent traits.

False and true positives in arthropod thermal adaptation candidate gene lists

Genome-wide studies are prone to false positives due to inherently low priors and statistical power. One approach to ameliorate this problem is to seek validation of reported candidate genes across independent studies: genes with repeatedly discovered effects are less likely to be false positives. Inversely, genes reported only as many times as expected by chance alone, while possibly representing novel discoveries, are also more likely to be false positives. We show that, across over 30 genome-wide studies that reported Drosophila and Daphnia genes with possible roles in thermal adaptation, the combined lists of candidate genes and orthologous groups are rapidly approaching the total number of genes and orthologous groups in the respective genomes. This is consistent with the expectation of high frequency of false positives. The majority of these spurious candidates have been identified by one or a few studies, as expected by chance alone. In contrast, a noticeable minority of genes have been identified by numerous studies with the probabilities of such discoveries occurring by chance alone being exceedingly small. For this subset of genes, different studies are in agreement with each other despite differences in the ecological settings, genomic tools and methodology, and reporting thresholds. We provide a reference set of presumed true positives among Drosophila candidate genes and orthologous groups involved in response to changes in temperature, suitable for cross-validation purposes. Despite this approach being prone to false negatives, this list of presumed true positives includes several hundred genes, consistent with the "omnigenic" concept of genetic architecture of complex traits.

Transcriptome-Wide Patterns of the Genetic and Expression Variations in Two Sympatric Schizothoracine Fishes in a Tibetan Plateau Glacier Lake

Sympatric speciation remains a central focus of evolutionary biology. Although some evidence shows speciation occurring in this way, little is known about the gene expression evolution and the characteristics of population genetics as species diverge. Two closely related Gymnocypris fish (Gymnocypris chui and Gymnocypris scleracanthus), which come from a small glacier lake in the Tibetan Plateau, Lake Langcuo, exist a possible incipient sympatric adaptive ecological speciation. We generated large amounts of RNA-Seq data from multiple individuals and tissues from each of the two species and compared gene expression patterns and genetic polymorphisms between them. Ordination analysis separated samples by organ rather than by species. The degree of expression difference between organs within and between species was different. Phylogenetic analyses indicated that the two closely related taxa formed a monophyletic complex. Population structure analysis displayed two distinctly divergent clusters of G. chui and G. scleracanthus populations. By contrast, G. scleracanthus population genetic diversity is higher than that of G. chui. Considerable sites of the two populations were differentiated with a coefficient of F_ST = 0.25–0.50, implying that a small proportion of loci nevertheless exhibited deep divergence in two comparisons. Concomitantly, putatively selected genes during speciation revealed functional categories are enriched in bone morphogenesis, cell growth, neurogenetics, enzyme activity, and binding activity in G. chui population. In contrast, nutrition and localization were highlighted in G. scleracanthus. Collectively, morphological traits and dietary preference combine with genetic variation and expression variation, probably contributed to the incipient speciation of two sympatric populations.

Sympatric speciation of wild emmer wheat driven by ecology and chromosomal rearrangements

In plants, the mechanism for ecological sympatric speciation (SS) is little known. Here, after ruling out the possibility of secondary contact, we show that wild emmer wheat, at the microclimatically divergent microsite of "Evolution Canyon" (EC), Mt. Carmel, Israel, underwent triple SS. Initially, it split following a bottleneck of an ancestral population, and further diversified to three isolated populations driven by disruptive ecological selection. Remarkably, two postzygotically isolated populations (SFS1 and SFS2) sympatrically branched within an area less than 30 m at the tropical hot and dry savannoid south-facing slope (SFS). A series of homozygous chromosomal rearrangements in the SFS1 population caused hybrid sterility with the SFS2 population. We demonstrate that these two populations developed divergent adaptive mechanisms against severe abiotic stresses on the tropical SFS. The SFS2 population evolved very early flowering, while the SFS1 population alternatively evolved a direct tolerance to irradiance by improved ROS scavenging activity that potentially accounts for its evolutionary fate with unstable chromosome status. Moreover, a third prezygotically isolated sympatric population adapted on the abutting temperate, humid, cool, and forested north-facing slope (NFS), separated by 250 m from the SFS wild emmer wheat populations. The NFS population evolved multiple resistant loci to fungal diseases, including powdery mildew and stripe rust. Our study illustrates how plants sympatrically adapt and speciate under disruptive ecological selection of abiotic and biotic stresses.

Aggression and courtship differences found in Drosophila melanogaster from two different microclimates at Evolution Canyon, Israel

Aggression and courtship behavior were examined of wild Drosophila melanogaster flies isolated from two contrasting microclimates found at Evolution Canyon in Mt. Carmel, Israel: an African-like dry tropical Slope (AS) and a European-like humid temperate Slope (ES), separated by 250 meters. Studies were carried out to ask whether behavioral differences existed between the two populations obtained from opposite slopes with divergent microclimates in Israel. First, we measured and compared intraslope aggression between same sex fly pairings collected from the same slope. Both male and female flies displayed similar fighting abilities from both slopes. ES males, however, from the humid biome, showed a tendency to lunge more per aggressive encounter, compared with AS males from the dry biome. Next, we tested interslope aggression by pairing flies from opposite slopes. ES males displayed higher numbers of lunges, and won more fights against their AS opponents. We also observed enhanced courtship performances in ES compared to AS males. The fighting and courtship superiority seen in ES males could reinforce fitness and pre-mating reproductive isolation mechanisms that underlie incipient sympatric speciation. This may support an evolutionary advantage of adaptively divergent fruit fly aggression phenotypes from different environments.

Rapid seasonal evolution in innate immunity of wild Drosophila melanogaster

Understanding the rate of evolutionary change and the genetic architecture that facilitates rapid adaptation is a current challenge in evolutionary biology. Comparative studies show that genes with immune function are among the most rapidly evolving genes across a range of taxa. Here, we use immune defence in natural populations of Drosophila melanogaster to understand the rate of evolution in natural populations and the genetics underlying rapid change. We probed the immune system using the natural pathogens Enterococcus faecalis and Providencia rettgeri to measure post-infection survival and bacterial load of wild D. melanogaster populations collected across seasonal time along a latitudinal transect along eastern North America (Massachusetts, Pennsylvania and Virginia). There are pronounced and repeatable changes in the immune response over the approximately 10 generations between spring and autumn collections, with a significant but less distinct difference observed among geographical locations. Genes with known immune function are not enriched among alleles that cycle with seasonal time, but the immune function of a subset of seasonally cycling alleles in immune genes was tested using reconstructed outbred populations. We find that flies containing seasonal alleles in Thioester-containing protein 3 (Tep3) have different functional responses to infection and that epistatic interactions among seasonal Tep3 and Drosomycin-like 6 (Dro6) alleles underlie the immune phenotypes observed in natural populations. This rapid, cyclic response to seasonal environmental pressure broadens our understanding of the complex ecological and genetic interactions determining the evolution of immune defence in natural populations.

A Large Panel of Drosophila simulans Reveals an Abundance of Common Variants

The rapidly expanding availability of large NGS data sets provides an opportunity to investigate population genetics at an unprecedented scale. Drosophila simulans is the sister species of the model organism Drosophila melanogaster, and is often presumed to share similar demographic history. However, previous population genetic and ecological work suggests very different signatures of selection and demography. Here, we sequence a new panel of 170 inbred genotypes of a North American population of D. simulans, a valuable complement to the DGRP and other D. melanogaster panels. We find some unexpected signatures of demography, in the form of excess intermediate frequency polymorphisms. Simulations suggest that this is possibly due to a recent population contraction and selection. We examine the outliers in the D. simulans genome determined by a haplotype test to attempt to parse the contribution of demography and selection to the patterns observed in this population. Untangling the relative contribution of demography and selection to genomic patterns of variation is challenging, however, it is clear that although D. melanogaster was thought to share demographic history with D. simulans different forces are at work in shaping genomic variation in this population of D. simulans.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Regulation of gene expression and RNA editing in Drosophila adapting to divergent microclimates

Determining the mechanisms by which a species adapts to its environment is a key endeavor in the study of evolution. In particular, relatively little is known about how transcriptional processes are fine-tuned to adjust to different environmental conditions. Here we study Drosophila melanogaster from 'Evolution Canyon' in Israel, which consists of two opposing slopes with divergent microclimates. We identify several hundred differentially expressed genes and dozens of differentially edited sites between flies from each slope, correlate these changes with genetic differences, and use CRISPR mutagenesis to validate that an intronic SNP in prominin regulates its editing levels. We also demonstrate that while temperature affects editing levels at more sites than genetic differences, genetically regulated sites tend to be less affected by temperature. This work shows the extent to which gene expression and RNA editing differ between flies from different microclimates, and provides insights into the regulation responsible for these differences.

Genomic signatures of local adaptation in the Drosophila immune response

As environments and pathogen landscapes shift, host defenses must evolve to remain effective. Due to this selection pressure, among-species comparisons of genetic sequence data often find immune genes to be among the fastest evolving genes across the genome. The full extent and nature of these immune adaptations, however, remain largely unexplored. In a recent study, we analyzed patterns of selection within distinct components of the Drosophila melanogaster immune pathway. While we found evidence of positive selection within some immune processes, immune genes were not universally characterized by signatures of strong selection. On the contrary, we even found that some immune functions show greater than expected constraint. Overall these results highlight 2 major factors that appear to play an outsize role in determining a gene's evolutionary rate: the type of pathogen the gene targets and the gene's position within the immune network. These results join a growing body of literature that highlight the complexity of immune adaptation. Rather than there being uniformly strong selection across all immune genes, a combination of pathogen-specificity and host genetic constraints appear to play key roles in determining each immune gene's individual evolutionary trajectory.

Adaptive differentiation coincides with local bioclimatic conditions along an elevational cline in populations of a lichen-forming fungus

Background

Many fungal species occur across a variety of habitats. Particularly lichens, fungi forming symbioses with photosynthetic partners, have evolved remarkable tolerances for environmental extremes. Despite their ecological importance and ubiquity, little is known about the genetic basis of adaption in lichen populations. Here we studied patterns of genome-wide differentiation in the lichen-forming fungus Lasallia pustulata along an altitudinal gradient in the Mediterranean region. We resequenced six populations as pools and identified highly differentiated genomic regions. We then detected gene-environment correlations while controlling for shared population history and pooled sequencing bias, and performed ecophysiological experiments to assess fitness differences of individuals from different environments.

Results

We detected two strongly differentiated genetic clusters linked to Mediterranean and temperate-oceanic climate, and an admixture zone, which coincided with the transition between the two bioclimates. High altitude individuals showed ecophysiological adaptations to wetter and more shaded conditions. Highly differentiated genome regions contained a number of genes associated with stress response, local environmental adaptation, and sexual reproduction.

Conclusions

Taken together our results provide evidence for a complex interplay between demographic history and spatially varying selection acting on a number of key biological processes, suggesting a scenario of ecological speciation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0929-8) contains supplementary material, which is available to authorized users.

Survey of Global Genetic Diversity Within the Drosophila Immune System

Numerous studies across a wide range of taxa have demonstrated that immune genes are routinely among the most rapidly evolving genes in the genome. This observation, however, does not address what proportion of immune genes undergo strong selection during adaptation to novel environments. Here, we determine the extent of very recent divergence in genes with immune function across five populations of Drosophila melanogaster and find that immune genes do not show an overall trend of recent rapid adaptation. Our population-based approach uses a set of carefully matched control genes to account for the effects of demography and local recombination rate, allowing us to identify whether specific immune functions are putative targets of strong selection. We find evidence that viral-defense genes are rapidly evolving in Drosophila at multiple timescales. Local adaptation to bacteria and fungi is less extreme and primarily occurs through changes in recognition and effector genes rather than large-scale changes to the regulation of the immune response. Surprisingly, genes in the Toll pathway, which show a high rate of adaptive substitution between the D. melanogaster and D. simulans lineages, show little population differentiation. Quantifying the flies for resistance to a generalist Gram-positive bacterial pathogen, we found that this genetic pattern of low population differentiation was recapitulated at the phenotypic level. In sum, our results highlight the complexity of immune evolution and suggest that Drosophila immune genes do not follow a uniform trajectory of strong directional selection as flies encounter new environments.

A Variable Genetic Architecture of Melanic Evolution in Drosophila melanogaster

Unraveling the genetic architecture of adaptive phenotypic divergence is a fundamental quest in evolutionary biology. In Drosophila melanogaster, high-altitude melanism has evolved in separate mountain ranges in sub-Saharan Africa, potentially as an adaptation to UV intensity. We investigated the genetic basis of this melanism in three populations using a new bulk segregant analysis mapping method. We identified 19 distinct QTL regions from nine mapping crosses, with several QTL peaks overlapping between two or all populations, and yet different crosses involving the same melanic population commonly yielded distinct QTL. The strongest QTL often overlapped well-known pigmentation genes, but we typically did not find wide signals of genetic differentiation (F_ST) between lightly and darkly pigmented populations at these genes. Instead, we found small numbers of highly differentiated SNPs at the probable causative genes. A simulation analysis showed that these patterns of polymorphism were consistent with selection on standing genetic variation. Overall, our results suggest that, even for potentially simpler traits like pigmentation, the complexity of adaptive trait evolution poses important challenges for QTL mapping and population genetic analysis.

Sympatric speciation of spiny mice, Acomys, unfolded transcriptomically at Evolution Canyon, Israel

Spiny mice, Acomys cahirinus, colonized Israel 30,000 y ago from dry tropical Africa and inhabited rocky habitats across Israel. Earlier, we had shown by mtDNA that A. cahirinus incipiently sympatrically speciates at Evolution Canyon I (EC I) in Mount Carmel, Israel because of microclimatic interslope divergence. The EC I microsite consists of a dry and hot savannoid "African" slope (AS) and an abutting humid and cool-forested "European" slope (ES). Here, we substantiate incipient SS in A. cahirinus at EC I based on the entire transcriptome, showing that multiple slope-specific adaptive complexes across the transcriptome result in two divergent clusters. Tajima's D distribution of the abutting Acomys interslope populations shows that the ES population is under stronger positive selection, whereas the AS population is under balancing selection, harboring higher genetic polymorphisms. Considerable sites of the two populations were differentiated with a coefficient of FST = 0.25-0.75. Remarkably, 24 and 37 putatively adaptively selected genes were detected in the AS and ES populations, respectively. The AS genes involved DNA repair, growth arrest, neural cell differentiation, and heat-shock proteins adapting to the local AS stresses of high solar radiation, drought, and high temperature. In contrast, the ES genes involved high ATP associated with energetics stress. The sharp ecological interslope divergence led to strong slope-specific selection overruling the interslope gene flow. Earlier tests suggested slope-specific mate choice. Habitat interslope-adaptive selection across the transcriptome and mate choice substantiate sympatric speciation (SS), suggesting its prevalence at EC I and commonality in nature.

Lack of population differentiation patterns of previously identified putatively adaptive transposable element insertions at microgeographic scales

Background

Transposable elements (TEs) play an important role in genome function and evolution. It has been shown that TEs are a considerable source of adaptive changes in the genome of Drosophila melanogaster. Specifically, footprints of selection at the DNA level, the presence of population differentiation patterns across environmental gradients, and detailed mechanistic and fitness analyses of a few candidate adaptive TEs pointed to the role of TEs in environmental adaptation. However, whether the population differentiation patterns observed at large geographic scales can be replicated at a microgeographic scale has never been assessed before.

Results

In this work, we explored the population patterns of putatively adaptive TEs at a micro-spatial scale level. We compared the frequencies of TEs, previously identified as putatively adaptive and putatively neutral, in populations collected in opposite slopes of the Evolution Canyon at Mt. Carmel in Israel separated by 200 m on average. However, the differentiation patterns previously observed across large geographic distances (2000–2200 km) were not replicated at the microscale level of the Evolution Canyon populations.

Conclusion

TE insertions previously associated with D. melanogaster environmental adaptation at a macro scale level do not play such a role at the microscale level of the Evolution Canyon populations. However, these results do not exclude a role of TEs in microgeographic adaptation because the dataset analyzed in this work is restricted to TEs identified in a single North American strain and as such is highly biased and incomplete.

Reviewers

This article was reviewed by Eugene Koonin, Limsoon Wong and Fyodor Kondrashov.

Putative adaptive inter-slope divergence of transposon frequency in fruit flies (Drosophila melanogaster) at "Evolution Canyon", Mount Carmel, Israel

Background

The current analysis of transposon elements (TE) in Drosophila melanogaster at Evolution Canyon, (EC), Israel, is based on data and analysis done by our collaborators (Drs. J. Gonzalez, J. Martinez and W. Makalowski, this issue). They estimated the frequencies of 28 TEs (transposon elements) in fruit flies (D. melanogaster) from the ecologically tropic, hot, and dry south-facing slope (SFS) or “African” slope (AS) of EC and compared it with the TE frequencies on the temperate-cool and humid north-facing slope (NFS) or “European” slope (ES), separated, on average, by 250 m. The flies were sampled from two stations on each slope. We received their results, including the frequencies of each TE on each slope, and the probabilities of the statistical analyses (G-tests) of each TE separately. We continued the analysis of the inter-slope differences of the frequencies of the TEs, and based our different conclusions on that analysis and on the difference between micro (=EC) and macro (2000 km.) comparisons [Gonzalez et al. 2015 doi:10.1186/s13062-015-0075-4].

Results

Our collaborators based all their conclusions on the non-significant results of each of the individual tests of the 28 TEs. We analysed also the distribution of the TE differences between the slopes, based on their results. Thirteen TEs were more frequent on the SFS, 11 were more frequent on the NFS, and four had equal frequencies. Because of the equalizing effect of the ongoing migration, only small and temporary differences between the slopes (0 – 0.06) were regarded by us as random fluctuations (drift). Three TEs were intermediate (0.08-0.09) and await additional research. The 11 TEs with large frequency differences (0.12 – 0.22) were regarded by us as putative adaptive TEs, because the equalizing power of ongoing migration will eliminate random large differences. Five of them were higher on the SFS and six were higher on the NFS. Gaps in the distribution of the differences distinguished between the large and small differences. The large gap among the 11 TEs favored on the NFS was significant and supports our rejection of drift as the only explanation of the distribution of the slope differences. The gaps in the distribution of the differences separated the putative TEs with strong enough selection from those TEs that couldn't overrule the migration.

The results are compared and contrasted with the directional effect of the frequencies of the same TEs in the study of global climatic comparisons across thousands of kilometers. From the 11 putative adaptive TEs in the local “Evolution Canyon,” six differentiate in the same direction as in the continental comparisons and four in the opposite direction. One TE, FBti0019144, differentiated in EC in the same direction as in Australia and in the opposite direction to that of North America.

Conclusions

We presume that the major divergent evolutionary driving force at the local EC microsite is natural selection overruling gene flow. Therefore, after we rejected drift as an explanation of all the large slope differences, we regarded them as putatively adaptive. In order to substantiate the individual TE adaptation, we need to increase the sample sizes and reveal the significant adaptive TEs.

The comparison of local and global studies show only partial similarity in the adaptation of the TEs, because of the dryness of the ecologically tropical climate in EC, in contrast to the wet tropical climate in the global compared climates. Moreover, adaptation of a TE may be expressed only in part of the time and specific localities.

Reviewers

Reviewed by Eugene Koonin, Limsoon Wong and Fyodor Kondrashov. For the full reviews, please go to the Reviewers’ comments section.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0074-5) contains supplementary material, which is available to authorized users.

Identifying signatures of sexual selection using genomewide selection components analysis

Sexual selection must affect the genome for it to have an evolutionary impact, yet signatures of selection remain elusive. Here we use an individual-based model to investigate the utility of genome-wide selection components analysis, which compares allele frequencies of individuals at different life history stages within a single population to detect selection without requiring a priori knowledge of traits under selection. We modeled a diploid, sexually reproducing population and introduced strong mate choice on a quantitative trait to simulate sexual selection. Genome-wide allele frequencies in adults and offspring were compared using weighted F_ST values. The average number of outlier peaks (i.e., those with significantly large F_ST values) with a quantitative trait locus in close proximity (“real” peaks) represented correct diagnoses of loci under selection, whereas peaks above the F_ST significance threshold without a quantitative trait locus reflected spurious peaks. We found that, even with moderate sample sizes, signatures of strong sexual selection were detectable, but larger sample sizes improved detection rates. The model was better able to detect selection with more neutral markers, and when quantitative trait loci and neutral markers were distributed across multiple chromosomes. Although environmental variation decreased detection rates, the identification of real peaks nevertheless remained feasible. We also found that detection rates can be improved by sampling multiple populations experiencing similar selection regimes. In short, genome-wide selection components analysis is a challenging but feasible approach for the identification of regions of the genome under selection.

A candidate multimodal functional genetic network for thermal adaptation

Vertebrate ectotherms such as reptiles provide ideal organisms for the study of adaptation to environmental thermal change. Comparative genomic and exomic studies can recover markers that diverge between warm and cold adapted lineages, but the genes that are functionally related to thermal adaptation may be difficult to identify. We here used a bioinformatics genome-mining approach to predict and identify functions for suitable candidate markers for thermal adaptation in the chicken. We first established a framework of candidate functions for such markers, and then compiled the literature on genes known to adapt to the thermal environment in different lineages of vertebrates. We then identified them in the genomes of human, chicken, and the lizard Anolis carolinensis, and established a functional genetic interaction network in the chicken. Surprisingly, markers initially identified from diverse lineages of vertebrates such as human and fish were all in close functional relationship with each other and more associated than expected by chance. This indicates that the general genetic functional network for thermoregulation and/or thermal adaptation to the environment might be regulated via similar evolutionarily conserved pathways in different vertebrate lineages. We were able to identify seven functions that were statistically overrepresented in this network, corresponding to four of our originally predicted functions plus three unpredicted functions. We describe this network as multimodal: central regulator genes with the function of relaying thermal signal (1), affect genes with different cellular functions, namely (2) lipoprotein metabolism, (3) membrane channels, (4) stress response, (5) response to oxidative stress, (6) muscle contraction and relaxation, and (7) vasodilation, vasoconstriction and regulation of blood pressure. This network constitutes a novel resource for the study of thermal adaptation in the closely related nonavian reptiles and other vertebrate ectotherms.

Divergence of Drosophila melanogaster repeatomes in response to a sharp microclimate contrast in Evolution Canyon, Israel

Repeat sequences, especially mobile elements, make up large portions of most eukaryotic genomes and provide enormous, albeit commonly underappreciated, evolutionary potential. We analyzed repeatomes of Drosophila melanogaster that have been diverging in response to a microclimate contrast in Evolution Canyon (Mount Carmel, Israel), a natural evolutionary laboratory with two abutting slopes at an average distance of only 200 m, which pose a constant ecological challenge to their local biotas. Flies inhabiting the colder and more humid north-facing slope carried about 6% more transposable elements than those from the hot and dry south-facing slope, in parallel to a suite of other genetic and phenotypic differences between the two populations. Nearly 50% of all mobile element insertions were slope unique, with many of them disrupting coding sequences of genes critical for cognition, olfaction, and thermotolerance, consistent with the observed patterns of thermotolerance differences and assortative mating.

Evolution in action: adaptation and incipient sympatric speciation with gene flow across life at “Evolution Canyon”, Israel

Various major evolutionary problems are still open, controversial or unsettled. These include even the basic evolutionary processes of adaptation and speciation. The “Evolution Canyon” model is a microscale natural laboratory that can highlight some of the basic problems requiring clarification (Nevo list of “Evolution Canyon” publications at http://evolution.haifa.ac.il). This is especially true if an interdisciplinary approach is practiced including ecological functional genomics, transcriptomics, proteomics, metabolomics and phenomics. Here I overview and reanalyze the incipient sympatric adaptive ecological speciation of five model organisms at “Evolution Canyon”, across life: the soil bacterium, Bacillus simplex; wild barley, the progenitor of cultivated barley, Hordeum spontaneum; the tiny beetle Oryzaephilus surinamensis; the cosmopolitan fruit-fly, Drosophila melanogaster, and the Africa-originated spiny mouse, Acomys cahirinus. All five models of organisms display evolution in action of microclimatic adaptation and incipient sympatric adaptive ecological speciation on the tropical and temperate abutting slopes, separated on average by only 250 meters. Some distant species converge in their micro-climatic adaptations to the hot and dry “African”, south-facing slope (SFS or AS) and to the cool and humid “European”, north-facing slope (NSF or ES). Natural selection overrules ongoing inter-slope gene-flow between the free interbreeding populations within and between slopes, and leads to adaptive incipient sympatric ecological speciation on the dramatically opposite abutting xeric savannoid and mesic forested slopes.