Search NASA⌕ Search

SEARCH · Search NASA

Results for “Effective population size”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Surveying Grassland Islands: the genetics and performance of Appalachian switchgrass (Panicum virgatum L.) collections

The interior Southeastern United States could contain novel germplasm for the bioenergy crop switchgrass due to its diverse habitats and geographic location between genetic subpopulations (Atlantic, Midwest, and Gulf). Collections from this region could accelerate breeding progress, contribute to conservation efforts, and improve understanding of isolated grasslands in the region. This study located 22 sites in the Midsouth region and obtained 1,521,210 single nucleotide polymorphism markers of 202 individuals through genotype-by-sequencing. Here, individuals were evaluated for flowering time, winter survival and tiller number. Comparison to a national diversity panel revealed that branches of two major subpopulations occur in the region with two levels of polyploidy: Atlantic tetraploids and Midwest octoploids. Two locations contained admixed octoploid individuals with Midwest and Gulf genetics. Field performance of the Midwest octoploids conformed with prior reported performance of the Midwest subpopulation, although three sites contained promising late flowering traits. The Atlantic tetraploids had moderate winter survival, short stature, and anomalously early flowering. Atlantic populations mostly occurred in marginal sites and their morphological and flowering time adaptations may be a resource conservation strategy. Demographic inference of historical effective population size variation in a subset of tetraploid locations indicated a widespread recent decline in effective population size. This pattern is consistent with isolation of these switchgrass communities from larger populations and is further supported by evidence of inbreeding within the populations (F I = 0.18). The populations documented in this study contain novel genetic diversity and adaptations to a range of marginal habitats. Therefore, this study provides a new source of germplasm for future breeding and conservation programs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Adaptive gene loss in the common bean pan-genome during range expansion and domestication

The common bean ( Phaseolus vulgaris L.) is a crucial legume crop and an ideal evolutionary model to study adaptive diversity in wild and domesticated populations. Here, we present a common bean pan-genome based on five high-quality genomes and whole-genome reads representing 339 genotypes. It reveals ~234 Mb of additional sequences containing 6,905 protein-coding genes missing from the reference, constituting 49% of all presence/absence variants (PAVs). More non-synonymous mutations are found in PAVs than core genes, probably reflecting the lower effective population size of PAVs and fitness advantages due to the purging effect of gene loss. Our results suggest pan-genome shrinkage occurred during wild range expansion. Selection signatures provide evidence that partial or complete gene loss was a key adaptive genetic change in common bean populations with major implications for plant adaptation. The pan-genome is a valuable resource for food legume research and breeding for climate change mitigation and sustainable agriculture.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Recombination smooths the time-signal disrupted by latency in within-host HIV phylogenies

Within-host HIV evolution involves several features that may disrupt standard phylogenetic reconstruction. One important feature is re-activation of latently integrated provirus, which has the potential to disrupt the temporal signal, leading to variation in the branch lengths and apparent evolutionary rates in a tree. Yet, real within-host HIV phylogenies tend to show clear, ladder-like trees structured by the time of sampling. Another important feature is recombination, which violates the fundamental assumption that evolutionary history can be represented by a single bifurcating tree. Thus, recombination complicates the within-host HIV dynamic by mixing genomes and creating evolutionary loop structures that cannot be represented in bifurcating trees. In this paper, we develop a coalescent-based simulator of within-host HIV evolution that includes latency, recombination, and effective population size dynamics that allows us to study the relationship between the true, complex genealogy of within-host HIV evolution, encoded as an Ancestral Recombination Graph (ARG), and the observed phylogenetic tree. To compare our ARG results to the familiar phylogeny format, we calculate the expected bifurcating tree after decomposing the ARG into all unique site trees, their combined distance matrix, and the overall corresponding bifurcating tree. While latency and recombination separately disrupt the phylogenetic signal, remarkably, we find that recombination recovers the temporal signal of within-host HIV evolution caused by latency by mixing fragments of old, latent genomes into the contemporary population. In effect, recombination averages over extant heterogeneity, whether it stems from mixed time-signals or population bottlenecks. Further, we establish that the signals of latency and recombination can be observed in phylogenetic trees despite being an incorrect representation of the true evolutionary history. Using an Approximate Bayesian Computation method, we develop a set of statistical probes to tune our simulation model to nine longitudinally-sampled within-host HIV phylogenies. Because ARGs are exceedingly difficult to infer from real HIV data, our simulation system allows investigating effects of latency, recombination, and population size bottlenecks by matching decomposed ARGs to real data as observed in standard phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Invasion and rapid adaptation of guppies ( Poecilia reticulata ) across the Hawaiian Archipelago

