Watershed boundaries and geographic isolation: patterns of diversification in cutthroat trout from western North America

Background For wide-ranging species, intraspecific variation can occur as a result of reproductive isolation from local adaptive differences or from physical barriers to movement. Cutthroat trout (Oncorhynchus clarkii), a widely distributed fish species from North America, has been divided into numerous putative subspecies largely based on its isolation in different watersheds. In this study, we examined mtDNA sequence variation of cutthroat trout to determine the major phylogenetic lineages of this polytypic species. We use these data as a means of testing whether geographic isolation by watershed boundaries can be a primary factor organizing intraspecific diversification. Results We collected cutthroat trout from locations spanning almost the entire geographic range of this species and included samples from all major subspecies of cutthroat trout. Based on our analyses, we reveal eight major lineages of cutthroat trout, six of which correspond to subspecific taxonomy commonly used to describe intraspecific variation in this species. The Bonneville cutthroat trout (O. c. utah) and Yellowstone cutthroat trout (O. c. bouvieri) did not form separate monophyletic lineages, but instead formed an intermixed clade. We also document the geographic distribution of a Great Basin lineage of cutthroat trout; a group typically defined as Bonneville cutthroat trout, but it appears more closely related to the Colorado River lineage of cutthroat trout. Conclusion Our study indicates that watershed boundaries can be an organizing factor isolating genetic diversity in fishes; however, historical connections between watersheds can also influence the template of isolation. Widely distributed species, like cutthroat trout, offer an opportunity to assess where historic watershed connections may have existed, and help explain the current distribution of biological diversity across a landscape.


Background
Species with wide geographic ranges often exhibit significant morphological, behavioral, or genetic variability across their range [1]. Such differences have been attributed to local adaptation and genetic drift as a result of isolation of populations by limited gene flow [2]. Although the importance of within-species diversity has been increasingly recognized by subspecific taxonomy, distinct population segments, and evolutionarily significant units [3,4], most studies have focused on understanding how phenotypic variability can influence the distributional limits of a species [5]. Intraspecific variability can represent significant adaptive potential allowing a species to occupy a broad range of ecological conditions or adapt to changing environmental conditions [6]. In addition to understanding the potential importance of phenotypic limits to a species distribution, past studies have also tried to identify physiographic features that may facilitate, organize, or constrain diversification of species that are distributed over large geographic areas [7,8].
Across a landscape, geographic isolation may arise as a consequence of barriers to movement from mountain building, glacial events, continental drift or a combination of factors [9]. While wide-ranging terrestrial species can exhibit geographic isolation and population structuring [10,11], aquatic species and freshwater fishes, in particular, often exhibit strong population structure likely resulting from their confinement to the network or 'ribbons' of water that flow across the landscape. The height of land that separates watershed boundaries and the relatively common occurrence of movement barriers within watersheds often isolate fish populations allowing for significant genetic differentiation [12][13][14][15].
In addition to the slow and continuous erosive changes to watersheds, major drainage patterns can be changed by events such as volcanic or glacial flows, or by climatic changes in precipitation that can connect or isolate watersheds [16][17][18]. Newly created watershed connections are thought to be important influences shaping community composition for organisms like fishes that are obligatorily confined to the network of streams and lakes by allowing species that were previously isolated to expand their geographic range [8]. In contrast, populations of the same species that become isolated by loss of connectivity may diversify into new forms representing early stages of speciation [7].
The cutthroat trout (Oncorhynchus clarkii; Teleostei: Salmonidae) is a widely distributed fish species native to western North America that occurs along a north-south axis from Alaska to New Mexico ( Figure 1). Populations of cutthroat trout are found in very diverse habitats ranging from coastal temperate rainforest watersheds to desert lakes and streams [19,20]. Perhaps not surprisingly, cutthroat trout exhibit dramatic morphological and life history diversity across these contrasting ecological conditions and have been divided into as many as 14 different subspecies to account for much of this diversity [19,21]. Despite the range of ecological diversity encompassed by cutthroat trout, different subspecies are primarily organized by major watershed boundaries that are often contiguous to each other and separated only by the height of land between adjacent areas (Figure 1). In addition to diversification created by isolation in different watersheds, cutthroat trout also occupy northern areas that were covered in ice during the last glacial maximum and re-invaded from southern ice-free watersheds by connections during pluvial periods. In western North America, the last period of glaciation reached a maximum about 18 000 years ago and covered much of present day Canada as well as northern portions of adjoining areas in the United States [22]. Previous comparisons of species occupying formerly glaciated and non-glaciated regions often reveal lower levels of genetic diversification in formerly glaciated portions of their range [17]; indicating the importance of climatic changes in the environment for species diversification.
Given their range and distribution, cutthroat trout populations offer an excellent opportunity to examine how geographic features such as major watershed boundaries can organize and shape genetic diversification in wide-ranging fish species. Similarly, their distribution across geographic regions that have experienced significant changes in climatic conditions also offers an opportunity to understand how physiographic changes in climate can influence the template of diversification by watershed isolation.
In this study we examine how the diversity and evolutionary history of a widely distributed fish species can be organized by primary geographic features such as watershed boundaries and how such organization may also be influenced by historic events that have connected watersheds but are no longer visible on the landscape. Although intraspecific taxonomic categories are often used in the management of cutthroat trout populations [23] there is remarkably little comparative data to assess whether subspecific differences capture or reliably organize the diversity present in this widely distributed species. In this study we sampled cutthroat trout from throughout their geographic range to compile a comprehensive database of mtDNA sequence data as a means of documenting the evolutionary history of cutthroat trout. By doing so, we test how well geographic features and past subspecific divisions of cutthroat trout represent the diversity within this species and how past connections between watersheds may have influenced this organization.

