USGS Science⌕ Search

SEARCH · USGS Science

Results for “Heredity”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources
USGS2023

Assembly of the largest squamate reference genome to date: The western fence lizard, Sceloporus occidentalis

Spiny lizards (genus Sceloporus ) have long served as important systems for studies of behavior, thermal physiology, dietary ecology, vector biology, speciation, and biogeography. The western fence lizard, Sceloporus occidentalis , is found across most of the major biogeographical regions in the western United States and northern Baja California, Mexico, inhabiting a wide range of habitats, from grassland to chaparral to open woodlands. As small ectotherms, Sceloporus lizards are particularly vulnerable to climate change, and S. occidentalis has also become an important system for studying the impacts of land use change and urbanization on small vertebrates. Here, we report a new reference genome assembly for S. occidentalis , as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomics strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 608 scaffolds spanning 2,856 Mb, has a contig N50 of 18.9 Mb, a scaffold N50 of 98.4 Mb, and BUSCO completeness score of 98.1% based on the tetrapod gene set. This reference genome will be valuable for understanding ecological and evolutionary dynamics in S. occidentalis , the species status of the California endemic island fence lizard ( S. becki ), and the spectacular radiation of Sceloporus lizards.

California↗
USGS2023

A genome assembly for the southern Pacific rattlesnake, Crotalus oreganus helleri, in the western rattlesnake species complex

Rattlesnakes play important roles in their ecosystems by regulating prey populations, are involved in complex coevolutionary dynamics with their prey, and exhibit a variety of unusual adaptations, including maternal care, heat-sensing pit organs, hinged fangs, and medically-significant venoms. The western rattlesnake ( Crotalus oreganus ) is one of the widest ranging rattlesnake species, with a distribution from British Columbia, where it is listed as threatened, to Baja California and east across the Great Basin to western Wyoming, Colorado and New Mexico. Here, we report a new reference genome assembly for one of six currently recognized subspecies, C. oreganus helleri , as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomic sequencing strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 698 scaffolds spanning 1,564,812,557 base pairs, has a contig N50 of 64.7 Mb, a scaffold N50 of 110.8 Mb, and BUSCO complete score of 90.5%. This reference genome will be valuable for studies on the genomic basis of venom evolution and variation within Crotalus , in resolving the taxonomy of C. oreganus and its relatives, and for the conservation and management of rattlesnakes in general.

California↗
USGS2023

A reference genome assembly for the continentally distributed ring-necked snake, Diadophis punctatus

Snakes in the family Colubridae include more than 2,000 currently recognized species, and comprise roughly 75% of the global snake species diversity on Earth. For such a spectacular radiation, colubrid snakes remain poorly understood ecologically and genetically. Two subfamilies, Colubrinae (788 species) and Dipsadinae (833 species), comprise the bulk of colubrid species richness. Dipsadines are a speciose and diverse group of snakes that largely inhabit Central and South America, with a handful of small-body-size genera that have invaded North America. Among them, the ring-necked snake, Diadophis punctatus , has an incredibly broad distribution with 14 subspecies. Given its continental distribution and high degree of variation in coloration, diet, feeding ecology, and behavior, the ring-necked snake is an excellent species for the study of genetic diversity and trait evolution. Within California, six subspecies form a continuously distributed “ring species” around the Central Valley, while a seventh, the regal ring-necked snake, Diadophis punctatus regalis is a disjunct outlier and Species of Special Concern in the state. Here, we report a new reference genome assembly for the San Diego ring-necked snake, D. p. similis , as part of the California Conservation Genomics Project. This assembly comprises a total of 444 scaffolds spanning 1,783 Mb and has a contig N50 of 8.0 Mb, scaffold N50 of 83 Mb, and BUSCO completeness score of 94.5%. This reference genome will be a valuable resource for studies of the taxonomy, conservation, and evolution of the ring-necked snake across its broad, continental distribution.

Journal of Heredity↗
USGS2023

Genomic insights into isolation of the threatened Florida crested caracara (Caracara plancus)

