An ancient viral epidemic involving host coronavirus interacting genes more than 20,000 years ago in East Asia.

The current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has emphasized the vulnerability of human populations to novel viral pressures, despite the vast array of epidemiological and biomedical tools now available. Notably, modern human genomes contain evolutionary information tracing back tens of thousands of years, which may help identify the viruses that have impacted our ancestors-pointing to which viruses have future pandemic potential. Here, we apply evolutionary analyses to human genomic datasets to recover selection events involving tens of human genes that interact with coronaviruses, including SARS-CoV-2, that likely started more than 20,000 years ago. These adaptive events were limited to the population ancestral to East Asian populations. Multiple lines of functional evidence support an ancient viral selective pressure, and East Asia is the geographical origin of several modern coronavirus epidemics. An arms race with an ancient coronavirus, or with a different virus that happened to use similar interactions as coronaviruses with human hosts, may thus have taken place in ancestral East Asian populations. By learning more about our ancient viral foes, our study highlights the promise of evolutionary information to better predict the pandemics of the future. Importantly, adaptation to ancient viral epidemics in specific human populations does not necessarily imply any difference in genetic susceptibility between different human populations, and the current evidence points toward an overwhelming impact of socioeconomic factors in the case of coronavirus disease 2019 (COVID-19).

Genozip: a universal extensible genomic data compressor.

We present Genozip, a universal and fully featured compression software for genomic data. Genozip is designed to be a general-purpose software and a development framework for genomic compression by providing five core capabilities-universality (support for all common genomic file formats), high compression ratios, speed, feature-richness and extensibility. Genozip delivers high-performance compression for widelyused genomic data formats in genomics research, namely FASTQ, SAM/BAM/CRAM, VCF, GVF, FASTA, PHYLIP and 23andMe formats. Our test results show that Genozip is fast and achieves greatly improved compression ratios, even when the files are already compressed. Further, Genozip is architected with a separation of the Genozip Framework from file-format-specific Segmenters and data-type-specific Codecs. With this, we intend for Genozip to be a general-purpose compression platform where researchers can implement compression for additional file formats, as well as new codecs for data types or fields within files, in the future. We anticipate that this will ultimately increase the visibility and adoption of these algorithms by the user community, thereby accelerating further innovation in this space. AVAILABILITY AND IMPLEMENTATION: Genozip is written in C. The code is open-source and available on http://www.genozip.com. The package is free for non-commercial use. It is distributed through the Conda package manager, github, and as a Docker container on DockerHub. Genozip is tested on Linux, Mac and Windows. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A simple genetic basis of adaptation to a novel thermal environment results in complex metabolic rewiring in Drosophila.

BACKGROUND: Population genetic theory predicts that rapid adaptation is largely driven by complex traits encoded by many loci of small effect. Because large-effect loci are quickly fixed in natural populations, they should not contribute much to rapid adaptation. RESULTS: To investigate the genetic architecture of thermal adaptation - a highly complex trait - we performed experimental evolution on a natural Drosophila simulans population. Transcriptome and respiration measurements reveal extensive metabolic rewiring after only approximately 60 generations in a hot environment. Analysis of genome-wide polymorphisms identifies two interacting selection targets, Sestrin and SNF4Aγ, pointing to AMPK, a central metabolic switch, as a key factor for thermal adaptation. CONCLUSIONS: Our results demonstrate that large-effect loci segregating at intermediate allele frequencies can allow natural populations to rapidly respond to selection. Because SNF4Aγ also exhibits clinal variation in various Drosophila species, we suggest that this large-effect polymorphism is maintained by temporal and spatial temperature variation in natural environments.

Molecular dissection of a natural transposable element invasion.

The first tracking of the dynamics of a natural invasion by a transposable element (TE) provides unprecedented details on the establishment of host defense mechanisms against TEs. We captured a D. simulans population at an early stage of a P-element invasion and studied the spread of the TE in replicated experimentally evolving populations kept under hot and cold conditions. We analyzed the factors controlling the invasion by NGS, RNA-FISH, and gonadal dysgenesis assays. Under hot conditions, the P-element spread rapidly for 20 generations, but no further spread was noted later on. This plateauing of the invasion was mediated by the rapid emergence of P-element-specific piRNAs. Under cold conditions, we observed a lower expression of the P-element and a slower emergence of the piRNA defense, resulting in a three times slower invasion that continued beyond 40 generations. We conclude that the environment is a major factor determining the evolution of TEs in their host.

High rate of translocation-based gene birth on the Drosophila Y chromosome.

