Samples and sequencing. To investigate possible immune-related genes under selection before
European contact, we sequenced the exomes of ancient and modern First Nations individuals of
the Prince Rupert Harbour (PRH) region of British Columbia, Canada (Supplementary Fig. 1).
We then performed genomic scans for positive selection and functional characterization of genes
exhibiting the strongest signals.
Exomes of 25 modern individuals from two Coast Tsimshian
communities, Metlakatla and Lax Kw’alaams (henceforth referred to as “Tsimshian”), were
sequenced to a mean depth of 9.66x. The exomes of 25 ancient individuals from archaeological
sites in the PRH region (henceforth referred to as “PRH Ancients”; Supplementary Fig. 1) were
sequenced to a mean depth of 7.97x. Contamination estimates using the exome-wide data
revealed a mean contamination of 0.94% with a 95% CI: 0.83-1.10 16 (Supplementary Table 3).
All 25 ancient individuals exhibited patterns of C→T and G→A transitions consistent with
deamination due to post-mortem DNA damage17,18 (Supplementary Fig. 3). Mitochondrial
haplogroups were determined for each ancient individual, all showing haplotypes previously
identified in Native Americans19 (Supplementary Table 1).