Genetic Structure of Remnant Populations and Cultivars of Switchgrass (Panicum virgatum) in the Context of Prairie Conservation and Restoration

Switchgrass (Panicum virgatum L.) is a dominant, perennial C₄ grass of North American tallgrass prairies with cultivars that are widely used in grassland restoration, pastures, and landscaping. However, these cultivars may be genetically dissimilar to small, remnant populations, raising concerns about altered genetic composition of native populations through gene flow. To address this issue on a local scale in Ohio and Illinois, we used microsatellite markers to characterize genetic diversity and differentiation of 10 remnant prairie populations (5 in each state) and 8 common cultivars. The bulk of genetic variation was found to reside within rather than among wild populations, consistent with the outcrossing breeding system of switchgrass. Genetic diversity was similar among the remnant populations despite large differences in area (approximately 2–2,590 ha), highlighting the importance of small native populations as reservoirs of variation and potential seed sources for prairie restoration. Cultivars generally had similar levels of variation to the wild populations, but we found clear genetic dissimilarity between wild and cultivated gene pools (especially for Kanlow, but also Trailblazer, Blackwell, Dacotah, Summer, and Sunburst cultivars). This suggests that using cultivars in local prairie restoration efforts may alter the genetic composition of wild populations. Whether such changes are deemed as negative depends on the cultivar under consideration and specific conservation goals for preserving native switchgrass populations. Patterns of genetic variation in remnant prairie populations and potential cultivar sources can be used to develop guidelines for restoration as well as future planting of cultivars for biofuels.     

Random amplified polymorphic DNA variation among remnant big bluestem (Andropogon gerardii vitman) populations from Arkansas' grand prairie

Random amplified polymorphic DNA (RAPD) analysis was used to characterize genetic diversity and genetic distinctiveness of Andropogon gerardii from remnant Arkansas prairies. Six oligonucleotide primers, which generated 37 RAPD bands, were used to analyse 30-32 plants from six Grand Prairie populations, Baker Prairie (Arkansas Ozarks), two Illinois prairies and two cultivars. Genetic diversity of the Arkansas remnants ranged from 82.7 to 99.3%, with 89% of the total genetic variation within and 11% among populations. The partitioning of genetic variation was consistent with that reported for other outcrossing perennial grasses, using the more conservative allozyme markers. Principal component analysis indicated a northern and southern association within Arkansas' Grand Prairie. Although there was no genetic structuring at the landscape level, the Illinois prairies and cultivars were different from all Arkansas prairies tested. There was significant within-population structuring in four of the seven Arkansas remnants, with a negative relationship between genetic similarity and geographical distance. The three nonstructured populations were from a linear railroad remnant, suggesting different population-level dynamics from nonlinear prairies. The results of this study indicated that small isolated remnant big bluestem populations were not genetically depauperate and that genetic relationships among populations could not be predicted solely on geographical proximity.     

High Plant Richness in Prairie Reconstructions Support Diverse Leafhopper Communities

A primary reason for restoring plant communities is to increase biodiversity to previous levels. It is expected that restoring land with greater plant diversity will increase biodiversity at higher trophic levels, but high diversity seed mixes are expensive. In this study, we used one insect family, leafhoppers (Hemiptera: Cicadellidae) to assess the difference in leafhopper communities that result from establishing high compared with low plant richness restorations. We tested the hypotheses that: (1) the added effort of a high richness restoration leads to measurable increases in both diversity and richness of leafhoppers; and (2) that leafhopper community composition is more similar to remnant prairies in high richness than in low plant richness restorations. We found that higher plant richness led to 3‐ to 7‐fold increases in leafhopper and prairie‐dependent leafhopper diversity and richness in restorations. Leafhopper communities in high richness restorations were not more similar to remnant prairies, rather they were distinct among high and low richness restorations and prairie interior. Leafhopper richness and diversity correlated with plant richness, and leafhopper community composition differed among plant community assemblages, but not with the occurrence of single plant species. For our sites, species‐rich restorations provided better quality habitat for leafhoppers that was comparable to remnant prairie. Our results suggest that restorations with high plant species richness better support animal food webs.     