Study sites
We sampled extant populations of cutthroat trout from 311 locations representing the majority of their geographic range (Figure 1, Additional file 1: Table S1) to compare mitochondrial (mtDNA) sequence divergence at the NADH dehydrogenase subunit 2 gene (ND2) and determine major phylogenetic branching among groups. Past taxonomic evaluations of cutthroat trout have listed up to 14 different subspecies of cutthroat trout [20,24]. In our study we included samples from eight major subspecies that are commonly recognized, including:  [19,20,25]. Two other subspecies, the Alvord cutthroat trout (O. c. alvordensis) and yellowfin cutthroat trout (O. c. macdonaldi), are thought to be extinct [20]. Because we were interested in examining geographic patterns of natural variation in sequence divergence, we sampled populations that are believed to be representative of native populations. In addition, we sampled nine rainbow trout (O. mykiss) populations to use as an outgroup for our phylogenetic analyses (Additional file 1: Table S1).
To collect cutthroat trout from as wide a geographic range as possible, we used a variety of sampling techniques. In streams, we used a backpack or boat-mounted electroshocker or angled for fish. For lake populations, we sampled fish by angling, or by using multi-panel monofilament gillnets, minnow traps, or trap nets. Figure 1 Historical range of cutthroat trout and sampling localities. Sampling locations of cutthroat trout from western North America (red dots) within native ranges of eight major subspecies of cutthroat trout (colored polygons). Historical distributions are based on Behnke [24] and McPhail [64].
Captured fish were held in a 20 L container with fresh water, and a 3-5 mm fin clip was taken from each fish and stored in a uniquely numbered vial with 95% ethanol. In four locations, supplementary genetic samples were provided by researchers using similar sampling methods.

GenBank sequences
In addition to the 311 locations we physically sampled or were sampled for us, we also searched the National Center for Biotechnology Information's GenBank database http://www.ncbi.nlm.nih.gov/genbank/ for ND2 sequences of cutthroat trout. We included samples from this database in geographic comparisons whenever we could identify locations either in the GenBank database or from corresponding publications that provided location information (Additional file 2: Table S2).

