Hybridization can profoundly affect the genomic composition and phenotypes of closely related species, and provides an opportunity to identify mechanisms that maintain reproductive isolation between species. Recent evidence suggests that hybridization outcomes within a species pair can vary across locations. However, we still do not know how variable outcomes of hybridization are across geographic replicates, and what mechanisms drive that variation. In this study, we described hybridization outcomes across 27 locations in the North Fork Shoshone River basin (Wyoming, USA) where native Yellowstone cutthroat trout and introduced rainbow trout co‐occur. We used genomic data and hierarchical Bayesian models to precisely identify ancestry of hybrid individuals. Hybridization outcomes varied across locations. In some locations, only rainbow trout and advanced backcrossed hybrids towards rainbow trout were present, while trout in other locations had a broader range of ancestry, including both parental species and first‐generation hybrids. Later‐generation intermediate hybrids were rare relative to backcrossed hybrids and rainbow trout individuals. Using an individual‐based simulation, we found that outcomes of hybridization in the North Fork Shoshone River basin deviate substantially from what we would expect under null expectations of random mating and no selection against hybrids. Since this deviation implies that some mechanisms of reproductive isolation function to maintain parental taxa and a diversity of hybrid types, we then modelled hybridization outcomes as a function of environmental variables and stocking history that are likely to affect prezygotic barriers to hybridization. Variables associated with history of fish stocking were the strongest predictors of hybridization outcomes, followed by environmental variables that might affect overlap in spawning time and location.

We used historical stocking and population survey records of Yellowstone Cutthroat Trout Oncorhynchus clarkii bouvieri and other salmonids in the North Fork Shoshone River drainage, Wyoming to summarize fish stocking history and population trends. Based on 98 years of historical records, we found that despite extensive stocking of Yellowstone Cutthroat Trout and minimal stocking of nonnative salmonids after about 1950, populations of wild Yellowstone Cutthroat Trout declined relative to those of nonnative salmonid species. The timing of increases in nonnative salmonids (1970s) did not coincide with their period of most intensive stocking (1935–1950). It is plausible that Yellowstone Cutthroat Trout populations persisted because of high levels of supplemental stocking from 1935 to 1965 and declined with reduced stocking efforts in the 1970s, thereby allowing the increase of introduced nonnative salmonids. The establishment of nonnative salmonids likely further reduced stocking success of Yellowstone Cutthroat Trout due to competition and hybridization. This study demonstrates that an understanding of long‐term stocking records and population survey data can be useful for developing and implementing successful management frameworks for the conservation of imperiled fish populations across the United States.