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

Analyses of 600+ insect genomes reveal repetitive element dynamics and highlight biodiversity-scale repeat annotation challenges

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal.

Satellite DNAs-From Localized to Highly Dispersed Genome Components

According to the established classical view, satellite DNAs are defined as abundant non-coding DNA sequences repeated in tandem that build long arrays located in heterochromatin. Advances in sequencing methodologies and development of specialized bioinformatics tools enabled defining a collection of all repetitive DNAs and satellite DNAs in a genome, the repeatome and the satellitome, respectively, as well as their reliable annotation on sequenced genomes. Supported by various non-model species included in recent studies, the patterns of satellite DNAs and satellitomes as a whole showed much more diversity and complexity than initially thought. Differences are not only in number and abundance of satellite DNAs but also in their distribution across the genome, array length, interspersion patterns, association with transposable elements, localization in heterochromatin and/or in euchromatin. In this review, we compare characteristic organizational features of satellite DNAs and satellitomes across different animal and plant species in order to summarize organizational forms and evolutionary processes that may lead to satellitomes' diversity and revisit some basic notions regarding repetitive DNA landscapes in genomes.

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.

Dissecting cricket genomes for the advancement of entomology and entomophagy

Significant advances in biophysical methods such as next-generation sequencing technologies have now opened the way to conduct evolutionary and applied research based on the genomic information of greatly diverse insects. Crickets belonging to Orthoptera (Insecta: Polyneoptera), one of the most flourishing groups of insects, have contributed to the development of multiple scientific fields including developmental biology and neuroscience and have been attractive targets in evolutionary ecology for their diverse ecological niches. In addition, crickets have recently gained recognition as food and feed. However, the genomic information underlying their biological basis and application research toward breeding is currently underrepresented. In this review, we summarize the progress of genomics of crickets. First, we outline the phylogenetic position of crickets in insects and then introduce recent studies on cricket genomics and transcriptomics in a variety of fields. Furthermore, we present findings from our analysis of polyneopteran genomes, with a particular focus on their large genome sizes, chromosome number, and repetitive sequences. Finally, how the cricket genome can be beneficial to the food industry is discussed. This review is expected to enhance greater recognition of how important the cricket genomes are to the multiple biological fields and how basic research based on cricket genome information can contribute to tackling global food security.

Genome of the pincer wasp Gonatopus flavifemur reveals unique venom evolution and a dual adaptation to parasitism and predation

Background

Hymenoptera comprise extremely diverse insect species with extensive variation in their life histories. The Dryinidae, a family of solitary wasps of Hymenoptera, have evolved innovations that allow them to hunt using venom and a pair of chelae developed from the fore legs that can grasp prey. Dryinidae larvae are also parasitoids of Auchenorrhyncha, a group including common pests such as planthoppers and leafhoppers. Both of these traits make them effective and valuable for pest control, but little is yet known about the genetic basis of its dual adaptation to parasitism and predation.

Results

We sequenced and assembled a high-quality genome of the dryinid wasp Gonatopus flavifemur, which at 636.5 Mb is larger than most hymenopterans. The expansion of transposable elements, especially DNA transposons, is a major contributor to the genome size enlargement. Our genome-wide screens reveal a number of positively selected genes and rapidly evolving proteins involved in energy production and motor activity, which may contribute to the predatory adaptation of dryinid wasp. We further show that three female-biased, reproductive-associated yellow genes, in response to the prey feeding behavior, are significantly elevated in adult females, which may facilitate the egg production. Venom is a powerful weapon for dryinid wasp during parasitism and predation. We therefore analyze the transcriptomes of venom glands and describe specific expansions in venom Idgf-like genes and neprilysin-like genes. Furthermore, we find the LWS2-opsin gene is exclusively expressed in male G. flavifemur, which may contribute to partner searching and mating.

Conclusions

Our results provide new insights into the genome evolution, predatory adaptation, venom evolution, and sex-biased genes in G. flavifemur, and present genomic resources for future in-depth comparative analyses of hymenopterans that may benefit pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01081-6.

Satellitome Analysis of the Pacific Oyster Crassostrea gigas Reveals New Pattern of Satellite DNA Organization, Highly Scattered across the Genome

