The easy road to genome‐wide medium density SNP screening in a non‐model species: development and application of a 10 K SNP‐chip for the house sparrow (Passer domesticus)

With the advent of next generation sequencing, new avenues have opened to study genomics in wild populations of non‐model species. Here, we describe a successful approach to a genome‐wide medium density Single Nucleotide Polymorphism (SNP) panel in a non‐model species, the house sparrow (Passer domesticus), through the development of a 10 K Illumina iSelect HD BeadChip. Genomic DNA and cDNA derived from six individuals were sequenced on a 454 GS FLX system and generated a total of 1.2 million sequences, in which SNPs were detected. As no reference genome exists for the house sparrow, we used the zebra finch (Taeniopygia guttata) reference genome to determine the most likely position of each SNP. The 10 000 SNPs on the SNP‐chip were selected to be distributed evenly across 31 chromosomes, giving on average one SNP per 100 000 bp. The SNP‐chip was screened across 1968 individual house sparrows from four island populations. Of the original 10 000 SNPs, 7413 were found to be variable, and 99% of these SNPs were successfully called in at least 93% of all individuals. We used the SNP‐chip to demonstrate the ability of such genome‐wide marker data to detect population sub‐division, and compared these results to similar analyses using microsatellites. The SNP‐chip will be used to map Quantitative Trait Loci (QTL) for fitness‐related phenotypic traits in natural populations.     

Seasonal variation in life history traits in two Drosophila species

Seasonal environmental heterogeneity is cyclic, persistent and geographically widespread. In species that reproduce multiple times annually, environmental changes across seasonal time may create different selection regimes that may shape the population ecology and life history adaptation in these species. Here, we investigate how two closely related species of Drosophila in a temperate orchard respond to environmental changes across seasonal time. Natural populations of Drosophila melanogaster and Drosophila simulans were sampled at four timepoints from June through November to assess seasonal change in fundamental aspects of population dynamics as well as life history traits. D. melanogaster exhibit pronounced change across seasonal time: early in the season, the population is inferred to be uniformly young and potentially represents the early generation following overwintering survivorship. D. melanogaster isofemale lines derived from the early population and reared in a common garden are characterized by high tolerance to a variety of stressors as well as a fast rate of development in the laboratory environment that declines across seasonal time. In contrast, wild D. simulans populations were inferred to be consistently heterogeneous in age distribution across seasonal collections; only starvation tolerance changed predictably over seasonal time in a parallel manner as in D. melanogaster. These results suggest fundamental differences in population and evolutionary dynamics between these two taxa associated with seasonal heterogeneity in environmental parameters and associated selection pressures.     

Development of an SNP‐based parentage verification panel for lovebirds

The genus Agapornis, or lovebirds, are popular pet parrots worldwide. Currently, breeders are dependent on pedigree records as a selection tool as no molecular parentage verification test is available for any of the nine species. The A. roseicollis reference genome was recently assembled. This was followed by the sequencing of the whole genomes of the parents of the reference genome individual at 30× coverage. The parents’ reads were mapped against the reference genome to identify SNPs. Over 1.6 million SNPs, shared between the parents, were discovered using the Genome Analysis Toolkit pipeline. SNPs were filtered to a panel of 480 SNPs based on Genome Analysis Toolkit parameters. The panel of 480 SNPs was genotyped in a population of 960 lovebirds across seven species. A panel of 262 SNPs was compiled that included SNPs successfully amplified across all species. The 262‐SNP panel was reduced based on the observed heterozygosity (H_O) and minor allele frequency (MAF) values per SNP to include the lowest number of SNPs with the highest exclusion power for parentage verification. Two smaller panels consisting of 195 SNPs with MAF and H_O values >0.1 and 40 SNPs with MAF and H_O values >0.3, were constructed. The panels were verified using 43 families from different species with known relationships to evaluate the exclusion power of each panel. The 195 SNP panel with an average exclusion probability of 99.9% and MAF and H_O values >0.1 was proposed as the routine Agapornis parentage verification panel.     

Are Drosophila preferences for yeasts stable or contextual?

