USGS Science⌕ Search

SEARCH · USGS Science

Results for “Genes”

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 1,027 records · Page 57Linked to original sources
USGS2017

Autotrophic microbial arsenotrophy in arsenic-rich soda lakes

A number of prokaryotes are capable of employing arsenic oxy-anions as either electron acceptors [arsenate; As(V)] or electron donors [arsenite; As(III)] to sustain arsenic-dependent growth (‘arsenotrophy’). A subset of these microorganisms function as either chemoautotrophs or photoautotrophs, whereby they gain sufficient energy from their redox metabolism of arsenic to completely satisfy their carbon needs for growth by autotrophy, that is the fixation of inorganic carbon (e.g. HCO 3 − ) into their biomass. Here we review what has been learned of these processes by investigations we have undertaken in three soda lakes of the western USA and from the physiological characterizations of the relevant bacteria, which include the critical genes involved, such as respiratory arsenate reductase ( arrA ) and the discovery of its arsenite-oxidizing counterpart ( arxA ). When possible, we refer to instances of similar process occurring in other, less extreme ecosystems and by microbes other than haloalkaliphiles.

FEMS Microbiology Letters↗
USGS2020

Nanopore amplicon sequencing reveals molecular convergence and local adaptation of rhodopsin in Great Lakes salmonids

Local adaptation can drive diversification of closely related species across environmental gradients and promote convergence of distantly related taxa that experience similar conditions. We examined a potential case of adaptation to novel visual environments in a species flock (Great Lakes salmonids, genus Coregonus ) using a new amplicon genotyping protocol on the Oxford Nanopore Flongle and MinION. We sequenced five visual opsin genes for individuals of C. artedi , C. hoyi , C. kiyi , and C. zenithicus. Comparisons revealed species-specific differences in a key spectral tuning amino acid in rhodopsin (Tyr261Phe substitution), suggesting local adaptation of C. kiyi to the blue-shifted depths of Lake Superior. Ancestral state reconstruction demonstrates that parallel evolution and “toggling” at this amino acid residue has occurred several times across the fish tree of life, resulting in identical changes to the visual systems of distantly related taxa across replicated environmental gradients. Our results suggest that ecological differences and local adaptation to distinct visual environments are strong drivers of both evolutionary parallelism and diversification.

Lake Superior↗
USGS2023

A highly-contiguous and annotated genome assembly of the Lesser Prairie-Chicken (Tympanuchus pallidicinctus).

The Lesser Prairie-Chicken ( Tympanuchus pallidicinctus ; LEPC) is an iconic North American prairie grouse, renowned for ornate and spectacular breeding season displays. Unfortunately, the species has disappeared across much of its historical range, with corresponding precipitous declines in contemporary population abundance, largely due to climatic and anthropogenic factors. These declines led to a 2022 US Fish and Wildlife decision to identify and list two distinct population segments (DPSs; i.e., northern and southern DPSs) as threatened or endangered under the 1973 Endangered Species Act. Herein, we describe an annotated reference genome that was generated from a LEPC sample collected from the southern DPS. We chose a representative from the southern DPS because of the potential for introgression in the northern DPS, where some populations hybridize with the Greater Prairie-Chicken ( Tympanuchus cupido ). This new LEPC reference assembly consists of 206 scaffolds, an N50 of 45 Mb, and 15,563 predicted protein-coding genes. We demonstrate the utility of this new genome assembly by estimating genome-wide heterozygosity in a representative LEPC and in related species. Heterozygosity in a LEPC sample was 0.0024, near the middle of the range (0.0003–0.0050) of related species. Overall, this new assembly provides a valuable resource that will enhance evolutionary and conservation genetic research in prairie grouse.

Genome Biology and Evolution↗
USGS2025

Speciation genomics in the tiger whiptail lizards ( Aspidoscelis tigris complex)

The transition from small genetic to genome-scale datasets for studying biodiversity has revealed that genetic exchange through introgressive hybridization is a widespread phenomenon in nature. Despite this, a lack of high-quality reference genomes for most non-model species limits our understanding of the impact of this process for many taxonomic groups. This restricts the range of insights that genomic tools can provide for conservation biologists, who often hope to employ genomic datasets to accurately identify historically isolated lineages to protect and to predict their evolutionary fate in the face of environmental change. Tiger whiptail lizards ( Aspidoscelis tigris complex) are an abundant and important ecological component of ecosystems across the southwestern United States. In this study, we assembled and annotated a chromosome-level reference genome for A. t. stejnegeri from coastal California. We then used this reference genome to reconstruct patterns of speciation and admixture within the larger species complex, finding evidence that gene flow is widespread both geographically and across the genome.

