Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

RapidRat: Development, validation and application of a genotyping-by-sequencing panel for rapid biosecurity and invasive species management

  • Bryson M. F. Sjodin,

    Roles Data curation, Formal analysis, Writing – original draft

    Affiliation Department of Biology, University of British Columbia, Kelowna, British Columbia, Canada

  • Robyn L. Irvine,

    Roles Data curation, Writing – review & editing

    Affiliation Gwaii Haanas National Park Reserve, National Marine Conservation Area Reserve and Haida Heritage Site, Parks Canada, Skidegate, British Columbia, Canada

  • Michael A. Russello

    Roles Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing

    michael.russello@ubc.ca

    Affiliation Department of Biology, University of British Columbia, Kelowna, British Columbia, Canada

Abstract

Invasive alien species (IAS) are among the main causes of global biodiversity loss. Invasive brown (Rattus norvegicus) and black (R. rattus) rats, in particular, are leading drivers of extinction on islands, especially in the case of seabirds where >50% of all extinctions have been attributed to rat predation. Eradication is the primary form of invasive rat management, yet this strategy has resulted in a ~10–38% failure rate on islands globally. Genetic tools can help inform IAS management, but such applications to date have been largely reactive, time-consuming, and costly. Here, we developed a Genotyping-in-Thousands by sequencing (GT-seq) panel for rapid species identification and population assignment of invasive brown and black rats (RapidRat) in Haida Gwaii, an archipelago comprising ~150 islands off the central coast of British Columbia, Canada. We constructed an optimized panel of 443 single nucleotide polymorphisms (SNPs) using previously generated double-digest restriction-site associated DNA (ddRAD) genotypic data (27,686 SNPs) from brown (n = 295) and black rats (n = 241) sampled throughout Haida Gwaii. The informativeness of this panel for identifying individuals to species and island of origin was validated relative to the ddRAD results; in all comparisons, admixture coefficients and population assignments estimated using RapidRat were consistent. To demonstrate application, 20 individuals from novel invasions of three islands (Agglomerate, Hotspring, Ramsay) were genotyped using RapidRat, all of which were confidently assigned (>98.5% probability) to Faraday and Murchison Islands as putative source populations. These results indicated that a previous eradication on Hotspring Island was conducted at an inappropriate geographic scale; future management should expand the eradication unit to include neighboring islands to prevent re-invasion. Overall, we demonstrated that RapidRat is an effective tool for managing invasive rat populations in Haida Gwaii and provided a clear framework for GT-seq panel development for informing biodiversity conservation in other systems.

Introduction

Human modification of the environment has greatly accelerated the rate of global biodiversity loss in recent decades [1,2]. A leading driver of these losses is the introduction of non-native species to new environments. Once established, these invasive alien species (IAS) can alter ecosystem dynamics, especially through predation, competition, and disease transmission, resulting in rapid declines to endemic species [3,4]. IAS effects on biodiversity are profound; 58% of extinctions listed in the International Union for Conservation of Nature (IUCN) Red List name IAS as main contributors [5]. IAS also persist as economically important pests, costing Canada between $16.6 and $34.5 billion dollars annually [6], with approximately $20 billion of that in the Forestry sector alone [7]. Likewise, the United States spends >US$220 billion annually in prevention efforts, damages, and habitat restoration [8]. These negative effects are amplified within island systems, evidenced by IAS being the leading cause of the majority (86%) of recent island extinctions [3]. Moreover, insular systems support disproportionately high levels of endemism relative to mainland areas [9], punctuating the need for immediate action towards improving IAS management on islands.

Invasive mammals are the main cause of animal extinctions on islands and constitute one of the most significant threats to insular biodiversity [10]. Brown rats (Rattus norvegicus), black rats (R. rattus), and Pacific rats (R. exulans), in particular, make up a triad of the most invasive mammals on the planet, found on every inhabited continent and >80% of all oceanic islands [11]. To date, invasive rats have caused >50 species extinctions and are a leading cause of extinction for insular species [1214]. Island seabirds, in particular, are adversely affected by rat invasions, as these species are poorly adapted to terrestrial predators due to their natural absence in many island systems [11,15,16]; in fact, approximately 60% of all seabird extinctions are linked to predation by invasive rats, 90% of these occurring on islands [14]. As part of seabird recovery efforts, management agencies often use whole-island eradications to manage invasive rat populations [14,17]. Yet, rodent eradication attempts on islands throughout the world have resulted in a ~10–38% failure rate [14,18]. In many cases, failure was due to knowledge gaps of key population parameters, including IAS population size, home range size, and dispersal capacity [14,18].

Genetic and genomic tools provide promising opportunities to inform vertebrate IAS management on islands. For example, genetic evidence was a large factor in the success of the largest rodent eradication in history on South Georgia Island (~390,000 ha) in the British Overseas Territory of the south Atlantic [19,20], providing critical information on population structure and demographic history that allowed for effective planning of a multi-stage eradication. In other cases, genetic and genomic tools can be used to explicitly evaluate the underlying cause of eradication failure, summarized in two prevailing hypotheses [14,21,22]. The survivor hypothesis suggests that the eradication effort failed to remove all individuals, leading to these individuals repopulating the island. On the other hand, the re-invader hypothesis posits that the eradication was successful in the complete removal of rats, but that the island was re-colonized from some connected source population. Genetic evidence was used to test these hypotheses following an unsuccessful rat eradication on Pearl Island in New Zealand; genetic analyses of population structure and individual assignment among pre- and post- eradication samples revealed evidence for re-colonization from a neighboring source population on Stewart Island (i.e., re-invader hypothesis), suggesting that a larger target area may be required to avoid eradication failure in this system in the future [23]. Although useful for informing IAS management, genetic analyses to date have been largely reactive (e.g., evaluating causes of eradication failure), time-consuming, and costly. Consequently, there is a need to harness leading-edge genomic technologies for providing genetic tools to aid managers in the decision-making process before, during and after IAS management to help improve efficiency, eradication durability, and cost-effectiveness.

Here, we developed a Genotyping-in-Thousands by sequencing (GT-seq) panel for rapid species identification and population assignment of invasive brown and black rats (RapidRat) in Haida Gwaii, an archipelago off the central coast of British Columbia (BC), Canada. GT-seq is a targeted, multiplexed amplicon sequencing approach that can simultaneously genotype sample sizes ranging from a single individual up to thousands of individuals at hundreds of single nucleotide polymorphisms (SNPs) [24]. We validated the RapidRat panel by comparing estimated population structure and assignment accuracy with results from a previously published double digest restriction site-associated DNA sequencing (ddRAD) analysis (n = 27,686 SNPs) of the same samples within this system [25]. Lastly, we deployed RapidRat to identify the source of recent, novel brown rat invasions on Agglomerate, Hotspring, and Ramsay Islands.

