Genome-wide patterns of adaptation to temperate environments associated with transposable elements in Drosophila

Evolution of Two Short Interspersed Elements in Callorhinchus milii (Chondrichthyes, Holocephali) and Related Elements in Sharks and the Coelacanth

Short interspersed elements (SINEs) are non-autonomous retrotransposons. Although they usually show fast evolutionary rates, in some instances highly conserved domains (HCDs) have been observed in elements with otherwise divergent sequences and from distantly related species. Here, we document the life history of two HCD-SINE families in the elephant shark Callorhinchus milii, one specific to the holocephalan lineage (CmiSINEs) and another one (SacSINE1-CM) with homologous elements in sharks and the coelacanth (SacSINE1s, LmeSINE1s). The analyses of their relationships indicated that these elements share the same 3′-tail, which would have allowed both elements to rise to high copy number by exploiting the C. milii L2-2_CM long interspersed element (LINE) enzymes. Molecular clock analysis on SINE activity in C. milii genome evidenced two replication bursts occurring right after two major events in the holocephalan evolution: the end-Permian mass extinction and the radiation of modern Holocephali. Accordingly, the same analysis on the coelacanth homologous elements, LmeSINE1, identified a replication wave close to the split age of the two extant Latimeria species. The genomic distribution of the studied SINEs pointed out contrasting results: some elements were preferentially sorted out from gene regions, but accumulated in flanking regions, while others appear more conserved within genes. Moreover, data from the C. milii transcriptome suggest that these SINEs could be involved in miRNA biogenesis and may be targets for miRNA-based regulation.

Genome-wide characterization of non-reference transposable element insertion polymorphisms reveals genetic diversity in tropical and temperate maize

Background

Maize was originally domesticated in a tropical environment but is now widely cultivated at temperate latitudes. Temperate and tropical maize populations have diverged both genotypically and phenotypically. Tropical maize lines grown in temperate environments usually exhibit delayed flowering, pollination, and seed set, which reduces their grain yield relative to temperate adapted maize lines. One potential mechanism by which temperate maize may have adapted to a new environment is novel transposable element insertions, which can influence gene regulation. Recent advances in sequencing technology have made it possible to study variation in transposon content and insertion location in large sets of maize lines.

Results

In total, 274,408 non-redundant TEs (NRTEs) were identified using resequencing data generated from 83 maize inbred lines. The locations of DNA TEs and copia-superfamily retrotransposons showed significant positive correlations with gene density and genetic recombination rates, whereas gypsy-superfamily retrotransposons showed a negative correlation with these two parameters. Compared to tropical maize, temperate maize had fewer unique NRTEs but higher insertion frequency, lower background recombination rates, and higher linkage disequilibrium, with more NRTEs close to flowering and stress-related genes in the genome. Association mapping demonstrated that the presence/absence of 48 NRTEs was associated with flowering time and that expression of neighboring genes differed between haplotypes where a NRTE was present or absent.

Conclusions

This study suggests that NRTEs may have played an important role in creating the variation in gene regulation that enabled the rapid adaptation of maize to diverse environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4103-x) contains supplementary material, which is available to authorized users.

Useful parasites: the evolutionary biology and biotechnology applications of transposable elements

Transposable elements usually comprise the most abundant nongenic fraction of eukaryotic genomes. Because of their capacity to selfreplicate and to induce a wide range of mutations, transposable elements have long been considered as 'parasitic' or 'selfish'. Today, we recognize that the findings about genomic changes affected by transposable elements have considerably altered our view of the ways in which genomes evolve and work. Numerous studies have provided evidences that mobile elements have the potential to act as agents of evolution by increasing, rearranging and diversifying the genetic repertoire of their hosts. With large-scale sequencing becoming increasingly available, more and more scientists come across transposable element sequences in their data. I will provide examples that transposable elements, although having signatures of 'selfish' DNA, play a significant biological role in the maintainance of genome integrity and providing novel regulatoty networks. These features, along with the transpositional and mutagenic capacity to produce a raw genetic diversity, make the genome mobile fraction, a key player in species adaptation and microevolution. The last but not least, transposable elements stand as informative DNA markers that may complement other conventional DNA markers. Altogether, transposable elements represent a promising, but still largely unexplored research niche and deserve to be included into the agenda of molecular ecologists, evolutionary geneticists, conservation biologists and plant breeders.

Estimation of long terminal repeat element content in the Helicoverpa zea genome from high-throughput sequencing of bacterial artificial chromosome pools

The lepidopteran pest insect Helicoverpa zea feeds on cultivated corn and cotton across the Americas where control remains challenging owing to the evolution of resistance to chemical and transgenic insecticidal toxins, yet genomic resources remain scarce for this species. A bacterial artificial chromosome (BAC) library having a mean genomic insert size of 145 ± 20 kbp was created from a laboratory strain of H. zea, which provides ∼12.9-fold coverage of a 362.8 ± 8.8 Mbp (0.37 ± 0.09 pg) flow cytometry estimated haploid genome size. Assembly of Illumina HiSeq 2000 reads generated from 14 pools that encompassed all BAC clones resulted in 165 485 genomic contigs (N₅₀ = 3262 bp; 324.6 Mbp total). Long terminal repeat (LTR) protein coding regions annotated from 181 contigs included 30 Ty1/copia, 78 Ty3/gypsy, and 73 BEL/Pao elements, of which 60 (33.1%) encoded all five functional polyprotein (pol) domains. Approximately 14% of LTR elements are distributed non-randomly across pools of BAC clones.

Rare horizontal transmission does not hide long-term inheritance of SINE highly conserved domains in the metazoan evolution