northern Mexico, southwestern United States↗
USGS2009

Comparative endocrinology in the 21st century

Hormones coordinate developmental, physiological, and behavioral processes within and between all living organisms. They orchestrate and shape organogenesis from early in development, regulate the acquisition, assimilation, and utilization of nutrients to support growth and metabolism, control gamete production and sexual behavior, mediate organismal responses to environmental change, and allow for communication of information between organisms. Genes that code for hormones; the enzymes that synthesize, metabolize, and transport hormones; and hormone receptors are important targets for natural selection, and variation in their expression and function is a major driving force for the evolution of morphology and life history. Hormones coordinate physiology and behavior of populations of organisms, and thus play key roles in determining the structure of populations, communities, and ecosystems. The field of endocrinology is concerned with the study of hormones and their actions. This field is rooted in the comparative study of hormones in diverse species, which has provided the foundation for the modern fields of evolutionary, environmental, and biomedical endocrinology. Comparative endocrinologists work at the cutting edge of the life sciences. They identify new hormones, hormone receptors and mechanisms of hormone action applicable to diverse species, including humans; study the impact of habitat destruction, pollution, and climatic change on populations of organisms; establish novel model systems for studying hormones and their functions; and develop new genetic strains and husbandry practices for efficient production of animal protein. While the model system approach has dominated biomedical research in recent years, and has provided extraordinary insight into many basic cellular and molecular processes, this approach is limited to investigating a small minority of organisms. Animals exhibit tremendous diversity in form and function, life-history strategies, and responses to the environment. A major challenge for life scientists in the 21st century is to understand how a changing environment impacts all life on earth. A full understanding of the capabilities of organisms to respond to environmental variation, and the resilience of organisms challenged by environmental changes and extremes, is necessary for understanding the impact of pollution and climatic change on the viability of populations. Comparative endocrinologists have a key role to play in these efforts.

Integrative and Comparative Biology↗
USGS2026

Cnidarian–algal partnerships structure bacterial communities during strobilation in Cassiopea xamachana

Cnidarian–algal (Symbiodiniaceae) symbioses rely on complex interactions among the cnidarian host, algal symbionts, and associated bacterial communities. In the upside-down jellyfish Cassiopea xamachana , the polyp-to-medusa transition (strobilation) requires the establishment of symbiosis with Symbiodiniaceae algal partners, yet bacterial community dynamics during this developmental process remain unknown. Here, we experimentally induced symbiosis in aposymbiotic polyps using four algal treatments: xenic Symbiodinium microadriaticum (native symbiont), xenic Breviolum minutum , antibiotic-treated B. minutum , and a photosynthetically impaired B. minutum mutant. We combined 16S rRNA gene sequencing with measurements of photosynthetic efficiency, asexual budding, and algal surface N-glycan profiles to characterize holobiont assembly during symbiosis onset and strobilation. Algal treatment structured bacterial communities in both algal cultures and polyp tissues. Our analyses identified a set of amplicon sequence variants that consistently distinguished strobilating polyps from non-strobilating aposymbiotic and mutant polyps, in addition to potential bacterial biomarkers associated with successful metamorphosis. Strobilation was associated with the enrichment of bacterial communities putatively involved in sulfur and nitrogen cycling, whereas non-strobilating aposymbiotic and mutant polyps were characterized by opportunistic bacteria and increased community variability. Together, these results reveal coordinated changes in algal physiology, surface glycan profiles, and bacterial community structure associated with successful strobilation in C. xamachana and support a model in which tripartite host–alga–bacteria interactions influence cnidarian life stage transitions.

ISME Communications↗
USGS2024

Transglobal spread of an ecologically relevant sea urchin parasite