DNA extraction and mtDNA sequencing
DNA was extracted from fin clips using the ZR Genomic DNA tissue extraction kit (Zymo Research) following the manufacturer's protocol. Mitochondrial DNA variation was analyzed using the entire ND2 gene and was amplified by polymerase chain reaction (PCR) using the sequencing primers NDintF6 (5' TAAGCTTTCGGGCC-CATACC 3') and NDVarR (5'GCT TTG AAG GCT CTT GGT CT 3') [26,27]. Our PCR reactions were performed in 25 μl volumes using 8 μL of 2X ReddyMix PCR Master Mix (ABgene), 1 μL (10 mM) of each primer, and 2 μL of genomic DNA. The thermal profile included an initial 94°C denaturation followed by 35 cycles at 94°for 30 s, annealing at 58°C for 45 s, and extension at 72°C for 75 s, with a final extension at 72°C for 10 min. PCR products were submitted to the University of Washington Highthroughput genomics unit for purification and DNA sequencing.

Sequence data analysis
We edited and aligned the sequences using Sequencher v4.9 (Gene Codes Corporation) and the online version of Muscle [28]. Unique sequences were submitted to Gen-Bank with accession numbers as reported in Additional file 1: Table S1. To illustrate the relationship between unique cutthroat trout haplotypes, we constructed a minimum spanning network using representatives of each unique ND2 haplotype and the program Arlequin v3.5 [29]. Mitochondrial DNA polymorphism was estimated as haplotype (h) and nucleotide (π) diversity, as well as the percent sequence divergence between groups using DNAsp v5 [30] and Mega v5.01 [31] software.

Phylogenetic analyses
For our phylogenetic analyses, we included one sequence representing each unique haplotype, as well as rainbow trout sequences, which were used as an outgroup. Phylogenies were estimated by maximum likelihood (ML) and Bayesian analyses. For our ML analysis, we used the Tamura-Nei substitution model with invariant sites and gamma-distributed rates among sites based on jModeltest [32] results and generated 1000 bootstrap replicates as implemented in the program Phyml v3.0 [33]. Our Bayesian analysis was conducted using MrBayes v3.1.2 [34,35] with the GTR + I + Γ model based on MrModeltest [36]. The analyses included two independent runs for 5.5 million generations each and were sampled every 100 th generation. The first 25% of samples were discarded as "burn-in" to ensure sampling from a stationary posterior distribution.
We used the primary lineages derived in phylogenetic analyses to compare the sequence diversity within and among groups and estimate their divergence times. The proportion of diversity within and among primary phylogenetic groups was estimated using an AMOVA as implemented in Arlequin v3.5. We estimated divergence times between primary cutthroat trout lineages by first testing for equal evolutionary rates between lineages using Tajima's relative rate test [37] and then calibrating and linearizing our phylogenetic tree using Mega 5.01 software. We applied the estimated divergence time of 6 million years ago (mya) to the node between rainbow trout and cutthroat trout to calibrate our tree. This divergence time is based on the combination of fossil data and DNA sequence data [38] and is within the range of proposed estimates of divergence between rainbow trout and cutthroat trout [39].

Results
Nucleotide sequences were generated for ND2 (1050 bp) for 384 trout representing 311 populations and 102 different haplotypes. We also compiled sequence data from an additional 95 cutthroat trout samples (Yellowstone cutthroat trout [n = 71], Lahontan cutthroat trout [n = 4], Rio Grande cutthroat trout [n = 10], Colorado River cutthroat trout [n = 6], and greenback cutthroat trout [n = 4]) deposited in the GenBank database [27,40,41]; providing an additional 38 haplotypes. Yellowstone cutthroat trout ND2 sequences from GenBank were 1050 bp, all other subspecies had ND2 sequences of 889 bp in length. Hence, when we combined our sequences with those from Genbank (both 889 bp and 1050 bp), a total of 140 different haplotypes of cutthroat trout were available for comparison (Table 1). A minimum spanning network of all haplotypes with complete sequences (1050 bp) revealed seven main clusters of haplotypes that corresponded to some of the subspecific classifications of cutthroat trout ( Figure 2). Haplotypes from the Bonneville and Yellowstone cutthroat trout did not form distinct clusters, but were intermixed. In contrast, haplotypes from populations of coastal, Colorado River, Lahontan, Rio Grande, and westslope cutthroat trout formed separate clusters. Interestingly, we also identified an additional cluster that included populations from within the Bonneville and Yellowstone cutthroat trout distribution, which we refer to as the Great Basin cutthroat trout lineage because of its geographic distribution within the central portion of the Great Basin region of North America (see below). The number of mutational steps separating each major group ranged from 8 to 17 ( Figure 2).

