Diets during critical brooding and winter periods likely influence the growth of Lesser Prairie-Chicken (Tympanuchus pallidicinctus) populations. During the brooding period, rapidly growing Lesser Prairie-Chicken chicks have high calorie demands and are restricted to foods within immediate surroundings. For adults and juveniles during cold winters, meeting thermoregulatory demands with available food items of limited nutrient content may be challenging. Our objective was to determine the primary animal and plant components of Lesser Prairie-Chicken diets among native prairie, cropland, and Conservation Reserve Program (CRP) fields in Kansas and Colorado, USA, during brooding and winter using a DNA metabarcoding approach. Lesser Prairie-Chicken fecal samples (n = 314) were collected during summer 2014 and winter 2014–2015, DNA was extracted, amplified, and sequenced. A region of the cytochrome oxidase I (COI) gene was sequenced to determine the arthropod component of the diet, and a portion of the trnL intron region was used to determine the plant component. Relying on fecal DNA to quantify dietary composition, as opposed to traditional visual identification of gut contents, revealed a greater proportion of soft-bodied arthropods than previously recorded. Among 80 fecal samples for which threshold arthropod DNA reads were obtained, 35% of the sequences were most likely from Lepidoptera, 26% from Orthoptera, 14% from Araneae, 13% from Hemiptera, and 12% from other orders. Plant sequences from 137 fecal samples were composed of species similar to Ambrosia (27%), followed by species similar to Lactuca or Taraxacum (10%), Medicago (6%), and Triticum (5%). Forbs were the predominant (>50% of reads) plant food consumed during both brood rearing and winter. The importance both of native forbs and of a broad array of arthropods that rely on forbs suggests that disturbance regimes that promote forbs may be crucial in providing food for Lesser Prairie-Chickens in the northern portion of their distribution.
Knowledge of how starvation, predation, and thermoregulation interact to regulate Lesser Prairie-Chicken populations (Tympanuchus pallidicinctus) is limited, in part, by a lack of knowledge of diets during critical ecological periods (McNamara and Houston 1987, Newton 1998, Patten et al. 2005, Haukos and Zavaleta 2016). Lesser Prairie-Chicken populations have experienced long-term declines and continue to decline in areas that appear to provide good-quality habitat at broad scales (Garton et al. 2016, Rodgers 2016, Spencer et al. 2017). Minimizing the degradation of remaining available habitat will require a comprehensive understanding of Lesser Prairie-Chicken biology, including dietary needs. Lesser Prairie-Chicken diets have not been well described but appear to be variable throughout the year (Olawsky 1987, Haukos and Zavaleta 2016). Most diet information is based on information from individuals collected in autumn over a small part of the species' range (Crawford and Bolen 1976, Smith 1979, Riley et al. 1993, Haukos and Zavaleta 2016). However, availability of food resources during brood rearing and winter may be most limiting for galliforms (Sedinger 1997, Sandercock et al. 2008, Hagen et al. 2009). Rapidly growing Lesser Prairie-Chicken, and other grouse (Phasianidae), chicks have high calorie demands and are restricted to foods within their immediate surroundings (Bergerud and Gratson 1988, Lautenbach 2015). For adults and juveniles, meeting thermoregulatory demands with available food items of limited nutrient content may be challenging during cold winters (Moss 1983, Olawsky 1987, Sedinger 1997).
During the brooding period, adult Lesser Prairie-Chickens and chicks consume an array of invertebrate taxa and are thought to specialize on grasshoppers (Orthoptera; Jones 1964, Suminski 1977, Davis et al. 1980). Yet this conclusion is based on only a few studies that assessed diets from crop and fecal contents and from sampling available invertebrates at locations visited by Lesser Prairie-Chickens (Haukos and Zavaleta 2016). Sampled plant and arthropod abundance may not always be a good estimator of food availability, and diets cannot always be assumed on the basis of association (Jones 1964, Davis et al. 1980, Litvaitis 2000). At feeding sites, the size, mobility, and phenology of invertebrates should constrain which arthropods are considered available prey for Lesser Prairie-Chicken chicks. Variation in arthropod prey vulnerability and availability at feeding sites, even within species, must be considered to identify optimal diets; a lack of accounting for this association may lead to erroneous conclusions (Sih and Christensen 2001).
Although arthropods are important food sources for Lesser Prairie-Chickens during summer and fall, Lesser Prairie-Chickens typically rely on plant matter to fulfill energetic demands during winter and spring (Haukos and Zavaleta 2016). Several research efforts have assessed winter diets in sand shinnery oak (Quercus havardii) prairie, where Lesser Prairie-Chickens readily use oak catkins and acorns when available (Jones 1964, Suminski 1977, Pettit 1986, Riley et al. 1993). Outside of periods when acorns are produced, and outside of the sand shinnery oak prairie, winter foods are less known (Salter et al. 2005, McDonald et al. 2014). The reliance on persistent woody vegetation during the winter months is well documented for grouse species, and Lesser Prairie-Chickens can make use of woody vegetation other than sand shinnery oak (Schmidt 1936, Schwilling 1955, Bergerud and Gratson 1988). For example, budding willows (Salix spp.) and cottonwoods (Populus deltoides) can be used during winter, as can portions of sand sagebrush (Artemisia filifolia) and skunkbrush sumac (Rhus aromatica; Schwilling 1955, Jones 1963). However, consumption of budding woody vegetation may be minimal in prairie-chickens in comparison to other grouse (Schmidt 1936).
Compared to other grouse, prairie-chickens may specialize on forb seeds and waste grain during winter (Schmidt 1936). Waste grain (e.g., Sorghum spp., Zea spp.) can provide an energy-rich food source for adult upland gamebirds (Evans and Dietz 1974, Bogenschutz et al. 1995, Guthery 2000). Use of grain fields by Lesser Prairie-Chickens has been reported during fall through early spring (Jamison et al. 2002); however, occurrence of Lesser Prairie-Chickens in cultivated fields has not been correlated with the amount of waste grain or related to increased body condition, survival, or reproductive output (Salter et al. 2005, Haukos and Zavaleta 2016). In addition to corn and sorghum, alfalfa (Medicago spp.) may be an important food resource in early spring (Jamison 2000, Larsson et al. 2013). It has been suggested that Lesser Prairie-Chickens use alfalfa fields primarily for the moisture content, and provision of moisture may make alfalfa fields more attractive than wheat (Triticum spp.; Larsson et al. 2013). Additionally, alfalfa may be used by prairie-chickens because it is richer in protein than other herbaceous foods (Mowat et al. 1965). In portions of their range removed from cultivation, broom snakeweed (Gutierrezia sarothrae), annual buckwheat (Eriogoum annum), and Johnny-jump-up (Viola spp.) may be primary winter food sources for Lesser Prairie-Chickens (Jones 1963).
True impacts on demography and contributions of food sources in the diet are difficult to estimate using traditional methods based on crop contents or scat dissection. For example, analysis of crop contents usually requires the harvesting of individuals and thus precludes any estimated impact on survival. Such post mortem analyses are not practical for species of concern. Microhistological analyses of feces are another option that can provide inference, and are noninvasive, but may underestimate easily digestible items (Bartolome et al. 1995, Litvaitis 2000). Additionally, not all contents in the crop are ultimately digested. Some of the material stored in the crop can be regurgitated (Jordan 2005). Therefore, DNA metabarcoding of fecal samples might be the best option for linking avian diets to fitness because it can identify prey items for species of conservation concern when collection of individuals is not practical (Pompanon et al. 2012). Instead of collecting individual crop samples, a standardized DNA region, or barcode, is identified that varies among, but is neutral within, taxa of interest. The DNA barcode region is amplified from fecal samples and compared to sequences from a reference database; then the relative contribution of food items can be estimated, based on the frequency of sequences (Ratnasingham and Hebert 2007, Zeale et al. 2011, Craine et al. 2015). DNA metabarcoding can be a particularly useful method for identifying soft-bodied arthropod prey items, which can be detected only by expert examination of gut contents or by histology of fecal samples (Burger et al. 1999, Zeale et al. 2011, Trevelline et al. 2016).
To estimate the effects of food availability on Lesser Prairie-Chicken populations, a stronger foundational understanding of diets used during critical life stages is needed, particularly in the northern extent of the species' range, which supports approximately two-thirds of the extant population (Garton et al. 2016, McDonald et al. 2016). Therefore, we used DNA metabarcoding of Lesser Prairie-Chicken fecal samples to quantify arthropod and plant taxa consumed by Lesser Prairie-Chickens during the brooding period and winter. We further used vegetation and arthropod survey data collected among 4 study sites in Kansas and Colorado, USA, to verify results.
The study area encompassed the northern extent of the Lesser Prairie-Chicken's distribution in Kansas and Colorado and included 4 study sites spread among the Mixed-Grass Prairie (Red Hills, Clark), Short-Grass Prairie/CRP Mosaic (Northwest), and Sand Sagebrush Prairie (Colorado, Clark) ecoregions (McDonald et al. 2014; Figure 1). Although the Colorado study site occurred within the Sand Sagebrush Prairie ecoregion, this site was predominantly composed of Conservation Reserve Program grassland (CRP) and cropland on the border of Prowers and Baca counties. Dominant grasses, forbs, subshrubs, shrubs, mean annual precipitation, and soil texture varied among study sites (Appendix Table 5). For example, subshrubs (e.g., Gutierrezia sarothrae and Amphiachyris dracunculoides) were more abundant than forbs in northwest Kansas and more abundant than shrubs at the Red Hills study site (Appendix Table 5). Forbs were predominantly Salsola tragus and Kochia scoparia, which were 2 of the top 3 most abundant forbs at all sites, excluding the Red Hills.
We collected fecal samples during the brooding period (May–September) and winter (November–March) from Lesser Prairie-Chickens captured at leks between early March and mid-May using walk-in funnel traps and drop nets (Haukos et al. 1990, Silvy et al. 1990). We sexed the birds on the basis of plumage coloration, length of pinnae, and tail pattern (Copelin 1963). We marked female Lesser Prairie-Chickens with either a 15 g VHF transmitter or a 22 g GPS satellite PTT transmitter. We obtained locations for each VHF-marked female 3–4 times wk−1, whereas females marked with GPS PTT transmitters accrued 8–10 locations day−1, contingent on available daily solar energy. GPS locations were recorded every 2 hr during the day, with a 6 hr gap between 2300 and 0500 hours.
During the brooding season, we collected fecal samples from marked hens and chicks (separate vials for each) during brood capture and weekly flush counts occurring within 1 hr of sunrise (2–98 days old). We classified fecal samples as either chick or adult samples on the basis of their relative size differences. During winter and early spring (December–March), we collected fecal samples (≥1 pellet) at roost sites. Fresh fecal samples that were still moist and appeared to have been dropped the previous night were placed in 20 mL vials using small plastic sampling spoons to minimize DNA contamination. Vials labeled with the date, unique bird ID, and coordinates of the collection location were stored in a freezer at field sites and at Kansas State University before being shipped frozen overnight for laboratory analyses.
We extracted Genomic DNA from fecal samples using the PowerSoil-htp 96-well Soil DNA Isolation Kit (MO BIO Laboratories, Carlsbad, California, USA). For arthropods, we amplified a fragment of the Folmer region of the cytochrome oxidase I (COI) gene using arthropod-specific primers (Bohmann et al. 2011, Zeale et al. 2011). To determine the contribution of plants to diets, a portion of the chloroplast trnL intron was PCR-amplified from each genomic DNA sample using the c and h trnL primers (Taberlet et al. 2007), but modified to include appropriate barcodes and adapter sequences for Illumina multiplexed sequencing. The barcodes used were 12 base pair (bp) error-correcting barcodes unique to each sample (Caporaso et al. 2012). Each 25 μL PCR reaction was mixed according to PCR Master Mix specifications (Promega, Madison, Wisconsin, USA), with 2 μL of genomic DNA template. For trnL, the thermocycling program used an initial step at 94°C for 1 min, a final extension at 72°C for 2 min, and the following steps cycled 36 times: 1 min at 94°C, 30 s at 55°C, and 30 s at 72°C. For COI, the thermocycling program used an initial step at 94°C for 5 min, a final extension at 72°C for 10 min, and the following steps cycled 45 times: 30 s at 94°C, 45 s at 45°C, and 45 s at 72°C. We cleaned amplicons from each sample and normalized them using SequalPrep Normalization Plates (Life Technologies, Carlsbad, California) before pooling them for sequencing on a MiSeq (Illumina, San Diego, California) running the 2 × 150 bp chemistry.
Assignment of Reads to Arthropod Genera
For COI reads indicating arthropod taxa, we demultiplexed sequences using “prep_fastq_for_uparse.py” (Leff 2018). Read 2s were used for downstream analysis, due to higher-quality scores. Sequences were filtered and operational taxonomic unit (OTU) picking was performed using the UPARSE pipeline (USEARCH 7). Quality filtering included trimming sequences to the expected amplicon length (158 bp—only for 250 bp reads), filtering by quality score (maxee value of 1.5), removing sequences below the minimum expected amplicon length (90 bp), and removing singletons. We clustered sequences de novo at 99% similarity for OTU picking. We performed taxonomy assignment in QIIME, using the hierarchical naive Bayesian classifer RDP, retrained with a custom reference database curated from the Barcode of Life Database (version 3). Taxonomy was assigned at 99% similarity, with a 50% confidence threshold. We further filtered sequences to remove non-arthropod sequences by removing sequences that were not resolved to at least the family level. All samples with <10 COI reads were excluded from analysis for arthropods in diet.
We calculated the percentages of all sequences assigned to a given OTU for each sample. This is referred to as RRA (relative read abundance; Kartzinel et al. 2015). For COI, an average of 9.67% of all sequences were matched to genera in the order Diptera, almost exclusively during summer. Due to observations of contact between fecal material and dipterans, we assumed that dipteran DNA entered fecal material through secondary contact after defecation and before collection. Therefore, we excluded all dipteran reads from analyses. We limited assignment of OTU to genera present among all study sites as estimated from arthropod sweep-net survey (see details below).
We constrained assignments to taxa available for consumption in western Kansas and eastern Colorado. We used sweep-net surveys at brood locations from May to August in 2013 and 2014 to sample available arthropod prey. Sweep netting is an efficient method for sampling a wide array of invertebrate species (Yi et al. 2012). However, sweep netting can be biased toward capture of Araneae, Orthoptera, Lepidoptera, and Thysanoptera (Doxon et al. 2011, Spafford and Lortie 2013). Therefore, we didn't compare biomass estimates from sweep-net surveys directly to items detected in diet using a resource-selection type analysis. Instead, we restricted DNA metabarcoding assignments to taxa detected among all sites including genera within Orthoptera, families within Hemiptera, families detected within Coleoptera, families within Araneae, and all other taxa to the order resolution.
To perform sweep-net surveys, three 100-sweep surveys were conducted at sites where fecal samples were collected and at nearby paired random locations. Survey sweeps moved north to south, passing along 3 parallel transects 10 m apart, with the center transect passing directly through the bird location (Hagen et al. 2005). We compared cumulative biomass (g) of arthropod orders (broader taxonomic resolution) at study sites to help explain relative differences in diets among sites.
Spatial and temporal influence on the consumption of arthropods
After RRA was estimated for all arthropod (COI) reads indicative of potential foods available in the study area, we summed genus-specific RRA to estimate RRA at the order level. Using RRA, we documented the relative contribution of all orders to Lesser Prairie-Chicken diets during the brood-rearing period and winter, and then assessed orders as dependent variables in separate beta regression model sets.
We used a regression based on a parameterization of the beta distribution to examine differences in RRA for orders that were predominant in fecal samples. We evaluated the relationships of RRA values among independent variables including period (brooding period and winter), chick (yes or no) during the brood-rearing period, and study sites (Northwest, Red Hills, Clark, and Colorado; Ferrari and Cribari-Neto 2004). We developed box plots to depict the median, first, and third quartiles, and maximum and minimum values of RRA for the 4 predominantly consumed orders at each site. After screening for differences among period, site, and age class, we used a multimodel inference approach to examine how spatially and temporally related covariates influenced the composition of arthropods in the diet during the brood-rearing and winter periods, separately. We examined periods separately because of the differences in available foods based on phenology and because Lesser Prairie-Chickens use a greater abundance of arthropods in the brood-rearing period than in winter, regardless of the composition of arthropods consumed (Jones 1963).
Spatial covariates were based on the location of the fecal sample and included binary covariates (occurred in cover type = 1, otherwise = 0) for native grassland, CRP, and cropland. Also included in the model set was land cover type as a categorical covariate with multiple levels, including native grassland, CRP, and cropland as separate factors and a study-site model with multiple levels (Northwest, Red Hills, Clark, and Colorado). “Native grassland” refers to grasslands occurring on soil never previously tilled and that were typically maintained for cattle production (but note that all CRP grasslands assessed were planted with native grasses and forbs). Temporally related covariates included day since start of period, chick age in days, and age class during the brood-rearing period (adult, juvenile). Day since start of period was set sequentially from 1, as the earliest date of bird use for a fecal sample collected, to the latest date of bird use for collected fecal samples in a period (brood rearing and winter). We conducted regression and performed multimodel inference using the packages “betareg” (Zeileis et al. 2016) and “AICmodavg” (Mazerolle 2016) in R (R Development Core Team 2016).
After fitting beta distribution regression models, we screened for period, age, and site effects based on informative beta coefficients. Beta coefficients were considered informative, or statistically meaningful, if not overlapping zero at the 85% confidence interval (CI; Arnold 2010). For multimodel inference, we ranked and selected the most parsimonious model based on Akaike's Information Criterion corrected for small sample sizes (AICc), for the 3 most abundant orders based on RRA. Models with ΔAICc ≤ 2 were considered equal in parsimony (Burnham and Anderson 2002, Arnold 2010).
Assignment of Reads to Plant Taxa and Functional Groups
Sequences were demultiplexed for trnL using a Python script (available from https://github.com/leffj/helper-code-for-uparse/blob/master/prep_fastq_for_uparse_paired.py). Paired end reads were then merged using “fastq_merge” pairs (Edgar 2010). We used “fastx_clipper” to trim primer and adaptor regions from both ends ( https://github.com/agordon/fastx_toolkit) because merged reads often extended beyond the amplicon region of the sequencing construct. Sequences lacking a primer region on both ends of the merged reads were discarded. Sequences were quality trimmed to have a maximum expected number of errors per read of <0.1, and only sequences with >3 identical replicates were included in downstream analyses. BLASTN 2.2.30+ was run locally, with a representative sequence for each OTU as the query and the current National Center for Biotechnology Information (NCBI) nucleotide and taxonomy database as the reference. The tabular BLAST hit tables for each OTU representative were then parsed so that only hits with >97% query coverage and identity were kept, using the “usearch7” approach (Edgar 2013, Craine et al. 2015). The NCBI genus names associated with each hit were used to populate the OTU taxonomy assignment lists. All samples with <50 trnL reads were excluded from analyses of trnL RRA (Kartzinel et al. 2015). We estimated OTU-specific RRA and defined a representative genus for each OTU to describe composition in diet. We used the representative genera when summarizing OTU composition in diets. For example, OTUs were from species in genera similar to Ambrosia. We limited plant genera within OTU to those detected during extensive vegetation surveys among sites (Appendix A).
For trnL, an average of 4% of sequences was from Pinus (range: 0–51%). Because of the unlikelihood of Pinus biomass being consumed and the presence of Pinus DNA in the blanks, the one OTU that matched with Pinus species was removed from the dataset. For trnL, among the top 10 OTUs, OTU 23 did not match at 97% levels for coverage and identity for any species in the NCBI database. However, OTU 23 matched at 100% coverage and 95% identity with a Chenopodium species in the NCBI database and was considered a species similar to Chenopodium for the purposes of this study.
Functional group assignments
Because OTUs often encompassed multiple genera, we grouped RRA from different plant genera into functional groups including forbs, shrubs, subshrubs (mostly Gutierrezia), legumes, grasses, crops (not including alfalfa), and alfalfa. Placing genera into each functional group presented challenges because the OTUs frequently encompassed genera indicative of multiple groups (see below). However, linking plant foods consumed to specific functional groups was necessary to allow for comparisons among sites and to make direct connections to the utility of landscapes with an agricultural component. In some instances, OTUs that included genera related to both grass and crop as well as shrub and subshrub functional groups included repeat values and, therefore, added values could surpass 100%. For example, 17 of 33 OTUs that identified either grass or crop foods included both crop and native grass genera (e.g., Triticum and Elymus); 2 of 45 OTUs of genera including shrub, subshrub, and forb species included representatives of >1 functional group (e.g., Artemisia and Ambrosia); and 1 of 5 OTUs for genera of legumes included both cultivated and native species (e.g., Medicago and Vicia). To overcome functional-group overlap within OTUs, we constrained the use of crop and shrub foods to instances when each land cover type occurred within 48 hr home ranges; and we used the Bayesian approach, similar to regional assignments in Royle and Rubenstein (2004), to estimate RRA for each functional groups using identity values as a prior probability:
We estimated an adjusted RRA for each functional group (RRAfg = k) by estimating the average identity value (Ig) among genera within an OTU and then dividing Ig by the sum of identity values for functional groups within each OTU. We then multiplied the quotient by the RRA estimated for each OTU (RRAOTU = j). The adjusted RRA accounts for the probability that each read is from a particular functional group based on the identity value. The identity value is a measure of the match between the OTU detected in the fecal sample and genus-specific reference sequence.
To limit plant forage possibilities to those available, and to minimize the overlap of certain OTUs encompassing multiple functional groups, we combined DNA metabarcoding inference with telemetry and extensive plant survey data. We limited native plant food availability to those genera detected during point-step transects among all study sites (Appendix A). At each study site, patches were delineated and digitized in ArcGIS 10.2 using aerial imagery from the Bing aerial basemap layer (product of ESRI, i-cubed, USDA FSA, USGS, AEX, GeoEye, Getmapping, Aerogrid, IGP) or the National Agriculture Imagery Program (NAIP) 2012 satellite imagery. Patches were identified as areas of homogeneous vegetation >2 ha in size, placed in categories (e.g., grassland, lowland, or CRP), and confirmed upon ground truthing. Within each patch, three 250 m point-step transects were conducted. Each point-step transect involved identifying the plant species for each pace (Evans and Love 1957). All delineated patches were surveyed during summer for each study site, and 20% of patches using a stratified random sample approach were surveyed during fall and winter.
To minimize overlap of certain OTUs that included multiple functional groups, we created home ranges encompassing the previous 48 hr period visited by each individual and identified the presence–absence of crop or shrub functional groups. We used minimum convex polygons for GPS-marked and buffered VHF-marked bird locations in ArcGIS 10.2 by maximum moved distance by GPS-marked birds during the 48 hr period. We used maximum distances to buffer sampled locations for VHF birds during each season. We excluded dispersing birds with straight-line movements >5 km from analyses. A 48 hr home range was used because it should encompass the spatiotemporal foraging extent incorporated into the fresh fecal sample. The 48 hr home interval encompassed a 9.9 hr fluid retention in Rock Ptarmigan (Lagopus muta), while providing enough locations to include foraging locations (Stevens and Hume 1998). We used occurrence of cultivated foods (row crops, alfalfa) and shrubs within an individual's home range to determine whether a bird had access to cultivated foods. We excluded cultivated crops as potential food items if there was no cropland in the 48 hr home range. After accounting for the availability of crop and shrub foods to each individual, we adjusted RRA to reflect availability by adding, or removing, functional-group possibilities. All home ranges included CRP or native grassland; therefore, forbs and grasses were included as possibilities for all individuals.
Spatial and temporal influence on the consumption of plants
After RRA was estimated for all plant functional groups (e.g., forbs, shrubs, subshrubs, legumes, grasses, and crops), we focused on univariate variation of specific functional groups among spatial and temporal independent covariates. Similar to methods described above, we used the package “betareg” in R to examine differences between periods (brooding period and winter) and among study sites (Northwest, Red Hills, Clark, and Colorado; Ferrari and Cribari-Neto 2004). Then we used a multimodel inference approach to test how differences in spatially and temporally related covariates influenced the composition of functional groups in the diet during the brood-rearing and winter periods separately.
We used the same spatially related covariates as we did for arthropods, including CRP, native grassland, crop, alfalfa, and land cover type. Temporally related covariates included day since start of period and the quadratic effect of day since start of period. We expected that the composition of functional plant groups may change later in the brood-rearing period and that plant composition of winter diets may change because only the most persistent shrub- and crop-based foods remain available during the coldest portions of winter. We followed the same multimodel inference protocol based on AICc and informative coefficients of beta regression models (85% CI) described above for arthropods (Burnham and Anderson 2002, Arnold 2010, Mazerolle 2016, Zeileis et al. 2016).
Evaluation of Sampled Taxonomic Richness
To examine whether sample sizes were sufficient to detect all arthropod and plant foods used by Lesser Prairie-Chickens at each study site, we used species accumulation curves depicting the relationship between number of OTUs and number of fecal samples. Species accumulation curves were generated in the R package “vegan” with the “specaccum” function, and the “Lomolino” function was used to describe the curves (Oksanen et al. 2015). From the function, we estimated an asymptote and the number of OTUs achieving a midpoint of the asymptote. We also estimated extrapolated species richness using the function “poolaccum” within package “vegan” following Chao (1987).
We collected a total of 314 fecal samples from Lesser Prairie-Chickens during the brood-rearing period (n = 211) and winter (n = 103) of 2014–2015. The number of samples collected varied by site and season (Table 1). Among all sites and seasons, arthropod DNA were obtained from 96 of the 314 samples, and readable plant DNA was sequenced in 152 of the 314 samples. A total of 334 plant and arthropod OTUs (unique DNA groupings) were identified among all fecal samples. Among the 80 samples that produced ≥10 COI sequences, there was an average of 376 sequences per sample. An average of 4,591 sequences per sample were present among the 150 samples that produced ≥50 trnL sequences (plant DNA). During the brood-rearing period, 6% (4) of the 48 hr home ranges included CRP, 22% (15) included cropland, and 72% (48) included native grassland. Of the winter 48 hr home ranges, 15% (21) included CRP, 27% (38) included cropland, and 57% (79) included native grassland.
Number of collected fecal samples and those with readable plant or animal DNA (in parentheses) at each study site in the northern portion of the Lesser Prairie-Chicken range in Kansas (KS) and Colorado, USA, during the brooding period and winter 2014–2015.
A total of 75 arthropod OTUs were identified in diets of Lesser Prairie-Chickens using COI analyses. Results from OTUs encompassed 4 classes: Insecta (63), Arachnida (9), Collembola (1), and Malacostraca (1). Among these 4 classes, 12 orders and 50 families were represented. Twenty-eight of the genera were Lepidoptera, 7 Araneae, and 6 Hemiptera (Appendix Table 6). On average, 35% of the RRA was from Lepidoptera, 26% from Orthoptera, 14% from Araneae, and 13% from Hemiptera (Appendix Figure 8 and Appendix Table 7).
Sweep-net transects indicated that arthropod communities varied among study sites. Orthoptera had the greatest percent biomass among taxa at each site (Clark = 90.2%, Red Hills = 71.5%, Northwest = 73.1%, and Colorado = 46.5%), followed by Lepidoptera, Phasmatodea, and Coleoptera (Appendix Figure 9). Lepidopterans comprised >4 times more of the arthropod community biomass in Northwest and Colorado sites than in the Red Hills site and 1.6 times more than in the Clark study site.
Beta regressions suggested no differences among Lepidoptera, Orthoptera, Hemiptera, and Araneae composition in diets between the brooding period and winter (winter β = 0.054 ± 0.303, 0.269 ± 0.293, 0.210 ± 0.265, −0.265 ± 0.279, respectively; brooding period as reference intercept). However, average reads per sample were fewer in the winter than in the brooding period for all sites except Colorado (Appendix Table 7). Given our sample size, the power of detecting a difference at an 85% CI was 0.24, 0.43, 0.47, and 0.56, respectively. Chick and adult diets during the brood-rearing period did not differ in consumption of Lepidoptera, Orthoptera, Araneae, and Hemiptera (chick β = 0.013 ± 0.403, 0.205 ± 0.386, 0.122 ± 0.388, −0.199 ± 0.370, respectively). Beta regressions also indicated no differential consumption of foods by age for Lepidoptera, Orthoptera, Araneae, and Hemiptera (age of chick days β = −0.004 ± 0.00779, 0.00732 ± 0.00788, −0.000999 ± 0.007839, −0.00218 ± 0.00700, respectively). There was an indication of more complicated nonlinear trends in the consumption of Lepidoptera and Orthoptera with minimal use of Lepidoptera after 40 days of age and greater consumption of Orthoptera when chicks surpassed 40 days of age (Figure 2).
The lack of variation among periods and ages is further indicated by stronger model support for land cover (Native Prairie, CRP, cropland) and site-based covariates for Lepidoptera and Araneae, which suggest that variation in arthropod diet consumption is more influenced by landscape characteristics than by temporal factors (Table 2). For Orthoptera during brood rearing, the model including date as a covariate was ranked highest but was equally parsimonious (ΔAICc < 2) with the native grassland, crop, and CRP models, and its beta coefficient overlapped zero at the 85% CI (Table 2). The combined effect of spatially related covariates in predicting the composition of each order during both brood rearing and winter carried an average model weight of 72% (Tables 2 and 3).
Results of beta regression model for the consumption of Lepidoptera, Orthoptera, and Araneae by Lesser Prairie-Chickens in Kansas and Colorado, USA, during the brood-rearing period (June–September) of 2014. K is the number of parameters, AICc is Akaike's Information Criterion adjusted for small sample size, ΔAICc is the difference in AICc compared to the smallest value, and wi is model weight. Models with beta coefficients not overlapping zero at the 85% confidence interval are in bold.
Beta regression model results for the consumption of Lepidoptera, Orthoptera, and Hymenoptera by Lesser Prairie-Chickens in Kansas and Colorado, USA, during winter 2014–2015. K is the number of parameters, AICc is Akaike's Information Criterion adjusted for small sample size, ΔAICc is the difference in AICc compared to the smallest value, and wi is model weight. Models with beta coefficients not overlapping zero at the 85% confidence interval are in bold.
Spatial variation in dietary composition was indicated by RRA among sites (Figure 3). During the brood-rearing period, presence of native grassland had the greatest influence on arthropod diet composition among Lepidoptera, Orthoptera, and Araneae but carried, on average, 30% of model weight (Table 2), which suggests that several variables were likely influential. The contribution of Lepidoptera in diets during the brood-rearing period decreased in native grassland (native grassland β = −0.657 ± 0.405; Table 2). Consumption of lepidopterans was 2.12× less in native grassland in comparison to cropland (23.2 ± 6.00% vs. 49.2 ± 11.8%; Figure 4). Similarly, the categorical native grassland covariate was the best predictor of the consumption of Araneae, based on AICc, and the beta coefficient did not overlap zero at the 85% CI (native grassland β = 0.559 ± 0.379). Araneae contributed 65× more to diets in native grassland than in other cover types and was rarely consumed in cropland (26.2 ± 7.02% vs. 0.04 ± 0.004%; Figure 4). For Orthoptera, the model including native grassland as a covariate was not informative (native grassland β = 0.154 ± 0.361). Despite not providing a statistically meaningful difference, point estimates for Orthoptera RRA was 21.7 ± 6.50% in native grassland vs. 12.7 ± 6.71% in other cover types. Hemiptera contributed relatively equally to diets among Lesser Prairie-Chickens using CRP grassland, native grassland, and cropland (Figure 4).
In winter, Lepidoptera, Orthoptera, and Hymenoptera (most likely galls) contributed most to arthropod-based food for Lesser Prairie-Chickens (Appendix Figure 8 and Appendix Table 7). Of the top 4 orders contributing to winter diets, Orthoptera was the only order that changed (decreased) as the winter progressed, which was significant at the 85% CI (day since start of period β = −0.035 ± 0.0131). Among sites, Clark birds had the greatest percentages of Orthoptera in their winter diet when compared to all other sites, and this was significant at the 85% CI (51.7 ± 12.6% in Clark vs. 18.3 ± 7.7% in Colorado vs. 0% in Red Hills and Northwest; Clark β = 1.86 ± 0.613).
Metabarcoding of fecal samples indicated that Lesser Prairie-Chickens consumed foods encompassing 2 classes (Magnoliopsida and Liliopsida), 19 orders (predominantly Asterales, Poales, and Fabales), 30 families, and 90 genera. A total of 235 OTUs were found to represent ≥1% of the plant diet for a given bird at a given time. In contrast to the assignment of OTU to specific arthropod taxa, trnL OTUs were not genus specific and, on average, comprised 4.15 ± 4.79 genera, ranging from 1 to 28 potential genera that were present at all study sites combined. Of the 235 recorded OTUs, 70 represented ≥10% of the diet for ≥1 of the samples. The most abundant OTUs were from species in genera similar to Ambrosia (27% OTU-specific RRA of all reads), followed by species in genera similar to Lactuca or Taraxacum (10%), Medicago (6%), and Triticum (5%).
For the brood-rearing period, the 10 most abundant OTUs included species similar to Ambrosia (16.2%), Lactuca (8.5%), Triticum (5.5%), Chenopodium (4.3%), Physalis (3.9%), Commelina (3.1%), Trifolium (1.8%), and Elymus (1.4%). Ambrosia and Triticum were represented by 2 separate OTUs as part of the top 10 most abundant summer OTU foods. During winter, the 10 most abundant OTUs consumed included species similar to Ambrosia (21.0%), Lactuca (5.6%), Medicago (4.8%), Triticum (4.4%), Bromus (1.1%), Oenothera (0.9%), Elymus (0.7%), Sorghum (0.6%), and Chenopodium (0.6%). Triticum was represented by 2 separate OTUs as part of the top 10 most abundant winter OTUs.
Home ranges (48 hr) averaged 45.06 ± 44.50 ha during the nonbreeding season and 11.17 ± 8.84 ha during brood rearing for GPS-marked birds. We then used the maximum-size home ranges of nondispersing GPS-marked individuals during each time period to estimate home ranges for VHF-marked Lesser Prairie-Chickens. Home ranges for VHF birds were derived from the higher-resolution GPS-marked bird data because GPS locations were obtained frequently enough to generate 48 hr home ranges. Maximum home range sizes during the nonbreeding and brooding periods were 191.52 ha and 32.83 ha, respectively, from which we derived 781 m and 323 m buffer distances around VHF fecal collection locations to account for all potentially used food sources.
In both the brood-rearing and winter periods, forbs were the predominant plant-based food source (winter 53.7 ± 3.7%, brooding 60.67 ± 5.5%; Appendix Figure 10). Differences in the overall use of functional groups among the winter and brood-rearing periods were minimal. However, subshrubs (e.g., Gutierrezia spp.) and grasses contributed 1.5 times (43.4 ± 3.7% vs. 29.8 ± 5.7%) more to Lesser Prairie-Chicken diets during winter than during brood rearing (winter β = 0.564 ± 0.220, 0.287 ± 0.195). By contrast, there was no difference in the consumption of forbs, legumes, shrubs, and crops between periods (brooding β = 0.198 ± 0.230, −0180 ± 0.209, 0.222 ± 0.175, −0.265 ± 0.185, respectively).
We assessed differences among all sites separately for each period. Within the brood-rearing period alone, foods in the forb, grass, and legume functional groups did not differ among sites. Shrub- and subshrub-based foods contributed more to diets during the brood-rearing period in the Red Hills and northwest Kansas compared to Clark and Colorado (Red Hills β = 1.82 ± 0.782, Northwest β = 0.769 ± 0.430, Clark β = 1.22 ± 0.779, Colorado β = 0.836 ± 0.444). Crop-based foods provided a greater contribution to brood-rearing diets in Colorado compared to other sites (β = 3.67 ± 0.509).
During winter, grass composition in diets varied among sites. More grasses were consumed during winter at the Northwest study than at the Clark study site (23.0 ± 2.6% vs. 11.0 ± 1.7%; β = 0.855 ± 0.289; Figure 5). Shrub foods contributed more in winter at the Red Hills study site than at Clark (β = 0.908 ± 0.391). Crop foods contributed more in winter to diets at the Northwest site than at Clark (β = 0.443 ± 0.288). Last, subshrub foods contributed more in winter to diets at the Northwest and Red Hills study sites than at Clark (β = 0.836 ± 0.445, 1.22 ± 0.779, respectively; Figure 5).
After screening for differences among periods and sites, we focused on winter diets, using a multimodel inference approach, because Lesser Prairie-Chickens predominantly consume plant material during winter (Jones 1963). Models including spatially related covariates carried, on average, 99% of model weight (AICc weight; Table 4). The top-ranking predictor for forb diet composition was occurrence in alfalfa and crop fields (Table 4). Forbs were consumed less in winter by Lesser Prairie-Chickens using alfalfa fields and crop fields in general (β = −1.57 ± 0.467; identical beta values for alfalfa and crop). Forbs were more readily consumed in native grassland and CRP (Figure 6). The proportion of grass in diets was best predicted by site (Table 4; see differences above), with use of native grassland ranking second among models (native grassland β = 0.386 ± 0.238). Birds using alfalfa and crop fields had the greatest relative proportion of legumes in their diet (β = 4.60 ± 0.507 for both alfalfa and crop). All fecal samples collected in cropland were collected in cultivated alfalfa, which confirms that birds can use alfalfa fields in winter as a food source. Shrubs contributed more to the diets of Lesser Prairie-Chickens using native grassland than to those in other cover types (native grassland β = 1.55 ± 0.254; Table 4). The relative diet composition of subshrub appears to be most strongly influenced by use of crop fields, with consumption of subshrub lower in cropland (β = −1.38 ± 0.454).
Akaike's Information Criterion adjusted for small sample size (AICc), difference in AICc compared to the smallest value (ΔAICc), and model weight (wi) for beta regression models explaining winter plant diets of Lesser Prairie-Chickens in Kansas and Colorado, USA, 2013–2014. K is the number of parameters. Models with beta coefficients not overlapping zero at the 85% confidence interval are in bold.
Evaluation of Sampled Taxonomic Richness
Among all sites, the arthropod species accumulation curve achieved an estimated asymptote at 156 OTUs, which suggests that we didn't sample all available forage; the midpoint for achieving an asymptote was estimated at 105 fecal samples (Figure 7). The extrapolated species richness at the OTU level (based on Chao 1987) was 101. The plant species accumulation curve achieved an estimated asymptote at 282 OTUs, which suggests that we sampled nearly all used plant forage at the OTU level. The midpoint for achieving the asymptote was estimated at 17 fecal samples (Figure 7). The extrapolated species richness at the OTU level (based on Chao 1987) was 262.
Using a combination of tools including DNA metabarcoding of fecal samples, telemetry data, and local plant and arthropod surveys, we identified foods consumed by Lesser Prairie-Chickens among 4 study sites. Lesser Prairie-Chickens that used native grassland maintained for cattle production consumed a greater diversity of arthropods and plant functional groups. In 48 hr home ranges that had a row-crop agriculture component, Lesser Prairie-Chickens largely used alfalfa when it was available during winter. Females and chicks, unexpectedly, preyed mostly on lepidopteran foods (likely larvae) during brood rearing. The use of shrub-based foods varied among sites but is likely not as important as in other regions (e.g., sand shinnery oak prairie) or in other grouse species (Schmidt 1936, Moss 1983, Olawsky 1987).
Arthropods in Lesser Prairie-Chicken Diets
The greater consumption of Lepidoptera in this study than was found in past research is likely a product of both the limited detection of soft-bodied prey using traditional methods and inclusion of study sites that have a strong row-crop agriculture component. Lesser Prairie-Chickens are known to consume lepidopteran larvae, yet the results of previous research suggest minimal consumption of Lepidoptera in comparison to Orthoptera (Davis et al. 1980). The traditional use of fecal dissection may not be effective in detecting lepidopteran larvae (e.g., butterfly and moth caterpillars). No study using fecal dissection identified Lepidoptera as a prey item for Lesser Prairie-Chickens (Jones 1963, Doerr and Guthery 1983). Only studies that examined crop contents have reported consumption of lepidopteran larvae (Crawford and Bolen 1976, Suminski 1977, Smith 1979, Davis et al. 1980, Riley et al. 1993). However, not all studies examining crop contents have explicitly identified Lepidoptera as a food item, and foods from this order may be clumped as “other insects” (Olawsky 1987), making comparisons among other studies challenging. Overall, the soft-bodied nature of caterpillars likely makes them easier to digest and subsequently harder to detect using traditional dissection approaches (Trevelline et al. 2016). DNA metabarcoding may be the least biased tool for comparing dietary composition among soft- and hard-bodied prey.
In addition to palatability, use of lepidopteran larvae during the brood-rearing period may be related to the ease of capture by a small, 13–35 g chick. Lepidopteran larvae would be easy for Lesser Prairie-Chicken chicks to obtain when occurring within reach on the ground or in shorter vegetation. It is possible that soft-bodied larvae from other orders (e.g., Coleoptera) could also be consumed when available. Although we didn't expect a greater consumption of Lepidoptera than of Orthoptera by Lesser Prairie-Chicken chicks, we predicted that chicks would be restricted to smaller arthropod prey of limited mobility (following optimal diet theory; Suminski 1977, Sih and Christensen 2001). The use of lepidopteran larvae by Lesser Prairie-Chicken chicks supports this prediction. The potential dietary selection of lepidopterans further identifies the necessity of matching life histories among predator and prey. The life-history strategies of arthropod species may largely determine their importance as a prey item.
Although Lepidoptera were used as a food source among all land cover types and sites, specific lepidopteran genera were used in agricultural landscapes. Diets of Lesser Prairie-Chickens during the brooding period were largely supported by the genera Euxoa and Dargida. These 2 genera comprise several known agricultural pest species, including army cutworms (Euxoa auxiliaris). Dietary use of cutworms by Lesser Prairie-Chickens was also detected in fall by Crawford and Bolen (1976) in fragmented sand shinnery oak prairie. Consumption of agricultural pests provides evidence of one ecological service provided by Lesser Prairie-Chickens that could be used to gain conservation support in private working landscapes throughout their distribution (Wenny et al. 2011).
In contrast to the prevalent consumption of Lepidoptera in their northern range, the predominant use of orthopteran foods by Lesser Prairie-Chickens is well supported by other published research (Jones 1964, Suminski 1977, Davis et al. 1980, Doerr and Guthery 1983). The difference in predominant foods (Orthoptera vs. Lepidoptera) may be a result of spatial variation among study areas, in addition to potential biases in detecting soft-bodied prey using traditional methods. Even within the present study, we detected substantial variation in diets among study sites. The greater consumption of orthopterans at the Clark study site could be driven by the limited availability of lepidopterans and an increased abundance of grasshoppers in the genus Melanoplus at the Clark site (Appendix Figure 9; D. A. Haukos et al. personal observation). Melanoplus was the main genus of orthopterans used as a food across all sites. At the Clark study site, Melanoplus sanguinipes was substantially more abundant, and the roosting and morning basking of this species on bare ground may make it an easily obtainable prey item for Lesser Prairie-Chickens (Pfadt 1994, D. A. Haukos et al. personal observation).
The similar consumption of Orthoptera by Lesser Prairie-Chickens using grassland compared to cropland or CRP also doesn't provide any indication of difference in use of Lepidoptera vs. Orthoptera in grassland. Although Orthoptera composition was greatest in grassland, the RRA of Orthoptera was nearly identical to that of Lepidoptera in native grassland. Because RRA data are proportional among arthropod orders, an estimate close to 25% (split among 4 main foods) within one cover type would suggest that individuals using that cover type have more diverse diets. Although the split among the 4 orders was not perfectly uniform, Lesser Prairie-Chickens that used native grassland consumed a more diverse arthropod diet, which contrasts with our hypothesis that Lesser Prairie-Chickens would specialize on Orthopteran prey. Lesser Prairie-Chicken broods using native grassland may be opportunistic predators when diets are assessed during 0–90 days of age (Davis et al. 1980).
Despite the fact that brood diets appeared to be opportunistic when examining the brooding period as a whole, there was some indication of a nonlinear transition from Lepidoptera- to Orthoptera-dominated diets as chicks surpassed ∼40 days of age. We were particularly interested in diet during the first few weeks of a Lesser Prairie-Chicken's life. Knowledge of factors influencing survival during the first 21 days can be crucial for understanding what drives overall population growth rates (Hagen et al. 2009, McNew et al. 2012, Lautenbach 2015). The finite rate of population growth (λ) among prairie grouse and other galliformes has consistently been shown to be sensitive to variation in the 0- to 21-day-old survival bottleneck (Tympanuchus spp.; Wisdom and Mills 1997, Sandercock et al. 2008, Hagen et al. 2009, McNew et al. 2012, Taylor et al. 2012). Food availability may be particularly important for survival through this life stage, as indicated by strong variation in the mass of chicks and by observations of dead, undepredated chicks that may have died from starvation or thermal stress (Lautenbach 2015). Knowledge of the effects of food availability on chick survival is largely limited to inference from a closely related species within the subfamily Tetraoninane, the Greater Sage-Grouse (Centrocercus urophasianus). Sage-grouse chick survival can increase with the availability of Lepidoptera, slender phlox (Phlox gracilis), and total forb cover (Gregg and Crawford 2009). The influence of food availability on chick survival may contrast with the remainder of a grouse's life when there is strong support that predation poses the greater survival risk (Bergerud and Gratson 1988). However, if food availability drives passage through the most influential life stage and survival bottleneck, even if only lasting up to 21 days (the first 7 days may be most influential; Lautenbach 2015), the influence of food availability may be paramount and materialize in population level trajectories at much broader scales.
Comparative Nutrient Values of Lepidopterans and Orthopterans
Lepidopteran and orthopteran foods both provide greater concentrations of protein than any plant-based foods at the nutrient level (Lassiter and Edwards 1982, Savory 1989, Rumpold and Schlüter 2013). Protein in arthropod foods are also likely more digestible than that in plants (Stiven 1961, Savory 1989). On average, orthopterans can provide a food source that is 61% protein and 13% fat, whereas lepidopterans are 45% protein and 27% fat (Sugimura et al. 1984). Among protein estimates, there is interspecific variation and differences in digestibility. Furthermore, assimilation of protein from chitin-rich orthopterans and soft-bodied lepidopterans may be similar amid differences in nutrient composition (Sugimura et al. 1984). Mineral and amino acid composition provided by the 2 families appears to be similar, with variation among prey species (Rumpold and Schlüter 2013).
The Need for Ancillary Data
The potential benefits of using DNA metabarcoding to understand diets of wildlife species are numerous, but the current utility of the method hinges on ancillary data that can be used to constrain and evaluate the completeness of reference databases. We were unable to distinguish among certain plant foods that were from grass and crop functional groups using the primers we selected. The addition of 48 hr home range data allowed for greater inference on the use of cultivated foods. Additionally, reference DNA sequences for species that did not occur at any of the field sites sometimes matched sequences in fecal samples. To avoid inaccurate predictions, we constrained possible food sources to those detected during vegetation and arthropod surveys. The amplification of plant and arthropod DNA in only a proportion of the samples may be a problem unique to Lesser Prairie-Chickens and, potentially, other grouse species. For example, DNA was successfully amplified in all fecal samples from Louisiana Waterthrush (Parkesia motacilla), in 100% of bison (Bison bison) fecal samples, and in 74% of fecal samples from bats (Bohmann et al. 2011, Craine et al. 2015, Trevelline et al. 2016).
The predominant use of forbs as a food source during both brood-rearing and winter periods highlights the need to maintain disturbance regimes that support healthy forb populations (Hagen et al. 2004). Forbs provided a critical habitat component for Lesser Prairie-Chickens as food resources, even though they often comprised <10% of the available vegetation.
We detected greater RRA of forbs during brood rearing and winter, with specific forbs showing greater use during specific periods. During the brood-rearing period, forbs consumed by Lesser Prairie-Chickens were largely from Chenopodium- and Abutilon-like species. Chenopodium album (lamb's quarters) was present at all field sites during summer. The leaves of C. album are known to be palatable and high in calcium, which may be particularly important for growing Lesser Prairie-Chicken chicks (Adedapo et al. 2011). The use of Abutilon-like species may indicate consumption of Callirhoe involucrata (purple poppy mallow) or Sphaeralcea coccinea (scarlet globemallow), both of which were present at all sites and actively growing during the brood-rearing period (D. A. Haukos et al. personal observation). Leaves of S. coccinea are high in vitamin A, calcium, and protein and can be selected as food by scaled quail (Callipepla squamata; Ault et al. 1983, Arthun et al. 1992). Although documentation of C. involucrata as food for grassland birds is limited, the plant has adequate phosphorus and crude protein content (Odocoileus virginianus; Everitt and Gonzalez 1981). It also functions as a known larval host and food source for several butterflies (Fernandez-Canero and Gonzalez-Redondo 2010, Scott 2014). Observations were made of several caterpillar larvae on the receptacles of C. involucrate flowers at the Clark study site during the brooding period (D. Sullins personal observation). The presence of Abutilon-like plants in Lesser Prairie-Chicken diets could be from either direct or indirect consumption mediated through lepidopteran herbivory. The presence of arthropod foods can be attained only by first providing necessary host plants.
Outside of the brooding period, plant matter becomes particularly important in Lesser Prairie-Chicken diets during winter and spring as available forage decreases, thermoregulatory needs are maximized, and stored energy becomes particularly important with approaching lekking and nesting seasons (Haukos and Zavaleta 2016). Winter diets of grouse are often limited to only a few items that can provide sustenance—typically high in fiber, low in nutrient content, and requiring longer digestive tracts to process (Moss 1983). In the present study, the greater consumption of forbs compared to all other functional groups suggests a reliance on noncultivated foods in the northern portion of the Lesser Prairie-Chicken range. Use of forbs by Lesser Prairie-Chickens contrasts with grouse of more ancestral Arctic and boreal origins that largely consume woody vegetation during winter (Schmidt 1936, Moss 1983, DeYoung and Williford 2016) but is consistent with comparatively greater predation of “weed seeds” by pinnated grouse (e.g., Greater Prairie-Chickens [Tympanuchus cupido]) in comparison to Sharp-tailed Grouse (T. phasianellus; Schmidt 1936). Forb DNA was nearly absent from fecal samples collected in cropland, which suggests that current use of herbicides may reduce the availability of forbs in cropland.
Although forbs were dominant plant foods used by Lesser Prairie-Chickens during brood rearing and in winter, the relative importance of crops, shrubs, legumes, and subshrubs as food sources increased from brood rearing to winter. The amount of grass consumed remained the same, in contrast to the results of Jones (1963), who documented a slight increase in grasses consumed during winter. The increased use of shrubs and subshrubs may be related to the persistence of shrub- and subshrub-based foods during winter. Broom snakeweed was present at all study sites. This subshrub maintains green basal leaves longer into the fall and winter compared to other plants in the region, thus providing a persistent source of leafy green vegetation (Ralphs and Wiedmeier 2004). Broom snakeweed is a known food for Lesser Prairie-Chickens and has protein and nutrient content similar to green grass, but numerous secondary metabolite compounds make broom snakeweed challenging to digest (Jones 1963, Davis et al. 1980, Ralphs and Wiedmeier 2004). Although subshrubs such as broom snakeweed may not be easy to digest, they may provide a food source, persistent throughout the winter, for which grouse have evolved advanced digestive systems to procure nutrients, as indicated by seasonal changes in gut morphology (Olawsky 1987, Sedinger 1997, Donaldson et al. 2006).
Shrub-based foods can be important for Lesser Prairie-Chickens (Jones 1964, Crawford and Bolen 1976, Suminski 1977, Olawsky 1987, Riley et al. 1993) and other grouse (Patterson 1952, Remington and Braun 1985). Most research indicating that shrubs are important for Lesser Prairie-Chickens has focused on the use of sand shinnery oak where available in Texas and New Mexico, USA (Suminski 1977, Olawsky 1987, Riley et al. 1993). Sand sagebrush, sumac (Rhus spp.), willow (Salix spp.), and cottonwood (Populus spp.) have also provided food for Lesser Prairie-Chickens (Schwilling 1955, Jones 1963, 1964). The increased use of shrub-based foods during winter corresponded with the increased consumption of sand sagebrush from December to February in northwest Oklahoma, USA (Jones 1963).
Outside of using persistent winter foods in the form of shrubs and subshrubs, cultivated crops can be used by Lesser Prairie-Chickens (Salter et al. 2005). Use of cultivated legumes during winter was largely restricted to the Clark study site, where the OTU containing alfalfa (Medicago spp., 100% identity and coverage) was consumed 1.95× more than the next leading OTU containing Triticum-like species. Cultivated alfalfa was available at the Clark study site and was consumed by Lesser Prairie-Chickens in distinct cropland areas. The use of alfalfa cropland at this site may explain differences in space use among regions (Robinson 2015).
Diversity and Food Stability
The greater diversity of forage in native working grassland may be key to food and nutrient stability in Lesser Prairie-Chickens. Lesser Prairie-Chickens occur in a region with the greatest variability of net primary productivity in the Great Plains (Sala et al 1988, Grisham et al. 2016). In such a variable environment, population viability may be more influenced by a stable presence of foods from year to year than by an abundance at any one time. Various arthropod and plant taxa can boom or bust in response to years of above-average precipitation or drought, and therefore food stability may be linked to a diversity of forage (Haglund 1980, Tilman and Downing 1994, Gutbrodt et al. 2011, Craine et al. 2013). Our results indicated that native working grassland provided forage for Lesser Prairie-Chickens, in addition to providing cover for reproduction and adult survival (Hagen et al. 2013). However, in some landscapes it is possible that the presence of small-scale row-crop agriculture adjacent to grassland could diversify food options (Rodgers 2016).
We thank all the landowners who allowed access during the research project and all our technicians, including B. Anderson, S. Baker, S. Bard, G. Brinkman, K. Broadfoot, R. Cooper, J. Danner, J. Decker, E. Entsminger, R. Galvin, N. Gilbert, A. Godar, G. Gould, B. Hardy, S. Hoffman, D. Holt, B. Irle, T. Karish, A. Klais, H. Kruckman, K. Kuechle, S. Lane, E. Leipold, J. Letlebo, E. Mangelinckx, L. McCall, A. Nichter, K. Phillips, J. Proescholdt, J. Rabon, T. Reed, A. Rhodes, B. Ross, D. Spencer, A. Steed, A. Swicegood, P. Waldron, B. Walter, I. Waters, W. White, E. Wiens, J. Yantachka, and A. Zarazua for assisting with data collection. We greatly appreciate the logistic and technical support provided by J. Pitman, C. Hagen, J. Prendergast, K. Fricke, D. Tacha, D. Kraft, J. Ungerer, D. Krehbiel, A. Elliott, D. Dahlgren, G. Kramos, and A. Flanders. Thanks to J. Zavaleta for thought-provoking conversations related to Lesser Prairie-Chicken diets. We thank C. Conway for reviewing an earlier version of this publication.
Funding statement: Funding for the project was provided by Kansas Wildlife, Parks, and Tourism (Federal Assistance Grant KS W-73-R-3); USDA Farm Services CRP Monitoring, Assessment, and Evaluation (12-IA-MRE CRP TA#7, KSCFWRU RWO 62); and USDA Natural Resources Conservation Service, Lesser Prairie-Chicken Initiative.
Ethics statement: We complied with the Guidelines to the Use of Wild Birds in Research for all components of this research. We prepared protocols and obtained collection permits to capture and handle through the Kansas State University Institutional Animal Care and Use Committee protocol nos. 3241 and 3703; Kansas Department of Wildlife, Parks, and Tourism scientific collection permits SC-042-2013, SC-079-2014, and SC-001-2015; and Colorado Parks and Wildlife scientific collection license numbers 13TRb2053, 14TRb2053, and 15TRb2053.
Author contributions: D.S.S. and D.A.H. conceived the idea, design, and experiment. J.M.C., N.F., D.S.S., J.M.L., S.G.R., J.D.L., J.D.K., R.T.P., and J.H.R. performed the experiments (collected data and conducted the research). D.S.S., J.M.C., and N.F. developed or designed the methods. D.S.S. and J.M.C. analyzed the data. D.S.S., D.A.H., J.M.C., and B.K.S. wrote the paper.
Plant Genera Detected (n = 257) during Vegetation Surveys at Study Sites in Western Kansas and Eastern Colorado, USA, 2013–2016
APPENDIX TABLE 5.
Proportional abundance of the most common grass, forb, subshrub, and shrub species estimated from point-step transects by study site for the northern portion of the Lesser Prairie-Chicken range in Kansas and Colorado, USA, 2014–2015. Mean annual precipitation (PRISM climate group averaged from a PRISM 800 m resolution raster file) and dominant soil textures (Soil Survey Staff 2016) are also included.
APPENDIX TABLE 6.
Families and genera of arthropods detected using DNA barcoding in fecal samples of Lesser Prairie-Chickens during brood rearing and winter at 4 study sites in Kansas and Colorado, USA, 2014–2015.
APPENDIX TABLE 7.
Relative read abundance (sample size, mean, and SD) of arthropod orders in the diets of Lesser Prairie-Chicken chicks and adults during the brooding period, and of adults during winter, from 4 study sites in Kansas and Colorado, USA, 2014–2015. Only one brood sample had readable DNA from Colorado.