Mass mortality of the dominant coral reef herbivore Diadema antillarum in the Caribbean in the early 1980s contributed to a persistent phase shift from coral- to algal-dominated reefs. In 2022, a scuticociliate most closely related to Philaster apodigitiformis caused further mass mortality of D. antillarum across the Caribbean, leading to >95% mortality at affected sites. Mortality was also reported in the related species Diadema setosum in the Mediterranean in 2022, though the causative agent of the Mediterranean outbreak has not yet been determined. In April 2023, mass mortality of Diadema setosum occurred along the Sultanate of Oman's coastline. Urchins displayed signs compatible with scuticociliatosis including abnormal behavior, drooping and loss of spines, followed by tissue necrosis and death. Here we report the detection of an 18S rRNA gene sequence in abnormal urchins from Muscat, Oman, that is identical to the Philaster strain responsible for D. antillarum mass mortality in the Caribbean. We also show that scuticociliatosis signs can be elicited in Diadema setosum by experimental challenge with the cultivated Philaster strain associated with Caribbean scuticociliatosis. These results demonstrate the Philaster sp. associated with D. antillarum mass mortality has rapidly spread to geographically distant coral reefs, compelling global-scale awareness and monitoring for this devastating condition through field surveys, microscopy, and molecular microbiological approaches, and prompting investigation of long-range transmission mechanisms.

ISME Journal↗
USGS2025

Manipulation of the Symbiodiniaceae microbiome confers multigenerational impacts on symbioses and reproductive ecology of its Exaiptasia diaphana host

Symbiodiniaceae-associated microbiota strongly influence cnidarian symbioses. We systematically reduced the bacterial and fungal communities associated with Symbiodiniaceae to study potential effects on the cnidarian holobiont Exaiptasia diaphana (Aiptasia). Clonal anemones were inoculated with xenic Breviolum minutum (SSB01) and microbiome-manipulated cultures after antibacterial or antifungal treatment. The asexual reproduction of pedal laceration allowed for three generations of clonal aposymbiotic Aiptasia to be utilised in this study, from the initial adult generation (G0), to the first (G1), and second (G2) generation. We inoculated small and large G1 Aiptasia with SSB01 algae and monitored onset of symbiosis, rate of algal proliferation, and holobiont characteristics. Sequencing the 16S and 18S rRNA gene regions identified significant differences in the bacterial and fungal communities of the G0 and G1 generations, alongside differences between the size classes of small and large G1 anemones. The microbiome of larger G1 individuals was distinct to the smaller G1 anemones, suggesting a microbiome maturation process. Control Breviolum minutum cultures exhibited a significantly greater proliferation rate in large G1 anemones when compared to antibacterial or antifungal treated cultures, whereas the opposite trend was documented in the small G1 anemones. Although no differences were observed between algal photochemical parameters, or the growth and polyp activity of G1 juveniles, we observed a significant influence in the production of G2 clones between treatments. Overall, we provide strong ecological implications of manipulating Symbiodiniaceae microbiome, not for the algae themselves, but for the maturation of the host Aiptasia, as well as for the cnidarian holobiont over multiple generations.

The ISME Journal↗
USGS2025

Sulfide stress tolerance as a controller of methane production in temperate wetlands

Wetlands are a major source of methane emissions and contribute to the observed increase in atmospheric methane over the last 20 years. Methane production in wetlands is the final step of carbon decomposition performed by anaerobic archaea. Although hydrogen/carbon dioxide and acetate are the substrates most often attributed to methanogenesis, other substrates—such as methylated compounds—may additionally play important roles in driving methane production in wetland systems. Here we conducted mesocosm experiments combined with genome-resolved metatranscriptomics to investigate the impact of diverse methanogenic substrate amendment on methanogenesis in two high methane-emitting wetlands with distinct geochemistry, termed P7 and P8. Methanol amendment resulted in high methane production at both sites, whereas acetate and formate amendment only stimulated methanogenesis in P7 mesocosms, where aqueous sulfide concentrations were lower. In P7 sediments, formate amendment fueled acetogenic microbes that produced acetate, which was subsequently utilized by acetoclastic methanogens. In contrast to expression profiles in P7 mesocosms, active methylotrophic methanogen genomes from P8 showed increased expression of genes related to membrane remodeling and DNA damage repair, indicative of stress tolerance mechanisms to counter sulfide toxicity. Methylotrophic methanogenesis generates higher free energy yields than acetoclastic methanogenesis, which likely enables allocation of more energy toward stress responses. These findings contribute to the growing body of literature highlighting methylotrophic methanogenesis as an important methane production pathway in wetlands. By using less competitive substrates like methanol that provide greater energy yields, methylotrophic methanogens may invest in physiological strategies that provide competitive advantages across a range of environmental stresses.

North Dakota↗
USGS2026

Diverse cyanopeptides follow distinct temporal succession patterns in freshwater harmful algal blooms