How much does natural selection, as opposed to genetic drift, admixture, and gene flow, contribute to the evolution of invasive species following introduction to a new environment? Here we assess how evolution can shape biological invasions by examining population genomic variation in non-native guppies (Poecilia reticulata) introduced to the Hawaiian Islands approximately a century ago. By examining 18 invasive populations from four Hawaiian islands and four populations from the native range in northern South America, we reconstructed the history of introductions and evaluated population structure as well as the extent of ongoing gene flow across watersheds and among islands. Patterns of differentiation indicate that guppies have developed significant population structure, with little natural or human-mediated gene flow having occurred among populations following introduction. Demographic modeling and admixture graph analyses together suggest that guppies were initially introduced to O‘ahu and Maui and then translocated to Hawai‘i and Kaua‘i. We detected evidence for only one introduction event from the native range, implying that any adaptive evolution in introduced populations likely utilized the genetic variation present in the founding population. Environmental association tests accounting for population structure identified loci exhibiting signatures of adaptive variation related to predators and landscape characteristics but not nutrient regimes. When paired with high estimates of effective population sizes and detectable population structure, the presence of environment-associated loci supports the role of natural selection in shaping contemporary evolution of Hawaiian guppy populations. Our findings indicate that local adaptation may engender invasion success, particularly in species with life histories that facilitate rapid evolution. Finally, evidence of low gene flow between populations suggests that removal could be an effective approach to control invasive guppies across the Hawaiian archipelago.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Giant kelp genetic monitoring before and after disturbance reveals stable genetic diversity in Southern California

Given the impacts of climate change and other anthropogenic stressors on marine systems, there is a need to accurately predict how species respond to changing environments and disturbance regimes. The use of genetic tools to monitor temporal trends in populations gives ecologists the ability to estimate changes in genetic diversity and effective population size that may be undetectable by traditional census methods. Although multiple studies have used temporal genetic analysis, they usually involve commercially important species, and rarely sample before and after disturbance. In this study, we run a temporal analysis of giant kelp, Macrocystis pyrifera , genetic diversity over the scope of 10 years (2008-2018) using the same microsatellite marker panel to assess the genetic consequences of disturbance in several populations of giant kelp ( Macrocystis pyrifera ) in the Southern California Bight. The study is a rare pre- and post-disturbance microsatellite analysis that included declines to giant kelp caused by the 2015/16 El Nino Southern Oscillation event. We used canopy biomass estimated by remote sensing (Landsat) to quantify the extent of disturbance to kelp beds, and sea surface temperature data to understand how kelp was pushed towards its temperature limits during this period. Despite prolonged periods with decreased canopy at several sites, no changes in genetic structure and allelic richness were observed. We argue that giant kelp in the region is best described as a “patchy population” system where true extinctions are rare. We discuss how deep refugia of subsurface sporophytes and cryptic microscopic life stages could have kept genetic diversity through disturbance. Given the increasing effects of climate change and uncertainty in modeling impacts of species with cryptic life history stages, we suggest further investigation to reveal the role such stages play in species resilience. Genetic monitoring studies of sites selected by remote census demographic and climate surveys should be continued in the future given the predicted impacts of climate change.

Klingbeil, III, William H.↗
DOE OSTI2022

Simulation of sugar kelp ( Saccharina latissima ) breeding guided by practices to accelerate genetic gains