Transposable elements (TEs) are self-replicating, mobile DNA sequences which constitute a significant fraction of eukaryotic genomes. They are generally considered selfish DNA, as their replication and random insertion may have deleterious effects on genome functionalities, although some beneficial effects and evolutionary potential have been recognized. Short interspersed elements (SINEs) are non-autonomous TEs with a modular structure: a small RNA-related head, a body, and a long interspersed element-related tail. Despite their high turnover rate and de novo emergence, the body may retain highly conserved domains (HCDs) shared among divergent SINE families: in metazoans, at least nine HCD-SINEs have been recognized. Data mining on public molecular databases allowed the retrieval of 16 new HCD-SINE families from cnidarian, molluscs, arthropods, and vertebrates. Tracking the ancestry of HCDs on the metazoan phylogeny revealed that some of them date back to the Radiata–Bilateria split. Moreover, phylogenetic and age versus divergence analyses of the most ancient HCDs suggested that long-term vertical inheritance is the rule, with few horizontal transfer events. We suggest that the evolutionary conservation of HCDs may be linked to their potential to serve as recombination hotspots. This indirectly affects host genomes by maintaining active and diverse SINE lineages, whose insertions may impact (either positively or negatively) on the evolution of the genome.

The industrial melanism mutation in British peppered moths is a transposable element

Discovering the mutational events that fuel adaptation to environmental change remains an important challenge for evolutionary biology. The classroom example of a visible evolutionary response is industrial melanism in the peppered moth (Biston betularia): the replacement, during the Industrial Revolution, of the common pale typica form by a previously unknown black (carbonaria) form, driven by the interaction between bird predation and coal pollution. The carbonaria locus has been coarsely localized to a 200-kilobase region, but the specific identity and nature of the sequence difference controlling the carbonaria-typica polymorphism, and the gene it influences, are unknown. Here we show that the mutation event giving rise to industrial melanism in Britain was the insertion of a large, tandemly repeated, transposable element into the first intron of the gene cortex. Statistical inference based on the distribution of recombined carbonaria haplotypes indicates that this transposition event occurred around 1819, consistent with the historical record. We have begun to dissect the mode of action of the carbonaria transposable element by showing that it increases the abundance of a cortex transcript, the protein product of which plays an important role in cell-cycle regulation, during early wing disc development. Our findings fill a substantial knowledge gap in the iconic example of microevolutionary change, adding a further layer of insight into the mechanism of adaptation in response to natural selection. The discovery that the mutation itself is a transposable element will stimulate further debate about the importance of 'jumping genes' as a source of major phenotypic novelty.

Zinc-Induced Transposition of Insertion Sequence Elements Contributes to Increased Adaptability of Cupriavidus metallidurans

Bacteria can respond to adverse environments by increasing their genomic variability and subsequently facilitating adaptive evolution. To demonstrate this, the contribution of Insertion Sequence (IS) elements to the genetic adaptation of Cupriavidus metallidurans AE126 to toxic zinc concentrations was determined. This derivative of type strain CH34, devoid of its main zinc resistance determinant, is still able to increase its zinc resistance level. Specifically, upon plating on medium supplemented with a toxic zinc concentration, resistant variants arose in which a compromised cnrYX regulatory locus caused derepression of CnrH sigma factor activity and concomitant induction of the corresponding RND-driven cnrCBA efflux system. Late-occurring zinc resistant variants likely arose in response to the selective conditions, as they were enriched in cnrYX disruptions caused by specific IS elements whose transposase expression was found to be zinc-responsive. Interestingly, deletion of cnrH, and consequently the CnrH-dependent adaptation potential, still enabled adaptation by transposition of IS elements (ISRme5 and IS1086) that provided outward-directed promoters driving cnrCBAT transcription. Finally, adaptation to zinc by IS reshuffling can also enhance the adaptation to subsequent environmental challenges. Thus, transposition of IS elements can be induced by stress conditions and play a multifaceted, pivotal role in the adaptation to these and subsequent stress conditions.

Genomic Analyses Reveal Demographic History and Temperate Adaptation of the Newly Discovered Honey Bee Subspecies Apis mellifera sinisxinyuan n. ssp

Studying the genetic signatures of climate-driven selection can produce insights into local adaptation and the potential impacts of climate change on populations. The honey bee (Apis mellifera) is an interesting species to study local adaptation because it originated in tropical/subtropical climatic regions and subsequently spread into temperate regions. However, little is known about the genetic basis of its adaptation to temperate climates. Here, we resequenced the whole genomes of ten individual bees from a newly discovered population in temperate China and downloaded resequenced data from 35 individuals from other populations. We found that the new population is an undescribed subspecies in the M-lineage of A. mellifera (Apis mellifera sinisxinyuan). Analyses of population history show that long-term global temperature has strongly influenced the demographic history of A. m. sinisxinyuan and its divergence from other subspecies. Further analyses comparing temperate and tropical populations identified several candidate genes related to fat body and the Hippo signaling pathway that are potentially involved in adaptation to temperate climates. Our results provide insights into the demographic history of the newly discovered A. m. sinisxinyuan, as well as the genetic basis of adaptation of A. mellifera to temperate climates at the genomic level. These findings will facilitate the selective breeding of A. mellifera to improve the survival of overwintering colonies.

Secondary contact and local adaptation contribute to genome-wide patterns of clinal variation in Drosophila melanogaster

Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome-wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome-wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an FST outlier approach. Our results provide novel insight into the well-studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.

Unique transposon landscapes are pervasive across Drosophila melanogaster genomes

To understand how transposon landscapes (TLs) vary across animal genomes, we describe a new method called the Transposon Insertion and Depletion AnaLyzer (TIDAL) and a database of >300 TLs in Drosophila melanogaster (TIDAL-Fly). Our analysis reveals pervasive TL diversity across cell lines and fly strains, even for identically named sub-strains from different laboratories such as the ISO1 strain used for the reference genome sequence. On average, >500 novel insertions exist in every lab strain, inbred strains of the Drosophila Genetic Reference Panel (DGRP), and fly isolates in the Drosophila Genome Nexus (DGN). A minority (<25%) of transposon families comprise the majority (>70%) of TL diversity across fly strains. A sharp contrast between insertion and depletion patterns indicates that many transposons are unique to the ISO1 reference genome sequence. Although TL diversity from fly strains reaches asymptotic limits with increasing sequencing depth, rampant TL diversity causes unsaturated detection of TLs in pools of flies. Finally, we show novel transposon insertions negatively correlate with Piwi-interacting RNA (piRNA) levels for most transposon families, except for the highly-abundant roo retrotransposon. Our study provides a useful resource for Drosophila geneticists to understand how transposons create extensive genomic diversity in fly cell lines and strains.

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.