Toxic cyanobacterial harmful algal blooms (cyanoHABs) threaten freshwater resources globally and are intensifying with increasing eutrophication. Bloom toxicity is strongly influenced by intraspecific variation in the biosynthetic repertoires of toxic cyanobacteria, yet few studies examine the diversity of cyanobacterial cyanopeptides beyond hepatotoxic microcystins . To understand the dynamics and drivers of cyanopeptide diversity in cyanoHABs, we analyzed temporal patterns of cyanobacteria, metabolites, and their biosynthetic gene clusters (BGCs) in western Lake Erie using a 7-year time series (2016–2022) of metagenomic and metabolomic data. Our findings demonstrate that shifts from Microcystis to Dolichospermum occur later in the bloom season, coinciding with lower temperatures. Modules of co-varying BGCs (biosynthesis modules) from these genera were identified with hierarchical clustering, with uncharacterized BGCs among the most abundant. Biosynthesis modules rich in nonribosomal peptide synthetases (NRPS) peaked in early August, coinciding with elevated levels of inorganic nitrogen, warmer temperatures, and high Microcystis abundance. In contrast, modules rich in polyketide synthases (PKS) and ribosomally synthesized and post-translationally modified peptides (RiPPs) peaked following the Microcystis maximum in mid-August. Metabolomic analyses confirmed that metabolites followed shared seasonal patterns with their associated biosynthesis modules, forming three phases characterized by (i) microcystins, (ii) anabaenopeptins and aeruginosins, and (iii) aerucyclamides. These phases co-varied with bottom-up and top-down pressures, with later phases coinciding with increased microbially processed organic nitrogen and reduced detection of grazers. This study demonstrates consistent seasonal patterns of cyanobacterial metabolite succession and co-occurrence beyond microcystins, suggesting tradeoffs between biosynthetic resource demands and ecological controls.

Ohio↗
USGS2025

Identical sequence types of Yersinia ruckeri associated with lethal disease in wild-caught invasive Blue Catfish and cultured hybrid catfish (Channel Catfish ♀ × Blue Catfish ♂) from disparate aquatic ecosystems

Objective The Blue Catfish Ictalurus furcatus is commonly raised in warmwater aquaculture in the United States to produce Channel Catfish I. punctatus × Blue Catfish hybrids. It is also a prominent aquatic invasive species of concern in the mid-Atlantic region of the United States. Here, Yersina ruckeri was isolated from moribund Blue Catfish and hybrid catfish from disparate regions of the USA. The goal of the research here was to compare these Y. ruckeri strains to each other and other known strains for which adequate sequence data was available. In addition, we sought to determine if the strain from Blue Catfish was pathogenic to Rainbow Trout Oncorhynchus mykiss . Methods Moribund hybrid catfish from culture ponds in Mississippi were processed for diagnostic evaluation in March 2016. In April 2022, a moribund Blue Catfish specimen was collected from a tributary of the Nanticoke River in Maryland. Bacterial isolates were identified and characterized using biochemical tests, antimicrobial sensitivity screening, serotyping, and complete or partial genome sequencing. Disease pathology was described via histology. The isolate from Blue Catfish was used in challenge experiments to determine if it was pathogenic to Rainbow Trout. Multilocus sequencing typing was conducted using the PubMLST database. Results Biochemical testing was consistent with Y. ruckeri . A draft genome of the Y. ruckeri isolate was assembled based on Oxford Nanopore Technology sequencing and identified a single genomic replicon (3,791,418 bp) consistent in size to other Y. ruckeri genomes and a pLT plasmid (60, 933 bp). The challenge study demonstrated no significant virulence of this isolate for Rainbow Trout ( Y. ruckeri ). This isolate was most similar to other strains isolated from ictalurids. Notably, the gyrase B gene from this isolate was identical to that of archived strains isolated from moribund Mississippi hybrid catfish aquaculture during 2016 and these isolates share identical PubMLST sequence type profiles. Similarly, they shared a pLT plasmid that differed by only 6 bp. This plasmid has never been reported from trout isolates and appears to be unique to ictalurids. Conclusions Analyses here provide preliminary genetic evidence that geographically distant (Maryland and Mississippi, USA) isolates of Y. ruckeri from ictalurids are genetically similar to each other and Y. ruckeri (strain SC09) that infects ictalurids in China. This strain is not a biothreat to Rainbow Trout at typical culture temperatures.

Journal of Aquatic Animal Health↗
USGS2001

Subspecific affinity of black bears in the White River National Wildlife Refuge