Abstract Though Saccharina japonica cultivation has been established for many decades in East Asian countries, the domestication process of sugar kelp (Saccharina latissima) in the Northeast United States is still at its infancy. In this study, by using data from our breeding experience, we will demonstrate how obstacles for accelerated genetic gain can be assessed using simulation approaches that inform resource allocation decisions. Thus far, we have used 140 wild sporophytes that were sampled in 2018 from the northern Gulf of Maine to southern New England. From these sporophytes, we sampled gametophytes and made and evaluated over 600 progeny sporophytes from crosses among the gametophytes in 2019 and 2020. The biphasic life cycle of kelp gives a great advantage in selective breeding as we can potentially select both on the sporophytes and gametophytes. However, several obstacles exist, such as the amount of time it takes to complete a breeding cycle, the number of gametophytes that can be maintained in the laboratory, and whether positive selection can be conducted on farm-tested sporophytes. Using the Gulf of Maine population characteristics for heritability and effective population size, we simulated a founder population of 1,000 individuals and evaluated the impact of overcoming these obstacles on rate of genetic gain. Our results showed that key factors to improve current genetic gain rely mainly on our ability to induce reproduction of the best farm-tested sporophytes, and to accelerate the clonal vegetative growth of released gametophytes so that enough gametophyte biomass is ready for making crosses by the next growing season. Overcoming these challenges could improve rates of genetic gain more than 2-fold. Future research should focus on conditions favorable for inducing spring reproduction, and on increasing the amount of gametophyte tissue available in time to make fall crosses in the same year.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Intra- and inter-subtype HIV diversity between 1994 and 2018 in southern Uganda: a longitudinal population-based study

There is limited data on human immunodeficiency virus (HIV) evolutionary trends in African populations. We evaluated changes in HIV viral diversity and genetic divergence in southern Uganda over a 24-year period spanning the introduction and scale-up of HIV prevention and treatment programs using HIV sequence and survey data from the Rakai Community Cohort Study, an open longitudinal population-based HIV surveillance cohort. Gag (p24) and env (gp41) HIV data were generated from people living with HIV (PLHIV) in 31 inland semi-urban trading and agrarian communities (1994–2018) and four hyperendemic Lake Victoria fishing communities (2011–2018) under continuous surveillance. HIV subtype was assigned using the Recombination Identification Program with phylogenetic confirmation. Inter-subtype diversity was evaluated using the Shannon diversity index, and intra-subtype diversity with the nucleotide diversity and pairwise TN93 genetic distance. Genetic divergence was measured using root-to-tip distance and pairwise TN93 genetic distance analyses. Demographic history of HIV was inferred using a coalescent-based Bayesian Skygrid model. Evolutionary dynamics were assessed among demographic and behavioral population subgroups, including by migration status. 9931 HIV sequences were available from 4999 PLHIV, including 3060 and 1939 persons residing in inland and fishing communities, respectively. In inland communities, subtype A1 viruses proportionately increased from 14.3% in 1995 to 25.9% in 2017 (P < .001), while those of subtype D declined from 73.2% in 1995 to 28.2% in 2017 (P < .001). The proportion of viruses classified as recombinants significantly increased by nearly four-fold from 12.2% in 1995 to 44.8% in 2017. Inter-subtype HIV diversity has generally increased. While intra-subtype p24 genetic diversity and divergence leveled off after 2014, intra-subtype gp41 diversity, effective population size, and divergence increased through 2017. Intra- and inter-subtype viral diversity increased across all demographic and behavioral population subgroups, including among individuals with no recent migration history or extra-community sexual partners. This study provides insights into population-level HIV evolutionary dynamics following the scale-up of HIV prevention and treatment programs. Continued molecular surveillance may provide a better understanding of the dynamics driving population HIV evolution and yield important insights for epidemic control and vaccine development.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

The western redcedar genome reveals low genetic diversity in a self-compatible conifer

We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Uniformly Ordered Binary Decision Algorithm for Benchmark Experiment Correlations in Whisper Validation

