Next-generation data filtering in the genomics era.

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

Opportunities and challenges of macrogenetic studies.

The rapidly emerging field of macrogenetics focuses on analysing publicly accessible genetic datasets from thousands of species to explore large-scale patterns and predictors of intraspecific genetic variation. Facilitated by advances in evolutionary biology, technology, data infrastructure, statistics and open science, macrogenetics addresses core evolutionary hypotheses (such as disentangling environmental and life-history effects on genetic variation) with a global focus. Yet, there are important, often overlooked, limitations to this approach and best practices need to be considered and adopted if macrogenetics is to continue its exciting trajectory and reach its full potential in fields such as biodiversity monitoring and conservation. Here, we review the history of this rapidly growing field, highlight knowledge gaps and future directions, and provide guidelines for further research.

Natural resistance to worms exacerbates bovine tuberculosis severity independently of worm coinfection.

Pathogen interactions arising during coinfection can exacerbate disease severity, for example when the immune response mounted against one pathogen negatively affects defense of another. It is also possible that host immune responses to a pathogen, shaped by historical evolutionary interactions between host and pathogen, may modify host immune defenses in ways that have repercussions for other pathogens. In this case, negative interactions between two pathogens could emerge even in the absence of concurrent infection. Parasitic worms and tuberculosis (TB) are involved in one of the most geographically extensive of pathogen interactions, and during coinfection worms can exacerbate TB disease outcomes. Here, we show that in a wild mammal natural resistance to worms affects bovine tuberculosis (BTB) severity independently of active worm infection. We found that worm-resistant individuals were more likely to die of BTB than were nonresistant individuals, and their disease progressed more quickly. Anthelmintic treatment moderated, but did not eliminate, the resistance effect, and the effects of resistance and treatment were opposite and additive, with untreated, resistant individuals experiencing the highest mortality. Furthermore, resistance and anthelmintic treatment had nonoverlapping effects on BTB pathology. The effects of resistance manifested in the lungs (the primary site of BTB infection), while the effects of treatment manifested almost entirely in the lymph nodes (the site of disseminated disease), suggesting that resistance and active worm infection affect BTB progression via distinct mechanisms. Our findings reveal that interactions between pathogens can occur as a consequence of processes arising on very different timescales.

Harnessing the power of RADseq for ecological and evolutionary genomics.

High-throughput techniques based on restriction site-associated DNA sequencing (RADseq) are enabling the low-cost discovery and genotyping of thousands of genetic markers for any species, including non-model organisms, which is revolutionizing ecological, evolutionary and conservation genetics. Technical differences among these methods lead to important considerations for all steps of genomics studies, from the specific scientific questions that can be addressed, and the costs of library preparation and sequencing, to the types of bias and error inherent in the resulting data. In this Review, we provide a comprehensive discussion of RADseq methods to aid researchers in choosing among the many different approaches and avoiding erroneous scientific conclusions from RADseq data, a problem that has plagued other genetic marker types in the past.

Landscape Connectivity and Genetic Structure in a Mainstem and a Tributary Stonefly (Plecoptera) Species Using a Novel Reference Genome.

Understanding how environmental variation influences population genetic structure can help predict how environmental change influences population connectivity, genetic diversity, and evolutionary potential. We used riverscape genomics modeling to investigate how climatic and habitat variables relate to patterns of genetic variation in 2 stonefly species, one from mainstem river habitats (Sweltsa coloradensis) and one from tributaries (Sweltsa fidelis) in 40 sites in northwest Montana, USA. We produced a draft genome assembly for S. coloradensis (N50 = 0.251 Mbp, BUSCO > 95% using "insecta_ob9" reference genes). We genotyped 1930 SNPs in 372 individuals for S. coloradensis and 520 SNPs in 153 individuals for S. fidelis. We found higher genetic diversity for S. coloradensis compared to S. fidelis, but nearly identical genetic differentiation among sites within each species (both had global loci median FST = 0.000), despite differences in stream network location. For landscape genomics and testing for selection, we produced a less stringently filtered data set (3454 and 1070 SNPs for S. coloradensis and S. fidelis, respectively). Environmental variables (mean summer precipitation, slope, aspect, mean June stream temperature, land cover type) were correlated with 19 putative adaptive loci for S. coloradensis, but there was only one putative adaptive locus for S. fidelis (correlated with aspect). Interestingly, we also detected potential hybridization between multiple Sweltsa species which has never been previously detected. Studies like ours, that test for adaptive variation in multiple related species are needed to help assess landscape connectivity and the vulnerability of populations and communities to environmental change.