The evolutionary dynamics of tRNA-gene copy number and codon-use in E. coli

Background

The introduction of foreign DNA by Lateral Gene Transfer (LGT) can quickly and drastically alter genome composition. Problems can arise if the genes introduced by LGT use codons that are not suited to the host’s translational machinery. Here we investigate compensatory adaptation of E. coli in response to the introduction of large volumes of codons that are rarely used by the host genome.

Results

We analyze genome sequences from the E. coli/Shigella complex, and find that certain tRNA genes are present in multiple copies in two pathogenic Shigella and O157:H7 subgroups of E. coli. Furthermore, we show that the codons that correspond to these multi-copy number tRNA genes are enriched in the high copy number Selfish Genetic Elements (SGE’s) in Shigella and laterally introduced genes in O157:H7. We analyze the duplicate copies and find evidence for the selective retention of tRNA genes introduced by LGT in response to the changed codon content of the genome.

Conclusion

These data support a model where the relatively rapid influx of LGT genes and SGE’s introduces a large number of genes maladapted to the host’s translational machinery. Under these conditions, it becomes advantageous for the host to retain tRNA genes that are required for the incorporation of amino acids at these codons. Subsequently, the increased number of copies of these specific tRNA genes adjusts the cellular tRNA pool to the demands set by global shifts in codon usage.

Electronic supplementary material

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

Revisiting classic clines in Drosophila melanogaster in the age of genomics

Adaptation to spatially varying environments has been studied for decades, but advances in sequencing technology are now enabling researchers to investigate the landscape of genetic variation underlying this adaptation genome wide. In this review we highlight some of the decades-long research on local adaptation in Drosophila melanogaster from well-studied clines in North America and Australia. We explore the evidence for parallel adaptation and identify commonalities in the genes responding to clinal selection across continents as well as discussing instances where patterns differ among clines. We also investigate recent studies utilizing whole-genome data to identify clines in D. melanogaster and several other systems. Although connecting segregating genomic variation to variation in phenotypes and fitness remains challenging, clinal genomics is poised to increase our understanding of local adaptation and the selective pressures that drive the extensive phenotypic diversity observed in nature.

Asian Citrus Psyllid Expression Profiles Suggest Candidatus Liberibacter Asiaticus-Mediated Alteration of Adult Nutrition and Metabolism, and of Nymphal Development and Immunity