Whether there are general mechanisms, driving interspecific chemical communication is uncertain. Saccharomycetaceae yeast and Drosophila fruit flies, both extensively studied research models, share the same fruit habitat, and it has been suggested their interaction comprises a facultative mutualism that is instigated and maintained by yeast volatiles. Using choice tests, experimental evolution, and volatile analyses, we investigate the maintenance of this relationship and reveal little consistency between behavioral responses of two isolates of sympatric Drosophila species. While D. melanogaster was attracted to a range of different Saccharomycetaceae yeasts and this was independent of fruit type, D. simulans preference appeared specific to a particular S. cerevisiae genotype isolated from a vineyard fly population. This response, however, was not consistent across fruit types and is therefore context‐dependent. In addition, D. simulans attraction to an individual S. cerevisiae isolate was pliable over ecological timescales. Volatile candidates were analyzed to identify a common signal for yeast attraction, and while D. melanogaster generally responded to fermentation profiles, D. simulans preference was more discerning and likely threshold‐dependent. Overall, there is no strong evidence to support the idea of bespoke interactions with specific yeasts for either of these Drosophila genotypes. Rather the data support the idea Drosophila are generally adapted to sense and locate fruits infested by a range of fungal microbes and/or that yeast–Drosophila interactions may evolve rapidly.     

Identification and validation of single nucleotide polymorphisms as tools to detect hybridization and population structure in freshwater stingrays

Single nucleotide polymorphism (SNP) markers were identified and validated for two stingrays species, Potamotrygon motoro and Potamotrygon falkneri, using double digest restriction‐site associated DNA (ddRAD) reads using 454‐Roche technology. A total of 226 774 reads (65.5 Mb) were obtained (mean read length 289 ± 183 bp) detecting a total of 5399 contigs (mean contig length: 396 ± 91 bp). Mining this data set, a panel of 143 in silico SNPs was selected. Eighty‐two of these SNPs were successfully validated and 61 were polymorphic: 14 in P. falkneri, 21 in P. motoro, 3 in both species and 26 fixed for alternative variants in both species, thus being useful for population analyses and hybrid detection.     

Development and evaluation of a genome‐wide Coffee 8.5K SNP array and its application for high‐density genetic mapping and for investigating the origin of Coffea arabica L.

Coffee species such as Coffea canephora P. (Robusta) and C. arabica L. (Arabica) are important cash crops in tropical regions around the world. C. arabica is an allotetraploid (2n = 4x = 44) originating from a hybridization event of the two diploid species C. canephora and C. eugenioides (2n = 2x = 22). Interestingly, these progenitor species harbour a greater level of genetic variability and are an important source of genes to broaden the narrow Arabica genetic base. Here, we describe the development, evaluation and use of a single‐nucleotide polymorphism (SNP) array for coffee trees. A total of 8580 unique and informative SNPs were selected from C. canephora and C. arabica sequencing data, with 40% of the SNP located in annotated genes. In particular, this array contains 227 markers associated to 149 genes and traits of agronomic importance. Among these, 7065 SNPs (~82.3%) were scorable and evenly distributed over the genome with a mean distance of 54.4 Kb between markers. With this array, we improved the Robusta high‐density genetic map by adding 1307 SNP markers, whereas 945 SNPs were found segregating in the Arabica mapping progeny. A panel of C. canephora accessions was successfully discriminated and over 70% of the SNP markers were transferable across the three species. Furthermore, the canephora‐derived subgenome of C. arabica was shown to be more closely related to C. canephora accessions from northern Uganda than to other current populations. These validated SNP markers and high‐density genetic maps will be useful to molecular genetics and for innovative approaches in coffee breeding.     

Linked‐read sequencing enables haplotype‐resolved resequencing at population scale

The feasibility to sequence entire genomes of virtually any organism provides unprecedented insights into the evolutionary history of populations and species. Nevertheless, many population genomic inferences – including the quantification and dating of admixture, introgression and demographic events, and inference of selective sweeps – are still limited by the lack of high‐quality haplotype information. The newest generation of sequencing technology now promises significant progress. To establish the feasibility of haplotype‐resolved genome resequencing at population scale, we investigated properties of linked‐read sequencing data of songbirds of the genus Oenanthe across a range of sequencing depths. Our results based on the comparison of downsampled (25×, 20×, 15×, 10×, 7×, and 5×) with high‐coverage data (46–68×) of seven bird genomes mapped to a reference suggest that phasing contiguities and accuracies adequate for most population genomic analyses can be reached already with moderate sequencing effort. At 15× coverage, phased haplotypes span about 90% of the genome assembly, with 50% and 90% of phased sequences located in phase blocks longer than 1.25–4.6 Mb (N50) and 0.27–0.72 Mb (N90). Phasing accuracy reaches beyond 99% starting from 15× coverage. Higher coverages yielded higher contiguities (up to about 7 Mb/1 Mb [N50/N90] at 25× coverage), but only marginally improved phasing accuracy. Phase block contiguity improved with input DNA molecule length; thus, higher‐quality DNA may help keeping sequencing costs at bay. In conclusion, even for organisms with gigabase‐sized genomes like birds, linked‐read sequencing at moderate depth opens an affordable avenue towards haplotype‐resolved genome resequencing at population scale.     