The Y chromosome is a unique genetic environment defined by a lack of recombination and male-limited inheritance. The Drosophila Y chromosome has been gradually acquiring genes from the rest of the genome, with only seven Y-linked genes being gained over the past 63 million years (0.12 gene gains per million years). Using a next-generation sequencing (NGS)-powered genomic scan, we show that gene transfers to the Y chromosome are much more common than previously suspected: at least 25 have arisen across three Drosophila species over the past 5.4 million years (1.67 per million years for each lineage). The gene transfer rate is significantly lower in Drosophila melanogaster than in the Drosophila simulans clade, primarily due to Y-linked retrotranspositions being significantly more common in the latter. Despite all Y-linked gene transfers being evolutionarily recent (<1 million years old), only three showed evidence for purifying selection (ω ≤ 0.14). Thus, although the resulting Y-linked functional gene acquisition rate (0.25 new genes per million years) is double the longer-term estimate, the fate of most new Y-linked genes is defined by rapid degeneration and pseudogenization. Our results show that Y-linked gene traffic, and the molecular mechanisms governing these transfers, can diverge rapidly between species, revealing the Drosophila Y chromosome to be more dynamic than previously appreciated. Our analytical method provides a powerful means to identify Y-linked gene transfers and will help illuminate the evolutionary dynamics of the Y chromosome in Drosophila and other species.

Aboriginal mitogenomes reveal 50,000 years of regionalism in Australia.

Aboriginal Australians represent one of the longest continuous cultural complexes known. Archaeological evidence indicates that Australia and New Guinea were initially settled approximately 50 thousand years ago (ka); however, little is known about the processes underlying the enormous linguistic and phenotypic diversity within Australia. Here we report 111 mitochondrial genomes (mitogenomes) from historical Aboriginal Australian hair samples, whose origins enable us to reconstruct Australian phylogeographic history before European settlement. Marked geographic patterns and deep splits across the major mitochondrial haplogroups imply that the settlement of Australia comprised a single, rapid migration along the east and west coasts that reached southern Australia by 49-45 ka. After continent-wide colonization, strong regional patterns developed and these have survived despite substantial climatic and cultural change during the late Pleistocene and Holocene epochs. Remarkably, we find evidence for the continuous presence of populations in discrete geographic areas dating back to around 50 ka, in agreement with the notable Aboriginal Australian cultural attachment to their country.

Ancestral population reconstitution from isofemale lines as a tool for experimental evolution.

Experimental evolution is a powerful tool to study adaptation under controlled conditions. Laboratory natural selection experiments mimic adaptation in the wild with better-adapted genotypes having more offspring. Because the selected traits are frequently not known, adaptation is typically measured as fitness increase by comparing evolved populations against an unselected reference population maintained in a laboratory environment. With adaptation to the laboratory conditions and genetic drift, however, it is not clear to what extent such comparisons provide unbiased estimates of adaptation. Alternatively, ancestral variation could be preserved in isofemale lines that can be combined to reconstitute the ancestral population. Here, we assess the impact of selection on alleles segregating in newly established Drosophila isofemale lines. We reconstituted two populations from isofemale lines and compared them to two original ancestral populations (AP) founded from the same lines shortly after collection. No significant allele frequency changes could be detected between both AP and simulations showed that drift had a low impact compared to Pool-Seq-associated sampling effects. We conclude that laboratory selection on segregating variation in isofemale lines is too weak to have detectable effects, which validates ancestral population reconstitution from isofemale lines as an unbiased approach for measuring adaptation in evolved populations.

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.

Parallel trait adaptation across opposing thermal environments in experimental Drosophila melanogaster populations.

Thermal stress is a pervasive selective agent in natural populations that impacts organismal growth, survival, and reproduction. Drosophila melanogaster exhibits a variety of putatively adaptive phenotypic responses to thermal stress in natural and experimental settings; however, accompanying assessments of fitness are typically lacking. Here, we quantify changes in fitness and known thermal tolerance traits in replicated experimental D. melanogaster populations following more than 40 generations of evolution to either cyclic cold or hot temperatures. By evaluating fitness for both evolved populations alongside a reconstituted starting population, we show that the evolved populations were the best adapted within their respective thermal environments. More strikingly, the evolved populations exhibited increased fitness in both environments and improved resistance to both acute heat and cold stress. This unexpected parallel response appeared to be an adaptation to the rapid temperature changes that drove the cycling thermal regimes, as parallel fitness changes were not observed when tested in a constant thermal environment. Our results add to a small, but growing group of studies that demonstrate the importance of fluctuating temperature changes for thermal adaptation and highlight the need for additional work in this area.

