DOI: 10.5061/dryad.4mt34
title: Data from: Phenotypic integration in the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus)

description: Selection can vary geographically across environments and temporally over
 the lifetime of an individual. Unlike geographic contexts, where different
 selective regimes can act on different alleles, age-specific selection is
 constrained to act on the same genome by altering age-specific expression.
 Snake venoms are exceptional traits for studying ontogeny because toxin
 expression variation directly changes the phenotype; relative amounts of
 venom components determine, in part, venom efficacy. Phenotypic
 integration is the dependent relationship between different traits that
 collectively produce a complex phenotype and, in venomous snakes, may
 include traits as diverse as venom, head shape and fang length. We
 examined the feeding system of the eastern diamondback rattlesnake
 (Crotalus adamanteus) across environments and over the lifetime of
 individuals and used a genotype–phenotype map approach, protein expression
 data and morphological data to demonstrate that: (i) ontogenetic effects
 explained more of the variation in toxin expression variation than
 geographic effects, (ii) both juveniles and adults varied geographically,
 (iii) toxin expression variation was a result of directional selection and
 (iv) different venom phenotypes covaried with morphological traits also
 associated with feeding in temporal (ontogenetic) and geographic
 (functional) contexts. These data are the first to demonstrate, to our
 knowledge, phenotypic integration between multiple morphological
 characters and a biochemical phenotype across populations and age classes.
 We identified copy number variation as the mechanism driving the
 difference in the venom phenotype associated with these morphological
 differences, and the parallel mitochondrial, venom and morphological
 divergence between northern and southern clades suggests that each clade
 may warrant classification as a separate evolutionarily significant unit.
description: Table_S1: Raw HPLC
 dataThis files contains the
 raw HPLC peak abundance data for all of our samples. The columns are as
 follows: Venom: unique venom sample identifier. Tissue: Unique animal
 identifier (one animal may have several venom samples; KW0719a is the same
 animal as KW0719b). Peak1a to Peak21: area under each HPLC peak (not a
 percentage). MaturityFactor: the age class each animal belongs to based on
 the 102 cm SVL cut-off. GeographyFactor: population of each animal. Eglin
 is Ca1, ANF is Ca2, ONF is Ca3, ENP is Ca4, LSG is Ca5, CAL is Ca6, SAP is
 Ca7. The last two columns indicate whether the sample was used in all
 analyses or only in the PCA and cluster analysis (i.e., removed prior to
 non-parametric MANOVA because it lacked size data and/or was a duplicate
 sample from an individual already present in the data
 set).Table_S1.xlsxTable_S2_raw_qPCR_dataThis file contains the raw qPCR data. For the “Results edited” and “Results” sheets: the sample name corresponds to the unique tissue identifier in Table S1. Each sample was run in duplicate at two different concentrations (e.g., 529 and 529-10). The CT is the cycle threshold, the CT mean is the cycle threshold mean, the CT SD is the cycle threshold standard deviation, and Tm corresponds to the melting temperature. The “pics” sheet shows the melting curve graph for northern individuals on the right and southern individuals on the left with the red line representing the standard (i.e., control). The “sheet1” sheet shows a draft of Figure 6 from the manuscript.Table_S2.xlsTable_S3_Raw_Morphological_DataThis table contains the raw morphological data. The “All Data” sheet contains the raw measurements for each individual. The columns are the following: “Museum ID” is the unique specimen identifier for the Florida Museum, SVL is snout-vent length in mm, HL is head length in mm, HW is head width in mm, IF is inter-fang distance in mm, and FL is fang length in mm. The “PCA_Kruskal_Wallis” sheet contains the data from the Kruskal-Wallis tests. The “Shapiro-Wilk” sheet contains the Shapiro-Wilk normality test data for each variable.Table_S3.xlsTable_S4_PCAThis file contains the loadings for the first 3 principal components broken down by RP-HPLC peak.Table_S4.pdfTable_S5_PCA_LoadingsThis file contains all PCA loadings.Table_S5.txtTable_S6This file contains the variance-covariance data used to identify the most variable RP-HPLC peaks in the entire data set, adult data set only, and juvenile data set only.Cada CytB AlignmentAlignment of cytochrome bCada ND5 AlignmentAlignment of ND5CadamanteusFinalTree.tre file for Crotalus adamanteus phylogeny