We conducted a population genomic study of the crested caracara ( Caracara plancus ) using samples ( n = 290) collected from individuals in Florida, Texas, and Arizona, United States. Crested caracaras are non-migratory raptors ranging from the southern tip of South America to the southern United States, including a federally protected relict population in Florida long thought to have been isolated since the last ice age. Our objectives were to evaluate genetic diversity and population structure of Florida’s apparently isolated population and to evaluate taxonomic relationships of crested caracaras at the northern edge of their range. Using DNA purified from blood samples, we conducted double-digest restriction site associated DNA sequencing and sequenced the mitochondrial ND2 gene. Analyses of population structure using over 9,000 SNPs suggest that two major clusters are best supported, one cluster including only Florida individuals and the other cluster including Arizona and Texas individuals. Both SNPs and mitochondrial haplotypes reveal the Florida population to be highly differentiated genetically from Arizona and Texas populations, whereas, Arizona and Texas populations are moderately differentiated from each other. The Florida population’s mitochondrial haplotypes form a separate monophyletic group, while Arizona and Texas populations share mitochondrial haplotypes. Results of this study provide substantial genetic evidence that Florida’s crested caracaras have experienced long-term isolation from caracaras in Arizona and Texas and thus, represent a distinct evolutionary lineage possibly warranting distinction as an Evolutionarily Significant Unit (ESU) or subspecies. This study will inform conservation strategies focused on long-term survival of Florida’s distinct, panmictic population.

Journal of Heredity↗
USGS2023

Reduction of genetic diversity in ‘Alalā (Hawaiian crow; Corvus hawaiiensis) between the late 1800s and the late 1900s

Genetic and genomic data are increasingly used to aid conservation management of endangered species by providing insights into evolutionary histories, factors associated with extinction risks, and potential for future adaptation. For the ‘Alalā, or Hawaiian crow ( Corvus hawaiiensis ), genetic concerns include negative correlations between inbreeding and hatching success. However, it is unclear if low genetic diversity and inbreeding depression are consequences of a historical population bottleneck, or if ‘Alalā had historically low genetic diversity that predated human influence, perhaps as a result of earlier declines or founding events. In this study, we applied a hybridization-based sequence capture to generate a genome-wide single nucleotide polymorphism (SNP) dataset for comparing historical specimens collected in the 1890s, when ‘Alalā were more numerous, to samples taken between 1973 and 1998, when ‘Alalā population densities were near the lowest documented levels in the wild, prior to all individuals being collected for captive rearing. We found low genome-wide diversity in both sample groups, however, the modern sample group (1973 to 1998 cohort) exhibited relatively fewer polymorphic alleles, a lower proportion of polymorphic loci, and lower observed heterozygosity, consistent with a population decline and potential bottleneck effects. These results combined with a current low population size highlight the importance of continued efforts by conservation managers to mitigate inbreeding and maintain founder representation to preserve what genetic diversity remains.

Hawaii↗
USGS2023

Glacial vicariance and secondary contact shape demographic histories in a freshwater mussel species complex

Characterizing the mechanisms influencing the distribution of genetic variation in aquatic species can be difficult due to the dynamic nature of hydrological landscapes. In North America’s Central Highlands, a complex history of glacial dynamics, long-term isolation, and secondary contact have shaped genetic variation in aquatic species. Although the effects of glacial history have been demonstrated in many taxa, responses are often lineage- or species-specific and driven by organismal ecology. In this study, we reconstruct the evolutionary history of a freshwater mussel species complex using a suite of mitochondrial and nuclear loci to resolve taxonomic and demographic uncertainties. Our findings do not support Pleurobema rubrum as a valid species, which is proposed for listing as threatened under the U.S. Endangered Species Act. We synonymize P. rubrum under Pleurobema sintoxia —a common and widespread species found throughout the Mississippi River Basin. Further investigation of patterns of genetic variation in P. sintoxia identified a complex demographic history, including ancestral vicariance and secondary contact, within the Eastern Highlands. We hypothesize these patterns were shaped by ancestral vicariance driven by the formation of Lake Green and subsequent secondary contact after the last glacial maximum. Our inference aligns with demographic histories observed in other aquatic taxa in the region and mirrors patterns of genetic variation of a freshwater fish species ( Erimystax dissimilis ) confirmed to serve as a parasitic larval host for P. sintoxia . Our findings directly link species ecology to observed patterns of genetic variation and may have significant implications for future conservation and recovery actions of freshwater mussels.