Patterns of linkage disequilibrium and long range hitchhiking in evolving experimental Drosophila melanogaster populations.

Whole-genome resequencing of experimental populations evolving under a specific selection regime has become a popular approach to determine genotype-phenotype maps and understand adaptation to new environments. Despite its conceptual appeal and success in identifying some causative genes, it has become apparent that many studies suffer from an excess of candidate loci. Several explanations have been proposed for this phenomenon, but it is clear that information about the linkage structure during such experiments is needed. Until now only Pool-Seq (whole-genome sequencing of pools of individuals) data were available, which do not provide sufficient information about the correlation between linked sites. We address this problem in two complementary analyses of three replicate Drosophila melanogaster populations evolving to a new hot temperature environment for almost 70 generations. In the first analysis, we sequenced 58 haploid genomes from the founder population and evolved flies at generation 67. We show that during the experiment linkage disequilibrium (LD) increased almost uniformly over much greater distances than typically seen in Drosophila. In the second analysis, Pool-Seq time series data of the three replicates were combined with haplotype information from the founder population to follow blocks of initial haplotypes over time. We identified 17 selected haplotype-blocks that started at low frequencies in the base population and increased in frequency during the experiment. The size of these haplotype-blocks ranged from 0.082 to 4.01 Mb. Moreover, between 42% and 46% of the top candidate single nucleotide polymorphisms from the comparison of founder and evolved populations fell into the genomic region covered by the haplotype-blocks. We conclude that LD in such rising haplotype-blocks results in long range hitchhiking over multiple kilobase-sized regions. LD in such haplotype-blocks is therefore a major factor contributing to an excess of candidate loci. Although modifications of the experimental design may help to reduce the hitchhiking effect and allow for more precise mapping of causative variants, we also note that such haplotype-blocks might be well suited to study the dynamics of selected genomic regions during experimental evolution studies.

Experimental evolution reveals habitat-specific fitness dynamics among Wolbachia clades in Drosophila melanogaster.

The diversity and infection dynamics of the endosymbiont Wolbachia can be influenced by many factors, such as transmission rate, cytoplasmic incompatibility, environment, selection and genetic drift. The interplay of these factors in natural populations can result in heterogeneous infection patterns with substantial differences between populations and strains. The causes of these heterogeneities are not yet understood, partly due to the complexity of natural environments. We present experimental evolution as a new approach to study Wolbachia infection dynamics in replicate populations exposed to a controlled environment. A natural Drosophila melanogaster population infected with strains of Wolbachia belonging to different clades evolved in two laboratory environments (hot and cold) for 1.5 years. In both treatments, the rate of Wolbachia infection increased until fixation. In the hot environment, the relative frequency of different Wolbachia clades remained stable over 37 generations. In the cold environment, however, we observed marked changes in the composition of the Wolbachia population: within 15 generations, one Wolbachia clade increased more than 50% in frequency, whereas the other two clades decreased in frequency, resulting in the loss of one clade. The frequency change was highly reproducible not only among replicates, but also when flies that evolved for 42 generations in the hot environment were transferred to the cold environment. These results document how environmental factors can affect the composition of Wolbachia in D. melanogaster. The high reproducibility of the pattern suggests that experimental evolution studies can efficiently determine the functional basis of habitat-specific fitness among Wolbachia strains.

Massive habitat-specific genomic response in D. melanogaster populations during experimental evolution in hot and cold environments.

Experimental evolution in combination with whole-genome sequencing (evolve and resequence [E&R]) is a promising approach to define the genotype-phenotype map and to understand adaptation in evolving populations. Many previous studies have identified a large number of putative selected sites (i.e., candidate loci), but it remains unclear to what extent these loci are genuine targets of selection or experimental noise. To address this question, we exposed the same founder population to two different selection regimes-a hot environment and a cold environment-and quantified the genomic response in each. We detected large numbers of putative selected loci in both environments, albeit with little overlap between the two sets of candidates, indicating that most resulted from habitat-specific selection. By quantifying changes across multiple independent biological replicates, we demonstrate that most of the candidate SNPs were false positives that were linked to selected sites over distances much larger than the typical linkage disequilibrium range of Drosophila melanogaster. We show that many of these mid- to long-range associations were attributable to large segregating inversions and confirm by computer simulations that such patterns could be readily replicated when strong selection acts on rare haplotypes. In light of our findings, we outline recommendations to improve the performance of future Drosophila E&R studies which include using species with negligible inversion loads, such as D. mauritiana and D. simulans, instead of D. melanogaster.