SSR_pipeline: a bioinformatic infrastructure for identifying microsatellites from paired-end Illumina high-throughput DNA sequencing data

Mark P. Miller; Brian J. Knaus; Thomas D. Mullins; Susan M. Haig

doi:10.1093/jhered/est056

Home/Archive/USGS

USGS · 70048358

SSR_pipeline: a bioinformatic infrastructure for identifying microsatellites from paired-end Illumina high-throughput DNA sequencing data

Mark P. Miller·Brian J. Knaus·Thomas D. Mullins·Susan M. Haig

View original source ↗Export metadata ↓

Abstract

SSR_pipeline is a flexible set of programs designed to efficiently identify simple sequence repeats (e.g., microsatellites) from paired-end high-throughput Illumina DNA sequencing data. The program suite contains 3 analysis modules along with a fourth control module that can automate analyses of large volumes of data. The modules are used to 1) identify the subset of paired-end sequences that pass Illumina quality standards, 2) align paired-end reads into a single composite DNA sequence, and 3) identify sequences that possess microsatellites (both simple and compound) conforming to user-specified parameters. The microsatellite search algorithm is extremely efficient, and we have used it to identify repeats with motifs from 2 to 25bp in length. Each of the 3 analysis modules can also be used independently to provide greater flexibility or to work with FASTQ or FASTA files generated from other sequencing platforms (Roche 454, Ion Torrent, etc.). We demonstrate use of the program with data from the brine fly Ephydra packardi (Diptera: Ephydridae) and provide empirical timing benchmarks to illustrate program performance on a common desktop computer environment. We further show that the Illumina platform is capable of identifying large numbers of microsatellites, even when using unenriched sample libraries and a very small percentage of the sequencing capacity from a single DNA sequencing run. All modules from SSR_pipeline are implemented in the Python programming language and can therefore be used from nearly any computer operating system (Linux, Macintosh, and Windows).

Explore related subjects

Journal of Heredity Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Mark P. Miller, Brian J. Knaus, Thomas D. Mullins, Susan M. Haig. 2013-09-19. SSR_pipeline: a bioinformatic infrastructure for identifying microsatellites from paired-end Illumina high-throughput DNA sequencing data. https://doi.org/10.1093/jhered%2Fest056

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

A chromosome-level genome assembly of a vernal pool specialist amphibian, the Western Spadefoot, Spea hammondii

We assembled and annotated a chromosome-level genome for the Western Spadefoot, Spea hammondii (Anura, Scaphiopodidae) representing one of only three amphibians included in the California Conservation Genomics Project (CCGP). Spea hammondii is a vernal pool breeding anuran native to California and northwestern Baja California which has undergone both range contractions and local extirpations across its distribution, primarily due to habitat loss and degradation and drought. The species is recognized by the state of California as a Species of Special Concern and is proposed for listing under the United States Endangered Species Act. Using the established CCGP pipeline, this S. hammondii genome was produced using Pacific Biosciences HiFi long-reads and Omni-C proximity ligation, resulting in a de novo genome assembly 1.14 Gb in length, distributed across 479 scaffolds (scaffold N50 = 120.8 Mb; largest scaffold = 183.6 Mb) with a BUSCO completeness score of 90.9% using a conserved tetrapod ortholog set. Our assembly shows high base accuracy (quality value [QV] = 63.7) and low frameshift error in coding regions (QV 50.42). Annotation of this genome yielded 20,434 genes with a BUSCO completeness score of 94.7%. This genome assembly, in combination with range-wide resequencing data from CCGP, will facilitate statewide population genomic assessments to delineate conservation units, quantify inbreeding and genomic load, and test for adaptive variation associated with vernal pool hydrology and drought tolerance, all of which are important considerations in the proposed federal listing.

Journal of Heredity↗

USGS2025

Bovine SNP array-based genetic assessment of American plains bison at American Prairie

American plains bison ( Bison bison bison , bison hereafter) experienced an extreme demographic bottleneck in the late 1800s. The species has since rebounded but is primarily managed as small and isolated herds due to habitat and sociopolitical limitations. Thus, reintroducing bison and allowing herds to achieve as much of their natural dynamics as possible is a major conservation goal. Concerns about genetic diversity loss in small, isolated herds and the persistence of cattle-origin variants from historical crossbreeding efforts have made genetic analysis an important part of bison conservation. The limitations of the current conservation genetic tools which are based on traditional markers such as microsatellites and mitochondrial DNA sequences, may be overcome with genome-wide genotyping panels commonly developed for agricultural species. Bison reintroduction in the grasslands at American Prairie began in 2005. Genetic analysis on these herds has yet to be conducted. We used the Illumina 777K Bovine genotyping panel to obtain data from 197 bison and 179 domestic cows to understand the current population genetic state of bison at American Prairie and gain insight on cattle ( Bos taurus ) introgression. Overall, bison at American Prairie currently have relatively high genetic diversity, low inbreeding, and no obvious signs of cattle introgression. A more comprehensive evaluation of introgression, likely including whole-genome sequence data, would clarify this finding. These results can serve as a baseline for future comparison as part of a genetic monitoring framework.

Journal of Heredity↗

USGS2024

Evidence of extensive home range sharing among mother–daughter bobcat pairs in the wildland–urban interface of the Tucson Mountains

Urbanization impacts the structure and viability of wildlife populations. Some habitat generalists, such as bobcats ( Lynx rufus ), maintain populations at the intersection of wild and urban habitats (wildland urban interface, WUI), but impacts of urbanization on bobcat social structure are not well understood. Although commonly thought to establish exclusive home ranges among females, instances of mother-daughter home range sharing have been documented. We combined GPS localities with genomic relatedness inferences from double-digest restriction site associated DNA sequencing (ddRADseq) to investigate mother-daughter home range sharing in bobcats ( n = 38) at the WUI in the Tucson Mountains, Arizona, USA. We found the highest relatedness among females, which showed stronger isolation by distance than males and the population as a whole. Using mother-daughter relationships inferred from pedigree reconstruction, we found extensive mother-daughter home range sharing, compared to other females. Every mother identified as having at least one daughter, shared home ranges with one daughter, while other confirmed daughters established adjacent home ranges. Our results provide substantial support for the mother-daughter home range sharing hypothesis, as well as evidence of spatiotemporal overlap between mothers and daughters, adding to the body of research complicating the solitary felid paradigm. These results additionally challenge the notion of home range partitioning by prior rights land tenure, suggesting a role of matrilineal land tenure in home range establishment of daughters. Habitat fragmentation due to human population growth and urbanization thus has the potential to alter landscape genetic structure and social dynamics in bobcats.

Arizona↗