DOI: 10.5061/dryad.8sm55
title: Data from: Genomic single-nucleotide polymorphisms confirm that Gunnison and Greater sage-grouse are genetically well differentiated and that the Bi-State population is distinct

description: Sage-grouse are iconic, declining inhabitants of sagebrush habitats in
 western North America, and their management depends on an understanding of
 genetic variation across the landscape. Two distinct species of
 sage-grouse have been recognized, Greater (Centrocercus urophasianus) and
 Gunnison sage-grouse (C. minimus), based on morphology, behavior, and
 variation at neutral genetic markers. A parapatric group of Greater
 Sage-Grouse along the border of California and Nevada (“Bi-State”) is also
 genetically distinct at the same neutral genetic markers, yet not
 different in behavior or morphology. Because delineating taxonomic
 boundaries and defining conservation units is often difficult in recently
 diverged taxa and can be further complicated by highly skewed mating
 systems, we took advantage of new genomic methods that improve our ability
 to characterize genetic variation at a much finer resolution. We
 identified thousands of single-nucleotide polymorphisms (SNPs) among
 Gunnison, Greater, and Bi-State sage-grouse and used them to
 comprehensively examine levels of genetic diversity and differentiation
 among these groups. The pairwise multilocus fixation index (FST) was high
 (0.49) between Gunnison and Greater sage-grouse, and both principal
 coordinates analysis and model-based clustering grouped samples
 unequivocally by species. Standing genetic variation was lower within the
 Gunnison Sage-Grouse. The Bi-State population was also significantly
 differentiated from Greater Sage-Grouse, albeit more weakly (FST = 0.09),
 and genetic clustering results were consistent with reduced gene flow with
 Greater Sage-Grouse. No comparable genetic divisions were found within the
 Greater Sage-Grouse sample, which spanned the southern half of the range.
 Thus, we provide much stronger genetic evidence supporting the recognition
 of Gunnison Sage-Grouse as a distinct species with low genetic diversity.
 Further, our work confirms that the Bi-State population is differentiated
 from other Greater Sage-Grouse. The level of differentiation is much less
 than the divergence between Greater and Gunnison sage-grouse, supporting
 the idea that the Bi-State represents a unique population within the
 Greater Sage-Grouse. New genomic methods like the
 restriction-site-associated DNA (RAD-tag) method used here illustrate how
 increasing the number of markers and coverage of the genome can better
 characterize patterns of genetic variation, particularly among recently
 diverged taxa, providing vital information for conservation and
 management.
description: sage grouse
 genepopInput genotype file for
 the program Genepopstructure_input_fileInput genotype file for the program Structure. DNA bases encoded as A=1, C=2, G=3, T=4.Genalex-distance-matrix-and-PCADistance matrix input file for the program GenAlEx, and output from PCA analysis using default settings.MAF10-final-lociExcel spreadsheet of SNPs with minor allele frequency of 0.10 or more, in variant call format (vcf). See samtools manual for detailed interpretation of this file.final-MAF10-tagsFASTA sequence file of "stacks" or restriction-site-associated tags. Stacks were created with the "stacks" pipeline as described in the text.Input-to-tblastx-searchFASTA file containing tags, or assembled contigs in lieu of tags, that had TBLASTX matches to the chicken reference genome at 1e-5. Contigs resulted from the assembly of paired reads plus an additional genomic sequencing run, see text for details. The tag that was replaced by a contig in this search is given in the header line.Sage-grouse-2-chicken-tblastx-summarySpreadsheet containing match information for Sage grouse sequences that had TBLASTX matches to the chicken genomeExpected-heterozygosity-FST-between-speciesExcel spreadsheet comparing expected heterozygosity and FST at SNP loci, along with actual base counts from reads. Comparison is between the Greater Sage-grouse and Gunnison Sage-grouse samples, ie it does not include the Bi-state population for this analysis.