Broad maternal geographic origin of domestic sheep in Anatolia and the Zagros.

We investigated the controversial origin of domestic sheep (Ovis aries) using large samples of contemporary and ancient domestic individuals and their closest wild relatives: the Asiatic mouflon (Ovis gmelini), the urial (Ovis vignei) and the argali (Ovis ammon). A phylogeny based on mitochondrial DNA, including 213 new cytochrome-b sequences of wild Ovism confirmed that O. gmelini is the maternal ancestor of sheep and precluded mtDNA contributions from O. vignei (and O. gmelini × O. vignei hybrids) to domestic lineages. We also produced 54 new control region sequences showing shared haplogroups (A, B, C and E) between domestic sheep and wild O. gmelini which localized the domestication center in eastern Anatolia and central Zagros, excluding regions further east where exclusively wild haplogroups were found. This overlaps with the geographic distribution of O. gmelini gmelini, further suggesting that the maternal origin of domestic sheep derives from this subspecies. Additionally, we produced 57 new CR sequences of Neolithic sheep remains from a large area covering Anatolia to Europe, showing the early presence of at least three mitochondrial haplogroups (A, B and D) in Western colonization routes. This confirmed that sheep domestication was a large-scale process that captured diverse maternal lineages (haplogroups).

Macrogenetic studies must not ignore limitations of genetic markers and scale.

Millette et al. (Ecology Letters, 2020, 23:55-67) reported no consistent worldwide anthropogenic effects on animal genetic diversity using repurposed mitochondrial DNA sequences. We reexamine data from this study, describe genetic marker and scale limitations which might lead to misinterpretations with conservation implications, and provide advice to improve future macrogenetic studies.

The Genetic Architecture of Fitness Drives Population Viability during Rapid Environmental Change.

AbstractThe rapid global loss of biodiversity calls for improved predictions of how populations will evolve and respond demographically to ongoing environmental change. The heritability (h2) of selected traits has long been known to affect evolutionary and demographic responses to environmental change. However, effects of the genetic architecture underlying the h2 of a selected trait on population responses to selection are less well understood. We use deterministic models and stochastic simulations to show that the genetic architecture underlying h2 can dramatically affect population viability during environmental change. Polygenic trait architectures (many loci, each with a small phenotypic effect) conferred higher population viability than genetic architectures with the same initial h2 and large-effect loci under a wide range of scenarios. Population viability also depended strongly on the initial frequency of large-effect beneficial alleles, with moderately low initial allele frequencies conferring higher viability than rare or already-frequent large-effect alleles. Greater population viability associated with polygenic architectures appears to be due to higher short-term evolutionary potential compared with architectures with large-effect loci. These results suggest that integrating information on the trait genetic architecture into quantitative genetic and population viability analysis will substantially improve our understanding and prediction of evolutionary and demographic responses following environmental change.

Evaluating the outcomes of genetic rescue attempts.

Augmenting gene flow is a powerful tool for the conservation of small, isolated populations. However, genetic rescue attempts have largely been limited to populations at the brink of extinction, in part due to concerns over negative outcomes (e.g., outbreeding depression). Increasing habitat fragmentation may necessitate more proactive genetic management. Broader application of augmented gene flow will, in turn, require rigorous evaluation to increase confidence and identify pitfalls in this approach. To date, there has been no assessment of best monitoring practices for genetic rescue attempts. We used genomically explicit, individual-based simulations to examine the effectiveness of common approaches (i.e., tests for increases in fitness, migrant ancestry, heterozygosity, and abundance) for determining whether genetic rescue or outbreeding depression occurred. Statistical power to detect the effects of gene flow on fitness was high (≥0.8) when effect sizes were large, a finding consistent with those from previous studies on severely inbred populations. However, smaller effects of gene flow on fitness can appreciably affect persistence probability but current evaluation approaches fail to provide results from which reliable inferences can be drawn. The power of the metrics we examined to evaluate genetic rescue attempts depended on the time since gene flow and whether gene flow was beneficial or deleterious. Encouragingly, the use of multiple metrics provided nonredundant information and improved inference reliability, highlighting the importance of intensive monitoring efforts. Further development of best practices for evaluating genetic rescue attempts will be crucial for a responsible transition to increased use of translocations to decrease extinction risk.