Methods

Study system

Haida Gwaii (Xaayda Gwaay in Haida) is an isolated archipelago consisting of approximately 150 smaller and two larger islands located ~80 km off the central coast of British Columbia (BC), Canada. This system contains more unique endemic sub-species than any other equal sized area in Canada [26] and is host to large numbers of seabirds and globally critical breeding habitat [2732]. Both brown and black rats have invaded Haida Gwaii, likely as stowaways on European ships using the archipelago as a fishing and whaling port in the 1700’s to early 1900s [29,30,33,34]. Since their introduction, invasive rats have adversely affected native seabird populations, resulting in the extirpation of several large multi-species colonies and population declines in at least half of the species present [29,31]. To meet seabird recovery goals and to protect ecological and cultural integrity [35], Parks Canada and other agencies started managing invasive rats in the 1990s primarily through whole-island eradications. Since 1997, invasive rats have been successfully eradicated from five islands throughout the archipelago (Fig 1); however, brown rats were re-detected on the Bischof Islands following two eradication attempts in 2003 and 2011, and brown rats were detected on Faraday and Murchison Islands four years after a successful black rat eradication in 2013 [29,30,36,37]. Brown rats have now also spread to Hotspring, Ramsay, and Agglomerate Islands, all of which have been historically rat-free (C. Bergman, pers. comm.). The origins of these invasions are unknown, though both Agglomerate and Hotspring Islands are situated ~1 km from the nearest brown rat population on Murchison Island (Fig 1), while the newly established population on Hotspring Island could have invaded Ramsay Island (900-1000m away); these spans are within the hypothesized swimming distance for brown rats in temperate waters [25,38,39]. Additionally, rats could have potentially dispersed from island to island by floating on debris (i.e., “rafting”) during high tide [40], and/or by commensal spread on ships.

thumbnail
Fig 1. Sampled islands and other locations within Haida Gwaii, British Columbia.