Several features already qualified the invasive bivalve species Crassostrea gigas as a valuable non-standard model organism in genome research. C. gigas is characterized by the low contribution of satellite DNAs (satDNAs) vs. mobile elements and has an extremely low amount of heterochromatin, predominantly built of DNA transposons. In this work, we have identified 52 satDNAs composing the satellitome of C. gigas and constituting about 6.33% of the genome. Satellitome analysis reveals unusual, highly scattered organization of relatively short satDNA arrays across the whole genome. However, peculiar chromosomal distribution and densities are specific for each satDNA. The inspection of the organizational forms of the 11 most abundant satDNAs shows association with constitutive parts of Helitron mobile elements. Nine of the inspected satDNAs are dominantly found in mobile element-associated form, two mostly appear standalone, and only one is present exclusively as Helitron-associated sequence. The Helitron-related satDNAs appear in more chromosomes than other satDNAs, indicating that these mobile elements could be leading satDNA propagation in C. gigas. No significant accumulation of satDNAs on certain chromosomal positions was detected in C. gigas, thus establishing a novel pattern of satDNA organization on the genome level.

Broad geographic sampling reveals the shared basis and environmental correlates of seasonal adaptation in Drosophila

To advance our understanding of adaptation to temporally varying selection pressures, we identified signatures of seasonal adaptation occurring in parallel among Drosophila melanogaster populations. Specifically, we estimated allele frequencies genome-wide from flies sampled early and late in the growing season from 20 widely dispersed populations. We identified parallel seasonal allele frequency shifts across North America and Europe, demonstrating that seasonal adaptation is a general phenomenon of temperate fly populations. Seasonally fluctuating polymorphisms are enriched in large chromosomal inversions, and we find a broad concordance between seasonal and spatial allele frequency change. The direction of allele frequency change at seasonally variable polymorphisms can be predicted by weather conditions in the weeks prior to sampling, linking the environment and the genomic response to selection. Our results suggest that fluctuating selection is an important evolutionary force affecting patterns of genetic variation in Drosophila.

Genomic Tackling of Human Satellite DNA: Breaking Barriers through Time

(Peri)centromeric repetitive sequences and, more specifically, satellite DNA (satDNA) sequences, constitute a major human genomic component. SatDNA sequences can vary on a large number of features, including nucleotide composition, complexity, and abundance. Several satDNA families have been identified and characterized in the human genome through time, albeit at different speeds. Human satDNA families present a high degree of sub-variability, leading to the definition of various subfamilies with different organization and clustered localization. Evolution of satDNA analysis has enabled the progressive characterization of satDNA features. Despite recent advances in the sequencing of centromeric arrays, comprehensive genomic studies to assess their variability are still required to provide accurate and proportional representation of satDNA (peri)centromeric/acrocentric short arm sequences. Approaches combining multiple techniques have been successfully applied and seem to be the path to follow for generating integrated knowledge in the promising field of human satDNA biology.

Satellite DNA-Mediated Gene Expression Regulation: Physiological and Evolutionary Implication

Satellite DNAs are tandemly repeated sequences organized in large clusters within (peri)centromeric and/or subtelomeric heterochromatin. However, in many species, satellite DNAs are not restricted to heterochromatin but are also dispersed as short arrays within euchromatin. Such genomic organization together with transcriptional activity seems to be a prerequisite for the gene-modulatory effect of satellite DNAs which was first demonstrated in the beetle Tribolium castaneum upon heat stress. Namely, enrichment of a silent histone mark at euchromatic repeats of a major beetle satellite DNA results in epigenetic silencing of neighboring genes. In addition, human satellite III transcripts induced by heat shock contribute to genome-wide gene silencing, providing protection against stress-induced cell death. Gene silencing mediated by satellite RNA was also shown to be fundamental for the early embryonic development of the mosquito Aedes aegypti. Apart from a physiological role during embryogenesis and heat stress response, activation of satellite DNAs in terms of transcription and proliferation can have an evolutionary impact. Spreading of satellite repeats throughout euchromatin promotes the variation of epigenetic landscapes and gene expression diversity, contributing to the evolution of gene regulatory networks and to genome adaptation in fluctuating environmental conditions.

Comparative analysis of morabine grasshopper genomes reveals highly abundant transposable elements and rapidly proliferating satellite DNA repeats

BACKGROUND

Repetitive DNA sequences, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), collectively called the "repeatome", are found in high proportion in organisms across the Tree of Life. Grasshoppers have large genomes, averaging 9 Gb, that contain a high proportion of repetitive DNA, which has hampered progress in assembling reference genomes. Here we combined linked-read genomics with transcriptomics to assemble, characterize, and compare the structure of repetitive DNA sequences in four chromosomal races of the morabine grasshopper Vandiemenella viatica species complex and determine their contribution to genome evolution.