Evaluación de los Resultados de los Intentos de Rescate Genético Resumen El aumento del flujo génico es una herramienta poderosa para la conservación de poblaciones pequeñas y aisladas. Sin embargo, los intentos de rescate genético en su mayoría se han limitado a las poblaciones que se encuentran al borde de la extinción, en parte debido a la preocupación que existe por los resultados negativos (es decir, la depresión exogámica). La creciente fragmentación del hábitat puede requerir un manejo genético más proactivo. La aplicación más extensa del flujo génico aumentado requerirá a su vez una evaluación rigurosa para incrementar la confianza e identificar las dificultades de esta estrategia. A la fecha, no ha habido una evaluación de las mejores prácticas de monitoreo para los intentos de rescate genético. Usamos simulaciones explícitas basadas en individuos para examinar la efectividad de las estrategias comunes (es decir, análisis del incremento en adaptabilidad, ascendencia migratoria, heterocigosidad y abundancia) para determinar si ocurrió el rescate genético o la depresión exogámica. El poder estadístico para detectar los efectos del flujo génico sobre la adaptabilidad fue elevado (≥0.8) cuando el tamaño de los efectos fue grande, un hallazgo consistente con aquellos realizados en estudios previos sobre poblaciones con una endogamia severa. Sin embargo, los efectos menores del flujo génico sobre la adaptabilidad pueden afectar de manera apreciable la probabilidad de persistencia, pero las estrategias actuales de evaluación no proporcionan resultados de los cuales se puedan hacer inferencias confiables. El poder de las medidas que examinamos para evaluar los intentos de rescate genético dependió del tiempo desde que inició el flujo génico y de si el flujo génico fue benéfico o perjudicial. De manera alentadora, el uso de múltiples medidas proporcionó información no redundante y mejoró la confiabilidad de la inferencia, resaltando así la importancia de los esfuerzos intensivos de monitoreo. El futuro desarrollo de mejores prácticas para la evaluación de los intentos de rescate genético será de suma importancia para la transición responsable hacia el mayor uso de reubicaciones para reducir el riesgo de extinción.

Big Data in Conservation Genomics: Boosting Skills, Hedging Bets, and Staying Current in the Field.

A current challenge in the fields of evolutionary, ecological, and conservation genomics is balancing production of large-scale datasets with additional training often required to handle such datasets. Thus, there is an increasing need for conservation geneticists to continually learn and train to stay up-to-date through avenues such as symposia, meetings, and workshops. The ConGen meeting is a near-annual workshop that strives to guide participants in understanding population genetics principles, study design, data processing, analysis, interpretation, and applications to real-world conservation issues. Each year of ConGen gathers a diverse set of instructors, students, and resulting lectures, hands-on sessions, and discussions. Here, we summarize key lessons learned from the 2019 meeting and more recent updates to the field with a focus on big data in conservation genomics. First, we highlight classical and contemporary issues in study design that are especially relevant to working with big datasets, including the intricacies of data filtering. We next emphasize the importance of building analytical skills and simulating data, and how these skills have applications within and outside of conservation genetics careers. We also highlight recent technological advances and novel applications to conservation of wild populations. Finally, we provide data and recommendations to support ongoing efforts by ConGen organizers and instructors-and beyond-to increase participation of underrepresented minorities in conservation and eco-evolutionary sciences. The future success of conservation genetics requires both continual training in handling big data and a diverse group of people and approaches to tackle key issues, including the global biodiversity-loss crisis.

Tolerance of aquifer stoneflies to repeated hypoxia exposure and oxygen dynamics in an alluvial aquifer.