Primary phylogenetic lineages
Our maximum likelihood analysis identified eight major lineages of cutthroat trout based on mtDNA sequence differences ( Figure 3a). Of these eight main groups, four, including coastal, westslope, Lahontan, and Colorado River cutthroat trout, occurred as monophyletic lineages and corresponded to subspecies designations previously used to describe intraspecific variation in cutthroat trout. Two groups (Rio Grande cutthroat trout and greenback cutthroat trout), while forming separate groups in our phylogeny, were not well supported. As indicated in the minimum spanning network, Bonneville and Yellowstone cutthroat trout did not form monophyletic lineages, but were mixed into two well-supported and very divergent clades. The Bonneville-Yellowstone clade included the majority of the sampling locations for these subspecies; however, a distinct and divergent clade is also present in the phylogeny (the Great Basin lineage), and is more closely related to populations from the Colorado River watershed despite its geographic proximity to the upper Snake River and Bear River ( Figure  3a). The Paiute and Whitehorse cutthroat trout were nested within the Lahontan branch of the phylogeny (Figure 3a), and fish representing fine-spotted Snake River cutthroat trout had the same haplotype of other Yellowstone cutthroat trout. When we repeated the   Table S1 and Additional file 2: Table S2. Hash marks (up to eight) indicate the number of mutational changes. One step is indicated by a line with no hash marks and numbers adjacent to lines indicate more than eight mutational changes. Only cutthroat trout haplotypes with 1050 bp are included in the network.
phylogenetic analysis using a Bayesian approach, we found the arrangement of branches with similar or greater support based on posterior probabilities ( Figure 3b). When we compared the variability in sequence data of cutthroat trout, we found that a predominant proportion was accounted for by differences between the major clades identified. Among the eight major clades, 71.84% of the variability in sequence data was accounted for by differences among clades; whereas 28.16% was accounted for by differences within clades (AMOVA, P < 0.0001). Excluding the four unique haplotypes deposited in GenBank for greenback cutthroat trout, haplotype diversity was highest in Bonneville-Yellowstone cutthroat trout and lowest in the westslope cutthroat trout; while nucleotide diversity was also highest in the Bonneville-Yellowstone group, it was lowest in Lahontan cutthroat trout group (Table 1) Across all major cutthroat clades, average sequence divergence was 2.4% (range: 1.3-3.2; Table 2). Average divergence of  cutthroat trout lineages with rainbow trout was 8% (range: 7.3-8.7%; Table 2).