Mississippi River basin↗
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↗
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↗
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↗
USGS2012

Genetic population substructure in bison at Yellowstone National Park

The Yellowstone National Park bison herd is 1 of only 2 populations known to have continually persisted on their current landscape since pre-Columbian times. Over the last century, the census size of this herd has fluctuated from around 100 individuals to over 3000 animals. Previous studies involving radiotelemetry, tooth wear, and parturition timing provide evidence of at least 2 distinct groups of bison within Yellowstone National Park. To better understand the biology of Yellowstone bison, we investigated the potential for limited gene flow across this population using multilocus Bayesian clustering analysis. Two genetically distinct and clearly defined subpopulations were identified based on both genotypic diversity and allelic distributions. Genetic cluster assignments were highly correlated with sampling locations for a subgroup of live capture individuals. Furthermore, a comparison of the cluster assignments to the 2 principle winter cull sites revealed critical differences in migration patterns across years. The 2 Yellowstone subpopulations display levels of differentiation that are only slightly less than that between populations which have been geographically and reproductively isolated for over 40 years. The identification of cryptic population subdivision and genetic differentiation of this magnitude highlights the importance of this biological phenomenon in the management of wildlife species.

Journal of Heredity↗
USGS2012

Variation in migratory behavior influences regional genetic diversity and structure among American kestrel populations (Falco sparverius) in North America

Birds employ numerous strategies to cope with seasonal fluctuations in high-quality habitat availability. Long distance migration is a common tactic; however, partial migration is especially common among broadly distributed species. Under partial migration systems, a portion of a species migrates, whereas the remainder inhabits breeding grounds year round. In this study, we identified effects of migratory behavior variation on genetic structure and diversity of American Kestrels (Falco sparverius), a widespread partial migrant in North America. American Kestrels generally migrate; however, a resident group inhabits the southeastern United States year round. The southeastern group is designated as a separate subspecies (F. s. paulus) from the migratory group (F. s. sparverius). Using mitochondrial DNA and microsatellites from 183 and 211 individuals, respectively, we illustrate that genetic structure is stronger among nonmigratory populations, with differentiation measures ranging from 0.060 to 0.189 depending on genetic marker and analysis approach. In contrast, measures from western North American populations ranged from 0 to 0.032. These findings suggest that seasonal migratory behavior is also associated with natal and breeding dispersal tendencies. We likewise detected significantly lower genetic diversity within nonmigratory populations, reflecting the greater influence of genetic drift in small populations. We identified the signal of population expansion among nonmigratory populations, consistent with the recent establishment of higher latitude breeding locations following Pleistocene glacial retreat. Differentiation of F. s. paulus and F. s. sparverius reflected subtle differences in allele frequencies. Because migratory behavior can evolve quickly, our analyses suggest recent origins of migratory American Kestrel populations in North America.

Journal of Heredity↗
USGS2012

Yellowstone bison genetics: let us move forward

White and Wallen (2012) disagree with the conclusions and suggestions made in our recent assessment of population structure among Yellowstone National Park (YNP) bison based on 46 autosomal microsatellite loci in 661 animals (Halbert et al. 2012). First, they suggest that "the existing genetic substructure (that we observed) was artificially created." Specifically, they suggest that the substructure observed between the northern and central populations is the result of human activities, both historical and recent. In fact, the genetic composition of all known existing bison herds was created by, or has been influenced by, anthropogenic activities, although this obviously does not reduce the value of these herds for genetic conservation (Dratch and Gogan 2010). As perspective, many, if not most, species of conservation concern have been influenced by human actions and as a result currently exist as isolated populations. However, it is quite difficult to distinguish between genetic differences caused by human actions and important ancestral variation contained in separate populations without data from early time periods. Therefore, to not lose genetic variation that may be significant or indicative of important genetic variation, the generally acceptable management approach is to attempt to retain this variation based on the observed population genetic subdivision (Hedrick et al. 1986).

Yellowstone National Park↗
USGS2013

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

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).