Aquatic insects cope with hypoxia and anoxia using a variety of behavioral and physiological responses. Most stoneflies (Plecoptera) occur in highly oxygenated surface waters, but some species live underground in alluvial aquifers containing heterogeneous oxygen concentrations. Aquifer stoneflies appear to be supported by methane-derived food resources, which they may exploit using anoxia-resistant behaviors. We documented dissolved oxygen dynamics and collected stoneflies over 5 years in floodplain wells of the Flathead River, Montana. Hypoxia regularly occurred in two wells, and nymphs of Paraperla frontalis were collected during hypoxic periods. We measured mass-specific metabolic rates (MSMRs) at different oxygen concentrations (12, 8, 6, 4, 2, 0.5 mg l-1, and during recovery) for 111 stonefly nymphs to determine whether aquifer and benthic taxa differed in hypoxia tolerance. Metabolic rates of aquifer taxa were similar across oxygen concentrations spanning 2 to 12 mg l-1 (P>0.437), but the MSMRs of benthic taxa dropped significantly with declining oxygen (P<0.0001; 2.9-times lower at 2 vs. 12 mg l-1). Aquifer taxa tolerated short-term repeated exposure to extreme hypoxia surprisingly well (100% survival), but repeated longer-term (>12 h) exposures resulted in lower survival (38-91%) and lower MSMRs during recovery. Our work suggests that aquifer stoneflies have evolved a remarkable set of behavioral and physiological adaptations that allow them to exploit the unique food resources available in hypoxic zones. These adaptations help to explain how large-bodied consumers might thrive in the underground aquifers of diverse and productive river floodplains.

Population Genomics Training for the Next Generation of Conservation Geneticists: ConGen 2018 Workshop.

The increasing availability and complexity of next-generation sequencing (NGS) data sets make ongoing training an essential component of conservation and population genetics research. A workshop entitled "ConGen 2018" was recently held to train researchers in conceptual and practical aspects of NGS data production and analysis for conservation and ecological applications. Sixteen instructors provided helpful lectures, discussions, and hands-on exercises regarding how to plan, produce, and analyze data for many important research questions. Lecture topics ranged from understanding probabilistic (e.g., Bayesian) genotype calling to the detection of local adaptation signatures from genomic, transcriptomic, and epigenomic data. We report on progress in addressing central questions of conservation genomics, advances in NGS data analysis, the potential for genomic tools to assess adaptive capacity, and strategies for training the next generation of conservation genomicists.

Evidence for adaptive introgression of exons across a hybrid swarm in deer.

BACKGROUND: Secondary contact between closely related lineages can result in a variety of outcomes, including hybridization, depending upon the strength of reproductive barriers. By examining the extent to which different parts of the genome introgress, it is possible to infer the strength of selection and gain insight into the evolutionary trajectory of lineages. Following secondary contact approximately 8000 years ago in the Pacific Northwest, mule deer (Odocoileus hemionus hemionus) and black-tailed deer (O. h. columbianus) formed a hybrid swarm along the Cascade mountain range despite substantial differences in body size (up to two times) and habitat preference. In this study, we examined genetic population structure, extent of introgression, and selection pressures in freely interbreeding populations of mule deer and black-tailed deer using mitochondrial DNA sequences, 9 microsatellite loci, and 95 SNPs from protein-coding genes. RESULTS: We observed bi-directional hybridization and classified approximately one third of the 172 individuals as hybrids, almost all of which were beyond the F1 generation. High genetic differentiation between black-tailed deer and mule deer at protein-coding genes suggests that there is positive divergent selection, though selection on these loci is relatively weak. Contrary to predictions, there was not greater selection on protein-coding genes thought to be associated with immune function and mate choice. Geographic cline analyses were consistent across genetic markers, suggesting long-term stability (over hundreds of generations), and indicated that the center of the hybrid swarm is 20-30 km to the east of the Cascades ridgeline, where there is a steep ecological transition from wet, forested habitat to dry, scrub habitat. CONCLUSIONS: Our data are consistent with a genetic boundary between mule deer and black-tailed deer that is porous but maintained by many loci under weak selection having a substantial cumulative effect. The absence of clear reproductive barriers and the consistent centering of geographic clines at a sharp ecotone suggests that ecology is a driver of hybrid swarm dynamics. Adaptive introgression in this study (and others) promotes gene flow and provides valuable insight into selection strength on specific genes and the evolutionary trajectory of hybridizing taxa.