The black bear population of the White River National Wildlife Refuge (NWR) is adjacent to populations of black bear in Louisiana ( Urusus americanus luteolus ) which are listed as threatened under the U.S. Endangered Species Act. Wildlife management plans can pose restrictions on bear harvests and timber extraction; therefore the management plan for the White River NWR is sensitive to subspecific classification of its bear population. The objective of this study was to analyze genetic variation in the White River NWR and seven adjacent populations of black bears to assess the subspecific affinity of the White River NWR population. Here we report the variation at seven microsatellite DNA loci among eight black bear populations. The patterns of genetic variation gave strong support for distinguishing a southern group of black bears comprised of the White River, Arkansas; Tensas River, Louisiana; Upper Atchafalaya, Louisiana; Lower Atchafalaya, Louisiana; and Alabama/Mississippi populations. Phylogenetic analysis of individual variation suggested that historical black bear introductions into Arkansas and Louisiana affected gene pools of certain southern receiving populations, but did not significantly change interpopulation relatedness. Phylogenetic inferences at both the population and individual levels support the hypothesis that the White River NWR population of black bears belongs to the U. a. luteolus subspecies.

Arkansas↗
USGS2022

Give and take: Effects of genetic admixture on mutation load in endangered Florida panthers

Genetic admixture is a biological event inherent to genetic rescue programs aimed at the long-term conservation of endangered wildlife. Although the success of such programs can be measured by the increase in genetic diversity and fitness of subsequent admixed individuals, predictions supporting admixture costs to fitness due to the introduction of novel deleterious alleles are necessary. Here, we analyzed nonsynonymous variation from conserved genes to quantify and compare levels of mutation load (i.e. proportion of deleterious alleles and genotypes carrying these alleles) among endangered Florida panthers and non-endangered Texas pumas. Specifically, we used canonical (i.e. non-admixed) Florida panthers, Texas pumas, and F 1 (canonical Florida × Texas) panthers dating from a genetic rescue program and Everglades National Park panthers with Central American ancestry resulting from an earlier admixture event. We found neither genetic drift nor selection significantly reduced overall proportions of deleterious alleles in the severely bottlenecked canonical Florida panthers. Nevertheless, the deleterious alleles identified were distributed into a disproportionately high number of homozygous genotypes due to close inbreeding in this group. Conversely, admixed Florida panthers (either with Texas or Central American ancestry) presented reduced levels of homozygous genotypes carrying deleterious alleles but increased levels of heterozygous genotypes carrying these variants relative to canonical Florida panthers. Although admixture is likely to alleviate the load of standing deleterious variation present in homozygous genotypes, our results suggest that introduced novel deleterious alleles (temporarily present in heterozygous state) in genetically rescued populations could potentially be expressed in subsequent generations if their effective sizes remain small.

Journal of Heredity↗
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

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↗
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↗
USGS2019

Elevated heterozygosity in adults relative to juveniles provides evidence of viability selection on eagles and falcons

Viability selection yields adult populations that are more genetically variable than those of juveniles, producing a positive correlation between heterozygosity and survival. Viability selection could be the result of decreased heterozygosity across many loci in inbred individuals and a subsequent decrease in survivorship resulting from the expression of the deleterious alleles. Alternatively, locus-specific differences in genetic variability between adults and juveniles may be driven by forms of balancing selection, including heterozygote advantage, frequency-dependent selection, or selection across temporal and spatial scales. We use a pooled-sequencing approach to compare genome-wide and locus-specific genetic variability between 74 golden eagle ( Aquila chrysaetos ), 62 imperial eagle ( Aquila heliaca ), and 69 prairie falcon ( Falco mexicanus ) juveniles and adults. Although genome-wide genetic variability is comparable between juvenile and adult golden eagles and prairie falcons, imperial eagle adults are significantly more heterozygous than juveniles. This evidence of viability selection may stem from a relatively smaller imperial eagle effective population size and potentially greater genetic load. We additionally identify ~2000 single-nucleotide polymorphisms across the 3 species with extreme differences in heterozygosity between juveniles and adults. Many of these markers are associated with genes implicated in immune function or olfaction. These loci represent potential targets for studies of how heterozygote advantage, frequency-dependent selection, and selection over spatial and temporal scales influence survivorship in avian species. Overall, our genome-wide data extend previous studies that used allozyme or microsatellite markers and indicate that viability selection may be a more common evolutionary phenomenon than often appreciated.

Journal of Heredity↗