RESULTS

We obtained linked-read genome assemblies of 2.73-3.27 Gb from estimated genome sizes of 4.26-5.07 Gb DNA per haploid genome of the four chromosomal races of V. viatica. These constitute the third largest insect genomes assembled so far. Combining complementary annotation tools and manual curation, we found a large diversity of TEs and satDNAs, constituting 66 to 75% per genome assembly. A comparison of sequence divergence within the TE classes revealed massive accumulation of recent TEs in all four races (314-463 Mb per assembly), indicating that their large genome sizes are likely due to similar rates of TE accumulation. Transcriptome sequencing showed more biased TE expression in reproductive tissues than somatic tissues, implying permissive transcription in gametogenesis. Out of 129 satDNA families, 102 satDNA families were shared among the four chromosomal races, which likely represent a diversity of satDNA families in the ancestor of the V. viatica chromosomal races. Notably, 50 of these shared satDNA families underwent differential proliferation since the recent diversification of the V. viatica species complex.

CONCLUSION

This in-depth annotation of the repeatome in morabine grasshoppers provided new insights into the genome evolution of Orthoptera. Our TEs analysis revealed a massive recent accumulation of TEs equivalent to the size of entire Drosophila genomes, which likely explains the large genome sizes in grasshoppers. Despite an overall high similarity of the TE and satDNA diversity between races, the patterns of TE expression and satDNA proliferation suggest rapid evolution of grasshopper genomes on recent timescales.

Decoding the Role of Satellite DNA in Genome Architecture and Plasticity-An Evolutionary and Clinical Affair

Repetitive DNA is a major organizational component of eukaryotic genomes, being intrinsically related with their architecture and evolution. Tandemly repeated satellite DNAs (satDNAs) can be found clustered in specific heterochromatin-rich chromosomal regions, building vital structures like functional centromeres and also dispersed within euchromatin. Interestingly, despite their association to critical chromosomal structures, satDNAs are widely variable among species due to their high turnover rates. This dynamic behavior has been associated with genome plasticity and chromosome rearrangements, leading to the reshaping of genomes. Here we present the current knowledge regarding satDNAs in the light of new genomic technologies, and the challenges in the study of these sequences. Furthermore, we discuss how these sequences, together with other repeats, influence genome architecture, impacting its evolution and association with disease.

Genomic divergence and adaptive convergence in Drosophila simulans from Evolution Canyon, Israel

Biodiversity refugia formed by unique features of the Mediterranean arid landscape, such as the dramatic ecological contrast of "Evolution Canyon," provide a natural laboratory in which local adaptations to divergent microclimate conditions can be investigated. Significant insights have been provided by studies of Drosophila melanogaster diversifying along the thermal gradient in Evolution Canyon, but a comparative framework to survey adaptive convergence across sister species at the site has been lacking. To fill this void, we present an analysis of genomic polymorphism and evolutionary divergence of Drosophila simulans, a close relative of Drosophila melanogaster with which it co-occurs on both slopes of the canyon. Our results show even deeper interslope divergence in D. simulans than in D. melanogaster, with extensive signatures of selective sweeps present in flies from both slopes but enhanced in the population from the hotter and drier south-facing slope. Interslope divergence was enriched for genes related to electrochemical balance and transmembrane transport, likely in response to increased selection for dehydration resistance on the hotter slope. Both species shared genomic regions that underwent major selective sweeps, but the overall level of adaptive convergence was low, demonstrating no shortage of alternative genomic solutions to cope with the challenges of the microclimate contrast. Mobile elements were a major source of genetic polymorphism and divergence, affecting all parts of the genome, including coding sequences of mating behavior-related genes.

Characterization of the satellitome in lower vascular plants: the case of the endangered fern Vandenboschia speciosa

BACKGROUND AND AIMS

Vandenboschia speciosa is a highly vulnerable fern species, with a large genome (10.5 Gb). Haploid gametophytes and diploid sporophytes are perennial, can reproduce vegetatively, and certain populations are composed only of independent gametophytes. These features make this fern a good model: (1) for high-throughput analysis of satellite DNA (satDNA) to investigate possible evolutionary trends in satDNA sequence features; (2) to determine the relative contribution of satDNA and other repetitive DNAs to its large genome; and (3) to analyse whether the reproduction mode or phase alternation between long-lasting haploid and diploid stages influences satDNA abundance or divergence.