Ecotypes of an ecologically dominant prairie grass (Andropogon gerardii) exhibit genetic divergence across the U.S. Midwest grasslands' environmental gradient

Big bluestem (Andropogon gerardii) is an ecologically dominant grass with wide distribution across the environmental gradient of U.S. Midwest grasslands. This system offers an ideal natural laboratory to study population divergence and adaptation in spatially varying climates. Objectives were to: (i) characterize neutral genetic diversity and structure within and among three regional ecotypes derived from 11 prairies across the U.S. Midwest environmental gradient, (ii) distinguish between the relative roles of isolation by distance (IBD) vs. isolation by environment (IBE) on ecotype divergence, (iii) identify outlier loci under selection and (iv) assess the association between outlier loci and climate. Using two primer sets, we genotyped 378 plants at 384 polymorphic AFLP loci across regional ecotypes from central and eastern Kansas and Illinois. Neighbour‐joining tree and PCoA revealed strong genetic differentiation between Kansas and Illinois ecotypes, which was better explained by IBE than IBD. We found high genetic variability within prairies (80%) and even fragmented Illinois prairies, surprisingly, contained high within‐prairie genetic diversity (92%). Using Bayenv 2, 14 top‐ranked outlier loci among ecotypes were associated with temperature and precipitation variables. Six of seven BayeScan F_ST outliers were in common with Bayenv 2 outliers. High genetic diversity may enable big bluestem populations to better withstand changing climates; however, population divergence supports the use of local ecotypes in grassland restoration. Knowledge of genetic variation in this ecological dominant and other grassland species will be critical to understanding grassland response and restoration challenges in the face of a changing climate.     

Responses of plant and bird communities to prescribed burning in tallgrass prairies

Historic losses and fragmentation of tallgrass prairie habitat to agriculture and urban development have led to declines in diversity and abundance of plants and birds associated with such habitat. Prescribed burning is a management strategy that has potential for restoring and rejuvenating prairies in fragmented landscapes, and through such restoration, might create habitat for birds dependent upon prairies. To provide improved data for management decision-making regarding the use of prescribed fire in tallgrass prairies, we compared responses of plant and bird communities on five burned and five unburned tallgrass prairie fragments at the DeSoto National Wildlife Refuge, Iowa, USA, from 1995 to 1997. Overall species richness and diversity were unaffected by burning, but individual species of plants and birds were affected by year-treatment interactions, including northern bobwhite (Colinus virginianus) and ring-necked pheasant (Phasianus colchicus), which showed time-delayed increases in density on burned sites. Analyses of species/area relationships indicated that, collectively, many small sites did make significant contributions to plant biodiversity at landscape levels, supporting the overall conservation value of prairie fragments. In contrast, most birds species were present on larger sites. Thus, higher biodiversity in bird communities which contain area-sensitive species might require larger sites able to support larger, more stable populations, greater habitat heterogeneity, and greater opportunity for niche separation.     

Patterns of Plant Species Diversity in Remnant and Restored Tallgrass Prairies