Geographic sub-structuring
Within several of the major clades of cutthroat trout, there was further geographically related sub-structuring among haplotypes. The Bonneville-Yellowstone clade, which included two subspecies of cutthroat trout, formed two primary subclades (Figure 4). Within this group, clade B included a mixture of haplotypes collected from geographic locations in the Bonneville and upper Snake River watersheds. In contrast, haplotypes from the clade A group were only found in the upper Snake River and upper Missouri River watersheds (Figure 5). Within the coastal cutthroat trout lineage three primary subdivisions were identified ( Figure 6). All three were found in areas below the southern end of Vancouver Island, British Columbia, but only haplotypes  Red circles represent sampling locations for the clade of Great Basin cutthroat trout (see Figure 8). Colored polygons represent the estimated boundary for the native range of Bonneville (blue) and Yellowstone (Yellow) subspecies of cutthroat trout [24]. The dashed line represents the extent of the Bear River watershed, excluding the Malad River sub-basin. associated with one of the subdivisions were found in more northern latitudes of British Columbia and southeast Alaska (Figure 7). The westslope cutthroat trout lineage also had two major subdivisions ( Figure 6). Both were found in the southern portions of its range in Idaho, but only haplotypes associated with one of the lineages were found in more northerly areas of Montana, British Columbia, and Alberta (Figure 7). Of the five different haplotypes observed in the Lahontan lineage of cutthroat trout, there was no large divergence among the geographic samples collected. The Whitehorse population did have a unique haplotype that was slightly different than populations from the Lahontan basin; whereas the Paiute cutthroat populations all had the same haplotype that was also found in the Humboldt and Reese River systems of Nevada (Additional file 1: Table S1). Within the Colorado River lineage we detected no well supported sub-structuring; haplotypes were distributed throughout the historic range of this subspecies (Additional file 1: Tables S1 and Additional file 2: Table S2). Some of the GenBank sequences also identified the Colorado River lineage across the Continental Divide into the range of greenback cutthroat trout [41]; and one greenback cutthroat trout was found east of the Continental Divide in the Colorado River range. As Metcalf and collaborators [41] note, the greatest divergence in greenback cutthroat trout occurs between Bear Creek, Colorado and the three remaining unique haplotypes (Figure 8). Similar to the findings of Pritchard and others [42], the Rio Grande lineage had two unique haplotypes that were most different from all other haplotypes (Figure 8), one of the two was found in   Figure 6 Phylogeny of Lahontan, westslope, and coastal cutthroat trout. Detailed maximum likelihood phylogeny of the Lahontan, westslope, and coastal lineages of cutthroat trout. Branch tips labeled as L represent haplotypes sampled within the historical distribution of Lahontan cutthroat trout; whereas those labeled as W and C represent haplotypes sampled within the range of westslope and coastal cutthroat trout, respectively (see Figures 1 and 7 for geographic ranges). Numbers in parentheses represent the number of sampling locations where the haplotype was detected followed by the haplotype number. Numbers on branches indicate percent bootstrap support based on 1000 replicates. Scale bar represents the proportion of sequence divergence in the ND2 gene used to construct the phylogeny. Figure 7 Geographic distribution of coastal and westslope cutthroat trout lineages. Geographic distribution of three major subclade lineages of coastal cutthroat trout (see Figure 6) depicted by red circles (clade A), white circles (clade B), and black circles (clade C); and for two major subclade lineages of westslope cutthroat trout depicted by white circles (clade A) and black circles (clade B). Colored polygons represent the estimates extent of the native range of coastal (blue) and westslope (brown) cutthroat trout [24,64].
the upper Canadian River watershed of New Mexico (Ricardo Creek), and one in Dalton Creek, a stream in the Pecos River watershed. All other haplotypes were found in the Rio Grande watershed. Within the Great Basin lineage only minor differences in sequence divergence were evident (Figure 8).

Divergence times
Results of the relative rate test indicated that we were able to apply a molecular clock to our data. Thus, based on the assumption of 6 million years between rainbow trout and cutthroat trout divergence, we estimated that all cutthroat trout shared a common ancestor approximately 1.97 mya, when there was a primary divergence between the lineage that includes the Bonneville-Yellowstone, Lahontan, westslope, and coastal groups versus the lineage including the greenback, Rio Grande, Great Basin, and Colorado River groups. Within each of these major cutthroat trout clades, the Bonneville-Yellowstone lineage was estimated to have diverged from other members of its group at 1.9 mya, followed by the coastal group at 1.85 mya, then between the westslope and Lahontan clades about 1.43 mya (Figure 3a). The Colorado River lineage of cutthroat trout was estimated to have diverged within its subclade at 1.41 mya, followed by the Great Basin lineage at 1.38 mya, and finally between the greenback and Rio Grande lineages at 0.94 mya (Figure 3a).  Figure 8 Phylogeny of greenback, Rio Grande, Great Basin, and Colorado River cutthroat trout lineages. Detailed maximum likelihood phylogeny of the greenback, Rio Grande, Great Basin, and Colorado River lineages of cutthroat trout. Branch tips labeled as G, RG, and CR represent haplotypes sampled within the historical distribution of greenback, Rio Grande, and Colorado River cutthroat trout; whereas those labeled as GB represent haplotypes sampled within the range of Bonneville or Yellowstone cutthroat trout, respectively (see Figures 1 and 5 for geographic ranges). Numbers in parentheses represent the number of sampling locations where the haplotype was detected followed by the haplotype number. Numbers on branches indicate percent bootstrap support based on 1000 replicates. Scale bar represents the proportion of sequence divergence in the ND2 gene used to construct the phylogeny.