METHODS

We analysed the repetitive fraction of the genome of this species in three different populations (one comprised only of independent gametophytes) using Illumina sequencing and bioinformatic analysis with RepeatExplorer and satMiner.

KEY RESULTS

The satellitome of V. speciosa is composed of 11 satDNA families, most of them showing a short repeat length and being A + T rich. Some satDNAs had complex repeats composed of sub-repeats, showing high similarity to shorter satDNAs. Three families had particular structural features and highly conserved motifs. SatDNA only amounts to approx. 0.4 % of its genome. Likewise, microsatellites do not represent more than 2 %, but transposable elements (TEs) represent approx. 50 % of the sporophytic genomes. We found high resemblance in satDNA abundance and divergence between both gametophyte and sporophyte samples from the same population and between populations.

CONCLUSIONS

(1) Longer (and older) satellites in V. speciosa have a higher A + T content and evolve from shorter ones and, in some cases, microsatellites were a source of new satDNAs; (2) the satellitome does not explain the huge genome size in this species while TEs are the major repetitive component of the V. speciosa genome and mostly contribute to its large genome; and (3) reproduction mode or phase alternation between gametophytes and sporophytes does not entail accumulation or divergence of satellites.

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.

Diversity and evolution of the transposable element repertoire in arthropods with particular reference to insects

BACKGROUND

Transposable elements (TEs) are a major component of metazoan genomes and are associated with a variety of mechanisms that shape genome architecture and evolution. Despite the ever-growing number of insect genomes sequenced to date, our understanding of the diversity and evolution of insect TEs remains poor.

RESULTS

Here, we present a standardized characterization and an order-level comparison of arthropod TE repertoires, encompassing 62 insect and 11 outgroup species. The insect TE repertoire contains TEs of almost every class previously described, and in some cases even TEs previously reported only from vertebrates and plants. Additionally, we identified a large fraction of unclassifiable TEs. We found high variation in TE content, ranging from less than 6% in the antarctic midge (Diptera), the honey bee and the turnip sawfly (Hymenoptera) to more than 58% in the malaria mosquito (Diptera) and the migratory locust (Orthoptera), and a possible relationship between the content and diversity of TEs and the genome size.

CONCLUSION

While most insect orders exhibit a characteristic TE composition, we also observed intraordinal differences, e.g., in Diptera, Hymenoptera, and Hemiptera. Our findings shed light on common patterns and reveal lineage-specific differences in content and evolution of TEs in insects. We anticipate our study to provide the basis for future comparative research on the insect TE repertoire.

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.

Rapid Increase in Genome Size as a Consequence of Transposable Element Hyperactivity in Wood-White (Leptidea) Butterflies

Characterizing and quantifying genome size variation among organisms and understanding if genome size evolves as a consequence of adaptive or stochastic processes have been long-standing goals in evolutionary biology. Here, we investigate genome size variation and association with transposable elements (TEs) across lepidopteran lineages using a novel genome assembly of the common wood-white (Leptidea sinapis) and population re-sequencing data from both L. sinapis and the closely related L. reali and L. juvernica together with 12 previously available lepidopteran genome assemblies. A phylogenetic analysis confirms established relationships among species, but identifies previously unknown intraspecific structure within Leptidea lineages. The genome assembly of L. sinapis is one of the largest of any lepidopteran taxon so far (643 Mb) and genome size is correlated with abundance of TEs, both in Lepidoptera in general and within Leptidea where L. juvernica from Kazakhstan has considerably larger genome size than any other Leptidea population. Specific TE subclasses have been active in different Lepidoptera lineages with a pronounced expansion of predominantly LINEs, DNA elements, and unclassified TEs in the Leptidea lineage after the split from other Pieridae. The rate of genome expansion in Leptidea in general has been in the range of four Mb/Million year (My), with an increase in a particular L. juvernica population to 72 Mb/My. The considerable differences in accumulation rates of specific TE classes in different lineages indicate that TE activity plays a major role in genome size evolution in butterflies and moths.

Satellite DNA: An Evolving Topic

