Jennifer A. Fike

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake.

Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct “Seq-to-SSR” approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus, and Clarks Nutcracker, Nucifraga columbiana), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as

Environmental DNA (eDNA) Sampling Improves Occurrence and Detection Estimates of Invasive Burmese Pythons

10 per sample – a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.

Landscape characteristics influencing the genetic structure of greater sage-grouse within the stronghold of their range: a holistic modeling approach.

Environmental DNA (eDNA) methods are used to detect DNA that is shed into the aquatic environment by cryptic or low density species. Applied in eDNA studies, occupancy models can be used to estimate occurrence and detection probabilities and thereby account for imperfect detection. However, occupancy terminology has been applied inconsistently in eDNA studies, and many have calculated occurrence probabilities while not considering the effects of imperfect detection. Low detection of invasive giant constrictors using visual surveys and traps has hampered the estimation of occupancy and detection estimates needed for population management in southern Florida, USA. Giant constrictor snakes pose a threat to native species and the ecological restoration of the Florida Everglades. To assist with detection, we developed species-specific eDNA assays using quantitative PCR (qPCR) for the Burmese python (Python molurus bivittatus), Northern African python (P. sebae), boa constrictor (Boa constrictor), and the green (Eunectes murinus) and yellow anaconda (E. notaeus). Burmese pythons, Northern African pythons, and boa constrictors are established and reproducing, while the green and yellow anaconda have the potential to become established. We validated the python and boa constrictor assays using laboratory trials and tested all species in 21 field locations distributed in eight southern Florida regions. Burmese python eDNA was detected in 37 of 63 field sampling events; however, the other species were not detected. Although eDNA was heterogeneously distributed in the environment, occupancy models were able to provide the first estimates of detection probabilities, which were greater than 91%. Burmese python eDNA was detected along the leading northern edge of the known population boundary. The development of informative detection tools and eDNA occupancy models can improve conservation efforts in southern Florida and support more extensive studies of invasive constrictors. Generic sampling design and terminology are proposed to standardize and clarify interpretations of eDNA-based occupancy models.

Two low coverage bird genomes and a comparison of reference-guided versus de novo genome assemblies

Given the significance of animal dispersal to population dynamics and geographic variability, understanding how dispersal is impacted by landscape patterns has major ecological and conservation importance. Speaking to the importance of dispersal, the use of linear mixed models to compare genetic differentiation with pairwise resistance derived from landscape resistance surfaces has presented new opportunities to disentangle the menagerie of factors behind effective dispersal across a given landscape. Here, we combine these approaches with novel resistance surface parameterization to determine how the distribution of high- and low-quality seasonal habitat and individual landscape components shape patterns of gene flow for the greater sage-grouse (Centrocercus urophasianus) across Wyoming. We found that pairwise resistance derived from the distribution of low-quality nesting and winter, but not summer, seasonal habitat had the strongest correlation with genetic differentiation. Although the patterns were not as strong as with habitat distribution, multivariate models with sagebrush cover and landscape ruggedness or forest cover and ruggedness similarly had a much stronger fit with genetic differentiation than an undifferentiated landscape. In most cases, landscape resistance surfaces transformed with 17.33-km-diameter moving windows were preferred, suggesting small-scale differences in habitat were unimportant at this large spatial extent. Despite the emergence of these overall patterns, there were differences in the selection of top models depending on the model selection criteria, suggesting research into the most appropriate criteria for landscape genetics is required. Overall, our results highlight the importance of differences in seasonal habitat preferences to patterns of gene flow and suggest the combination of habitat suitability modeling and linear mixed models with our resistance parameterization is a powerful approach to discerning the effects of landscape on gene flow.

Latent spatial models and sampling design for landscape genetics

As a greater number and diversity of high-quality vertebrate reference genomes become available, it is increasingly feasible to use these references to guide new draft assemblies for related species. Reference-guided assembly approaches may substantially increase the contiguity and completeness of a new genome using only low levels of genome coverage that might otherwise be insufficient for de novo genome assembly. We used low-coverage (∼3.5–5.5x) Illumina paired-end sequencing to assemble draft genomes of two bird species (the Gunnison Sage-Grouse, Centrocercus minimus, and the Clarks Nutcracker, Nucifraga columbiana). We used these data to estimate de novo genome assemblies and reference-guided assemblies, and compared the information content and completeness of these assemblies by comparing CEGMA gene set representation, repeat element content, simple sequence repeat content, and GC isochore structure among assemblies. Our results demonstrate that even lower-coverage genome sequencing projects are capable of producing informative and useful genomic resources, particularly through the use of reference-guided assemblies.

Development of 13 microsatellites for Gunnison Sage-grouse (Centrocercus minimus) using next-generation shotgun sequencing and their utility in Greater Sage-grouse (Centrocercus urophasianus)