Assessing thermal adaptation using family-based association and FST outlier tests in a threatened trout species.

Discovering genetic markers associated with phenotypic or ecological characteristics can improve our understanding of adaptation and guide conservation of key evolutionary traits. The Lahontan cutthroat trout (Oncorhynchus clarkii henshawi) of the northern Great Basin Desert, USA, demonstrated exceptional tolerance to high temperatures in the desert lakes where it resided historically. This trait is central to a conservation hatchery effort to protect the genetic legacy of the nearly extinct lake ecotype. We genotyped full-sibling families from this conservation broodstock and samples from the only two remaining, thermally distinct, native lake populations at 4,644 new single nucleotide polymorphisms (SNPs). Family-based genome-wide association testing of the broodstock identified nine and 26 SNPs associated with thermal tolerance (p < 0.05 and p < 0.1), measured in a previous thermal challenge experiment. Genes near the associated SNPs had complex functions related to immunity, growth, metabolism and ion homeostasis. Principal component analysis using the thermotolerance-related SNPs showed unexpected divergence between the conservation broodstock and the native lake populations at these loci. FST outlier tests on the native lake populations identified 18 loci shared between two or more of the tests, with two SNPs identified by all three tests (p < 0.01); none overlapped with loci identified by association testing in the broodstock. A recent history of isolation and the complex genetic and demographic backgrounds of Lahontan cutthroat trout probably limited our ability to find shared thermal tolerance loci. Our study extends the still relatively rare application of genomic tools testing for markers associated with important phenotypic or environmental characteristics in species of conservation concern.

Signatures of polygenic adaptation associated with climate across the range of a threatened fish species with high genetic connectivity.

Adaptive differences across species' ranges can have important implications for population persistence and conservation management decisions. Despite advances in genomic technologies, detecting adaptive variation in natural populations remains challenging. Key challenges in gene-environment association studies involve distinguishing the effects of drift from those of selection and identifying subtle signatures of polygenic adaptation. We used paired-end restriction site-associated DNA sequencing data (6,605 biallelic single nucleotide polymorphisms; SNPs) to examine population structure and test for signatures of adaptation across the geographic range of an iconic Australian endemic freshwater fish species, the Murray cod Maccullochella peelii. Two univariate gene-association methods identified 61 genomic regions associated with climate variation. We also tested for subtle signatures of polygenic adaptation using a multivariate method (redundancy analysis; RDA). The RDA analysis suggested that climate (temperature- and precipitation-related variables) and geography had similar magnitudes of effect in shaping the distribution of SNP genotypes across the sampled range of Murray cod. Although there was poor agreement among the candidate SNPs identified by the univariate methods, the top 5% of SNPs contributing to significant RDA axes included 67% of the SNPs identified by univariate methods. We discuss the potential implications of our findings for the management of Murray cod and other species generally, particularly in relation to informing conservation actions such as translocations to improve evolutionary resilience of natural populations. Our results highlight the value of using a combination of different approaches, including polygenic methods, when testing for signatures of adaptation in landscape genomic studies.

Legacy introductions and climatic variation explain spatiotemporal patterns of invasive hybridization in a native trout.

Hybridization between invasive and native species, a significant threat to worldwide biodiversity, is predicted to increase due to climate-induced expansions of invasive species. Long-term research and monitoring are crucial for understanding the ecological and evolutionary processes that modulate the effects of invasive species. Using a large, multidecade genetics dataset (N = 582 sites, 12,878 individuals) with high-resolution climate predictions and extensive stocking records, we evaluate the spatiotemporal dynamics of hybridization between native cutthroat trout and invasive rainbow trout, the world's most widely introduced invasive fish, across the Northern Rocky Mountains of the United States. Historical effects of stocking and contemporary patterns of climatic variation were strongly related to the spread of hybridization across space and time. The probability of occurrence, extent of, and temporal changes in hybridization increased at sites in close proximity to historical stocking locations with greater rainbow trout propagule pressure, warmer water temperatures, and lower spring precipitation. Although locations with warmer water temperatures were more prone to hybridization, cold sites were not protected from invasion; 58% of hybridized sites had cold mean summer water temperatures (<11°C). Despite cessation of stocking over 40 years ago, hybridization increased over time at half (50%) of the locations with long-term data, the vast majority of which (74%) were initially nonhybridized, emphasizing the chronic, negative impacts of human-mediated hybridization. These results show that effects of climate change on biodiversity must be analyzed in the context of historical human impacts that set ecological and evolutionary trajectories.