When performing a validation exercise for determining the upper subcritical limit of a nuclear criticality safety application, an analyst should select and perform a statistical analysis on a population of benchmark experiments that are neutronically similar to the application. The size of this population should be sufficiently large such that the statistical analysis has a high degree of confidence that the bias plus bias uncertainty (calculational margin) has been accurately quantified. A complication arises because many benchmark experiments share common components, leading to correlations in their measured effective multiplication factors. Correlations between benchmark experiments within the population reduces its predictive power. This motivates the need for methods that consider benchmark experiment correlations and ensure adequate statistical significance of results. The Whisper code is a statistical analysis pack- age that incorporates nuclear data sensitivity coefficients from MCNP to assess benchmark experiment similarity and then performs an extreme-value analysis to estimate the bias plus bias uncertainty. The original methodology in Whisper does not consider the effect of benchmark experiment correlations when making this estimation, and this summary proposes the uniformly ordered binary decision algorithm to address this shortcoming. The original methodology in Whisper computes similarity coefficients ck for an application compared to all benchmark experiments in its library and develops weighting factors for a selected population proportional to the ck values. The methodology can be interpreted as statistically emulating a validation exercise for a particular application where the weighting factors may be viewed as the likelihood that an analyst would include a particular benchmark experiment within the population. The effective sample size of the population is the expected or mean number of benchmark experiments in the population. The uniformly ordered binary decision algorithm identifies clusters of correlated benchmark experiments within the population and then computes adjusted weighting factors based on the magnitude of the correlation coefficients within the cluster to compute a reduced effective sample size accounting for the lower information content because of correlations. Benchmark experiments within the cluster are ordered randomly with equal probability and probabilistic decisions are made as to whether a benchmark. experiment within the cluster should treated as redundant with a previous one; if two redundant benchmark experiments are included, then the conservative worst case bias plus bias uncertainty is used and the pair is counted as a single benchmark experiment in the population. Results are provided for HEU solutions in a research version of the Whisper software using benchmark experiment correlations provided by DICE, the Database for the International Criticality Safety Benchmark Evaluation Project (ICSBEP). These show that there can be a significant increase in the bias plus bias uncertainty because the effective sample size is reduced, and therefore the algorithm, needing to meet sample size requirements, expands the benchmark experiment population by accepting less similar benchmark experiments that would have otherwise not been included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Origin and evolution of HIV-1 subtype A6

Background: HIV outbreaks in the Former Soviet Union (FSU) countries were characterized by repeated transmission of the HIV variant AFSU, which is now classified as a distinct subtype A sub-subtype called A6. The current study used phylogenetic/phylodynamic and signature mutation analyses to determine likely evolutionary relationship between subtype A6 and other subtype A sub-subtypes. Methods: For this study, an initial Maximum Likelihood phylogenetic analysis was performed using a total of 553 full-length, publicly available, reverse transcriptase sequences, from A1, A2, A3, A4, A5, and A6 sub-subtypes of subtype A. For phylogenetic clustering and signature mutation analysis, a total of 5961 and 3959 pol and env sequences, respectively, were used. Results: Phylogenetic and signature mutation analysis showed that HIV-1 sub-subtype A6 likely originated from sub-subtype A1 of African origin. A6 and A1 pol and env genes shared several signature mutations that indicate genetic similarity between the two subtypes. For A6, tMRCA dated to 1975, 15 years later than that of A1. Conclusion: The current study provides insights into the evolution and diversification of A6 in the backdrop of FSU countries and indicates that A6 in FSU countries evolved from A1 of African origin and is getting bridged outside the FSU region.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Lineage frequency time series reveal elevated levels of genetic drift in SARS-CoV-2 transmission in England

Genetic drift in infectious disease transmission results from randomness of transmission and host recovery or death. The strength of genetic drift for SARS-CoV-2 transmission is expected to be high due to high levels of superspreading, and this is expected to substantially impact disease epidemiology and evolution. However, we don’t yet have an understanding of how genetic drift changes over time or across locations. Furthermore, noise that results from data collection can potentially confound estimates of genetic drift. To address this challenge, we develop and validate a method to jointly infer genetic drift and measurement noise from time-series lineage frequency data. Our method is highly scalable to increasingly large genomic datasets, which overcomes a limitation in commonly used phylogenetic methods. We apply this method to over 490,000 SARS-CoV-2 genomic sequences from England collected between March 2020 and December 2021 by the COVID-19 Genomics UK (COG-UK) consortium and separately infer the strength of genetic drift for pre-B.1.177, B.1.177, Alpha, and Delta. We find that even after correcting for measurement noise, the strength of genetic drift is consistently, throughout time, higher than that expected from the observed number of COVID-19 positive individuals in England by 1 to 3 orders of magnitude, which cannot be explained by literature values of superspreading. Our estimates of genetic drift suggest low and time-varying establishment probabilities for new mutations, inform the parametrization of SARS-CoV-2 evolutionary models, and motivate future studies of the potential mechanisms for increased stochasticity in this system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Thermal transport exceeding bulk heat conduction due to nonthermal micro/nanoscale phonon populations