To restore diversity of native vegetation, we must understand factors responsible for diversity in targeted communities. These factors operate at different spatial scales and may affect the number and relative abundances of species differently. We measured diversity of plant species and functional groups of species in replicated plots within paired restored and remnant (relic) tallgrass prairies at three locations in central Texas, U.S.A. To determine the contributions of species abundances and of spatial patterns of diversity to differences between prairie types, we separated diversity into richness and evenness (relative biomass) and into within‐plot (α), among‐plot (β), and prairie (γ) components. Species diversity was greater in remnant than in restored prairies at all spatial scales. At the γ scale, both species richness and species evenness were greater in remnants because of greater spatial variation in species composition. At the α scale, remnants were more diverse because of greater richness alone. Mean α richness correlated positively with the size of the species pool in restored prairies only, implying that in remnants, α richness was influenced more by colonization dynamics than by the number of species available for colonization. Plots in remnant prairies contained more functional groups and fewer species per group than did plots in restored prairies, suggesting that resource partitioning was greater in relic prairies. Our results are consistent with the interpretation that local ecological processes, like resource partitioning and limitations on seed dispersal, contribute to the greater diversity of remnant than restored prairies in central Texas. Restoration practices that limit abundances of competitive dominants, increase the number of species in seed mixtures, and increase the proximity of plants of different functional groups thus may be required to better simulate the plant diversity of tallgrass prairies.     

Inoculation with remnant prairie soils increased the growth of three native prairie legumes but not necessarily their associations with beneficial soil microbes

Restoring the diversity of plant species found in remnant communities is a challenge for restoration practitioners, in part because many reintroduced plant species fail to establish in restored sites. Legumes establish particularly poorly, perhaps because they depend on two guilds of soil microbial mutualists, rhizobial bacteria and arbuscular mycorrhizal (AM) fungi, that may be absent from restored sites. We tested the effect of soil microorganisms from remnant and restored prairies on legume growth by inoculating seedlings of Lespedeza capitata, Amorpha canescens, and Dalea purpurea with soil from 10 restored prairies and 6 remnant (untilled) prairies from southwest Michigan. We generally found support for the hypothesis that restored prairie soils lack microbes that enhance prairie plant growth, although there was variation across species and mutualist guilds. All three legumes grew larger and two legumes (Lespedeza and Amorpha) produced more nodules when inoculated with soil from remnant prairies, suggesting that low quantity and/or quality of rhizobial partners may limit the establishment of those species in restored prairies. In contrast, no legume experienced greater root colonization by AM fungi in remnant prairie soils, suggesting equivalent quantity (but not necessarily quality) of fungal partners in remnant and restored prairie soils. We detected no evidence of spontaneous recovery of the community of beneficial soil microbes in restorations. These results suggest that the absence of rhizobia, a largely overlooked component of prairie soils, could play a strong role in limiting restored prairie diversity by hindering legume establishment. Active reintroduction of appropriate rhizobial strains could enhance prairie restoration outcomes.     

Biochar soil amendments in prairie restorations do not interfere with benefits from inoculation with native arbuscular mycorrhizal fungi

Little of the historical extent of tallgrass prairie ecosystems remains in North America, and therefore there is strong interest in restoring prairies. However, slow‐growing prairie plants are initially weak competitors with the fast‐growing yet short‐lived weedy plant species that are typically abundant in recently established prairie restorations. One way to aid establishment of slow‐growing plant species is through adding soil amendments to prairie restorations before planting. Arbuscular mycorrhizal (AM) fungi form mutualisms with the roots of most terrestrial plants and are particularly important for the growth of slow‐growing prairie plant species. As prairie ecosystems are adapted to fires that leave biochar (charred organic material) in the soil, adding biochar as well as AM fungal strains from undisturbed remnant prairies into the soil of prairie restorations may improve restoration outcomes. Here, we test this prediction during the first four growing seasons of a prairie restoration. When prairie plant seedlings were inoculated prior to planting into the field with AM fungi derived from remnant prairies, that one‐time inoculation significantly increased growth of five of the nine tested plant species through at least two growing seasons. This long‐term benefit of AM fungal inoculation was unaffected by biochar addition to the soil. Biochar application rates of at least 10 tons/ha significantly decreased Coreopsis tripteris growth but acted synergistically with AM fungal inoculation to significantly improve survival of Schizachyrium scoparium. Overall, inoculation with native AM fungi can help promote prairie plant establishment, but concomitant use of biochar soil amendments had relatively little effect.     

Dominance of an invasive earthworm in native and non-native grassland ecosystems