Journal of Heredity↗
USGS2013

The population history of endogenous retroviruses in mule deer (Odocoileus heminous)

Mobile elements are powerful agents of genomic evolution and can be exceptionally informative markers for investigating species and population-level evolutionary history. While several studies have utilized retrotransposon-based insertional polymorphisms to resolve phylogenies, few population studies exist outside of humans. Endogenous retroviruses are LTR-retrotransposons derived from retroviruses that have become stably integrated in the host genome during past infections and transmitted vertically to subsequent generations. They offer valuable insight into host-virus co-evolution and a unique perspective on host evolutionary history because they integrate into the genome at a discrete point in time. We examined the evolutionary history of a cervid endogenous gammaretrovirus (CrERVγ) in mule deer ( Odocoileus hemionus ). We sequenced 14 CrERV proviruses (CrERV-in1 to -in14), and examined the prevalence and distribution of 13 proviruses in 262 deer among 15 populations from Montana, Wyoming, and Utah. CrERV absence in white-tailed deer ( O. virginianus ), identical 5′ and 3′ long terminal repeat (LTR) sequences, insertional polymorphism, and CrERV divergence time estimates indicated that most endogenization events occurred within the last 200000 years. Population structure inferred from CrERVs (F ST = 0.008) and microsatellites (θ = 0.01) was low, but significant, with Utah, northwestern Montana, and a Helena herd being particularly differentiated. Clustering analyses indicated regional structuring, and non-contiguous clustering could often be explained by known translocations. Cluster ensemble results indicated spatial localization of viruses, specifically in deer from northeastern and western Montana. This study demonstrates the utility of endogenous retroviruses to elucidate and provide novel insight into both ERV evolutionary history and the history of contemporary host populations.

Montana;Utah;Wyoming↗
USGS2017

Evolutionary and functional mitogenomics associated with the genetic restoration of the Florida panther

Florida panthers are endangered pumas that currently persist in reduced patches of habitat in South Florida, USA. We performed mitogenome reference-based assemblies for most parental lines of the admixed Florida panthers that resulted from the introduction of female Texas pumas into South Florida in 1995. With the addition of 2 puma mitogenomes, we characterized 174 single nucleotide polymorphisms (SNPs) across 12 individuals. We defined 5 haplotypes (Pco1–Pco5), one of which (Pco1) had a geographic origin exclusive to Costa Rica and Panama and was possibly introduced into the Everglades National Park, Florida, prior to 1995. Haplotype Pco2 was native to Florida. Haplotypes Pco3 and Pco4 were exclusive to Texas, whereas haplotype Pco5 had an undetermined geographic origin. Phylogenetic inference suggests that haplotypes Pco1–Pco4 diverged ~202000 (95% HPDI = 83000–345000) years ago and that haplotypes Pco2–Pco4 diverged ~61000 (95% HPDI = 9000–127000) years ago. These results are congruent with a south-to-north continental expansion and with a recent North American colonization by pumas. Furthermore, pumas may have migrated from Texas to Florida no earlier than ~44000 (95% HPDI = 2000–98000) years ago. Synonymous mutations presented a greater mean substitution rate than other mitochondrial functional regions: nonsynonymous mutations, tRNAs, rRNAs, and control region. Similarly, all protein-coding genes were under predominant negative selection constraints. We directly and indirectly assessed the presence of potential deleterious SNPs in the ND2 and ND5 genes in Florida panthers prior to and as a consequence of the introduction of Texas pumas. Screenings for such variants are recommended in extant Florida panthers.

Journal of Heredity↗
USGS2017

Estimating inbreeding rates in natural populations: Addressing the problem of incomplete pedigrees