Accounting for adaptive capacity and uncertainty in assessments of species' climate-change vulnerability.

Climate-change vulnerability assessments (CCVAs) are valuable tools for assessing species' vulnerability to climatic changes, yet failure to include measures of adaptive capacity and to account for sources of uncertainty may limit their effectiveness. We took a more comprehensive approach that incorporates exposure, sensitivity, and capacity to adapt to climate change. We applied our approach to anadromous steelhead trout (Oncorhynchus mykiss) and nonanadromous bull trout (Salvelinus confluentus), threatened salmonids within the Columbia River Basin (U.S.A.). We quantified exposure on the basis of scenarios of future stream temperature and flow, and we represented sensitivity and capacity to adapt to climate change with metrics of habitat quality, demographic condition, and genetic diversity. Both species were found to be highly vulnerable to climate change at low elevations and in their southernmost habitats. However, vulnerability rankings varied widely depending on the factors (climate, habitat, demographic, and genetic) included in the CCVA and often differed for the 2 species at locations where they were sympatric. Our findings illustrate that CCVA results are highly sensitive to data inputs and that spatial differences can complicate multispecies conservation. Based on our results, we suggest that CCVAs be considered within a broader conceptual and computational framework and be used to refine hypotheses, guide research, and compare plausible scenarios of species' vulnerability to climate change.

Vive la résistance: genome-wide selection against introduced alleles in invasive hybrid zones.

Evolutionary and ecological consequences of hybridization between native and invasive species are notoriously complicated because patterns of selection acting on non-native alleles can vary throughout the genome and across environments. Rapid advances in genomics now make it feasible to assess locus-specific and genome-wide patterns of natural selection acting on invasive introgression within and among natural populations occupying diverse environments. We quantified genome-wide patterns of admixture across multiple independent hybrid zones of native westslope cutthroat trout and invasive rainbow trout, the world's most widely introduced fish, by genotyping 339 individuals from 21 populations using 9380 species-diagnostic loci. A significantly greater proportion of the genome appeared to be under selection favouring native cutthroat trout (rather than rainbow trout), and this pattern was pervasive across the genome (detected on most chromosomes). Furthermore, selection against invasive alleles was consistent across populations and environments, even in those where rainbow trout were predicted to have a selective advantage (warm environments). These data corroborate field studies showing that hybrids between these species have lower fitness than the native taxa, and show that these fitness differences are due to selection favouring many native genes distributed widely throughout the genome.

Conservation genomics of natural and managed populations: building a conceptual and practical framework.

The boom of massive parallel sequencing (MPS) technology and its applications in conservation of natural and managed populations brings new opportunities and challenges to meet the scientific questions that can be addressed. Genomic conservation offers a wide range of approaches and analytical techniques, with their respective strengths and weaknesses that rely on several implicit assumptions. However, finding the most suitable approaches and analysis regarding our scientific question are often difficult and time-consuming. To address this gap, a recent workshop entitled 'ConGen 2015' was held at Montana University in order to bring together the knowledge accumulated in this field and to provide training in conceptual and practical aspects of data analysis applied to the field of conservation and evolutionary genomics. Here, we summarize the expertise yield by each instructor that has led us to consider the importance of keeping in mind the scientific question from sampling to management practices along with the selection of appropriate genomics tools and bioinformatics challenges.

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics.

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population-specific and pairwise FST ) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate-related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and FST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin-wide to the metapopulation scale). Sensitivity analysis (leave-one-population-out) revealed consistent relationships between climate variables and FST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.