More attention is currently being focused on earthworm invasions; however, in many ecosystems the relative abundance of native and invasive earthworm species is unknown. We characterized earthworm populations of two grassland types within the Palouse region: native prairie remnants and Conservation Reserve Program (CRP) set asides planted with exotic grasses. The earthworm community in both grassland types was completely dominated by the exotic-invasive Aporrectodea trapezoides. Only one individual of a native species, Driloleirus americanus (the giant Palouse earthworm), was found in a prairie remnant. No differences were found between prairie remnants and CRP sites for mean earthworm density (24–106 individuals m⁻²) or fresh weight (12–45 g m⁻²). Our results suggest that the combined effects of land-use change, habitat fragmentation and competitive interactions have resulted in the decimation of native earthworm populations and dominance of invasive earthworms in native and non-native grasslands of the Palouse region.     

Intraspecific Chromosome Number Variation and Prairie Restoration—A Case Study in Northeast Iowa,U.S.A

Genome duplication has played an important role in plant evolution. Variation in genome size within species is common, particularly in grasses, but rarely considered when planning restorations. We surveyed ploidy variation within one habitat (tall grass prairie) in one region (northeast Iowa, U.S.A.) to assess the risk of ploidy mismatch in restoration plantings. Genome sizes were estimated using flow cytometry for samples from 19 remnant prairies, 5 restoration plantings, and 2 seed sources. Intraspecific ploidy variation in remnant prairie populations was found for two species of grasses, Andropogon gerardii (big bluestem) and Panicum virgatum (switchgrass). Restoration seeds differed from remnants in ploidy for three species of grass, P. virgatum, Sorghastrum nutans (Indian grass), and Spartina pectinata (prairie cordgrass), and for one species of forb, Amorpha canescens (lead plant). In the case of S. pectinata, local ecotype seeds were found to consist of two different ploidy levels. Restorations in grasslands in the United States and elsewhere are likely to create mixed ploidy populations, probably resulting in lower reproductive success for the remnant population. Prevention of mixed ploidy populations will require the screening of restoration seed sources and regional surveys for ploidy variation.     

Patterns of Species Richness and Composition in Re-Created Grassland

The success of many prairie restorations is not well documented. A restoration begun in 1975 at the Fermi National Accelerator Laboratory near Chicago, Illinois allows assessment of restoration efforts as well as changes through time. Data are presented on species richness and composition for 13 restorations planted in successive years between 1975 and 1990 and two remnant prairies. Presence of species was recorded using a stratified random design. Species richness at several scales and non‐metric multidimensional scaling ordination were used to assess trends in the vegetation. Species richness declined through time at all scales examined and was always less in the restored prairies than that found in the remnant prairies. Species composition changed with time but not in the direction of the composition found in the remnants. Our understanding of the maintenance of species richness is not sufficient to allow the re‐creation of patterns of species found in remnant grassland communities.     

Legacy effects of afforestation on prairie plant and seed bank communities in a northern Canadian prairie

Afforestation resulting from fire suppression, modified grazing, plantation establishment and climate change poses a threat to northern prairie ecosystems. Trees alter the composition and function of plant and soil fauna communities and can compromise the restoration of afforested prairies. To evaluate the hypothesis that legacies of afforestation persist in restored prairie communities and decrease the potential for restoration, we examined the composition, structure, and diversity of plant and seed bank communities along a 20 year chronosequence of plantation tree removal from a northern fescue prairie in Riding Mountain National Park, Manitoba, Canada. Tree removal increased the abundance of weedy species in the plant and seed bank communities of restored prairies and plant diversity peaked and declined over the 20 year period of passive restoration. As a result, time since tree removal and the encroachment of invasive species were key in explaining the composition of restored prairie communities. Low correlation between the species composition of plant and seed bank communities, including the complete absence of Festuca hallii in restored treatments, demonstrated that legacies of afforestation compromised the potential of seed banks to facilitate prairie restoration. We conclude that tree removal alone is insufficient for the restoration of northern fescue prairies and that, in the absence of active management, the persistence of low-diversity plant and seed bank communities constitutes an important legacy of afforestation and an important barrier to future restoration.     

