| Popis: |
Additional file 1: FigureS1. A Results of GEMMA analysis utilizing workers from colonies fromall regions, visualized via Manhattan plot. We used a linear mixed model withcolony social form as the independent variable. Each point represents anindividual SNP, with the corresponding chromosome on the x-axis and thenegative logarithm of the SNP p-value on the y-axis. Only one SNP, from chromosome 3, exceeds the significance level (Bonferroni corrected significance threshold: 1.52E−04;p-value: 7.26E−06). B Results of GEMMA analysis onworkers from colonies in Alberta only, to reduce the effect of underlyingpopulation structure on the GWAS, visualized via Manhattan plot. We used alinear mixed model with colony social form as the independent variable. Eachpoint represents an individual SNP, with the corresponding chromosome on thex-axis and the negative logarithm of the SNP p-value on the y-axis. Five SNPs,all from chromosome 3, are above the significance threshold (Bonferroni corrected significance threshold: 1.89E−04; p-values: 7.79E−07, three at 4.43E−07, and 3.79E−07). Figure S2. Stacked bar plot displaying genotypes of samples from colonies labeled as ambiguous in social form (three from Alberta, one from southern British Columbia). Each bar represents all samples from an individual colony. Genotype in relation to each individual worker, is indicated by color: green = Sm/Sm andbrown = Sm/Sp. Figure S3. Line graph showing the effects of missingness on the Huang estimator. Individual (red) and whole sample (black) relatedness values are represented as lines. Figure S4. Density plot showing pairwise relatedness (determined by PolyRelatedness) among all sequencedworkers of colonies from the Evan Thomas trailhead in Alberta, Canada. On average, workers from monogyne colonies are the most related, whereas intermediate colonies (“ambiguous” in Fig. 1) have a bimodal distribution offull siblings and individuals with intermediate relatedness values, andpolygyne colonies have the lowest overall pairwise relatedness. Figure S5. A principal component analysis utilizing all markers except those on chromosome 3 in individuals fromcolony samples, with the year of the respective batches color coded. Table S1. Overview of consensus social form variables. Table S2. Sampling regions and their associated coordinate values. Table S3. PC axis weightings of variants on chromosome 3 (see also Fig. 2). Weightings greater than 1 or less than − 1 are highlighted in grey for each PC axis. The outliers from GWAS analyses are highlighted in green (Alberta only) and orange (GWAS of the full dataset). |