LATENT SPATIAL MODELS AND SAMPLING DESIGN FOR LANDSCAPE GENETICS1 BY EPHRAIM M. HANKS∗, MEVIN B. HOOTEN†,‡, STEVEN T. KNICK§, SARA J. OYLER-MCCANCE¶, JENNIFER A. FIKE¶, TODD B. CROSS‖,∗∗ AND MICHAEL K. SCHWARTZ‖ Pennsylvania State University∗, U.S. Geological Survey, Colorado Cooperative Fish and Wildlife Research Unit†, Colorado State University‡, U.S. Geological Survey, Forest and Rangeland Ecosystem Science Center§, U.S. Geological Survey, Fort Collins Science Center¶, U.S. Forest Service, Rocky Mountain Research Station‖ and University of Montana∗∗

Development and characterization of thirteen microsatellite loci in Clark's nutcracker (Nucifraga columbiana)

Gunnison Sage-grouse are an obligate sagebrush species that has experienced significant population declines and has been proposed for listing under the U.S. Endangered Species Act. In order to examine levels of connectivity among Gunnison Sage-grouse leks, we identified 13 novel microsatellite loci though next-generation shotgun sequencing, and tested them on the closely related Greater Sage-grouse. The number of alleles per locus ranged from 2 to 12. No loci were found to be linked, although 2 loci revealed significant departures from Hardy–Weinberg equilibrium or evidence of null alleles. While these microsatellites were designed for Gunnison Sage-grouse, they also work well for Greater Sage-grouse and could be used for numerous genetic questions including landscape and population genetics.

Genetic Mark–Recapture Improves Estimates of Maternity Colony Size for Indiana Bats

Clark’s nutcrackers are important seed dispersers for two widely-distributed western North American conifers, whitebark pine and limber pine, which are declining due to outbreaks of mountain pine beetle and white pine blister rust. Because nutcracker seed dispersal services are key to maintaining viable populations of these imperiled pines, knowledge of movement patterns of Clark’s nutcrackers helps managers understand local extinction risks for these trees. To investigate population structure within Clark’s nutcracker, we developed primers for and characterized 13 polymorphic microsatellite loci. In a screen of 22 individuals from one population, levels of variability ranged from 6 to 15 alleles. No loci were found to be linked, although 4 loci revealed significant departures from Hardy–Weinberg equilibrium and evidence of null alleles. These microsatellite loci will enable population genetic analyses of Clark’s nutcrackers, which could provide insights into the spatial relationships between nutcrackers and the trees they help disperse.

Genetic characterization of the Pacific sheath-tailed bat (Emballonura semicaudata rotensis) using mitochondrial DNA sequence data

Abstract Genetic mark–recapture methods are increasingly being used to estimate demographic parameters in species where traditional techniques are problematic or imprecise. The federally endangered Indiana bat Myotis sodalis has declined dramatically and threats such as white-nose syndrome continue to afflict this species. To date, important demographic information for Indiana bats has been difficult to estimate precisely using traditional techniques such as emergence counts. Successful management and protection of Indiana bats requires better methods to estimate population sizes and survival rates throughout the year, particularly during summer when these bats reproduce and are widely dispersed away from their winter hibernacula. In addition, the familial makeup of maternity colonies is unknown, yet important for understanding local and regional population dynamics. We had four objectives in this study. For the first two objectives we investigated the potential use of DNA from fecal samples (fecal DNA) c...

Long-term persistence of horse fecal DNA in the environment makes equids particularly good candidates for noninvasive sampling

Abstract Emballonura semicaudata occurs in the southwestern Pacific and populations on many islands have declined or disappeared. One subspecies (E. semicaudata rotensis) occurs in the Northern Mariana Islands, where it has been extirpated from all but 1 island (Aguiguan). We assessed genetic similarity between the last population of E. s. rotensis and 2 other subspecies, and examined genetic diversity on Aguiguan. We sampled 12 E. s. rotensis, sequenced them at 3 mitochondrial loci, and compared them with published sequences from 2 other subspecies. All 12 E. s. rotensis had identical sequences in each of the 3 regions. Using cytochrome-b (Cytb) data E. s. rotensis was sister to E. s. palauensis in a clade separate from E. s. semicaudata. 12S ribosomal RNA (12S) sequences grouped all E. s. semicaudata in 1 clade with E. s. rotensis in a clade by itself. Genetic distances among the 3 subspecies at Cytb were smallest between E. s. palauensis and E. s. rotensis. Distance between E. s. semicaudata and the other 2 subspecies was not different from the distance between E. s. semicaudata and the full species E. raffrayana. A similar relationship was found using the 12S data. These distances are larger than those typically reported for mammalian subspecies using Cytb sequence and within the range of sister species.