Understanding and estimating inbreeding is essential for managing threatened and endangered wildlife populations. However, determination of inbreeding rates in natural populations is confounded by incomplete parentage information. We present an approach for quantifying inbreeding rates for populations with incomplete parentage information. The approach exploits knowledge of pedigree configurations that lead to inbreeding coefficients of F = 0.25 and F = 0.125, allowing for quantification of Pr( I | k ): the probability of observing pedigree I given the fraction of known parents ( k ). We developed analytical expressions under simplifying assumptions that define properties and behavior of inbreeding rate estimators for varying values of k . We demonstrated that inbreeding is overestimated if Pr( I | k ) is not taken into consideration and that bias is primarily influenced by k . By contrast, our new estimator, incorporating Pr( I | k ), is unbiased over a wide range of values of k that may be observed in empirical studies. Stochastic computer simulations that allowed complex inter- and intragenerational inbreeding produced similar results. We illustrate the effects that accounting for Pr( I | k ) can have in empirical data by revisiting published analyses of Arabian oryx ( Oryx leucoryx ) and Red deer ( Cervus elaphus ). Our results demonstrate that incomplete pedigrees are not barriers for quantifying inbreeding in wild populations. Application of our approach will permit a better understanding of the role that inbreeding plays in the dynamics of populations of threatened and endangered species and may help refine our understanding of inbreeding avoidance mechanisms in the wild.

Journal of Heredity↗
USGS2018

Genome-wide analysis of SNPs is consistent with no domestic dog ancestry in the endangered Mexican Wolf (Canis lupus baileyi)

The Mexican gray wolf ( Canis lupus baileyi ) was historically distributed throughout the southwestern United States and northern Mexico. Extensive predator removal campaigns during the early 20th century, however, resulted in its eventual extirpation by the mid 1980s. At this time, the Mexican wolf existed only in 3 separate captive lineages (McBride, Ghost Ranch, and Aragón) descended from 3, 2, and 2 founders, respectively. These lineages were merged in 1995 to increase the available genetic variation, and Mexican wolves were reintroduced into Arizona and New Mexico in 1998. Despite the ongoing management of the Mexican wolf population, it has been suggested that a proportion of the Mexican wolf ancestry may be recently derived from hybridization with domestic dogs. In this study, we genotyped 87 Mexican wolves, including individuals from all 3 captive lineages and cross-lineage wolves, for more than 172000 single nucleotide polymorphisms. We identified levels of genetic variation consistent with the pedigree record and effects of genetic rescue. To identify the potential to detect hybridization with domestic dogs, we compared our Mexican wolf genotypes with those from studies of domestic dogs and other gray wolves. The proportion of Mexican wolf ancestry assigned to domestic dogs was only between 0.06% (SD 0.23%) and 7.8% (SD 1.0%) for global and local ancestry estimates, respectively; and was consistent with simulated levels of incomplete lineage sorting. Overall, our results suggested that Mexican wolves lack biologically significant ancestry with dogs and have useful implications for the conservation and management of this endangered wolf subspecies.

Journal of Heredity↗
USGS2018

Hatchery strain contributions to emerging wild lake trout populations in Lake Huron

Recent assessments indicate the emergence of naturally produced lake trout ( Salvelinus namaycush ) recruitment throughout Lake Huron in the North American Laurentian Great Lakes (>50% of fish <7 years). Because naturally produced fish derived from different stocked hatchery strains are unmarked, managers cannot distinguish strains contributing to natural recruitment. We used 15 microsatellite loci to identify strains of naturally produced lake trout ( N = 1567) collected in assessment fisheries during early (2002–2004) and late (2009–2012) sampling periods. Individuals from 13 American and Canadian hatchery strains ( N = 1143) were genotyped to develop standardized baseline information. Strain contributions were estimated using a Bayesian inferential approach. Deviance information criteria were used to compare models evaluating strain contributions at different spatial and temporal scales. The best performing models were the most complex models, suggesting that hatchery strain contributions to naturally produced lake trout varied spatially among management districts and temporally between time periods. Contributions of Seneca strain lake trout were consistently high across most management districts, with contributions increasing from early to late time periods (estimates ranged from 52% to 94% for the late period across 8 of 9 districts). Strain contributions deviated from expectations based on historical stocking levels, indicating strains differed with respect to survival, reproductive success, and/or dispersal. Knowledge of recruitment levels of strains stocked in different management districts, and how strain-specific recruitment varies temporally, spatially, and as a function of local or regional stocking is important to prioritize strains for future stocking and management of the transition process from primarily hatchery to naturally produced stocks.

Journal of Heredity↗