Colonization and Recovery of Invertebrate Ecosystem Engineers during Prairie Restoration

Ants (Hymenoptera: Formicidae) and earthworms (Oligochaeta) are considered ecosystem engineers because they form biogenic structures in the soil that influence resource supply. The objectives of this study were to quantify recovery dynamics of these invertebrate groups across a chronosequence of restored prairies and elucidate whether changes in the abundance and biomass of ants and earthworms were related to key plant and ecosystem properties. We sampled ants and earthworms from cultivated fields, grasslands restored from 1 to 21 years, and native prairie. Ant abundance and biomass peaked between 5 and 8 years of restoration and abundance was 198 times greater than cultivated fields. Earthworm abundance increased linearly across the chronosequence and became representative of native prairie, but all earthworm populations were dominated by European species. Ant abundance and biomass were positively correlated with plant diversity and plant richness, whereas earthworm abundance biomass was only related to surface litter. These results demonstrate that earthworm abundance increases with time since cessation of cultivation and concomitant with prairie establishment, whereas the abundance and biomass of ants are more related to the structure of restored plant communities than time. The dominance of exotic earthworms in these restorations, coupled with their capacity to alter soil properties and processes may represent novel conditions for grassland development.     

A paired study of prairie carbon stocks, fluxes, and phenology: comparing the world's oldest prairie restoration with an adjacent remnant

We measured carbon (C) stocks and fluxes and vegetation phenology in the world's oldest prairie restoration (∼65 years) and an adjacent prairie remnant in southern Wisconsin from 2001–2004 to quantify structural and functional differences. While the species distributions and frequency differed, the number of species measured per 1 m² quadrat were not significantly different (15.8±4.4 and 14.1±2.1 for remnant and planted [order for all reported values in abstract]; P=0.29), and the annual average aboveground net primary productivity (271±51 and 330±55 g C m⁻²) and peak leaf area index (2.9–4.9 m² m⁻²) were comparable under similar fire management. Total root biomass was not significantly different in 2002 (1736±1062 and 1690±459 g dry matter m⁻²) or 2003 (3029±2081 and 2146±898 g m⁻²), but annual average soil respiration (1229±77 and 1428±24 g C m⁻² yr⁻¹) was significantly higher in the restoration (P<0.0001). However, the prairie remnant contained 37% greater soil C (P<0.0001) in the top 25 cm. Soil respiration response to 10 cm soil temperature (Q₁₀) varied with respect to prairie and soil moisture conditions as annual Q₁₀ values ranged from 2.5 to 3.6. We calculated a range of net ecosystem production (NEP) values using estimated heterotrophic respiration and three root turnover values. Average NEP varied from −1.4 to 1.9 and −2.3 to 1.3 Mg C ha⁻¹ yr⁻¹ for the remnant and planted prairies, respectively. While these two prairies share similar structural components and functional attributes, the large uncertainty in NEP casts doubt as to whether we can verify these prairies as C sources or sinks without direct measures of heterotrophic respiration and root turnover. We argue that quantitative studies of C exchange in prairies, which differ in restoration methodology, management intensity, and fire frequency, are needed to solidify the relationship between prairie structure and potentially desired functions such as C sequestration.     

Microsatellite analysis reveals interpopulation differentiation and gene flow in the endangered tree Changiostyrax dolichocarpa (Styracaceae) with fragmented distribution in central China