While classical size effects usually lead to a reduced effective thermal conductivity, we report here that nonthermal phonon populations produced by a micro/nanoscale heat source can lead to enhanced heat conduction, exceeding the prediction from Fourier's law. We study nondiffusive thermal transport by phonons at small distances within the framework of the Boltzmann transport equation (BTE) and demonstrate that the transport is significantly affected by the distribution of phonons emitted by the source. We discuss analytical solutions of the steady-state BTE for a source with a sinusoidal spatial profile, as well as for a three-dimensional Gaussian “hot spot,” and provide numerical results for single crystal silicon at room temperature. If a micro/nanoscale heat source produces a thermal phonon distribution, it gets hotter than that predicted by the heat diffusion equation; however, if the source predominantly produces low-frequency acoustic phonons with long mean free paths, it may get significantly cooler than that predicted by the heat equation, yielding an enhanced heat transport beyond bulk heat conduction.

42 ENGINEERING↗
DOE OSTI2022

Identification of quantitative trait loci for sorghum leaf blight resistance

Sorghum leaf blight and northern corn leaf blight, both caused by Exserohilum turcicum {(Pass.) K. J. Leonard and Suggs [syn. Setosphaeria turcica (Luttr.) K. J. Leonard and Suggs.]}, are major diseases of sorghum [Sorghum bicolor (L.) Moench] and maize (Zea mays L.), respectively. Examining the genetic architecture of resistance in sorghum will lead to a better understanding of the relationship between resistance in sorghum and maize, which can ultimately enhance management options in both crops. In 2018 and 2019, we evaluated two sorghum recombinant inbred line (RIL) populations for resistance to E. turcicum. The BTx623 × IS3620C and BTx623 × SC155 populations consisted of 235 and 81 RILs, respectively. Resistance in both populations was moderately to highly heritable. We identified a total of six quantitative trait loci (QTL) across the two populations. Three QTL with small- to moderate-effect sizes were identified in the BTx623 × IS3620C population. Three QTL, including a large-effect QTL on chromosome three that explained 24% of the variation, were identified in the BTx623 × SC155 population. We compared the identified QTL with the position of northern corn leaf blight candidate genes and found eight candidate resistance gene orthologs that colocalize with the sorghum leaf blight QTL. There were also several nucleotide-binding leucine-rich repeat encoding genes within the candidate intervals. Understanding host resistance in multiple species furthers our understanding of the Exserohilum turcicum patho-system.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Group size mediates effects of intraspecific competition and forest structure on productivity in a recovering social woodpecker population