Island names with superscripts indicate past invasive rat (Rattus spp.) eradications. Eradications were performed in: (a) 1997; (b) 1998; (c) 2003; (d) 2011; and (e) 2013. Shaded are indicates the Gwaii Haanas National Park Reserve, National Marine Conservation Area Reserve and Haida Heritage site. Samples from islands written in red were used for GT-seq panel validation. The map was generated using data licensed under the Open Government Licence–British Columbia (https://www2.gov.bc.ca/gov/content/data/open-data/open-government-licence-bc).

https://doi.org/10.1371/journal.pone.0234694.g001

We previously characterized population connectivity among brown and black rat populations in Haida Gwaii using ddRAD sequencing (n = 27,686 SNPs) across the archipelago [25]. We found discrete genetic units largely corresponding to island populations, and we characterized within-island dynamics on the larger islands of Graham, Lyell, and Kunghit. We also assigned individuals from recent invasions on Bischof, Faraday, and Murchison Islands as invaders from neighbouring Lyell Island, demonstrating the utility of genetic analysis in this system, and providing an excellent source of reference data for the current study [25].

Study design

To identify candidate loci for our GT-seq panel, we used ddRAD data previously generated from brown rats (n = 295) and black rats (n = 241) collected throughout the Haida Gwaii archipelago [25]. Of these samples, brown rat (n = 59) and black rat (n = 37) DNA extractions (used in [25]) were re-sequenced and genotyped using GT-seq to optimize panel loci (see below). In addition to these samples, we sequenced and genotyped brown rats using GT-seq from novel invasions on Agglomerate Island (n = 18), Ramsay Island (n = 1), and Hotspring Island (n = 1). See Table 1 for detailed sample distribution. All sample collection was performed under Parks Canada Agency Animal Care Committee protocol GHNPR11-5. Whole genomic DNA (gDNA) was extracted from 10–20 mg of dried ear tissue using the Qiagen DNeasy® Blood and Tissue Kit and treated with RNase A (5PRIME) following the manufacturer’s protocol and stored in distilled water until GT-seq library preparation.

thumbnail
Table 1. Sample size and distribution for brown (Rattus norvegicus) and black rats (R. rattus) collected across Haida Gwaii, BC.

Samples were either genotyped using double-digest restriction site-associated DNA sequencing (ddRAD) or genotyping-in-thousands by sequencing (GT-seq). All brown rats were genotyped at 315 SNPs, while all black rats were genotyped at 429 SNPs. Numbers in brackets indicate sample size less replicate samples. (*) These samples are referred to as “Kunghit” for the black rats. (†) These samples are referred to as “Lyell-SW” for the black rats. Haida island names are indicated in parentheses.

https://doi.org/10.1371/journal.pone.0234694.t001

SNP discovery, quality control, and panel selection

Raw sequencing reads were first demultiplexed from the ddRAD samples using process_radtags as implemented in Stacks v2.0b [41]. Reads were trimmed to 100 bp to remove low-quality base calls at the 3’ ends. Individual reads were aligned to the most-current brown rat reference genome (Rnor_6.0, GenBank assembly accession: GCA_000001895.4) using the default settings in Bowtie 2 v2.2.9 [42]. Single nucleotide polymorphisms were identified and genotyped using gstacks and populations in Stacks v2.0b. To ensure only high-quality SNPs were retained, all loci were filtered such that they were genotyped in >90% of total individuals, had a minimum minor allele frequency of 5%, and had a maximum observed heterozygosity of 50%. Additionally, only a single SNP per RADtag was retained to reduce the likelihood of linkage disequilibrium among loci.

To ensure maximum quality and functionality of the panel, we additionally filtered loci such that they were variable in at least one species, genotyped in >90% of individuals in each species separately, as well as with a mean depth of coverage of 6x in each species separately using VCFtools v0.1.15 [43]. Of these, only SNPs found within 40–60 bp of the start of the RADtag sequence were retained to ensure sufficient flanking sequence remained for primer design. We then estimated population differentiation for each locus both within and between species using Weir and Cockerham’s unbiased estimator θ [44]. Estimates were calculated for all pairwise observations as implemented in Genetix v4.05.2 [45]. To ensure linkage equilibrium, locus pairs were positioned at least 1 Mb apart within the genome; in the event of a conflict, the locus with the lowest θ value was removed. We retained the top 350 loci with the highest discriminatory power (i.e., highest θ) within each species as well as an additional eight loci for differentiating between species for a total putative panel of 708 SNPs. The associated RADtag sequences for these loci were sent to GTseek LLC (https://gtseek.com/), where custom PCR1 primers were designed for each locus incorporating Illumina primer binding sites.

Genotyping-in-thousands by sequencing and panel optimization

To test the efficacy of our GT-seq panel, we sequenced and genotyped 92 individuals across both species (brown rats, n = 59; black rats, n = 37), preferentially selecting individuals with the highest quality ddRAD data and incorporating the entire geographic sampling distribution (Table 1). Of these, four individuals were replicated to explicitly assess genotyping error rate (see below), resulting in a total 96 samples in our initial library. GT-seq library preparation followed the original protocol [24] with some modifications [46]. DNA extractions for all individuals were normalized to 20 ng/uL to ensure even amplification across samples and diluted the PCR1 products 1:10 before use in PCR2. We quantified PCR2 products using PicoGreen (Molecular Probes, Inc.) and manually normalized these products to a concentration of 10 ng/μL. We pooled 2.5 μL into a final sequencing library and purified this library using the MinElute PCR Purification kit (Qiagen®), eluting into a final volume of 50 μL of nuclease-free water. The library was sequenced using a single lane of Illumina MiSeq paired-end 100 bp sequencing at the McGill University and Génome Québec Innovation Centre.

We genotyped raw sequencing data using the GT-seq pipeline found on GitHub (https://github.com/GT-seq/GT-seq-Pipeline). We used the GT-seq_SeqTest.pl and GT-seq_Primer-Interaction-Test.pl scripts to identify over-represented primer sequences and PCR artifacts and primer-dimers between loci. We removed loci accounting for >1% of total forward primer reads, loci with total counts >1% of PCR artifacts, as well as one locus per hetero-dimer pair. The remaining loci were combined into a new primer pool, and the library was re-sequenced and optimized following the above protocol. Loci that did not genotype in the final library, as well as individuals with >50% missing data, were removed from downstream analysis.

Assessment of population structure and panel validation

To assess genotyping accuracy of the panel within each species, we compared individual genotypes from GT-seq to ddRAD at species-specific loci and measured the percent discordance across methods. We likewise calculated genotyping error rates within each sequencing method by comparing individual genotypes between replicate samples, also measured as percent discordance. Loci with missing data in either one or both of each pair of samples were not included in error calculations.

We examined the ability of our panel to differentiate between species across all panel loci using discriminant analysis of principle components (DAPC) as implemented in the R package adegenet v2.1.1 [47]. This procedure first applies principle component analysis (PCA) to identify genetic clusters and then uses discriminant analysis to maximize the variation among these clusters while minimizing within-cluster variation [48]. These results were compared to species-assignments from Sjodin et al. [25] to assess the accuracy of assignment. To ensure our subset of loci resulted in similar patterns of population structure, we ran projected PCA using the ddRAD samples to infer eigenvectors and compared to Sjodin et al. [25] for accuracy using the smartpca function in EIGENSOFT v7.2.1 [4951]. We then projected GT-seq samples onto these eigenvectors to compare clustering across methods. Additionally, we estimated admixture coefficients for ddRAD using sparse non-negative matrix factorization (SNMF) as implemented by the snmf() function in the R-package LEA v2.6.0 [52]. We averaged admixture coefficients over 20 iterations at the previously identified optimal k = 9 for each species and compared these results to Sjodin et al. [25]. We also conducted population assignment following the method of Rannala and Mountain [53] as implemented in Geneclass2.0 [54]. Reference populations were defined using only the ddRAD samples, and both ddRAD and GT-seq samples were posteriorly assigned to these populations.

Assignment of novel invasions

We evaluated the utility of the panel for assigning individuals of unknown origin to known source populations. We examined samples collected from novel invasions on Agglomerate Island (n = 18), Ramsay Island (n = 1), and Hotspring Island (n = 1) genotyped using GT-seq. To identify the source of these invasions, we conducted projected PCA using the smartpca function from EIGENSOFT v7.2.1 [4951] and the ddRAD data to define eigenvectors. As the brown rats in Haida Gwaii are strongly clustered to northern, central, and southern regional populations [25], we ran the projected PCA using only the regional population to which the unknown samples most closely clustered during initial evaluation of population structure (see above). We also assigned individuals to reference populations using the method outlined in Rannala and Mountain [53] as implemented in Geneclass2.0 [54], including all potential reference populations. For this analysis, the Faraday and Murchison Island populations were grouped into a single reference population based on previous results [25].

Results

GT-seq panel quality and species determination

From the initial 708 putative loci, we were able to design primers targeting 526 loci. Following optimization, we retained 443 loci across both species, including 3 of 8 SNPs that were diagnostically distinct between brown and black rats. Overall, this panel had high discriminatory power to the species-level, resulting in 100% concordance with ddRAD species assignment results based on 27,686 SNPs. In the optimized GT-seq panel, all 443 loci were significantly differentiated between brown and black rats, with 150 loci exhibiting pairwise θ > 0.90. From this optimized set of 443 loci, 315 and 429 were variable within brown and black rats, respectively, and used as species-specific datasets.

Mean GT-seq genotyping rates for brown rats was 88.6%; five individuals were removed due to a genotyping rate of <50%. Mean GT-seq genotyping rates for black rats was 78.2%; two individuals were removed due to a genotyping rate of <50%. Mean genotyping discordance across methods (ddRAD v. GT-seq) was 2.8% and 3.5% in brown and black rats, respectively. No genotyping error was observed between GT-seq replicates in brown rats; black rats revealed only a 0.2% genotyping error between GT-seq replicates. All GT-seq individuals accurately assigned to species, and all individuals from the novel invasions on Agglomerate, Hotspring, and Ramsay Islands were confirmed as brown rats (Fig 2).

thumbnail
Fig 2. Posterior probability of species assignment for brown (Rattus norvegicus; n = 53; brown) and black rats (R. rattus; n = 35; black) collected in Haida Gwaii, BC.

Rats were genotyped using a GT-seq panel of 443 SNPs. Posterior probabilities were calculated by discriminant analysis of principle components as implemented in the R-package adegenet v2.1.1 [47]. “Unknown” individuals were from novel invasions on Agglomerate Island (n = 18), Ramsay Island (n = 1) and Hotspring Island (n = 1) and were not previously assigned to species.

https://doi.org/10.1371/journal.pone.0234694.g002

Assignment accuracy of brown rats

We detected the same three regional clusters in the brown rats as previously identified [25] via PCA (Fig 3a). Furthermore, GT-seq samples accurately projected to their sampled ddRAD clusters. Ancestry coefficients estimated at the reduced panel were consistent with previous results, though there was a loss of some finer-scale structure relative to the previous ddRAD analysis [25], especially among centrally located islands (Fig 4). Population assignment resulted in one individual sampled on the Bischof Islands in the GT-seq dataset assigning to a population other than the one in which it was sampled (Tanu Island), though this individual was previously identified as a first-generation migrant from Richardson Island in the full ddRAD analysis [25] (S1 Table).

thumbnail
Fig 3.

(a) Principle component analysis (PCA) for brown rats (Rattus norvegicus, n = 338) collected throughout Haida Gwaii, BC. Eigenvectors were inferred using ddRAD samples (filled circles; n = 265), with GT-seq samples (hollow circles, triangles; n = 73) posteriorly projected onto inferred eigenvectors. All analyses were performed using the same 315 loci, and analysis was implemented using EIGENSOFT v7.2.1 [4951]. PCA was re-run using only populations centrally-located in Haida Gwaii (b) for assignment of individuals from novel island invasions (triangles, n = 20).

https://doi.org/10.1371/journal.pone.0234694.g003

thumbnail
Fig 4. Admixture coefficients for n = 281 brown rats (Rattus norvegicus) sampled in Haida Gwaii, BC.

Admixture coefficients were calculated using the full ddRAD dataset of 27,686 SNPs (top) or with the reduce GT-seq panel of 315 SNPs (bottom). Estimates were obtained using the snmf() function in the R-package LEA v2.6.0 [52] at an optimal k = 9 and averaged over 20 iterations.

https://doi.org/10.1371/journal.pone.0234694.g004

Assignment accuracy of black rats

Black rats also clustered similar to previous PCA results based on ddRAD [25] (Fig 5). As with the brown rats, admixture coefficients estimated using the reduced SNP panel were consistent with previous results, with all individuals accurately assigned to their sampled population (S2 Table). However, there was a loss of some finer-scale structure among the Graham Island and Sandspit populations, as well as within the Lyell Island populations relative to the ddRAD dataset [25] (Fig 6).

thumbnail
Fig 5. Principle component analysis (PCA) for black rats (Rattus rattus, n = 249) collected throughout Haida Gwaii, BC.

Eigenvectors were inferred using ddRAD samples (filled circles; n = 214), with GT-seq samples (hollow circles; n = 35) posteriorly projected onto inferred eigenvectors. All analyses were performed using the same 429 loci, and implemented using EIGENSOFT v7.2.1 [4951].

https://doi.org/10.1371/journal.pone.0234694.g005

thumbnail
Fig 6. Admixture coefficients for 238 black rats (Rattus rattus) sampled in Haida Gwaii, BC.

Admixture coefficients were calculated using the full ddRAD dataset of 10,770 SNPs (top) or with a reduce GT-seq panel of 428 SNPs (bottom). Estimates were obtained using the snmf() function in the R-package LEA v2.6.0 [52] at an optimal k = 9 and averaged over 20 iterations.

https://doi.org/10.1371/journal.pone.0234694.g006

Identification of origin for novel invasions

All samples from novel invasions clustered with the central brown rat populations when projected onto eigenvectors inferred using all reference populations (Fig 3a); as such, we only considered central populations as putative sources. The projected PCA for the Agglomerate, Hotspring, and Ramsay Island samples indicated a Faraday or Murchison Island origin, with these samples closely clustering with Faraday and Murchison Island projected GT-seq samples (Fig 3b). Population assignment of the Agglomerate, Hotspring, and Ramsay Island individuals strongly indicated a Faraday/Murchison Island origin (Table 2).

thumbnail
Table 2. Population assignment probabilities for 20 recently invaded brown rats (Rattus norvegicus) of unknown ancestry sampled on three islands in Haida Gwaii, BC.

Individuals were assigned following the methods of Rannala and Mountain [53] as implemented in Geneclass2 [54]. (Far/Mur) These samples were collected from Faraday and Murchison Island and grouped into a single unit based on a previous study [25]. Only reference populations with >0 probability of assignment in at least one individual are shown.

https://doi.org/10.1371/journal.pone.0234694.t002

Discussion

Invasive rats are having devastating impacts on biodiversity globally, necessitating the development of innovative tools for informing IAS management. Native species, especially seabirds, are unlikely to adapt to the invasive threat in time to prevent extirpation [55]. Without intervention, these losses could result in large scale changes at the ecosystem level [11,5658]. Given the isolation of the Haida Gwaii archipelago, seabirds play a vital role in nutrient cycling as significant transporters of nutrients from off-island via guano deposition [58,59]. Disruption of this transport chain could lead to increased depletion of resources on the islands and whole-scale regime shifts, a consequence of rat invasions observed in other island systems around the globe [16,5860]. The effects of removing seabirds from Haida Gwaii are not purely ecological; seabirds, such as the rat-threatened ancient murrelet (Synthliboramphus antiquus), also have a significant cultural importance to the Haida Nation [31]. To promote seabird recovery in Haida Gwaii, invasive rats must be eradicated.

Recent work has shown the utility of genetic data for informing invasive rat management in Haida Gwaii [25]. By harnessing the power of next-generation sequencing and ddRAD, previous work characterized population structure and connectivity among both brown and black rats throughout the archipelago, laying the foundation for robust inference of population origin for nearly any novel invasion [25]. To demonstrate application, this earlier study revealed that failed eradications on the Bischof Islands were due to re-invasion from nearby Lyell Island and not from incomplete removal [i.e., the survivor hypothesis; 25]. The mechanism of reinvasion was likely via swimming given the close geographic proximity (<500m) to Lyell Island. “Rafting” during high tide and/or commensal spread via ships remain possibilities. Additionally, the origin of novel invasions on Faraday and Murchison Islands was also identified as neighboring Lyell Island [25]. These results not only highlighted the need for potentially increased biosecurity between islands in Haida Gwaii, but also provided insights into the dispersal capacity of brown rats through temperate waters.

This initial population genomic analysis was achieved using ddRAD, a reduced representation genome sequencing approach that allows for the genotyping of dozens or more individuals at a large number of SNPs (>10,000) in a single lane of Illumina sequencing [61]. Yet, this method requires labour-intensive library preparation, where speed and cost-effectiveness do not scale linearly for smaller sample sizes. These considerations are of the utmost importance, as novel rat invasions often occur with relatively few individuals [62,63]; in some cases, only a single individual may be captured for analysis when populations are at low densities (for example, on Hotspring and Ramsay Islands). In addition, ddRAD requires substantial bioinformatic processing and filtering of data post-sequencing to ensure genotype calls are accurate and not the result of sequencing error, sample contamination, or other confounding factors [64]. While ddRAD does afford connectivity between datasets, these bioinformatic processing steps must be repeated for inclusion of any new individuals. Taken together, ddRAD and similar approaches are not ideal when applied to a single or few samples, especially when the required information is time-sensitive. Rat invasions, in particular, have been likened to other acute impacts to biodiversity such as wildfires or oil spills that require rapid and decisive responses based on timely information [65].

GT-seq overcomes these limitations, providing a more expedient and cost-effective approach for generating needed genetic information, especially when involving a single or few high priority samples. To start, GT-seq library preparation is highly simplified compared to ddRAD [24,61]. While ddRAD has many time-consuming purification and size selection steps along with multiple ligations and PCR steps, GT-seq only requires two PCRs and a final library purification. Consequently, GT-seq libraries can be prepared in a matter of hours rather than days or weeks as may be the case with ddRAD. Library preparation for GT-seq uses few reagents and a commercially available PCR kit, whereas ddRAD library preparation requires a fuller suite of more costly reagents including restriction enzymes and ligases. Moreover, bioinformatic processing and analyses of raw data are also much more streamlined in GT-seq, which uses open-source software (https://github.com/GTseq/GTseq-Pipeline) to simultaneously identify locus-specific primers and genotype the corresponding SNP using minimal computational resources (e.g., standard laptop). These genotypes require no additional filtering and are immediately suitable for use in population genetic analyses, further positioning GT-seq as an effective rapid response tool.

While GT-seq does offer many advantages over other approaches (e.g., reduced representation sequencing, including ddRAD), panel development does require existing genomic data for locus identification and primer design [24,46]. For this study, we had previously genotyped brown and black rats across their distributions in Haida Gwaii, as well as characterized population structure among island populations [25]. This previous work not only provided substantial genomic data for GT-seq panel design, but also allowed for an explicit evaluation of the ability of the panel to identify known population structure based on ddRAD results. While we were fortunate in this regard, many systems are lacking these foundational data. In such systems, prior genetic work will be needed to best inform GT-seq panel development; however, once completed, the full advantages of GT-seq over more costly and labour-intensive methods can be exploited [24,46].

Moving forward, GT-seq and other targeted amplicon sequencing approaches hold great promise for informing conservation and management of biodiversity [66]. Originally demonstrated in steelhead trout [24], GT-seq has been subsequently applied to the study of other fish species, including redband trout [66], brook trout [67], Pacific lamprey [68], and walleye [69]. More recently, the efficacy of using GT-seq for genotyping minimally-invasive samples has been demonstrated [46], and successfully applied to further understanding of the molecular ecology and conservation status of the at-risk Western rattlesnake [70]. To the best of our knowledge, our study is the first to demonstrate use of GT-seq to inform invasive species management to help mitigate biodiversity loss.

Panel efficacy

Across both species-specific panels, we were able to accurately and consistently detect previously identified population structure [25] among invasive rat populations in Haida Gwaii. Principle component analysis resulted in near identical clustering when compared to the previous ddRAD study [25] (Figs 3 and 5). We also identified genetic structure consistent with previous results, though there was a slight loss in resolution of fine-scale structure among highly connected populations (Figs 4 and 6) likely due to the reduced number of SNPs in the GT-seq dataset. Both brown and black rats are highly vagile species; significant connectivity and gene flow can occur among island populations separated by as much as 1 km of ocean, and terrestrial populations can remain connected at distances of >10 km [38], resulting in less differentiation and more admixture. Using large SNP datasets (≤ 25,686 loci), Sjodin et al. [25] were able to detect this fine-scale structure even among minimally differentiated populations. To detect this same structure using our reduced GT-seq panel, targeted SNPs likely will need to be added to further maximize among-cluster variation for these admixed clusters, especially for brown rats sampled on Lyell, Tanu, and Richardson Islands, and for both species with sampling distributions across larger islands like Graham and Kunghit Islands. Nevertheless, these panels still adequately captured among-population variation across both species capable of correctly assigning individuals to their true populations. Furthermore, as GT-seq is a targeted amplicon sequencing approach, this method is highly replicable and completely connectible, so additional reference populations can always be added to the database to allow for even more robust inferences. The ability to cost-effectively genotype individuals is also conducive to applications such as long-term monitoring to better estimate contemporary gene flow and inform biosecurity measures. These estimates would be particularly useful for monitoring brown rat populations on Lyell and surrounding islands, as these populations appear to be rapidly expanding and are currently the highest biosecurity threat in terms of seabird conservation [25,31].

Novel invasions

We found strong support for a Faraday and/or Murchison origin for novel invasions on Agglomerate, Hotspring, and Ramsay Islands (Fig 3b, Table 2). The Ramsay Island invasion is particularly concerning as this island hosts the most significant remaining seabird colony in the entire archipelago [31]. Rat invasions have already led to the complete, or near complete, extirpation of large seabird colonies in Haida Gwaii, such as on Langara, Kunghit and Lyell Islands [37,71]; establishment of a rat population on Ramsay Island will likely lead to the same devasting consequences for extant seabird populations. Prior to this study, we investigated the origin of the incursion onto Hotspring Island using ddRAD (Sjodin & Russello, unpublished report). These results also indicated a Faraday/Murchison Island origin and provided information to support an eradication on Hotspring Island in November 2018 with the goal of preventing spread to Ramsay Island. But for reasons discussed above, ddRAD, while accurate, is not suitable for a single or few individuals as resources and processing time do not scale down linearly for small sample sizes. For island eradications and control of incursions, rapid response is critical and targeting the response at an appropriate scale is crucial to the durability of the eradication. Eradicating the rats on Hotspring Island kept Ramsay Island rat-free, but only for a year; rats subsequently re-invaded Hotspring Island, which likely acted as a stepping-stone leading to recent detections on Ramsay Island (Wojtaszek, pers. comm.). This outcome indicates that the management response to the initial incursions on Hotspring Island was conducted at an inappropriate geographic scale; the distance from Ramsay Island that needed to remain rat free to provide some eradication durability was greater than just the most adjacent islands (Wojtaszek, pers. comm.). Armed with information from the RapidRat GT-seq panel, future initiatives could determine whether all islands within a certain area should be considered a single eradication unit or if a more localized response is appropriate. These data can also point to mechanisms of invasion, which can guide how management agencies deal with the incursion. For example, if the rat on Hotspring Island had come from another location within the archipelago via a vessel moored on the north side of the island (e.g., human-mediated), rather than via natural dispersal from adjacent islands, the response would be scaled, and biosecurity measures enhanced.

Conclusions

Here, we developed, validated and deployed an effective and efficient genotyping tool (RapidRat) for informing invasive rat management in Haida Gwaii. It is important to highlight, however, that GT-seq panel development was greatly facilitated by the existence of an archipelago-wide SNP dataset previously collected via ddRAD [25]. When no such data are available, an initial population genomic analysis would most likely need to be conducted. Moreover, we initially designed RapidRat to be informative for both invasive brown and black rats, likely limiting the resolution at the single species level. Although, this does not appear to be an issue for informing invasive rat management in Haida Gwaii, independent, species-specific panels could be developed using a similar workflow as described here. The potential multi-species ascertainment bias also likely translates to a finer-level; for example, we do not expect RapidRat to be informative in other systems given that locus selection was solely based on among-population genetic variation in Haida Gwaii. Nevertheless, invasive rats constitute a global threat to biodiversity and innovative tools are required to inform management action and minimize impact on novel ecosystems. RapidRat demonstrates the applicability of GT-seq as one such IAS rapid response tool. Importantly, the framework demonstrated here for panel development and validation can be directly applied for informing invasive rat management and biodiversity conservation in other systems.

Supporting information

S1 Table. Population assignment probabilites for n = 53 brown rats (Rattus norvegicus) sampled in Haida Gwaii, BC.

Individuals were assigned to ddRAD reference populations using the methods outlined in Rannala and Mountain (1997) as implemented in the program Geneclass2 (Piry et al. 2004). Highlighted lines indicate samples that assigned to population other than their sampled population.

https://doi.org/10.1371/journal.pone.0234694.s001

(XLSX)

S2 Table. Population assignment probabilites for n = 34 black rats (Rattus rattus) sampled in Haida Gwaii, BC.

Individuals were assigned to ddRAD reference populations using the methods outlined in Rannala and Mountain (1997) as implemented in the program Geneclass2 (Piry et al. 2004).

https://doi.org/10.1371/journal.pone.0234694.s002

(XLSX)

Acknowledgments

We thank Parks Canada Agency employees from Gwaii Haanas National Park Reserve, National Marine Conservation Area Reserve and Haida Heritage site for the collection of the field samples. We are particularly grateful to Danielle Schmidt and Nathan Campbell for providing support with GT-seq panel development and optimization.

References

  1. 1. De Vos JM, Joppa LN, Gittleman JL, Stephens PR, Pimm SL. Estimating the normal background rate of species extinction. Conserv Biol. 2015;29(2):452–62. pmid:25159086
  2. 2. Otto SP. Adaptation, speciation and extinction in the Anthropocene. Proc R Soc B Biol Sci. 2018 Nov 21;285(1891):20182047.
  3. 3. Bellard C, Cassey P, Blackburn TM. Alien species as a driver of recent extinctions. Biol Lett. 2016 Feb 29;12(2):20150623. pmid:26888913
  4. 4. Clavero M, García-Berthou E. Invasive species are a leading cause of animal extinctions. Trends Ecol Evol. 2005 Mar 1;20(3):110. pmid:16701353
  5. 5. Doherty TS, Glen AS, Nimmo DG, Ritchie EG, Dickman CR. Invasive predators and global biodiversity loss. Proc Natl Acad Sci. 2016 Oct 4;113(40):11261–5. pmid:27638204
  6. 6. Colautti RI, Bailey SA, van Overdijk CDA, Amundsen K, MacIsaac HJ. Characterised and projected costs of nonindigenous species in Canada. Biol Invasions. 2006 Jan 1;8(1):45–59.
  7. 7. Government of Canada CFIA. 2013–2014 Departmental Performance Report [Internet]. [cited 2019 Dec 31]. https://www.inspection.gc.ca/about-the-cfia/accountability/reports-to-parliament/2013-2014-dpr/eng/1409769354767/1409769355486?chap=0#c32s3c
  8. 8. Pimentel D. Environmental and economic costs associated with alien invasive species in the United States. In: Biological Invasions: Economic and Environmental Costs of Alien Plant, Animal, and Microbe Species, Second Edition. CRC Press; 2011. p. 411–30.
  9. 9. Kier G, Kreft H, Lee TM, Jetz W, Ibisch PL, Nowicki C, et al. A global assessment of endemism and species richness across island and mainland regions. Proc Natl Acad Sci. 2009 Jun 9;106(23):9322–7. pmid:19470638
  10. 10. Jones HP, Holmes ND, Butchart SHM, Tershy BR, Kappes PJ, Corkery I, et al. Invasive mammal eradication on islands results in substantial conservation gains. Proc Natl Acad Sci. 2016 Apr 12;113(15):4033–8. pmid:27001852
  11. 11. Atkinson IA. The spread of commensal species of Rattus to oceanic islands and their effects on island avifaunas. In: Conservation of island birds: case studies for the management of threatened island species: proceedings of a symposium held at the XVIII ICBP World Conference in Cambridge, England, in August 1982 under the chairmanship of Sir Peter Scott. Cambridge, England: International Council for Bird Preservation; 1985. p. 35–81. (ICBP technical publication).
  12. 12. Gurevitch J, Padilla DK. Are invasive species a major cause of extinctions? Trends Ecol Evol. 2004 Sep 1;19(9):470–4. pmid:16701309
  13. 13. Towns DR, Atkinson IAE, Daugherty CH. Have the harmful effects of introduced rats on islands been exaggerated? Biol Invasions. 2006 Jun 1;8(4):863–91.
  14. 14. Howald G, Donlan CJ, Galván JP, Russell JC, Parkes J, Samaniego A, et al. Invasive rodent eradication on islands. Conserv Biol. 2007 Oct;21(5):1258–68. pmid:17883491
  15. 15. Jones HP, Tershy BR, Zavaleta ES, Croll DA, Keitt BS, Finkelstein ME, et al. Severity of the effects of invasive rats on seabirds: a global review: Effects of Rats on Seabirds. Conserv Biol. 2008 Feb;22(1):16–26. pmid:18254849
  16. 16. Harper GA, Bunbury N. Invasive rats on tropical islands: Their population biology and impacts on native species. Glob Ecol Conserv. 2015 Jan;3:607–27.
  17. 17. Damin-Pernik M, Espana B, Lefebvre S, Fourel I, Caruel H, Benoit E, et al. Management of rodent populations by anticoagulant rodenticides: toward third-generation anticoagulant rodenticides. Drug Metab Dispos. 2017 Feb 1;45(2):160–5. pmid:27934637
  18. 18. MacKay JWB, Russell JC, Murphy EC. Eradicating house mice from islands: successes, failures and the way forward. Manag Vertebr Invasive Species Proc Internatial Symp. 2007 Aug 1;294–304.
  19. 19. Piertney SB, Black A, Watt L, Christie D, Poncet S, Collins MA. Resolving patterns of population genetic and phylogeographic structure to inform control and eradication initiatives for brown rats Rattus norvegicus on South Georgia. J Appl Ecol. 2016;53(2):332–9.
  20. 20. Russell JC, Broome KG. Fifty years of rodent eradications in New Zealand: another decade of advances. N Z J Ecol Christch. 2016;40(2):197–204.
  21. 21. Holmes ND, Griffiths R, Pott M, Alifano A, Will D, Wegmann AS, et al. Factors associated with rodent eradication failure. Biol Conserv. 2015 May;185:8–16.
  22. 22. Amos W, Nichols HJ, Churchyard T, Brooke M de L. Rat eradication comes within a whisker! A case study of a failed project from the South Pacific. R Soc Open Sci. 2016 Apr;3(4):160110. pmid:27152226
  23. 23. Russell JC, Miller SD, Harper GA, MacInnes HE, Wylie MJ, Fewster RM. Survivors or reinvaders? Using genetic assignment to identify invasive pests following eradication. Biol Invasions. 2010 Jun;12(6):1747–57.
  24. 24. Campbell NR, Harmon SA, Narum SR. Genotyping-in-Thousands by sequencing (GT-seq): A cost effective SNP genotyping method based on custom amplicon sequencing. Mol Ecol Resour. 2015;15(4):855–67. pmid:25476721
  25. 25. Sjodin B, Irvine R, Ford A, Howald GR, Russello M. Rattus population genomics across the Haida Gwaii archipelago provides a framework for guiding invasive species management. Evol Appl. 2019 Dec 16;eva.12907.
  26. 26. Parks Canada Agency G of C. Wildlife—Gwaii Haanas National Park Reserve, National Marine Conservation Area Reserve, and Haida Heritage Site [Internet]. 2019 [cited 2019 Oct 30]. https://www.pc.gc.ca/en/pn-np/bc/gwaiihaanas/nature/faune-wildlife
  27. 27. Calder JA, Taylor RL, Mulligan GA. Flora of the Queen Charlotte Islands. Ottawa: Queen’s Printer; 1968. (Canada. Dept. of Agriculture. Research Branch. Monograph).
  28. 28. Foster JB. Evolution of the mammals of Haida Gwaii, British Columbia. Occas Pap Br Columbia Prov Mus. 1965;14:130.
  29. 29. Gaston AJ, Golumbia TE, Martin J-L, Sharpe ST. Lessons from the islands: introduced species and what they tell us about how ecosystems work. In: Proc Res Group Introd Species 2002 Symp Queen Charlotte City Queen Charlotte Isl Br Columbia Can Wildl Serv Spec Publ Environ Can Ott. 2008;
  30. 30. Golumbia TE. Introduced species management in Haida Gwaii (Queen Charlotte Islands). In: Proc Conf Biol Manag Species Habitats Risk Kamloops BC. 1999;1:327–32.
  31. 31. Harfenist A. Seabird Colonies Background Report for the Haida Gwaii/Queen Charlotte Islands Land Use Plan. Ministry of Sustainable Resource Management; 2003.
  32. 32. Moodie GEE, Reimchen TE. Endemism and conservation of sticklebacks in the Queen Charlotte islands. Can Field-Nat. 1973;87:173–5.
  33. 33. Bertram DF, Nagorsen DW. Introduced rats, Rattus spp., on the Queen Charlotte Islands: implications for seabird conservation. Can Field-Nat. 1995;109(1):6–10.
  34. 34. Harrison C. Ancient warriors of the north Pacific: the Haidas, their laws, customs and legends, with some historical account of the Queen Charlotte Islands. London: H. F. & G. Witherby; 1925. 222 p.
  35. 35. Parks Canada Agency. A natural priority–A report on Parks Canada’s Conservation and Restoration Program. Ottawa: Parks Canada Agency; 2018 p. 47.
  36. 36. Gill C, Wein L, Howald G, McClelland P. Challenges associated with eradicating invasive rodents from islands: lessons learned. Proc Vertebr Pest Conf [Internet]. 2014 [cited 2019 Dec 21];26. https://escholarship.org/uc/item/0xr266wk
  37. 37. Kaiser GW. The Langara Island seabird habitat recovery project: eradication of Norway rats, 1993–97. Delta, BC: Environment Canada, Canadian Wildlife Service, Pacific & Yukon Region, Environmental Conservation Branch; 1997. 85 p. (Technical report series / Canadian Wildlife Service; no. 304 ISSN 0831-6481).
  38. 38. Tabak MA, Poncet S, Passfield K, Carling MD, Martinez del Rio C. The relationship between distance and genetic similarity among invasive rat populations in the Falkland Islands. Conserv Genet. 2015 Feb 1;16(1):125–35.
  39. 39. Tabak MA, Poncet S, Passfield K, del Rio CM. Modeling the distribution of Norway rats (Rattus norvegicus) on offshore islands in the Falkland Islands. NeoBiota. 2015 Jan 16;24:33–48.
  40. 40. Spennemann DHR. Distribution of rat species (Rattus spp.) on the atolls of the Marshall Islands: past and present dispersal. Atoll Res Bull. 1997;446:1–18.
  41. 41. Catchen J, Hohenlohe PA, Bassham S, Amores A, Cresko WA. Stacks: an analysis tool set for population genomics. Mol Ecol. 2013 Jun 1;22(11):3124–40. pmid:23701397
  42. 42. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012 Apr;9(4):357–9. pmid:22388286
  43. 43. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011 Aug 1;27(15):2156–8. pmid:21653522
  44. 44. Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984 Nov 1;38(6):1358–70. pmid:28563791
  45. 45. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, Windows TM software for population genetics. Laboratoire génome, populations, interactions, CNRS UMR, 5000; 2004.
  46. 46. Schmidt DA, Campbell NR, Govindarajulu P, Larsen KW, Russello MA. Genotyping‐in‐Thousands by sequencing (GT‐seq) panel development and application to minimally invasive DNA samples to support studies in molecular ecology. Mol Ecol Resour. 2019 Sep 24;1755–0998.13090.
  47. 47. Jombart T, Ahmed I. adegenet 1.3–1: new tools for the analysis of genome-wide SNP data. Bioinformatics. 2011 Nov 1;27(21):3070–1. pmid:21926124
  48. 48. Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010 Oct 15;11(1):94.
  49. 49. Galinsky KJ, Bhatia G, Loh P-R, Georgiev S, Mukherjee S, Patterson NJ, et al. Fast principal-component analysis reveals convergent evolution of ADH1B in Europe and East Asia. Am J Hum Genet. 2016 Mar;98(3):456–72. pmid:26924531
  50. 50. Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLOS Genet. 2006 Dec 22;2(12):e190. pmid:17194218
  51. 51. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006 Jul 23;38:904. pmid:16862161
  52. 52. Frichot E, François O. LEA: An R package for landscape and ecological association studies. Methods Ecol Evol. 2015 Aug 1;6(8):925–9.
  53. 53. Rannala B, Mountain JL. Detecting immigration by using multilocus genotypes. Proc Natl Acad Sci. 1997 Aug 19;94(17):9197–201. pmid:9256459
  54. 54. Piry S, Alapetite A, Cornuet J-M, Paetkau D, Baudouin L, Estoup A. GENECLASS2: A software for genetic assignment and first-generation migrant detection. J Hered. 2004 Nov 1;95(6):536–9. pmid:15475402
  55. 55. Strauss SY, Lau JA, Carroll SP. Evolutionary responses of natives to introduced species: what do introductions tell us about natural communities? Ecol Lett. 2006;9(3):357–74. pmid:16958902
  56. 56. Harris DB, Macdonald DW. Interference competition between introduced black rats and endemic Galápagos rice rats. Ecology. 2007;88(9):2330–44. pmid:17918410
  57. 57. Clout MN. Biodiversity loss caused by invasive alien vertebrates. Z Für Jagdwiss. 2002 Dec 1;48(1):51–8.
  58. 58. Mulder CPH, Grant-Hoffman MN, Towns DR, Bellingham PJ, Wardle DA, Durrett MS, et al. Direct and indirect effects of rats: does rat eradication restore ecosystem functioning of New Zealand seabird islands? Biol Invasions. 2009 Aug;11(7):1671–88.
  59. 59. Fukami T, Wardle DA, Bellingham PJ, Mulder CPH, Towns DR, Yeates GW, et al. Above- and below-ground impacts of introduced predators in seabird-dominated island ecosystems. Ecol Lett. 2006 Dec;9(12):1299–307. pmid:17118004
  60. 60. Towns DR, Wardle DA, Mulder CPH, Yeates GW, Fitzgerald BM, Parrish GR, et al. Predation of seabirds by invasive rats: multiple indirect consequences for invertebrate communities. Oikos. 2009 Mar 1;118(3):420–30.
  61. 61. Peterson BK, Weber JN, Kay EH, Fisher HS, Hoekstra HE. Double digest RADseq: An inexpensive method for de novo SNP discovery and genotyping in model and non-model species. Orlando L, editor. PLoS ONE. 2012 May 31;7(5):e37135. pmid:22675423
  62. 62. Russell JC, Beaven BM, MacKay JWB, Towns DR, Clout MN. Testing island biosecurity systems for invasive rats. Wildl Res. 2008;35(3):215.
  63. 63. Russell JC, Towns DR, Anderson SH, Clout MN. Intercepting the first rat ashore. Nature. 2005 Oct;437(7062):1107–1107. pmid:16237434
  64. 64. Shafer ABA, Peart CR, Tusso S, Maayan I, Brelsford A, Wheat CW, et al. Bioinformatic processing of RAD-seq data dramatically impacts downstream population genetic inference. Methods Ecol Evol. 2018 May 1;907–17.
  65. 65. Ramsey D, Will D. Framework for undertaking eradication programs on insular populations of vertebrate pests. Melbourne: Arthur Rylah Institute for Enviornmental Research; 2012 Nov
  66. 66. Meek MH, Larson WA. The future is now: Amplicon sequencing and sequence capture usher in the conservation genomics era. Mol Ecol Resour. 2019;19(4):795–803. pmid:30681776
  67. 67. Delomas TA. Differentiating diploid and triploid individuals using single nucleotide polymorphisms genotyped by amplicon sequencing. Mol Ecol Resour. 2019;19(6):1545–51. pmid:31448561
  68. 68. Parker KA, Hess JE, Narum SR, Kinziger AP. Evidence for the genetic basis and epistatic interactions underlying ocean- and river-maturing ecotypes of Pacific Lamprey (Entosphenus tridentatus) returning to the Klamath River, California. Mol Ecol. 2019;28(13):3171–85. pmid:31132186
  69. 69. Bootsma ML, Gruenthal KM, McKinney GJ, Simmons L, Miller L, Sass GG, et al. A GT-seq panel for walleye (Sander vitreus) provides a generalized workflow for efficient development and implementation of amplicon panels in non-model organisms. bioRxiv. 2020 Feb 14;2020.02.13.948331.
  70. 70. Schmidt DA, Govindarajulu P, Larsen KW, Russello MA. Genotyping-in-Thousands by sequencing reveals marked population structure in Western Rattlesnakes to inform conservation status. Ecol Evol. in press;
  71. 71. Taylor RH, Kaiser GW, Drever MC. Eradication of Norway rats for recovery of seabird habitat on Langara island, British Columbia. Restor Ecol. 2000 Jun 1;8(2):151–60.