Complexity of the genetic basis of ageing in nature revealed by a clinal study of lifespan and methuselah,a gene for ageing,in Drosophila from eastern Australia

Clinal studies are a powerful tool for understanding the genetic basis of climatic adaptation. However, while clines in quantitative traits and genetic polymorphisms have been observed within and across continents, few studies have attempted to demonstrate direct links between them. The gene methuselah in Drosophila has been shown to have a major effect on stress response and longevity phenotypes based largely on laboratory studies of induced mutations in the mth gene. Clinal patterns in the most common mth haplotype and for lifespan (both increasing with latitude) have been observed in North American populations of D. melanogaster, implicating climatic selection. While these clinal patterns have led some to suggest that mth influences ageing in natural populations, limited evidence on the association between the two has so far been collected. Here, we describe a significant cline in the mth haplotype in eastern Australian D. melanogaster populations that parallel the cline in North America. We also describe a cline in mth gene expression. These findings further support the idea that mth is itself under selection. In contrast, we show that lifespan has a strong nonlinear clinal pattern, increasing southwards from the tropics, but then decreasing again from mid‐latitudes. Furthermore, in association studies, we find no evidence for a direct link between mth haplotype and lifespan. Thus, while our data support a role for mth variation being under natural selection, we found no link to naturally occurring variation in lifespan and ageing in Australian populations of D. melanogaster. Our results indicate that the mth locus likely has genetic background and environment‐specific effects.     

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high‐coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single‐nucleotide polymorphisms (SNPs), of which 1 590 653 high‐confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype‐restricted, non‐synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near‐reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up‐CNVs), 112 630 genes present in fewer copies than in the reference genome (down‐CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland‐specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.     

Spatial dynamics and mixing of bluefin tuna in the Atlantic Ocean and Mediterranean Sea revealed using next‐generation sequencing

The Atlantic bluefin tuna is a highly migratory species emblematic of the challenges associated with shared fisheries management. In an effort to resolve the species’ stock dynamics, a genomewide search for spatially informative single nucleotide polymorphisms (SNPs) was undertaken, by way of sequencing reduced representation libraries. An allele frequency approach to SNP discovery was used, combining the data of 555 larvae and young‐of‐the‐year (LYOY) into pools representing major geographical areas and mapping against a newly assembled genomic reference. From a set of 184,895 candidate loci, 384 were selected for validation using 167 LYOY. A highly discriminatory genotyping panel of 95 SNPs was ultimately developed by selecting loci with the most pronounced differences between western Atlantic and Mediterranean Sea LYOY. The panel was evaluated by genotyping a different set of LYOY (n = 326), and from these, 77.8% and 82.1% were correctly assigned to western Atlantic and Mediterranean Sea origins, respectively. The panel revealed temporally persistent differentiation among LYOY from the western Atlantic and Mediterranean Sea (F_ST = 0.008, p = .034). The composition of six mixed feeding aggregations in the Atlantic Ocean and Mediterranean Sea was characterized using genotypes from medium (n = 184) and large (n = 48) adults, applying population assignment and mixture analyses. The results provide evidence of persistent population structuring across broad geographic areas and extensive mixing in the Atlantic Ocean, particularly in the mid‐Atlantic Bight and Gulf of St. Lawrence. The genomic reference and genotyping tools presented here constitute novel resources useful for future research and conservation efforts.     

Development of Anncaliia algerae (Microsporidia) in Drosophila melanogaster

Representatives of the genus Anncaliia are known as natural parasites of dipteran and coleopteran insects, amphipod crustaceans, but also humans, primarily with immunodeficiency. Anncaliia algerae‐caused fatal myositis is considered as an emergent infectious disease in humans. A. (=Nosema, Brachiola) algerae, the best studied species of the genus, demonstrates the broadest among microsporidia range of natural and experimental hosts, but it has never been propagated in Drosophila. We present ultrastructural analysis of development of A. algerae in visceral muscles and adipocytes of Drosophila melanogaster 2 wk after per oral experimental infection. We observed typical to Anncaliia spp. features of ultrastructure and cell pathology including spore morphology, characteristic extensions of the plasma membrane, and presence of “ridges” and appendages of tubular material at proliferative stages. Anncaliia algerae development in D. melanogaster was particularly similar to one of A. algerae and A.(Brachiola) vesicularum in humans with acute myositis. Given D. melanogaster is currently the most established genetic model, with a fully sequenced genome and easily available transgenic forms and genomic markers, a novel host–parasite system might provide new genetic tools to investigate host–pathogen interactions of A. algerae, as well to test antimicrosporidia drugs.     

Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy

Increased inbreeding is an inevitable consequence of selection in livestock populations. The analysis of high‐density single nucleotide polymorphisms (SNPs) facilitates the identification of long and uninterrupted runs of homozygosity (ROH) that can be used to identify chromosomal regions that are identical by descent. In this work, the distribution of ROH of different lengths in five Italian cattle breeds is described. A total of 4095 bulls from five cattle breeds (2093 Italian Holstein, 749 Italian Brown, 364 Piedmontese, 410 Marchigiana and 479 Italian Simmental) were genotyped at 54K SNP loci. ROH were identified and used to estimate molecular inbreeding coefficients (F_ROH), which were compared with inbreeding coefficients estimated from pedigree information (F_PED) and using the genomic relationship matrix (F_GRM). The average number of ROH per animal ranged from 54 ± 7.2 in Piedmontese to 94.6 ± 11.6 in Italian Brown. The highest number of short ROH (related to ancient consanguinity) was found in Piedmontese, followed by Simmental. The Italian Brown and Holstein had a higher proportion of longer ROH distributed across the whole genome, revealing recent inbreeding. The F_PED were moderately correlated with F_ROH > 1 Mb (0.662, 0.700 and 0.669 in Italian Brown, Italian Holstein and Italian Simmental respectively) but poorly correlated with F_GRM (0.134, 0.128 and 0.448 for Italian Brown, Italian Holstein and Italian Simmental respectively). The inclusion of ROH > 8 Mb in the inbreeding calculation improved the correlation of F_ROH with F_PED and F_GRM. ROH are a direct measure of autozygosity at the DNA level and can overcome approximations and errors resulting from incomplete pedigree data. In populations with high linkage disequilibrium (LD) and recent inbreeding (e.g. Italian Holstein and Italian Brown), a medium‐density marker panel, such as the one used here, may provide a good estimate of inbreeding. However, in populations with low LD and ancient inbreeding, marker density would have to be increased to identify short ROH that are identical by descent more precisely.     

Morphological convergence in a Mexican garter snake associated with the ingestion of a novel prey

Morphological convergence is expected when organisms which differ in phenotype experience similar functional demands, which lead to similar associations between resource utilization and performance. To consume prey with hard exoskeletons, snakes require either specialized head morphology, or to deal with them when they are vulnerable, for example, during molting. Such attributes may in turn reduce the efficiency with which they prey on soft‐bodied, slippery animals such as fish. Snakes which consume a range of prey may present intermediate morphology, such as that of Thamnophiine (Natricinae), which may be classified morphometrically across the soft–hard prey dietary boundary. In this study, we compared the dentition and head structure of populations of Thamnophis melanogaster that have entered the arthropod–crustacean (crayfish)‐eating niche and those that have not, and tested for convergence between the former and two distantly related crayfish specialists of the genus Regina (R. septemvittata and R. grahamii). As a control, we included the congener T. eques. Multivariate analysis of jaw length, head length, head width, and number of maxillary teeth yielded three significant canonical variables that together explained 98.8% of the variance in the size‐corrected morphological data. The first canonical variable significantly discriminated between the three species. The results show that head dimensions and number of teeth of the two Regina species are more similar to those of crayfish‐eating T. melanogaster than to non‐crayfish‐eating snakes or of T. eques. It is unclear how particular head proportions or teeth number facilitates capture of crayfish, but our results and the rarity of soft crayfish ingestion by T. melanogaster may reflect the novelty of this niche expansion, and are consistent with the hypothesis that some populations of T. melanogaster have converged in their head morphology with the two soft crayfish‐eating Regina species, although we cannot rule out the possibility of a morphological pre‐adaptation to ingest crayfish.     

New resources for genetic studies in Populus nigra: genome‐wide SNP discovery and development of a 12k Infinium array

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water‐use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead‐Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5–7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural‐population based genetic association studies in P. nigra.     

Multilocus nuclear DNA markers reveal population structure and demography of Anopheles minimus

Utilization of multiple putatively neutral DNA markers for inferring evolutionary history of species population is considered to be the most robust approach. Molecular population genetic studies have been conducted in many species of Anopheles genus, but studies based on single nucleotide polymorphism (SNP) data are still very scarce. Anopheles minimus is one of the principal malaria vectors of Southeast (SE) Asia including the Northeastern (NE) India. Although population genetic studies with mitochondrial genetic variation data have been utilized to infer phylogeography of the SE Asian populations of this species, limited information on the population structure and demography of Indian An. minimus is available. We herewith have developed multilocus nuclear genetic approach with SNP markers located in X chromosome of An. minimus in eight Indian and two SE Asian population samples (121 individual mosquitoes in total) to infer population history and test several hypotheses on the phylogeography of this species. While the Thai population sample of An. minimus presented the highest nucleotide diversity, majority of the Indian samples were also fairly diverse. In general, An. minimus populations were moderately substructured in the distribution range covering SE Asia and NE India, largely falling under three distinct genetic clusters. Moreover, demographic expansion events could be detected in the majority of the presently studied populations of An. minimus. Additional DNA sequencing of the mitochondrial COII region in a subset of the samples (40 individual mosquitoes) corroborated the existing hypothesis of Indian An. minimus falling under the earlier reported mitochondrial lineage B.     

Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single‐nucleotide polymorphism array

With the access to draft genome sequence assemblies and whole‐genome resequencing data from population samples, molecular ecology studies will be able to take truly genome‐wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single‐nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10‐fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F₁ hybrids but no later‐generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as <1 Ma, this suggests strong selection against F₁ hybrids and unusually rapid evolution of reproductive incompatibility in an avian system.     

Ingestible insecticides for spotted wing Drosophila control: a polyol,Erythritol, and an insect growth regulator,Lufenuron

Bioassays tested insecticidal activity of Erythritol from the nutritive sweetener, Truvia, and an insect growth regulator, Lufenuron, against life stages (eggs, larvae, pupae, adults) of Drosophila melanogaster (Meigen) and Drosophila suzukii (Matsumura), the spotted wing Drosophila (SWD). These compounds were chosen for their demonstrated acute toxicity to adult and larval Drosophila and potential use on organic fruit farms. D. melanogaster fed on standard Drosophila diet media moistened with water containing known concentrations of Erythritol. Likewise, SWD consumed standard diet media as well as thawed host fruit (blackberries and blueberries) treated with solutions of Erythritol, Lufenuron or both. During the first bioassay, Erythritol at lower concentrations between 0 and 500 mm (~61 000 ppm) in water and mixed with instant diet media increased adult survival from ~80% to 97% for D. melanogaster and SWD. However, from aqueous concentrations ranging from 1750 (~414 000 ppm) to 2000 mm (~244 000 ppm), Erythritol killed 100% of adult Drosophila in culture vials. One hundred per cent mortality for SWD and D. melanogaster occurred at ≥0.5 m (~61 000 ppm) Erythritol added to diet media or topically applied to host fruit. In a second bioassay, 0.013–1.000 ppm of aqueous Lufenuron, a chitin synthase inhibitor, when added to dry diet media prevented 90–99% of SWD from reaching the pupal stage. In another assay, ~67% of SWD eggs or neonates (early first instars) died inside blackberries pre‐treated with (dipped in) a soapy solution of 10 ppm Lufenuron. Pre‐treating blackberry fruit with an Erythritol–Lufenuron mixture reduced SWD brood survival by 99%. Likewise, during our last fruit‐based bioassay, 98% of eggs and neonates died inside blueberries similarly pre‐treated. During the last experiment, Lufenuron in diet media also rendered adult females sterile. Sterility, however, dissipated over 7 days once females began feeding on a Lufenuron‐free diet media.     

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.     

Patterns of transposable element variation and clinality in Drosophila

Natural populations often exist in spatially diverse environments and may experience variation in the strength and targets of natural selection across their ranges. Drosophila provides an excellent opportunity to study the effects of spatially varying selection in natural populations, as both Drosophila melanogaster and Drosophila simulans live across a wide range of environments in North America. Here, we characterize patterns of variation in transposable elements (TEs) from six populations of D. melanogaster and nine populations of D. simulans sampled from multiple latitudes across North America. We find a nearly twofold excess of TEs in D. melanogaster relative to D. simulans, with this difference largely driven by TEs segregating at the lowest and highest allele frequencies. We find no effect of latitude on either total TE abundance or average TE allele frequencies in either species. Moreover, we show that, as a class of mutations, the most common patterns of TE variation do not coincide with the sampled latitudinal gradient, nor are they consistent with local adaptation acting on environmental differences found in the most extreme latitudes. We also do not find a cline in ancestry for North American D. melanogaster—for either TEs or single nucleotide polymorphisms—suggesting a limited role for demography in shaping patterns of TE variation. Though we find little evidence for widespread clinality among TEs in Drosophila, this does not necessarily imply a limited role for TEs in adaptation. We discuss the need for improved models of adaptation to large‐scale environmental heterogeneity, and how these might be applied to TEs.