Discussion
In this study we examined how cutthroat trout diversity may be organized by isolation into major geographic boundaries. Our data indicate that watershed boundaries do organize primary phylogenetic divisions in cutthroat trout, but historical connections that are no longer evident also appear to have influenced the pattern of diversification. Our range-wide analysis of cutthroat trout diversification provides an additional line of evidence for historic watershed connections that have long been suggested by geologists and are of interest to biologists as a means of explaining the current distribution of fish fauna across a landscape [8,16,18,43].
Intraspecific taxonomic evaluations of cutthroat trout have described up to 14 different subspecies of cutthroat trout, two of which are now considered to be extinct [20,24]. Of the subspecies we evaluated, the Bonneville-Yellowstone, Colorado River, Rio Grande, greenback, and Lahontan lineages are delimited at least in part by the height of land defining major portions of watersheds. Based on the mtDNA sequence data we compared and the partitioning of genetic variance, cutthroat trout can be divided into eight major lineages that correspond with six of the primary subspecies of cutthroat trout commonly recognized. Hence, watershed boundaries and the intraspecific taxonomic categories that have been used to describe this diversity have been relatively successful in capturing a significant component of the evolutionary diversification in cutthroat trout.
If a molecular clock is applied to the sequence divergence, we estimate an initial divergence from a common ancestor approximately 2 mya. This primary divergence led to one major clade that includes the Bonneville-Yellowstone, Lahontan, westslope, and coastal lineages. The second branch includes the greenback, Rio Grande, Great Basin, and Colorado River lineages. The pattern of haplotype diversification indicates that the lineage leading to the Bonneville-Yellowstone clade was first to diverge from the common ancestor of all cutthroat trout, which may have first colonized the ancestral Snake River as well as areas that are now part of Colorado River watershed. Although the pattern of colonization can never be known with exact certainty, drainage patterns have changed significantly over the period of time estimated from the divergence of the ancestral cutthroat trout. The greatest divergence of cutthroat trout between the Snake and Colorado River drainages indicates a former connection between the two drainages that has since isolated these two main lineages for the longest period of time. Geologic evidence does point to connections between areas that are now in the upper Colorado River watershed but were once part of the Bear River of Idaho, Utah, and Wyoming [16].
Within each of the main lineages of cutthroat trout, the pattern of relatedness also indicates how several sub-basins may have been connected to each other. The next closest associations within the Bonneville-Yellowstone group are those from coastal cutthroat trout populations followed by a split between populations from the Lahontan basin (and adjacent areas) and westslope cutthroat trout from the Columbia River as well as populations from eastern drainages across the Continental Divide. Downstream colonization and origin of coastal cutthroat trout could have occurred if the ancestral Snake River was connected to the Columbia watershed or other rivers draining to the Pacific Ocean, but the association of the coastal cutthroat trout with westslope cutthroat trout found in Columbia drainages suggests that the Snake River was part of the Columbia watershed when the lineages emerged. Despite the proximity of some parts of the Lahontan basin to areas of the Colorado River watershed, Hubbs and Miller [43] noted the similarity of the Lahontan fish fauna with that of Columbia, Sacramento, and Death Valley watersheds, indicating that the Lahontan basin was at one time connected to one of these systems and not the Colorado River watershed. Our sequence data on cutthroat trout again point to a connection with the Columbia watershed based on similarity of lineages also found between the watersheds. Finally, the associations of the greenback, Rio Grande, and Great Basin lineages with that of the Colorado lineage all indicate that they were derived in an earlier form of the Colorado River basin. Based on the evolutionary history of the ND2 gene, our data indicate a complex history of isolation and evolutionary divergence of the ancestral Colorado lineage of cutthroat trout. This divergence led to the current Colorado lineage as well as the Rio Grande and Great Basin lineages. The lineage that led to the Rio Grande also appears to have given rise to the greenback lineage, which subsequently colonized its current distribution. The presence of greenback cutthroat trout in eastern draining rivers, or rivers draining to the south in the case of Rio Grande cutthroat, indicate transfers of populations through headwater connections and subsequent isolation. As Minckley and others [16] note, fish species that commonly penetrate into headwater streams often have representative populations across drainage divides presumably by finding connections, even if they occur infrequently, and are further isolated over longer time periods by the elevation of land.
Two subspecies that did not form monophyletic lineages were the Bonneville and Yellowstone cutthroat trout. The core distribution of Bonneville cutthroat trout is thought to occur in the Bear River, a watershed that originates in the Uinta Mountains of northeastern Utah and flows northward into Wyoming and Idaho before turning sharply back southward towards Utah and the Great Salt Lake ( Figure 5). Such a dramatic change of direction is probably representative of stream capture, and in fact the Bear River is thought to have been a tributary to the upper Snake River [44]. Between 600 000 and 50 000 years ago, lava flows crossed the Soda Springs area of the Bear River valley and forced the path of the river to change direction southward [44]. A past connection between the Bear River and upper Snake River would explain the intermixing of the primary lineages of cutthroat trout found in these two watersheds, as well as the morphological similarities observed between these two subspecies [24,45]. The presence of the two main cutthroat trout lineages in the upper Snake River, but only one in the Bear River would suggest that it originated in the Bear River and subsequently colonized downstream to the Snake River, whereas the second lineage was unable to colonize upstream to the Bear River.
A new perspective generated from our range wide analysis was the geographic distribution and position of what we have referred to as the Great Basin clade of cutthroat trout. Although some previous studies have detected and noted the dramatic genetic divergence of representative samples from this lineage [27,39], our study provides a clearer picture of its geographic distribution and its position within the phylogeny of cutthroat trout. Our data indicate that despite its proximity to the Bear River and upper Snake River watersheds, it is more closely related to the major diversification of cutthroat trout that includes the Colorado River, greenback, and Rio Grande clades of cutthroat trout. Such a relationship again points to historic connections between watersheds now isolated from each other. As Smith and collaborators [39] note, previous studies have documented possible hydrological connections between the Bonneville basin and the Colorado River watershed. Connections between these watersheds would be congruent with the Great Basin lineage of cutthroat trout arising in an ancient Colorado River area that either invaded the Bonneville basin by a past connection or evolved within the basin when some parts drained toward the Colorado River watershed. The presence of the Great Basin cutthroat lineage in the upper Snake River may also represent past hydrological connections between the Snake River and the Bonneville basin, but were probably more recent, corresponding to well-documented Pleistocene connections [44]. Approximately 15 000 years ago, Red Rock pass in southeast Idaho was the path of overflow of Lake Bonneville into the upper Snake River travelling through upper Marsh Creek and the Portneuf River; the same location where we detected the Great Basin cutthroat trout lineage [46]. The second location where we observed the Great Basin lineage in the upper Snake River watershed was in the Raft River drainage in southern Idaho, about 150 km to the west of the Portneuf River ( Figure 5). This area borders the Bonneville basin and our data suggest that part of the current Raft River drainage was previously in the Bonneville basin and has been captured by the Snake River watershed, isolating the Great Basin cutthroat trout in the headwater tributary streams observed today. Such captures may not be too unexpected because of the gradual subsidence of the Snake River plain that may capture streams bordering the region [47]. As noted in phylogeographic studies of galaxiid fish in New Zealand, geographic patterns of genetic diversity can provide an additional line of evidence to infer historical changes to drainage patterns on the landscape [48].
Watershed boundaries and barriers to movement can organize major phylogenetic lineages for freshwater fishes in some cases but appear to be most important when populations have been isolated within them for extended periods of time. Glaciated areas often appear to have been more recently connected, allowing major lineages to disperse from glacial refugia over large geographic areas, as illustrated by examples of lake trout (Salvelinus namaycush) and lake whitefish (Coregonus clupeaformis). Each of these species is thought to have dispersed through large proglacial lakes that submerged or connected a number of contemporary watersheds in North America [49,50]. Populations of freshwater fish that have not established themselves by dispersal into recently de-glaciated areas appear much more likely to exhibit significantly greater evolutionary diversification as a consequence of isolation within watershed boundaries [13,15,51]. As a widely distributed species spanning both glaciated and non-glaciated regions of North America, cutthroat trout exhibit reduced haplotype diversity at the northern periphery of its range and higher levels to the south. A reduction in genetic diversity in northern regions is concordant with a pattern of dispersal from a past southern glacial refuge, illustrating the importance of past geologic events on the phylogenetic structure of a freshwater fish species [17].
As a widely distributed trout species in North America, cutthroat trout have received a great deal of interest as a sport fish, but have also suffered severe declines in abundance when populations have been affected by human plans for consumption of water, land, or introductions of non-native species (see Trotter [19] for a review). Despite early and more recent attempts to describe the diversity of cutthroat trout [21,24] it is somewhat surprising that there is remarkably little comparative data to support subspecific designations commonly used to manage different cutthroat trout populations [52]. Perhaps as a result of a lack of such data, only one taxonomic ranking of cutthroat trout is currently recognized by the joint Committee on Names of Fishes, sponsored by the American Fisheries Society and the American Society of Ichthyologists and Herpetologists [53], and possibly because of the difficulty in diagnosing subspecies of cutthroat trout with traditional meristic characters [54]. A number of previous studies have provided evidence for significant genetic diversification of cutthroat trout [55][56][57][58] and even phylogenetic comparisons of representative populations from different parts of the species range [27,39,59]. However, our study provides the first range-wide comparison attempting to identify major phylogenetic lineages over much of the geographic area occupied by this polytypic species.
Our analysis provides a better understanding of the main evolutionary lineages of cutthroat trout and their geographic distribution. Similar components have been proposed as a basis for organizing management and conservation units of species [3,60]. Although primary evolutionary lineages of cutthroat trout may be a logical starting point for management units, even with the addition of genetic diversity at nuclear loci, they may not capture important ecological variation that can occur has as result of local adaption to specific conditions [61,62]. Hence, in addition to measures of mtDNA and nuclear DNA divergence, ecological diversification not always apparent from molecular phylogenies should also be considered in attempts to preserve the diversity present within a species. Indeed, salmonid fishes often exhibit significant ecological variation [63] and cutthroat trout in particular have been observed to exhibit significant morphological variability that is associated with specific ecological conditions [45]. Fortunately, biologists have often used the precautionary principle in arguing for protecting populations with unusual life history or ecological relationships even in the absence of direct comparative data. By combining comprehensive range-wide phylogenetic comparisons to identify the distribution of major species lineages, with an understanding of the main ecological factors shaping the phenotype of a species, biologists are likely to provide the best chance of protecting the evolutionary potential of a species in a changing landscape.

Conclusions
Cutthroat trout populations have diversified into major phylogenetic lineages in western North America. While geographic isolation within major watershed boundaries is a primary organizing factor of genetic diversity in cutthroat trout, past hydrological connections have also influenced the evolutionary history of cutthroat trout. Our study illustrates how genetic diversity can be used to identify where historic watershed connections may have existed and how that information can help explain the current distribution of biological diversity across a landscape.

Additional material
Additional file 1: Table S1. Population, phylogenetic grouping, geographic location, and accession number of haplotypes for trout sampled in this study.
Additional file 2: Table S2. Phylogenetic clade, haplotype, and accession number for trout used from Genbank database.