The Asian citrus psyllid (ACP) Diaphorina citri Kuwayama (Hemiptera: Psyllidae) is the insect vector of the fastidious bacterium Candidatus Liberibacter asiaticus (CLas), the causal agent of citrus greening disease, or Huanglongbing (HLB). The widespread invasiveness of the psyllid vector and HLB in citrus trees worldwide has underscored the need for non-traditional approaches to manage the disease. One tenable solution is through the deployment of RNA interference technology to silence protein-protein interactions essential for ACP-mediated CLas invasion and transmission. To identify psyllid interactor-bacterial effector combinations associated with psyllid-CLas interactions, cDNA libraries were constructed from CLas-infected and CLas-free ACP adults and nymphs, and analyzed for differential expression. Library assemblies comprised 24,039,255 reads and yielded 45,976 consensus contigs. They were annotated (UniProt), classified using Gene Ontology, and subjected to in silico expression analyses using the Transcriptome Computational Workbench (TCW) (http://www.sohomoptera.org/ACPPoP/). Functional-biological pathway interpretations were carried out using the Kyoto Encyclopedia of Genes and Genomes databases. Differentially expressed contigs in adults and/or nymphs represented genes and/or metabolic/pathogenesis pathways involved in adhesion, biofilm formation, development-related, immunity, nutrition, stress, and virulence. Notably, contigs involved in gene silencing and transposon-related responses were documented in a psyllid for the first time. This is the first comparative transcriptomic analysis of ACP adults and nymphs infected and uninfected with CLas. The results provide key initial insights into host-parasite interactions involving CLas effectors that contribute to invasion-virulence, and to host nutritional exploitation and immune-related responses that appear to be essential for successful ACP-mediated circulative, propagative CLas transmission.

Exploring the Phenotypic Space and the Evolutionary History of a Natural Mutation in Drosophila melanogaster

A major challenge of modern Biology is elucidating the functional consequences of natural mutations. Although we have a good understanding of the effects of laboratory-induced mutations on the molecular- and organismal-level phenotypes, the study of natural mutations has lagged behind. In this work, we explore the phenotypic space and the evolutionary history of a previously identified adaptive transposable element insertion. We first combined several tests that capture different signatures of selection to show that there is evidence of positive selection in the regions flanking FBti0019386 insertion. We then explored several phenotypes related to known phenotypic effects of nearby genes, and having plausible connections to fitness variation in nature. We found that flies with FBti0019386 insertion had a shorter developmental time and were more sensitive to stress, which are likely to be the adaptive effect and the cost of selection of this mutation, respectively. Interestingly, these phenotypic effects are not consistent with a role of FBti0019386 in temperate adaptation as has been previously suggested. Indeed, a global analysis of the population frequency of FBti0019386 showed that climatic variables explain well the FBti0019386 frequency patterns only in Australia. Finally, although FBti0019386 insertion could be inducing the formation of heterochromatin by recruiting HP1a (Heterochromatin Protein 1a) protein, the insertion is associated with upregulation of sra in adult females. Overall, our integrative approach allowed us to shed light on the evolutionary history, the relevant fitness effects, and the likely molecular mechanisms of an adaptive mutation and highlights the complexity of natural genetic variants.

Transposable elements as agents of rapid adaptation may explain the genetic paradox of invasive species

Rapid adaptation of invasive species to novel habitats has puzzled evolutionary biologists for decades, especially as this often occurs in the face of limited genetic variability. Although some ecological traits common to invasive species have been identified, little is known about the possible genomic/genetic mechanisms that may underlie their success. A common scenario in many introductions is that small founder population sizes will often lead to reduced genetic diversity, but that invading populations experience large environmental perturbations, such as changes in habitat and environmental stress. Although sudden and intense stress is usually considered in a negative context, these perturbations may actually facilitate rapid adaptation by affecting genome structure, organization and function via interactions with transposable elements (TEs), especially in populations with low genetic diversity. Stress-induced changes in TE activity can alter gene action and can promote structural variation that may facilitate the rapid adaptation observed in new environments. We focus here on the adaptive potential of TEs in relation to invasive species and highlight their role as powerful mutational forces that can rapidly create genetic diversity. We hypothesize that activity of transposable elements can explain rapid adaptation despite low genetic variation (the genetic paradox of invasive species), and provide a framework under which this hypothesis can be tested using recently developed and emerging genomic technologies.

Transposable element islands facilitate adaptation to novel environments in an invasive species

Adaptation requires genetic variation, but founder populations are generally genetically depleted. Here we sequence two populations of an inbred ant that diverge in phenotype to determine how variability is generated. Cardiocondyla obscurior has the smallest of the sequenced ant genomes and its structure suggests a fundamental role of transposable elements (TEs) in adaptive evolution. Accumulations of TEs (TE islands) comprising 7.18% of the genome evolve faster than other regions with regard to single-nucleotide variants, gene/exon duplications and deletions and gene homology. A non-random distribution of gene families, larvae/adult specific gene expression and signs of differential methylation in TE islands indicate intragenomic differences in regulation, evolutionary rates and coalescent effective population size. Our study reveals a tripartite interplay between TEs, life history and adaptation in an invasive species.

Genetic variation is key to species evolution. Here the authors sequence two phenotypically distinct populations of the ant Cardiocondyla obscurior, and find accumulations of transposable elements correlating with genetic variation that may have a role in differentiation, adaptation and speciation.

T-lex2: genotyping, frequency estimation and re-annotation of transposable elements using single or pooled next-generation sequencing data

Transposable elements (TEs) constitute the most active, diverse and ancient component in a broad range of genomes. Complete understanding of genome function and evolution cannot be achieved without a thorough understanding of TE impact and biology. However, in-depth analysis of TEs still represents a challenge due to the repetitive nature of these genomic entities. In this work, we present a broadly applicable and flexible tool: T-lex2. T-lex2 is the only available software that allows routine, automatic and accurate genotyping of individual TE insertions and estimation of their population frequencies both using individual strain and pooled next-generation sequencing data. Furthermore, T-lex2 also assesses the quality of the calls allowing the identification of miss-annotated TEs and providing the necessary information to re-annotate them. The flexible and customizable design of T-lex2 allows running it in any genome and for any type of TE insertion. Here, we tested the fidelity of T-lex2 using the fly and human genomes. Overall, T-lex2 represents a significant improvement in our ability to analyze the contribution of TEs to genome function and evolution as well as learning about the biology of TEs. T-lex2 is freely available online at http://sourceforge.net/projects/tlex.

Adaptations to climate-mediated selective pressures in sheep

Following domestication, sheep (Ovis aries) have become essential farmed animals across the world through adaptation to a diverse range of environments and varied production systems. Climate-mediated selective pressure has shaped phenotypic variation and has left genetic "footprints" in the genome of breeds raised in different agroecological zones. Unlike numerous studies that have searched for evidence of selection using only population genetics data, here, we conducted an integrated coanalysis of environmental data with single nucleotide polymorphism (SNP) variation. By examining 49,034 SNPs from 32 old, autochthonous sheep breeds that are adapted to a spectrum of different regional climates, we identified 230 SNPs with evidence for selection that is likely due to climate-mediated pressure. Among them, 189 (82%) showed significant correlation (P ≤ 0.05) between allele frequency and climatic variables in a larger set of native populations from a worldwide range of geographic areas and climates. Gene ontology analysis of genes colocated with significant SNPs identified 17 candidates related to GTPase regulator and peptide receptor activities in the biological processes of energy metabolism and endocrine and autoimmune regulation. We also observed high linkage disequilibrium and significant extended haplotype homozygosity for the core haplotype TBC1D12-CH1 of TBC1D12. The global frequency distribution of the core haplotype and allele OAR22_18929579-A showed an apparent geographic pattern and significant (P ≤ 0.05) correlations with climatic variation. Our results imply that adaptations to local climates have shaped the spatial distribution of some variants that are candidates to underpin adaptive variation in sheep.

Population genomics of transposable elements in Drosophila

Studies of the population dynamics of transposable elements (TEs) in Drosophila melanogaster indicate that consistent forces are affecting TEs independently of their modes of transposition and regulation. New sequencing technologies enable biologists to sample genomes at an unprecedented scale in order to quantify genome-wide polymorphism for annotated and novel TE insertions. In this review, we first present new insights gleaned from high-throughput data for population genomics studies of D. melanogaster. We then consider the latest population genomics models for TE evolution and present examples of functional evidence revealed by genome-wide studies of TE population dynamics in D. melanogaster. Although most of the TE insertions are deleterious or neutral, some TE insertions increase the fitness of the individual that carries them and play a role in genome adaptation.

A transposable element insertion confers xenobiotic resistance in Drosophila

The increase in availability of whole genome sequences makes it possible to search for evidence of adaptation at an unprecedented scale. Despite recent progress, our understanding of the adaptive process is still very limited due to the difficulties in linking adaptive mutations to their phenotypic effects. In this study, we integrated different levels of biological information to pinpoint the ecologically relevant fitness effects and the underlying molecular and biochemical mechanisms of a putatively adaptive TE insertion in Drosophila melanogaster: the pogo transposon FBti0019627. We showed that other than being incorporated into Kmn1 transcript, FBti0019627 insertion also affects the polyadenylation signal choice of CG11699 gene. Consequently, only the short 3'UTR transcript of CG11699 gene is produced and the expression level of this gene is higher in flies with the insertion. Our results indicated that increased CG11699 expression leads to xenobiotic stress resistance through increased ALDH-III activity: flies with FBti0019627 insertion showed increased survival rate in response to benzaldehyde, a natural xenobiotic, and to carbofuran, a synthetic insecticide. Although differences in survival rate between flies with and without the insertion were not always significant, when they were, they were consistent with FBti0019627 mediating resistance to xenobiotics. Taken together, our results provide a plausible explanation for the increase in frequency of FBti0019627 in natural populations of D. melanogaster and add to the limited number of examples in which a natural genetic mutation has been linked to its ecologically relevant phenotype. Furthermore, the widespread distribution of TEs across the tree of life and conservation of stress response pathways across organisms make our results relevant not only for Drosophila, but for other organisms as well.

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.

Transposable elements: powerful contributors to angiosperm evolution and diversity

Transposable elements (TEs) are a dominant feature of most flowering plant genomes. Together with other accepted facilitators of evolution, accumulating data indicate that TEs can explain much about their rapid evolution and diversification. Genome size in angiosperms is highly correlated with TE content and the overwhelming bulk (>80%) of large genomes can be composed of TEs. Among retro-TEs, long terminal repeats (LTRs) are abundant, whereas DNA-TEs, which are often less abundant than retro-TEs, are more active. Much adaptive or evolutionary potential in angiosperms is due to the activity of TEs (active TE-Thrust), resulting in an extraordinary array of genetic changes, including gene modifications, duplications, altered expression patterns, and exaptation to create novel genes, with occasional gene disruption. TEs implicated in the earliest origins of the angiosperms include the exapted Mustang, Sleeper, and Fhy3/Far1 gene families. Passive TE-Thrust can create a high degree of adaptive or evolutionary potential by engendering ectopic recombination events resulting in deletions, duplications, and karyotypic changes. TE activity can also alter epigenetic patterning, including that governing endosperm development, thus promoting reproductive isolation. Continuing evolution of long-lived resprouter angiosperms, together with genetic variation in their multiple meristems, indicates that TEs can facilitate somatic evolution in addition to germ line evolution. Critical to their success, angiosperms have a high frequency of polyploidy and hybridization, with resultant increased TE activity and introgression, and beneficial gene duplication. Together with traditional explanations, the enhanced genomic plasticity facilitated by TE-Thrust, suggests a more complete and satisfactory explanation for Darwin's "abominable mystery": the spectacular success of the angiosperms.

The transposable element Bari-Jheh mediates oxidative stress response in Drosophila

Elucidating the fitness effects of natural genetic variants is one of the current major challenges in evolutionary biology. Understanding the interplay between genotype, phenotype and environment is necessary to make accurate predictions of important biological outcomes such as stress resistance or yield in economically important plants and animals, and disease in humans. Based on population frequency patterns and footprints of selection at the DNA level, the transposable element Bari-Jheh, inserted in the intergenic region of Juvenile Hormone Epoxy Hydrolase (Jheh) genes, was previously identified as putatively adaptive. However, the adaptive effect of this mutation remained elusive. In this work, we integrate information on transcription factor binding sites, available ChIP-Seq data, gene expression analyses and phenotypic assays to identify the functional and the mechanistic underpinnings of Bari-Jheh. We show that Bari-Jheh adds extra antioxidant response elements upstream of Jheh1 and Jheh2 genes. Accordingly, we find that Bari-Jheh is associated with upregulation of Jheh1 and Jheh2 and with resistance to oxidative stress induced by two different compounds relevant for natural D. melanogaster populations. We further show that TEs other than Bari-Jheh might be playing a role in the D. melanogaster response to oxidative stress. Overall our results contribute to the understanding of resistance to oxidative stress in natural populations and highlight the role of transposable elements in environmental adaptation. The replicability of fitness effects on different genetic backgrounds also suggests that epistatic interactions do not seem to dominate the genetic architecture of oxidative stress resistance.

The bright side of transposons in crop evolution

The past decades have revealed an unexpected yet prominent role of so-called 'junk DNA' in the regulation of gene expression, thereby challenging our view of the mechanisms underlying phenotypic evolution. In particular, several mechanisms through which transposable elements (TEs) participate in functional genome diversity have been depicted, bringing to light the 'TEs bright side'. However, the relative contribution of those mechanisms and, more generally, the importance of TE-based polymorphisms on past and present phenotypic variation in crops species remain poorly understood. Here, we review current knowledge on both issues, and discuss how analyses of massively parallel sequencing data combined with statistical methodologies and functional validations will help unravelling the impact of TEs on crop evolution in a near future.

P element activity and molecular structure in Drosophila melanogaster populations from Firtina Valley, Turkey

In order to study P element dynamics in natural populations of Drosophila melanogaster, 88 isofemale lines were examined from the Firtina Valley, Turkey. The P-M gonadal dysgenesis characteristics and the molecular patterns of P and KP elements were analyzed. Gonadal dysgenesis tests showed a slight variation both for P activity and P susceptibility, however the results showed a predominant M' phenotype for this region. The P and KP element were also characterized by polymerase chain reaction. The molecular analyses showed that all the populations examined had the entire 1.15 kb KP element. The molecular patterns of KP elements were the same for the populations studied. No clear relationship was found between phenotype and genomic P element composition. The correlations between the level of gonadal dysgenesis percentage (as an index for P activity and P susceptibility) and several geoclimatic factors were tested, and no general effects of altitude, temperature, rainfall, or humidity were found. The theoretical P' strain, which is very rare in natural populations, was also recorded for this region.

A population genetic model for the maintenance of R2 retrotransposons in rRNA gene loci

R2 retrotransposable elements exclusively insert into the tandemly repeated rRNA genes, the rDNA loci, of their animal hosts. R2 elements form stable long-term associations with their host, in which all individuals in a population contain many potentially active copies, but only a fraction of these individuals show active R2 retrotransposition. Previous studies have found that R2 RNA transcripts are processed from a 28S co-transcript and that the likelihood of R2-inserted units being transcribed is dependent upon their distribution within the rDNA locus. Here we analyze the rDNA locus and R2 elements from nearly 100 R2-active and R2-inactive individuals from natural populations of Drosophila simulans. Along with previous findings concerning the structure and expression of the rDNA loci, these data were incorporated into computer simulations to model the crossover events that give rise to the concerted evolution of the rRNA genes. The simulations that best reproduce the population data assume that only about 40 rDNA units out of the over 200 total units are actively transcribed and that these transcribed units are clustered in a single region of the locus. In the model, the host establishes this transcription domain at each generation in the region with the fewest R2 insertions. Only if the host cannot avoid R2 insertions within this 40-unit domain are R2 elements active in that generation. The simulations also require that most crossover events in the locus occur in the transcription domain in order to explain the empirical observation that R2 elements are seldom duplicated by crossover events. Thus the key to the long-term stability of R2 elements is the stochastic nature of the crossover events within the rDNA locus, and the inevitable expansions and contractions that introduce and remove R2-inserted units from the transcriptionally active domain.

Selfish transposable elements survive in eukaryotic genomes despite the elaborate mechanisms developed by the hosts to limit their activity. One accessible system that simplifies the complex interactions between element and host involves the R2 elements, which exclusively insert in the tandemly arranged rRNA genes. R2 exhibits remarkable stability in animal lineages even though each insertion inactivates one rRNA gene. Here we determine the size of the rDNA locus and R2 number in natural isolates of Drosophila simulans. Combined with previous data concerning the expression and regulation of R2, we develop a detailed population genetic model for rRNA gene and R2 evolution that duplicates all properties of the rRNA loci in natural populations. Critical components of the model are that only a contiguous 40 unit array of rRNA gene units are needed for transcription, that R2 elements are active only when present in this transcription domain, and that most of the crossovers in the rDNA loci occur in this domain. These results suggest that the key to the long-term survival of R2 is the redistribution of rDNA units in the locus brought about by the crossovers that maintain sequence identity in all rDNA units.

Genomic dynamics of transposable elements in the western clawed frog (Silurana tropicalis)

Transposable elements (TEs) are repetitive DNA sequences that can make new copies of themselves that are inserted elsewhere in a host genome. The abundance and distributions of TEs vary considerably among phylogenetically diverse hosts. With the aim of exploring the basis of this variation, we evaluated correlations between several genomic variables and the presence of TEs and non-TE repeats in the complete genome sequence of the Western clawed frog (Silurana tropicalis). This analysis reveals patterns of TE insertion consistent with gene disruption but not with the insertional preference model. Analysis of non-TE repeats recovered unique features of their genome-wide distribution when compared with TE repeats, including no strong correlation with exons and a particularly strong negative correlation with GC content. We also collected polymorphism data from 25 TE insertion sites in 19 wild-caught S. tropicalis individuals. DNA transposon insertions were fixed at eight of nine sites and at a high frequency at one of nine, whereas insertions of long terminal repeat (LTR) and non-LTR retrotransposons were fixed at only 4 of 16 sites and at low frequency at 12 of 16. A maximum likelihood model failed to attribute these differences in insertion frequencies to variation in selection pressure on different classes of TE, opening the possibility that other phenomena such as variation in rates of replication or duration of residence in the genome could play a role. Taken together, these results identify factors that sculpt heterogeneity in TE distribution in S. tropicalis and illustrate that genomic dynamics differ markedly among TE classes and between TE and non-TE repeats.

Population and sex differences in Drosophila melanogaster brain gene expression

Background

Changes in gene regulation are thought to be crucial for the adaptation of organisms to their environment. Transcriptome analyses can be used to identify candidate genes for ecological adaptation, but can be complicated by variation in gene expression between tissues, sexes, or individuals. Here we use high-throughput RNA sequencing of a single Drosophila melanogaster tissue to detect brain-specific differences in gene expression between the sexes and between two populations, one from the ancestral species range in sub-Saharan Africa and one from the recently colonized species range in Europe.

Results

Relatively few genes (<100) displayed sexually dimorphic expression in the brain, but there was an enrichment of sex-biased genes, especially male-biased genes, on the X chromosome. Over 340 genes differed in brain expression between flies from the African and European populations, with the inter-population divergence being highly correlated between males and females. The differentially expressed genes included those involved in stress response, olfaction, and detoxification. Expression differences were associated with transposable element insertions at two genes implicated in insecticide resistance (Cyp6g1 and CHKov1).

Conclusions

Analysis of the brain transcriptome revealed many genes differing in expression between populations that were not detected in previous studies using whole flies. There was little evidence for sex-specific regulatory adaptation in the brain, as most expression differences between populations were observed in both males and females. The enrichment of genes with sexually dimorphic expression on the X chromosome is consistent with dosage compensation mechanisms affecting sex-biased expression in somatic tissues.

Transposable elements and viruses as factors in adaptation and evolution: an expansion and strengthening of the TE-Thrust hypothesis

In addition to the strong divergent evolution and significant and episodic evolutionary transitions and speciation we previously attributed to TE-Thrust, we have expanded the hypothesis to more fully account for the contribution of viruses to TE-Thrust and evolution. The concept of symbiosis and holobiontic genomes is acknowledged, with particular emphasis placed on the creativity potential of the union of retroviral genomes with vertebrate genomes. Further expansions of the TE-Thrust hypothesis are proposed regarding a fuller account of horizontal transfer of TEs, the life cycle of TEs, and also, in the case of a mammalian innovation, the contributions of retroviruses to the functions of the placenta. The possibility of drift by TE families within isolated demes or disjunct populations, is acknowledged, and in addition, we suggest the possibility of horizontal transposon transfer into such subpopulations. “Adaptive potential” and “evolutionary potential” are proposed as the extremes of a continuum of “intra-genomic potential” due to TE-Thrust. Specific data is given, indicating “adaptive potential” being realized with regard to insecticide resistance, and other insect adaptations. In this regard, there is agreement between TE-Thrust and the concept of adaptation by a change in allele frequencies. Evidence on the realization of “evolutionary potential” is also presented, which is compatible with the known differential survivals, and radiations of lineages. Collectively, these data further suggest the possibility, or likelihood, of punctuated episodes of speciation events and evolutionary transitions, coinciding with, and heavily underpinned by, intermittent bursts of TE activity.

Keeping pace with climate change: what is wrong with the evolutionary potential of upper thermal limits?

The potential of populations to evolve in response to ongoing climate change is partly conditioned by the presence of heritable genetic variation in relevant physiological traits. Recent research suggests that Drosophila melanogaster exhibits negligible heritability, hence little evolutionary potential in heat tolerance when measured under slow heating rates that presumably mimic conditions in nature. Here, we study the effects of directional selection for increased heat tolerance using Drosophila as a model system. We combine a physiological model to simulate thermal tolerance assays with multilocus models for quantitative traits. Our simulations show that, whereas the evolutionary response of the genetically determined upper thermal limit (CTmax) is independent of methodological context, the response in knockdown temperatures varies with measurement protocol and is substantially (up to 50%) lower than for CTmax. Realized heritabilities of knockdown temperature may grossly underestimate the true heritability of CTmax. For instance, assuming that the true heritability of CTmax in the base population is h(2) = 0.25, realized heritabilities of knockdown temperature are around 0.08-0.16 depending on heating rate. These effects are higher in slow heating assays, suggesting that flawed methodology might explain the apparently limited evolutionary potential of cosmopolitan D. melanogaster.

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.

Longitudinal trends in climate drive flowering time clines in North American Arabidopsis thaliana

Introduced species frequently show geographic differentiation, and when differentiation mirrors the ancestral range, it is often taken as evidence of adaptive evolution. The mouse-ear cress (Arabidopsis thaliana) was introduced to North America from Eurasia 150–200 years ago, providing an opportunity to study parallel adaptation in a genetic model organism. Here, we test for clinal variation in flowering time using 199 North American (NA) accessions of A. thaliana, and evaluate the contributions of major flowering time genes FRI, FLC, and PHYC as well as potential ecological mechanisms underlying differentiation. We find evidence for substantial within population genetic variation in quantitative traits and flowering time, and putatively adaptive longitudinal differentiation, despite low levels of variation at FRI, FLC, and PHYC and genome-wide reductions in population structure relative to Eurasian (EA) samples. The observed longitudinal cline in flowering time in North America is parallel to an EA cline, robust to the effects of population structure, and associated with geographic variation in winter precipitation and temperature. We detected major effects of FRI on quantitative traits associated with reproductive fitness, although the haplotype associated with higher fitness remains rare in North America. Collectively, our results suggest the evolution of parallel flowering time clines through novel genetic mechanisms.

Distribution of genes and repetitive elements in the Diabrotica virgifera virgifera genome estimated using BAC sequencing

Feeding damage caused by the western corn rootworm, Diabrotica virgifera virgifera, is destructive to corn plants in North America and Europe where control remains challenging due to evolution of resistance to chemical and transgenic toxins. A BAC library, DvvBAC1, containing 109,486 clones with 104 ± 34.5 kb inserts was created, which has an ~4.56X genome coverage based upon a 2.58 Gb (2.80 pg) flow cytometry-estimated haploid genome size. Paired end sequencing of 1037 BAC inserts produced 1.17 Mb of data (~0.05% genome coverage) and indicated ~9.4 and 16.0% of reads encode, respectively, endogenous genes and transposable elements (TEs). Sequencing genes within BAC full inserts demonstrated that TE densities are high within intergenic and intron regions and contribute to the increased gene size. Comparison of homologous genome regions cloned within different BAC clones indicated that TE movement may cause haplotype variation within the inbred strain. The data presented here indicate that the D. virgifera virgifera genome is large in size and contains a high proportion of repetitive sequence. These BAC sequencing methods that are applicable for characterization of genomes prior to sequencing may likely be valuable resources for genome annotation as well as scaffolding.

Evolution of genome content: population dynamics of transposable elements in flies and humans

Recent research is starting to shed light on the factors that influence the population and evolutionary dynamics of transposable elements (TEs) and TE life cycles. Genomes differ sharply in the number of TE copies, in the level of TE activity, in the diversity of TE families and types, and in the proportion of old and young TEs. In this chapter, we focus on two well-studied genomes with strikingly different architectures, humans and Drosophila, which represent two extremes in terms of TE diversity and population dynamics. We argue that some of the answers might lie in (1) the larger population size and consequently more effective selection against new TE insertions due to ectopic recombination in flies compared to humans; and (2) in the faster rate of DNA loss in flies compared to humans leading to much faster removal of fixed TE copies from the fly genome.

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.

Mobilizing the genome of Lepidoptera through novel sequence gains and end creation by non-autonomous Lep1 Helitrons

Transposable elements (TEs) can affect the structure of genomes through their acquisition and transposition of novel DNA sequences. The 134-bp repetitive elements, Lep1, are conserved non-autonomous Helitrons in lepidopteran genomes that have characteristic 5′-CT and 3′-CTAY nucleotide termini, a 3′-terminal hairpin structure, a 5′- and 3′-subterminal inverted repeat (SIR), and integrations that occur between AT or TT nucleotides. Lep1 Helitrons have acquired and propagated sequences downstream of their 3′-CTAY termini that are 57–344-bp in length and have termini composed of a 3′-CTRR preceded by a 3′-hairpin structure and a region complementary to the 5′-SIR (3′-SIRb). Features of both the Lep1 Helitron and multiple acquired sequences indicate that secondary structures at the 3′-terminus may have a role in rolling circle replication or genome integration mechanisms, and are a prerequisite for novel end creation by Helitron-like TEs. The preferential integration of Lep1 Helitrons in proximity to gene-coding regions results in the creation of genetic novelty that is shown to impact gene structure and function through the introduction of novel exon sequence (exon shuffling). These findings are important in understanding the structural requirements of genomic DNA sequences that are acquired and transposed by Helitron-like TEs.

Identification of a functional transposon insertion in the maize domestication gene tb1

Genetic diversity created by transposable elements is an important source of functional variation upon which selection acts during evolution. Transposable elements are associated with adaptation to temperate climates in Drosophila, a SINE element is associated with the domestication of small dog breeds from the gray wolf and there is evidence that transposable elements were targets of selection during human evolution. Although the list of examples of transposable elements associated with host gene function continues to grow, proof that transposable elements are causative and not just correlated with functional variation is limited. Here we show that a transposable element (Hopscotch) inserted in a regulatory region of the maize domestication gene, teosinte branched1 (tb1), acts as an enhancer of gene expression and partially explains the increased apical dominance in maize compared to its progenitor, teosinte. Molecular dating indicates that the Hopscotch insertion predates maize domestication by at least 10,000 years, indicating that selection acted on standing variation rather than new mutation.

Parallel adaptive divergence among geographically diverse human populations

Few genetic differences between human populations conform to the classic model of positive selection, in which a newly arisen mutation rapidly approaches fixation in one lineage, suggesting that adaptation more commonly occurs via moderate changes in standing variation at many loci. Detecting and characterizing this type of complex selection requires integrating individually ambiguous signatures across genomically and geographically extensive data. Here, we develop a novel approach to test the hypothesis that selection has favored modest divergence at particular loci multiple times in independent human populations. We find an excess of SNPs showing non-neutral parallel divergence, enriched for genic and nonsynonymous polymorphisms in genes encompassing diverse and often disease related functions. Repeated parallel evolution in the same direction suggests common selective pressures in disparate habitats. We test our method with extensive coalescent simulations and show that it is robust to a wide range of demographic events. Our results demonstrate phylogenetically orthogonal patterns of local adaptation caused by subtle shifts at many widespread polymorphisms that likely underlie substantial phenotypic diversity.

Identifying regions of the human genome that differ among populations because of natural selection is both essential for understanding evolutionary history and a powerful method for finding functionally important variants that contribute to phenotypic diversity and disease. Adaptive events on timescales corresponding to the human diaspora may often manifest as relatively small changes in allele frequencies at numerous loci that are difficult to distinguish from stochastic changes due to genetic drift, rather than the more dramatic selective sweeps described by classic models of natural selection. In order to test whether a substantial proportion of interpopulation genetic differences are indeed adaptive, we identify loci that have undergone moderate allele frequency changes in multiple independent human lineages, and we test whether these parallel divergence events are more frequent than expected by chance. We report a significant excess of polymorphisms showing parallel divergence, especially within genes, a pattern that is best explained by geographically varying natural selection. Our results indicate that local adaptation in humans has occurred by subtle, repeated changes at particular genes that are likely to be associated with important morphological and physiological differences among human populations.

A brief history of the status of transposable elements: from junk DNA to major players in evolution

The idea that some genetic factors are able to move around chromosomes emerged more than 60 years ago when Barbara McClintock first suggested that such elements existed and had a major role in controlling gene expression and that they also have had a major influence in reshaping genomes in evolution. It was many years, however, before the accumulation of data and theories showed that this latter revolutionary idea was correct although, understandably, it fell far short of our present view of the significant influence of what are now known as "transposable elements" in evolution. In this article, I summarize the main events that influenced my thinking about transposable elements as a young scientist and the influence and role of these specific genomic elements in evolution over subsequent years. Today, we recognize that the findings about genomic changes affected by transposable elements have considerably altered our view of the ways in which genomes evolve and work.

Adaptations to climate-mediated selective pressures in humans

Humans inhabit a remarkably diverse range of environments, and adaptation through natural selection has likely played a central role in the capacity to survive and thrive in extreme climates. Unlike numerous studies that used only population genetic data to search for evidence of selection, here we scan the human genome for selection signals by identifying the SNPs with the strongest correlations between allele frequencies and climate across 61 worldwide populations. We find a striking enrichment of genic and nonsynonymous SNPs relative to non-genic SNPs among those that are strongly correlated with these climate variables. Among the most extreme signals, several overlap with those from GWAS, including SNPs associated with pigmentation and autoimmune diseases. Further, we find an enrichment of strong signals in gene sets related to UV radiation, infection and immunity, and cancer. Our results imply that adaptations to climate shaped the spatial distribution of variation in humans.

Classical studies that examined the global distributions of human physiological traits such as pigmentation, basal metabolic rate, and body shape and size suggested that natural selection related to climate has been important during recent human evolutionary history. We scanned the human genome using data for about 650,000 variants in 61 worldwide populations to look for correlations between allele frequencies and 9 climate variables and found evidence for adaptations to climate at the genome-wide level. In addition, we detected compelling signals for individual SNPs involved in pigmentation and immune response, as well as for pathways related to UV radiation, infection and immunity, and cancer. A particularly appealing aspect of this approach is that we identify a set of candidate advantageous SNPs associated with specific biological hypotheses, which will be useful for follow-up testing. We developed an online resource to browse the results of our data analyses, allowing researchers to quickly assess evidence for selection in a particular genomic region and to compare it across several studies.