Conservation of endangered social wildlife in disturbance-prone forests is challenging because direct and indirect effects of management strategies developed at the time of species’ listing when population density is low may change under high-density conditions in locally recovered populations. Here, we used piecewise structural equation modeling to evaluate direct and indirect drivers of productivity in the federally endangered cooperatively breeding red-cockaded woodpecker Dryobates borealis (RCW) on Savannah River Site, South Carolina, USA. We estimated direct and indirect relationships among group size, neighboring group sizes, fledgling production, density of cavity tree clusters occupied by RCWs, area satisfying threshold criteria of ≥22 stems ha –1 of pines ≥35.6 cm diameter at breast height (dbh), <1.4 m 2 ha –1 basal area (BA) of hardwoods 7.6–22.9 cm dbh, and <6% hardwood canopy cover, and area treated with prescribed fire, and tested whether group size mediated indirect effects of area satisfying threshold criteria on fledgling production. Increases in area with ≥22 stems ha –1 of pines ≥35.6 cm dbh and <1.4 m 2 ha –1 BA of hardwoods 7.6–22.9 cm dbh, and area treated with prescribed fire, but not area with <6% hardwood canopy cover, had direct positive effects on group size. Group size and area treated with prescribed fire, but not area satisfying threshold criteria, had direct positive effects on fledgling production. The direct effect of neighboring group sizes on fledgling production was negative and smaller relative to the direct positive effect of group size on fledgling production. Overall, our results indicate positive direct effects of group size on fledgling production outweighed negative direct effects of neighboring group sizes, and that group size mediated positive indirect effects of area satisfying structural threshold criteria on fledgling production. These findings indicate that ongoing forest management aimed to increase area with ≥22 pines ha –1 ≥35.6 cm dbh and <1.4 m 2 ha –1 BA of hardwoods 7.6–22.9 cm dbh will promote large group sizes, which in turn improve fledgling production and offset costs of heightened competition with neighboring groups under high-density conditions. Additionally, positive effects of area treated with prescribed fire on RCW group size and fledgling production indicate prescribed fire has unique contributions to woodpecker productivity, likely via direct effects on forest structure and potentially indirect effects on arthropod prey available to foraging RCWs. Furthermore, by simultaneously accounting for multiple drivers of productivity in social wildlife, our study contributes to the understanding of how increases in social wildlife population sizes can alter previously documented habitat-fitness relationships.

prescribed fire↗
DOE OSTI2024

Body size and early marine conditions drive changes in Chinook salmon productivity across northern latitude ecosystems

Disentangling the influences of climate change from other stressors affecting the population dynamics of aquatic species is particularly pressing for northern latitude ecosystems, where climate‐driven warming is occurring faster than the global average. Chinook salmon (Oncorhynchus tshawytscha) in the Yukon‐Kuskokwim (YK) region occupy the northern extent of their species' range and are experiencing prolonged declines in abundance resulting in fisheries closures and impacts to the well‐being of Indigenous people and local communities. These declines have been associated with physical (e.g., temperature, streamflow) and biological (e.g., body size, competition) conditions, but uncertainty remains about the relative influence of these drivers on productivity across populations and how salmon–environment relationships vary across watersheds. To fill these knowledge gaps, we estimated the effects of marine and freshwater environmental indicators, body size, and indices of competition, on the productivity (adult returns‐per‐spawner) of 26 Chinook salmon populations in the YK region using a Bayesian hierarchical stock‐recruitment model. Across most populations, productivity declined with smaller spawner body size and sea surface temperatures that were colder in the winter and warmer in the summer during the first year at sea. Decreased productivity was also associated with above average fall maximum daily streamflow, increased sea ice cover prior to juvenile outmigration, and abundance of marine competitors, but the strength of these effects varied among populations. Maximum daily stream temperature during spawning migration had a nonlinear relationship with productivity, with reduced productivity in years when temperatures exceeded thresholds in main stem rivers. These results demonstrate for the first time that well‐documented declines in body size of YK Chinook salmon were associated with declining population productivity, while taking climate into account.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Spatially Explicit Individual-Based Modeling Approach to Evaluate the Cumulative Effects of Wind Energy Development on the Greater Sage-Grouse: Sensitivity Analysis and Validation of Key Parameters