Satellite DNA represents one of the most fascinating parts of the repetitive fraction of the eukaryotic genome. Since the discovery of highly repetitive tandem DNA in the 1960s, a lot of literature has extensively covered various topics related to the structure, organization, function, and evolution of such sequences. Today, with the advent of genomic tools, the study of satellite DNA has regained a great interest. Thus, Next-Generation Sequencing (NGS), together with high-throughput in silico analysis of the information contained in NGS reads, has revolutionized the analysis of the repetitive fraction of the eukaryotic genomes. The whole of the historical and current approaches to the topic gives us a broad view of the function and evolution of satellite DNA and its role in chromosomal evolution. Currently, we have extensive information on the molecular, chromosomal, biological, and population factors that affect the evolutionary fate of satellite DNA, knowledge that gives rise to a series of hypotheses that get on well with each other about the origin, spreading, and evolution of satellite DNA. In this paper, I review these hypotheses from a methodological, conceptual, and historical perspective and frame them in the context of chromosomal organization and evolution.

Quantitative sequence characterization for repetitive DNA content in the supernumerary chromosome of the migratory locust

Repetitive DNA is a major component in most eukaryotic genomes but is ignored in most genome sequencing projects. Here, we report the quantitative composition in repetitive DNA for a supernumerary (B) chromosome, in the migratory locust (Locusta migratoria), by Illumina sequencing of genomic DNA from B-carrying and B-lacking individuals and DNA obtained from a microdissected B chromosome, as well as the physical mapping of some elements. B chromosome DNA of 94.9% was repetitive, in high contrast with the 64.1% of standard (A) chromosomes. B chromosomes are enriched in satellite DNA (satDNA) (65.2% of B-DNA), with a single satellite (LmiSat02-176) comprising 55% of the B. Six satDNAs were visualized by FISH on the B chromosome, and the only A chromosome carrying all these satellites was autosome 9, pointing to this chromosome, along with autosome 8 (which shares histone genes with the B) as putative ancestors of the B chromosome. We found several transposable elements (TEs) showing nucleotidic variation specific to B-carrying individuals, which was also present in B-carrying transcriptomes. Remarkably, an interstitial region of the B chromosome included a 17 kb chimera composed of 29 different TEs, suggesting reiterative TE insertion in this B chromosome region.

High-throughput analysis of the satellitome illuminates satellite DNA evolution

Satellite DNA (satDNA) is a major component yet the great unknown of eukaryote genomes and clearly underrepresented in genome sequencing projects. Here we show the high-throughput analysis of satellite DNA content in the migratory locust by means of the bioinformatic analysis of Illumina reads with the RepeatExplorer and RepeatMasker programs. This unveiled 62 satDNA families and we propose the term "satellitome" for the whole collection of different satDNA families in a genome. The finding that satDNAs were present in many contigs of the migratory locust draft genome indicates that they show many genomic locations invisible by fluorescent in situ hybridization (FISH). The cytological pattern of five satellites showing common descent (belonging to the SF3 superfamily) suggests that non-clustered satDNAs can become into clustered through local amplification at any of the many genomic loci resulting from previous dissemination of short satDNA arrays. The fact that all kinds of satDNA (micro- mini- and satellites) can show the non-clustered and clustered states suggests that all these elements are mostly similar, except for repeat length. Finally, the presence of VNTRs in bacteria, showing similar properties to non-clustered satDNAs in eukaryotes, suggests that this kind of tandem repeats show common properties in all living beings.

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.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

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.

Multilevel Selection Theory and the Evolutionary Functions of Transposable Elements

One of several issues at play in the renewed debate over “junk DNA” is the organizational level at which genomic features might be seen as selected, and thus to exhibit function, as etiologically defined. The intuition frequently expressed by molecular geneticists that junk DNA is functional because it serves to “speed evolution” or as an “evolutionary repository” could be recast as a claim about selection between species (or clades) rather than within them, but this is not often done. Here, we review general arguments for the importance of selection at levels above that of organisms in evolution, and develop them further for a common genomic feature: the carriage of transposable elements (TEs). In many species, not least our own, TEs comprise a large fraction of all nuclear DNA, and whether they individually or collectively contribute to fitness—or are instead junk— is a subject of ongoing contestation. Even if TEs generally owe their origin to selfish selection at the lowest level (that of genomes), their prevalence in extant organisms and the prevalence of extant organisms bearing them must also respond to selection within species (on organismal fitness) and between species (on rates of speciation and extinction). At an even higher level, the persistence of clades may be affected (positively or negatively) by TE carriage. If indeed TEs speed evolution, it is at these higher levels of selection that such a function might best be attributed to them as a class.

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.