Polymorphic simple sequence repeat (SSR) markers were used to investigate the impact of habitat fragmentation on the population structure and gene flow of Changiostyrax dolichocarpa, a critically endangered tree in central China. Intrapopulation genetic diversity, population structure and gene flow in the five extant populations of this species were analysed by eight SSR markers. Intrapopulation genetic diversity results suggest that C. dolichocarpa remnants maintained a relatively high degree of genetic diversity despite severe fragmentation. Low genetic differentiation among populations was found based on Wright's F(ST) and amova analysis. Both the F(ST)-based estimate and private allele method revealed high historical gene flow among the remnant populations. Recent immigrants, detected by assignment tests, tend to decrease from the grandparent generation to the current generation. The potentially highly restricted current gene flow among fragments may render the fragmented populations of C. dolichocarpa at a higher risk of local extinction several generations after fragmentation. Both in situ and ex situ conservation management for the remnant populations of C. dolichocarpa are therefore urgently needed to rescue remaining genetic diversity.     

Insect Herbivory of Ovules and Seeds in Native and Restored Prairies

Insect herbivory of ovules and seeds is known to have a negative impact on individual plant reproductive success by reducing overall seed set. Although this seed set reduction has been well documented in plant populations found in native habitats, little work has been done on populations found in restored habitats. In a 4-year study (2001–2004), I investigated the herbivory of the Gelechiid moth ( Coleotechnites eryngiella ) of Eryngium yuccifolium populations in native and restored prairies. Data were collected from 20 E. yuccifolium populations (10 native prairies/10 restored prairies) in Illinois. Percent herbivory and percent seed set were determined for each population. From 2001 to 2004, percent herbivory in prairie restorations ranged from 0 to 93 and percent seed set ranged from 2 to 82. In native prairies, percent herbivory ranged from 0 to 98 and percent seed set ranged from 0 to 71. No differences were found between native and restored prairies for percent herbivory or percent seed set. However, significant differences were found among years for percent seed set. This study shows that the antagonistic effects of insect herbivory can reach similar levels in restored and native prairies.     

Population genetic consequences of geographic disjunction: a prairie plant isolated on Great Lakes alvars

Species may often exhibit geographic variation in population genetic structure due to contemporary and historical variation in population size and gene flow. Here, we test the predictions that populations on the margins of a species' distribution contain less genetic variation and are more differentiated than populations towards the core of the range by comparing patterns of genetic variation at five microsatellite loci between disjunct and core populations of the perennial, allohexaploid herb Geum triflorum. We sampled nine populations isolated on alvar habitat within the eastern Great Lakes region in North America, habitats that include disjunct populations of several plant species, and compared these to 16 populations sampled from prairie habitat throughout the core of the species' distribution in midwestern Canada and the USA. Alvar populations exhibited much lower within-population diversity and contained only a subset of alleles found in prairie populations. We detected isolation by distance across the species' range and within alvar and prairie regions separately. As predicted, genetic differentiation was higher among alvar populations than among prairie populations, even after controlling for the geographic distance between sampled populations. Low diversity and high differentiation can be accounted for by the greater contemporary spatial isolation of alvar populations. However, the genetic structure of alvar populations may also have been influenced by postglacial range expansion and contraction. Our results are consistent with alvar populations being founded during an expansion of prairie habitat during the warmer, hypsithermal period approximately 5000 bp and subsequently becoming stranded on isolated alvar habitat as the climate grew cooler and wetter.     

Wild bee community change over a 26‐year chronosequence of restored tallgrass prairie

Restoration efforts often focus on plants, but additionally require the establishment and long‐term persistence of diverse groups of nontarget organisms, such as bees, for important ecosystem functions and meeting restoration goals. We investigated long‐term patterns in the response of bees to habitat restoration by sampling bee communities along a 26‐year chronosequence of restored tallgrass prairie in north‐central Illinois, U.S.A. Specifically, we examined how bee communities changed over time since restoration in terms of (1) abundance and richness, (2) community composition, and (3) the two components of beta diversity, one‐to‐one species replacement, and changes in species richness. Bee abundance and raw richness increased with restoration age from the low level of the pre‐restoration (agricultural) sites to the target level of the remnant prairie within the first 2–3 years after restoration, and these high levels were maintained throughout the entire restoration chronosequence. Bee community composition of the youngest restored sites differed from that of prairie remnants, but 5–7 years post‐restoration the community composition of restored prairie converged with that of remnants. Landscape context, particularly nearby wooded land, was found to affect abundance, rarefied richness, and community composition. Partitioning overall beta diversity between sites into species replacement and richness effects revealed that the main driver of community change over time was the gradual accumulation of species, rather than one‐to‐one species replacement. At the spatial and temporal scales we studied, we conclude that prairie restoration efforts targeting plants also successfully restore bee communities.     

