Unique genetic signatures of local adaptation over space and time for diapause, an ecologically relevant complex trait, in Drosophila melanogaster
Data files
Sep 15, 2020 version files 6.34 GB
-
6d_data.Rdata
188.09 MB
-
accuracy.summary.dat.tar.gz
276.63 KB
-
accuracy.txt
1.01 MB
-
all_mapping_environmental_data.txt
319.46 MB
-
all_reconstructed_haplotypes_melted.txt
5.33 MB
-
all_simulated_reconstructed_haplotypes_melted.txt
4.91 MB
-
avg_qqplot_data_thinned.Rdat
20.34 MB
-
bergland_2014_cline_clean.txt
55.51 MB
-
bergland_2014_seasonal_clean.txt
53.86 MB
-
bergland2014_sign_universal_threshold_dropmissing_bottom_quantiles.Rdata
801.53 KB
-
bergland2014_sign_universal_threshold_dropmissing_droptlk.Rdata
803.32 KB
-
bergland2014_sign_universal_threshold_dropmissing.Rdata
808.61 KB
-
bergland2019_sign_universal_threshold_dropmissing.Rdata
828.60 KB
-
BME_LN_LNPA.ihs.txt
163.23 MB
-
CM_dmel_dsim_2017.2018.csv
14.81 KB
-
core20delta.rdat
690.60 MB
-
east_coast_cline_V2_allele_freqs.txt
79.94 MB
-
east_coast_cline_V2_clean.txt
131.11 MB
-
final2.vcf.gz
2.55 GB
-
fst_ZI_AT_gr_12_fall.txt
168.88 MB
-
gcta_hsq_1000perms_st7.txt
211.18 KB
-
gcta_hsq_1000perms.txt
211.41 KB
-
GIFbychr.Rdata
240.93 KB
-
haplotypesA.vcf.gz
10.58 MB
-
haplotypesB.vcf.gz
8.72 MB
-
hs.hc.snp.99.9.both.linenames.ann.vcf
10.26 MB
-
hs.hc.snp.99.9.both.linenames.ann.vcf.seq.gds
116.97 MB
-
hwe_missing_maf_filters.txt
847.34 MB
-
hybrid_parent_PCA_dropmissing.Rdata
1.31 MB
-
hybrid_parent_PCA.Rdata
1.35 MB
-
ihs_dgrp_universal_threshold_dropmissing_bottomquantiles.Rdat
215.29 KB
-
ihs_dgrp_universal_threshold_dropmissing_droptlk.Rdat
444.09 KB
-
ihs_dgrp_universal_threshold_dropmissing.Rdat
445.95 KB
-
ihs_north_universal_threshold_dropmissing_bottomquantiles.Rdat
213.36 KB
-
ihs_north_universal_threshold_dropmissing_droptlk.Rdat
441.80 KB
-
ihs_north_universal_threshold_dropmissing.Rdat
443.72 KB
-
ihs.txt
381.35 MB
-
interchromosomal_ld_dropmissing.txt
5.80 MB
-
intrachromosomal_ld_dropmissing.txt
3.83 MB
-
inv_fixed_coord.txt
9.23 KB
-
ld_decay_dgrp_mafover0.05_downsampled_revised34.txt
7.51 MB
-
ld_decay_dgrp_mafover0.05_revised.txt
7.71 MB
-
ld_decay_mafover0.05_parents_revised_popA.txt
7.92 MB
-
ld_decay_mafover0.05_parents_revised_popB.txt
7.92 MB
-
ld_decay_mafover0.05_parents_revised_popboth.txt
8.24 MB
-
library_info.txt
283.26 KB
-
mesocosm_F0.csv
20.56 KB
-
mesocosm_F2.csv
183.02 KB
-
PAE_AOB_library_metadata_all.txt
72.72 KB
-
parents_northsouth_PCA.txt
96.21 KB
-
parents_PCA.txt
112.63 KB
-
parents_season_PCA.txt
30.30 KB
-
phenos_wolbachia_missinggeno.txt
1.74 MB
-
popinfo_toplot.csv
588 B
-
pro_diapause_alleles_africa_clump200kb.txt
300.10 MB
-
recombination.bed
38.12 KB
-
reconstructA.vcf.gz
64.38 MB
-
seas_glm_NOTswitch_clean.txt
85.58 MB
-
simulated_actual_haplotypes.txt
1.93 MB
-
snp_annotations.Rdat
16.36 MB
-
strain_geography.csv
3.28 KB
-
tim_genotypes.csv
181.10 KB
-
universal_input.txt
229.33 KB
-
weather_data.csv
5.12 KB
-
yeast_supplementation.csv
347.26 KB
-
ZI_admix_universal_threshold.Rdata
116.64 KB
Oct 26, 2020 version files 6.34 GB
-
6d_data.Rdata
188.09 MB
-
accuracy.summary.dat.tar.gz
276.63 KB
-
accuracy.txt
1.01 MB
-
all_mapping_environmental_data.txt
319.46 MB
-
all_reconstructed_haplotypes_melted.txt
5.33 MB
-
all_simulated_reconstructed_haplotypes_melted.txt
4.91 MB
-
avg_qqplot_data_thinned.Rdat
20.34 MB
-
bergland_2014_cline_clean.txt
55.51 MB
-
bergland_2014_seasonal_clean.txt
53.86 MB
-
bergland2014_sign_universal_threshold_dropmissing_bottom_quantiles.Rdata
801.53 KB
-
bergland2014_sign_universal_threshold_dropmissing_droptlk.Rdata
803.32 KB
-
bergland2014_sign_universal_threshold_dropmissing.Rdata
808.61 KB
-
bergland2019_sign_universal_threshold_dropmissing.Rdata
828.60 KB
-
BME_LN_LNPA.ihs.txt
163.23 MB
-
CM_dmel_dsim_2017.2018.csv
14.81 KB
-
core20delta.rdat
690.60 MB
-
east_coast_cline_V2_allele_freqs.txt
79.94 MB
-
east_coast_cline_V2_clean.txt
131.11 MB
-
final2.vcf.gz
2.55 GB
-
fst_ZI_AT_gr_12_fall.txt
168.88 MB
-
gcta_hsq_1000perms_st7.txt
211.18 KB
-
gcta_hsq_1000perms.txt
211.41 KB
-
GIFbychr.Rdata
240.93 KB
-
haplotypesA.vcf.gz
10.58 MB
-
haplotypesB.vcf.gz
8.72 MB
-
hs.hc.snp.99.9.both.linenames.ann.vcf
10.26 MB
-
hs.hc.snp.99.9.both.linenames.ann.vcf.seq.gds
116.97 MB
-
hwe_missing_maf_filters.txt
847.34 MB
-
hybrid_parent_PCA_dropmissing.Rdata
1.31 MB
-
hybrid_parent_PCA.Rdata
1.35 MB
-
ihs_dgrp_universal_threshold_dropmissing_bottomquantiles.Rdat
215.29 KB
-
ihs_dgrp_universal_threshold_dropmissing_droptlk.Rdat
444.09 KB
-
ihs_dgrp_universal_threshold_dropmissing.Rdat
445.95 KB
-
ihs_north_universal_threshold_dropmissing_bottomquantiles.Rdat
213.36 KB
-
ihs_north_universal_threshold_dropmissing_droptlk.Rdat
441.80 KB
-
ihs_north_universal_threshold_dropmissing.Rdat
443.72 KB
-
ihs.txt
381.35 MB
-
interchromosomal_ld_dropmissing.txt
5.80 MB
-
intrachromosomal_ld_dropmissing.txt
3.83 MB
-
inv_fixed_coord.txt
9.23 KB
-
ld_decay_dgrp_mafover0.05_downsampled_revised34.txt
7.51 MB
-
ld_decay_dgrp_mafover0.05_revised.txt
7.71 MB
-
ld_decay_mafover0.05_parents_revised_popA.txt
7.92 MB
-
ld_decay_mafover0.05_parents_revised_popB.txt
7.92 MB
-
ld_decay_mafover0.05_parents_revised_popboth.txt
8.24 MB
-
library_info.txt
283.26 KB
-
mesocosm_F0.csv
20.56 KB
-
mesocosm_F2.csv
183.02 KB
-
PAE_AOB_library_metadata_all.txt
72.72 KB
-
parents_northsouth_PCA.txt
96.21 KB
-
parents_PCA.txt
112.63 KB
-
parents_season_PCA.txt
30.30 KB
-
phenos_wolbachia_missinggeno.txt
1.74 MB
-
popinfo_toplot-2.csv
590 B
-
popinfo_toplot.csv
588 B
-
pro_diapause_alleles_africa_clump200kb.txt
300.10 MB
-
recombination.bed
38.12 KB
-
reconstructA.vcf.gz
64.38 MB
-
seas_glm_NOTswitch_clean.txt
85.58 MB
-
simulated_actual_haplotypes.txt
1.93 MB
-
snp_annotations.Rdat
16.36 MB
-
strain_geography.csv
3.28 KB
-
tim_genotypes.csv
181.10 KB
-
universal_input.txt
229.33 KB
-
weather_data.csv
5.12 KB
-
yeast_supplementation.csv
347.26 KB
-
ZI_admix_universal_threshold.Rdata
116.64 KB
Oct 15, 2021 version files 7.62 GB
-
6d_data.Rdata
188.09 MB
-
accuracy.summary.dat.tar.gz
276.63 KB
-
accuracy.txt
1.01 MB
-
admix_universal_lasso_dropmissing.Rdata
143 B
-
all_GIF.Rdata
96.72 KB
-
all_mapping_environmental_data.txt
319.46 MB
-
all_reconstructed_haplotypes_melted.txt
5.33 MB
-
all_simulated_reconstructed_haplotypes_melted.txt
4.91 MB
-
avg_qqplot_data_thinned.Rdat
20.34 MB
-
bergland_2014_cline_clean.txt
55.51 MB
-
bergland_2014_seasonal_clean.txt
53.86 MB
-
bergland2014_sign_universal_threshold_dropmissing_bottom_quantiles.Rdata
801.53 KB
-
bergland2014_sign_universal_threshold_dropmissing_droptlk.Rdata
803.32 KB
-
bergland2014_sign_universal_threshold_dropmissing.Rdata
808.61 KB
-
bergland2019_sign_universal_threshold_dropmissing.Rdata
828.60 KB
-
BME_LN_LNPA.ihs.txt
163.23 MB
-
CM_dmel_dsim_2017.2018.csv
14.81 KB
-
code_by_figure_R2.R
130.17 KB
-
core20delta.rdat
690.60 MB
-
east_coast_cline_V2_allele_freqs.txt
79.94 MB
-
east_coast_cline_V2_clean.txt
131.11 MB
-
final2.vcf.gz
2.55 GB
-
fst_ZI_AT_gr_12_fall.txt
168.88 MB
-
gcta_hsq_1000perms_st7.txt
211.18 KB
-
gcta_hsq_1000perms.txt
211.41 KB
-
GIFbychr.Rdata
240.93 KB
-
gwas_top1percent_dropmissing.Rdat
1.06 GB
-
haplotypesA.vcf.gz
10.58 MB
-
haplotypesB.vcf.gz
8.72 MB
-
hs.hc.snp.99.9.both.linenames.ann.vcf
10.26 MB
-
hs.hc.snp.99.9.both.linenames.ann.vcf.seq.gds
116.97 MB
-
hwe_missing_maf_filters.txt
847.34 MB
-
hybrid_parent_PCA_dropmissing.Rdata
1.31 MB
-
hybrid_parent_PCA.Rdata
1.35 MB
-
ihs_dgrp_universal_threshold_dropmissing_bottomquantiles.Rdat
215.29 KB
-
ihs_dgrp_universal_threshold_dropmissing_droptlk.Rdat
444.09 KB
-
ihs_dgrp_universal_threshold_dropmissing.Rdat
445.95 KB
-
ihs_north_universal_threshold_dropmissing_bottomquantiles.Rdat
213.36 KB
-
ihs_north_universal_threshold_dropmissing_droptlk.Rdat
441.80 KB
-
ihs_north_universal_threshold_dropmissing.Rdat
443.72 KB
-
ihs.txt
381.35 MB
-
interchromosomal_ld_dropmissing.txt
5.80 MB
-
intrachromosomal_ld_dropmissing.txt
3.83 MB
-
inv_fixed_coord.txt
9.23 KB
-
lasso_snps_dropmissing.Rdata
23.90 MB
-
lassoPred.dropmissing.Rdata
181.05 MB
-
ld_decay_dgrp_mafover0.05_downsampled_revised34.txt
7.51 MB
-
ld_decay_dgrp_mafover0.05_revised.txt
7.71 MB
-
ld_decay_maf0.05_popA.txt
7.41 MB
-
ld_decay_maf0.05_popB.txt
7.39 MB
-
ld_decay_mafover0.05_parents_revised_popA.txt
7.92 MB
-
ld_decay_mafover0.05_parents_revised_popB.txt
7.92 MB
-
ld_decay_mafover0.05_parents_revised_popboth.txt
8.24 MB
-
library_info.txt
283.26 KB
-
mesocosm_F0.csv
20.56 KB
-
mesocosm_F2.csv
183.02 KB
-
PAE_AOB_library_metadata_all.txt
72.72 KB
-
parents_hybrids_karyotype_calls.Rdata
10.87 KB
-
parents_northsouth_PCA.txt
96.21 KB
-
parents_PCA.txt
112.63 KB
-
parents_season_PCA.txt
30.30 KB
-
phenos_wolbachia_missinggeno.txt
1.74 MB
-
popinfo_toplot-2.csv
590 B
-
popinfo_toplot.csv
588 B
-
pro_diapause_alleles_africa_clump200kb.txt
300.10 MB
-
recombination.bed
38.12 KB
-
reconstructA.vcf.gz
64.38 MB
-
seas_glm_NOTswitch_clean.txt
85.58 MB
-
simulated_actual_haplotypes.txt
1.93 MB
-
single_population_sign_universal_threshold_dropmissing.Rdata
4.97 MB
-
snp_annotations.Rdat
16.36 MB
-
strain_geography.csv
3.28 KB
-
tim_genotypes.csv
181.10 KB
-
universal_input.txt
229.33 KB
-
weather_data.csv
5.12 KB
-
yeast_supplementation.csv
347.26 KB
-
zambia_universal_threshold_dropmissing.Rdata
496.92 KB
-
ZI_admix_universal_threshold.Rdata
116.64 KB
Abstract
Organisms living in seasonally variable environments utilize cues such as light and temperature to induce plastic responses, enabling them to exploit favorable seasons and avoid unfavorable ones. Local adapation can result in variation in seasonal responses, but the genetic basis and evolutionary history of this variation remains elusive. Many insects, including Drosophila melanogaster, are able to undergo an arrest of reproductive development (diapause) in response to unfavorable conditions. In D. melanogaster, the ability to diapause is more common in high latitude populations, where flies endure harsher winters, and in the spring, reflecting differential survivorship of overwintering populations. Using a novel hybrid swarm-based genome wide association study, we examined the genetic basis and evolutionary history of ovarian diapause. We exposed outbred females to different temperatures and day lengths, characterized ovarian development for over 2800 flies, and reconstructed their full phased genomes. We found that diapause scored at two different developmental cutoffs has modest heritability, and we identified hundreds of SNPs associated with each of the two phenotypes. Alleles associated with one of the diapause phenotypes tend to be more common at higher latitudes, but these alleles do not show predictable seasonal variation. The collective signal of many small-effect, clinally varying SNPs can plausibly explain latitudinal phenotypic variation seen in North America. SNPs associated with diapause do not exhibit signs of recent selective sweeps, but most are segregating at relatively high frequencies in Africa, suggesting that variation in diapause relies on ancestral polymorphisms. Finally, we utilized outdoor mesocosms to track diapause under natural conditions. We found that hybrid swarms reared outdoors evolved increased propensity for diapause in late fall, whereas indoor control populations experienced no such change. Our results indicate that diapause is a complex, quantitative trait with different evolutionary patterns across time and space.
Usage notes
All file names should match the file names found in the code for this analysis, located at: https://github.com/ericksonp/diapause-scripts-clean