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
-
accuracy.summary.dat.tar.gz
-
accuracy.txt
-
all_mapping_environmental_data.txt
-
all_reconstructed_haplotypes_melted.txt
-
all_simulated_reconstructed_haplotypes_melted.txt
-
avg_qqplot_data_thinned.Rdat
-
bergland_2014_cline_clean.txt
-
bergland_2014_seasonal_clean.txt
-
bergland2014_sign_universal_threshold_dropmissing_bottom_quantiles.Rdata
-
bergland2014_sign_universal_threshold_dropmissing_droptlk.Rdata
-
bergland2014_sign_universal_threshold_dropmissing.Rdata
-
bergland2019_sign_universal_threshold_dropmissing.Rdata
-
BME_LN_LNPA.ihs.txt
-
CM_dmel_dsim_2017.2018.csv
-
core20delta.rdat
-
east_coast_cline_V2_allele_freqs.txt
-
east_coast_cline_V2_clean.txt
-
final2.vcf.gz
-
fst_ZI_AT_gr_12_fall.txt
-
gcta_hsq_1000perms_st7.txt
-
gcta_hsq_1000perms.txt
-
GIFbychr.Rdata
-
haplotypesA.vcf.gz
-
haplotypesB.vcf.gz
-
hs.hc.snp.99.9.both.linenames.ann.vcf
-
hs.hc.snp.99.9.both.linenames.ann.vcf.seq.gds
-
hwe_missing_maf_filters.txt
-
hybrid_parent_PCA_dropmissing.Rdata
-
hybrid_parent_PCA.Rdata
-
ihs_dgrp_universal_threshold_dropmissing_bottomquantiles.Rdat
-
ihs_dgrp_universal_threshold_dropmissing_droptlk.Rdat
-
ihs_dgrp_universal_threshold_dropmissing.Rdat
-
ihs_north_universal_threshold_dropmissing_bottomquantiles.Rdat
-
ihs_north_universal_threshold_dropmissing_droptlk.Rdat
-
ihs_north_universal_threshold_dropmissing.Rdat
-
ihs.txt
-
interchromosomal_ld_dropmissing.txt
-
intrachromosomal_ld_dropmissing.txt
-
inv_fixed_coord.txt
-
ld_decay_dgrp_mafover0.05_downsampled_revised34.txt
-
ld_decay_dgrp_mafover0.05_revised.txt
-
ld_decay_mafover0.05_parents_revised_popA.txt
-
ld_decay_mafover0.05_parents_revised_popB.txt
-
ld_decay_mafover0.05_parents_revised_popboth.txt
-
library_info.txt
-
mesocosm_F0.csv
-
mesocosm_F2.csv
-
PAE_AOB_library_metadata_all.txt
-
parents_northsouth_PCA.txt
-
parents_PCA.txt
-
parents_season_PCA.txt
-
phenos_wolbachia_missinggeno.txt
-
popinfo_toplot.csv
-
pro_diapause_alleles_africa_clump200kb.txt
-
recombination.bed
-
reconstructA.vcf.gz
-
seas_glm_NOTswitch_clean.txt
-
simulated_actual_haplotypes.txt
-
snp_annotations.Rdat
-
strain_geography.csv
-
tim_genotypes.csv
-
universal_input.txt
-
weather_data.csv
-
yeast_supplementation.csv
-
ZI_admix_universal_threshold.Rdata
Oct 26, 2020 version files 6.34 GB
-
6d_data.Rdata
-
accuracy.summary.dat.tar.gz
-
accuracy.txt
-
all_mapping_environmental_data.txt
-
all_reconstructed_haplotypes_melted.txt
-
all_simulated_reconstructed_haplotypes_melted.txt
-
avg_qqplot_data_thinned.Rdat
-
bergland_2014_cline_clean.txt
-
bergland_2014_seasonal_clean.txt
-
bergland2014_sign_universal_threshold_dropmissing_bottom_quantiles.Rdata
-
bergland2014_sign_universal_threshold_dropmissing_droptlk.Rdata
-
bergland2014_sign_universal_threshold_dropmissing.Rdata
-
bergland2019_sign_universal_threshold_dropmissing.Rdata
-
BME_LN_LNPA.ihs.txt
-
CM_dmel_dsim_2017.2018.csv
-
core20delta.rdat
-
east_coast_cline_V2_allele_freqs.txt
-
east_coast_cline_V2_clean.txt
-
final2.vcf.gz
-
fst_ZI_AT_gr_12_fall.txt
-
gcta_hsq_1000perms_st7.txt
-
gcta_hsq_1000perms.txt
-
GIFbychr.Rdata
-
haplotypesA.vcf.gz
-
haplotypesB.vcf.gz
-
hs.hc.snp.99.9.both.linenames.ann.vcf
-
hs.hc.snp.99.9.both.linenames.ann.vcf.seq.gds
-
hwe_missing_maf_filters.txt
-
hybrid_parent_PCA_dropmissing.Rdata
-
hybrid_parent_PCA.Rdata
-
ihs_dgrp_universal_threshold_dropmissing_bottomquantiles.Rdat
-
ihs_dgrp_universal_threshold_dropmissing_droptlk.Rdat
-
ihs_dgrp_universal_threshold_dropmissing.Rdat
-
ihs_north_universal_threshold_dropmissing_bottomquantiles.Rdat
-
ihs_north_universal_threshold_dropmissing_droptlk.Rdat
-
ihs_north_universal_threshold_dropmissing.Rdat
-
ihs.txt
-
interchromosomal_ld_dropmissing.txt
-
intrachromosomal_ld_dropmissing.txt
-
inv_fixed_coord.txt
-
ld_decay_dgrp_mafover0.05_downsampled_revised34.txt
-
ld_decay_dgrp_mafover0.05_revised.txt
-
ld_decay_mafover0.05_parents_revised_popA.txt
-
ld_decay_mafover0.05_parents_revised_popB.txt
-
ld_decay_mafover0.05_parents_revised_popboth.txt
-
library_info.txt
-
mesocosm_F0.csv
-
mesocosm_F2.csv
-
PAE_AOB_library_metadata_all.txt
-
parents_northsouth_PCA.txt
-
parents_PCA.txt
-
parents_season_PCA.txt
-
phenos_wolbachia_missinggeno.txt
-
popinfo_toplot-2.csv
-
popinfo_toplot.csv
-
pro_diapause_alleles_africa_clump200kb.txt
-
recombination.bed
-
reconstructA.vcf.gz
-
seas_glm_NOTswitch_clean.txt
-
simulated_actual_haplotypes.txt
-
snp_annotations.Rdat
-
strain_geography.csv
-
tim_genotypes.csv
-
universal_input.txt
-
weather_data.csv
-
yeast_supplementation.csv
-
ZI_admix_universal_threshold.Rdata
Oct 15, 2021 version files 7.62 GB
-
6d_data.Rdata
-
accuracy.summary.dat.tar.gz
-
accuracy.txt
-
admix_universal_lasso_dropmissing.Rdata
-
all_GIF.Rdata
-
all_mapping_environmental_data.txt
-
all_reconstructed_haplotypes_melted.txt
-
all_simulated_reconstructed_haplotypes_melted.txt
-
avg_qqplot_data_thinned.Rdat
-
bergland_2014_cline_clean.txt
-
bergland_2014_seasonal_clean.txt
-
bergland2014_sign_universal_threshold_dropmissing_bottom_quantiles.Rdata
-
bergland2014_sign_universal_threshold_dropmissing_droptlk.Rdata
-
bergland2014_sign_universal_threshold_dropmissing.Rdata
-
bergland2019_sign_universal_threshold_dropmissing.Rdata
-
BME_LN_LNPA.ihs.txt
-
CM_dmel_dsim_2017.2018.csv
-
code_by_figure_R2.R
-
core20delta.rdat
-
east_coast_cline_V2_allele_freqs.txt
-
east_coast_cline_V2_clean.txt
-
final2.vcf.gz
-
fst_ZI_AT_gr_12_fall.txt
-
gcta_hsq_1000perms_st7.txt
-
gcta_hsq_1000perms.txt
-
GIFbychr.Rdata
-
gwas_top1percent_dropmissing.Rdat
-
haplotypesA.vcf.gz
-
haplotypesB.vcf.gz
-
hs.hc.snp.99.9.both.linenames.ann.vcf
-
hs.hc.snp.99.9.both.linenames.ann.vcf.seq.gds
-
hwe_missing_maf_filters.txt
-
hybrid_parent_PCA_dropmissing.Rdata
-
hybrid_parent_PCA.Rdata
-
ihs_dgrp_universal_threshold_dropmissing_bottomquantiles.Rdat
-
ihs_dgrp_universal_threshold_dropmissing_droptlk.Rdat
-
ihs_dgrp_universal_threshold_dropmissing.Rdat
-
ihs_north_universal_threshold_dropmissing_bottomquantiles.Rdat
-
ihs_north_universal_threshold_dropmissing_droptlk.Rdat
-
ihs_north_universal_threshold_dropmissing.Rdat
-
ihs.txt
-
interchromosomal_ld_dropmissing.txt
-
intrachromosomal_ld_dropmissing.txt
-
inv_fixed_coord.txt
-
lasso_snps_dropmissing.Rdata
-
lassoPred.dropmissing.Rdata
-
ld_decay_dgrp_mafover0.05_downsampled_revised34.txt
-
ld_decay_dgrp_mafover0.05_revised.txt
-
ld_decay_maf0.05_popA.txt
-
ld_decay_maf0.05_popB.txt
-
ld_decay_mafover0.05_parents_revised_popA.txt
-
ld_decay_mafover0.05_parents_revised_popB.txt
-
ld_decay_mafover0.05_parents_revised_popboth.txt
-
library_info.txt
-
mesocosm_F0.csv
-
mesocosm_F2.csv
-
PAE_AOB_library_metadata_all.txt
-
parents_hybrids_karyotype_calls.Rdata
-
parents_northsouth_PCA.txt
-
parents_PCA.txt
-
parents_season_PCA.txt
-
phenos_wolbachia_missinggeno.txt
-
popinfo_toplot-2.csv
-
popinfo_toplot.csv
-
pro_diapause_alleles_africa_clump200kb.txt
-
recombination.bed
-
reconstructA.vcf.gz
-
seas_glm_NOTswitch_clean.txt
-
simulated_actual_haplotypes.txt
-
single_population_sign_universal_threshold_dropmissing.Rdata
-
snp_annotations.Rdat
-
strain_geography.csv
-
tim_genotypes.csv
-
universal_input.txt
-
weather_data.csv
-
yeast_supplementation.csv
-
zambia_universal_threshold_dropmissing.Rdata
-
ZI_admix_universal_threshold.Rdata
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