Flower visitor communities are similar on remnant and reconstructed tallgrass prairies despite forb community differences

One common goal of habitat restoration and reconstruction is to reinstate the biodiversity found at intact reference sites. However, few researchers have examined whether these practices reinstate communities of flower‐visiting insects. This is unfortunate, as anthropogenically mediated declines in flower visitors, including bees (the primary pollinators for most terrestrial ecosystems), beetles, flies, and butterflies, have been reported worldwide. Biodiversity declines may be especially severe in North America's tallgrass prairie, a once‐vast grassland that has experienced severe destruction and degradation due to agricultural conversion. As such, we assessed the structure of forb and flower‐visiting insect communities as a whole and two subsets of the flower visitor community—bees and phytophagous beetles—across five tallgrass prairie remnants and five reconstructed prairies (former crop fields) in Kansas from 2013 to 2015. Remnant prairies had significantly higher forb diversity and differed significantly in forb composition, compared to reconstructed prairies. Despite the dissimilarities in forb community structure, there were no differences in flower visitor diversity or abundance between remnants and reconstructed prairies. However, when considered separately, bee communities exhibited significantly greater variability in composition on reconstructed prairies, likely due to the abundance of generalist bee species visiting non‐native legumes at two reconstructed prairies. Our work provides evidence that prairie habitat reconstruction is a valuable tool for reestablishing flower‐visiting insect communities and also emphasizes the considerable role that non‐native species may play in structuring grassland plant–bee interactions.     

Restoration Potential of Bighorn Sheep in a Prairie Region

Efforts to recover Rocky Mountain bighorn sheep (Ovis canadensis canadensis) throughout western North America have had limited success with the majority of current populations remaining in small and isolated areas on a fraction of their historical range. Prairie environments with rugged topography throughout the Northern Great Plains ecoregion were historically occupied by relatively robust bighorn sheep populations. We predicted there is likely unrealized potential habitat for restoring bighorn sheep to these areas; however, relatively little attention has been devoted to identifying habitat in unoccupied prairie regions. We used global positioning system (GPS)-collar data collected from 43 female bighorn sheep in 2 populations located in the eastern Montana, USA, portion of the Northern Great Plains during 2014–2018 to estimate a population-level annual resource selection model and identify the important factors affecting bighorn sheep resource selection. We extrapolated model predictions across eastern Montana's prairie region and identified potential habitat to understand restoration potential and assist with future translocations of bighorn sheep. Resource selection of bighorn sheep was most strongly associated with terrain slope and ruggedness, tree canopy cover, and a normalized difference vegetation index metric. Within currently unoccupied areas of the historical range, the model extrapolation predicted 7,211 km² of habitat, with most owned and managed by private landowners (44%), Bureau of Land Management (33%), and the United States Fish and Wildlife Service (15%). Our results provide an empirical evaluation of landscape covariates influencing resource selection of bighorn sheep occupying prairie environments and provide a habitat model that may be generalizable to other areas in the Northern Great Plains ecoregion. We demonstrate substantial potential for restoration opportunities of bighorn sheep in the Northern Great Plains ecoregion. Broad restoration of bighorn sheep across the ecoregion would likely require strong collaboration among and between public resource managers, private landowners, and livestock producers given the heterogeneous land ownership patterns, management strategies, and domestic sheep distributions. © 2020 The Wildlife Society.