To address the need for an approach to evaluate cumulative ecological impacts of wind energy development, particularly those on critical wildlife habitats and species, the authors developed a proof-of-concept, landscape-based, spatially explicit individual-based modeling (IBM) framework for populations of the greater sage-grouse (Centrocercus urophasianus) in Albany County, Wyoming, based on published and other available information on the life history of the species. This sage-grouse IBM estimates the spatial and temporal movements of sage-grouse based on modeled effects of habitat characteristics and developments on sage-grouse condition, survivorship, and reproduction. After completing the initial version of the model, we conducted a series of sensitivity analyses of key model parameters and validation of lek occurrences and spatial distributions to better understand performance of the model relative to observed values of important life history characteristics. The parameter sensitivity analysis was conducted as a series of one-factor-at-a-time tests in which the values of selected parameters were individually varied to evaluate the corresponding effects on model output. We examined the model’s response to changes in parameter values to gain insight into the overall robustness of the model, identify key driver parameters, and initiate qualitative validation. Five functions selected for sensitivity analyses were (1) the rate of change in body condition in relationship to habitat suitability index (HSI) of occupied habitat (condition change rate), (2) the effect of competition between individuals with overlapping home ranges expressed as a percentage reduction in HSI (competition factor), (3) value of condition at which survivorship probability is equal to zero (zero survivorship condition), (4) the relationship between female body condition and clutch size, and (5) the relationship between female body condition and nest success. We varied the parameter values to be tested for the sensitivity analysis around the original value used in the model (i.e., tested values were higher and lower than the model’s original value) to examine the effect of this variation on the predicted output of the model. Model outputs evaluated in the sensitivity analysis included (1) predicted age-class distribution (as measured by the percentage of immature individuals within the population); (2) predicted sex ratio; (3) predicted age-sex class distribution (percent of population in each of six age-sex classes, i.e., male and female juveniles, yearlings, and adults); (4) lifespan; (5) population size; and (6) spatial distribution. Of the five parameters analyzed, condition change rate, survivorship, clutch size, and nest success did not substantially affect any of the evaluated model outputs. For the predicted age-class distribution, the percentage of immature individuals predicted by the model for these four parameters ranged from 61.9% to 65.0%, which was reasonable with respect to the reference values of 51.4 to 57.8% in published studies. The ratio of females to males ranged from 1.09 to 1.14, which was comparable to the reference values of 1.2 to 3.0 presented in an existing study. The predicted mean lifespan ranged from 1.62 to 1.74 years, which appears to be reasonable with respect to the reference values of 0.9 to 1.1 years for sharp-tail grouse and greater prairie chickens presented by published studies. The predicted population size ranged from 1,877 to 3,140. This prediction would be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county. The predicted spatial distribution was not substantially affected by the parameter values examined. In our sensitivity analyses, the only tested parameter that had a noticeable effect on model output was competition factor. Over the range of values tested for this parameter, predicted population size exhibited the largest range (from 1,833 to 4,737) of all parameters tested. As the competition factor increased, the predicted population size decreased. This reflects the effect competition had on the number of individuals that could co occur in high-quality habitat patches. Despite this apparently significant effect, the range of values appears to be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county.

17 WIND ENERGY↗
DOE OSTI2023

Genetic variation associated with adult migration timing in lineages of Steelhead and Chinook Salmon in the Columbia River

Abstract With the discovery of a major effect region (GREB1L, ROCK1) for adult migration timing in genomes of both Chinook Salmon and Steelhead, several subsequent studies have investigated the effect size and distribution of early and late migration alleles among populations in the Columbia River. Here, we synthesize the results of these studies for the major lineages of Chinook Salmon and Steelhead that include highly distinct groups in the interior Columbia River that exhibit atypical life histories from most coastal lineage populations of these two species. Whole‐genome studies with high marker density have provided extensive insight into SNPs most associated with adult migration timing, and suites of markers for each species have been genotyped in large numbers of individuals to further validate phenotypic effects. For Steelhead, the largest phenotypic effect sizes have been observed in the coastal lineage (36% of variation for passage timing at Bonneville Dam; 43% of variation for tributary arrival timing) compared to the inland lineage (7.5% of variation for passage timing at Bonneville Dam; 8.4% of variation for tributary arrival timing) that overwinter in freshwater prior to spawning. For Chinook Salmon, large effect sizes have been observed in all three lineages for multiple adult migration phenotypes (Coastal lineage: percentage of variation of 27.9% for passage timing at Bonneville Dam, 28.7% for arrival timing for spawning; Interior ocean type: percentage of variation of 47.6% for passage timing at Bonneville Dam, 39.6% for tributary arrival timing, 77.9% for arrival timing for spawning; Interior stream type: percentage of variation of 35.3% for passage at Bonneville Dam, 9.8% for tributary arrival timing, 4.7% for arrival timing for spawning). Together, these results have extended our understanding of genetic variation associated with life history diversity in distinct populations of the Columbia River, however, much research remains necessary to determine the causal mechanism for this major effect region on migration timing in these species.

Narum, Shawn R.↗