Distinct patterns of inheritance shape life-history traits in steelhead trout.

Life-history variation is the raw material of adaptation, and understanding its genetic and environmental underpinnings is key to designing effective conservation strategies. We used large-scale genetic pedigree reconstruction of anadromous steelhead trout (Oncorhynchus mykiss) from the Russian River, CA, USA, to elucidate sex-specific patterns of life-history traits and their heritability. SNP data from adults returning from sea over a 14-year period were used to identify 13,474 parent-offspring trios. These pedigrees were used to determine age structure, size distributions and family sizes for these fish, as well as to estimate the heritability of two key life-history traits, spawn date and age at maturity (first reproduction). Spawn date was highly heritable (h2 = 0.73) and had a cross-sex genetic correlation near unity. We provide the first estimate of heritability for age at maturity in ocean-going fish from this species and found it to be highly heritable (h2 from 0.29 to 0.62, depending on sex and method), with a much lower genetic correlation across sexes. We also evaluated genotypes at a migration-associated inversion polymorphism and found sex-specific correlations with age at maturity. The significant heritability of these two key reproductive traits in these imperiled fish, and their patterns of inheritance in the two sexes, is consistent with predictions of both natural and sexually antagonistic selection (sexes experience opposing selection pressures). This emphasizes the importance of anthropogenic factors, including hatchery practices and ecosystem modifications, in shaping the fitness of this species, thus providing important guidance for management and conservation efforts.

Low-coverage whole genome sequencing for highly accurate population assignment: Mapping migratory connectivity in the American Redstart (Setophaga ruticilla).

Understanding the geographic linkages among populations across the annual cycle is an essential component for understanding the ecology and evolution of migratory species and for facilitating their effective conservation. While genetic markers have been widely applied to describe migratory connections, the rapid development of new sequencing methods, such as low-coverage whole genome sequencing (lcWGS), provides new opportunities for improved estimates of migratory connectivity. Here, we use lcWGS to identify fine-scale population structure in a widespread songbird, the American Redstart (Setophaga ruticilla), and accurately assign individuals to genetically distinct breeding populations. Assignment of individuals from the nonbreeding range reveals population-specific patterns of varying migratory connectivity. By combining migratory connectivity results with demographic analysis of population abundance and trends, we consider full annual cycle conservation strategies for preserving numbers of individuals and genetic diversity. Notably, we highlight the importance of the Northern Temperate-Greater Antilles migratory population as containing the largest proportion of individuals in the species. Finally, we highlight valuable considerations for other population assignment studies aimed at using lcWGS. Our results have broad implications for improving our understanding of the ecology and evolution of migratory species through conservation genomics approaches.

Limited, asymmetric hybridization between coastal cutthroat trout and steelhead in a Northern California river.

Hybridization between coastal cutthroat trout (Oncorhynchus clarkii clarkii) and steelhead (O. mykiss) was assessed in the Smith River, California. Individuals were categorized as pure or as 1 of 10 hybrid classes using 30 "diagnostic" single-nucleotide polymorphisms positioned on 26 separate chromosomes. Most of the individuals examined (n = 876), were pure coastal cutthroat trout (n = 634) or pure steelhead (n = 213), and 29 individuals were identified as having hybrid ancestry. Among hybrids, first generation hybrids (n = 15) and coastal cutthroat trout backcrosses (n = 12) were the most common. No individuals were identified as backcrosses to SH, suggesting the presence of genetic or behavioral mechanisms constraining such backcrosses, or the growth and survival of their progeny. Mitochondrial DNA of 14 of 15 F1 hybrids was of steelhead origin, suggesting that hybridization was driven primarily by sneak-mating of male coastal cutthroat trout with female steelhead. Evaluation of classical phenotypic characters for coastal cutthroat trout and steelhead (i.e. jaw slash, maxillary length, and hyoid teeth) were not reliable by themselves for identification of either pure parental fish or hybrids. In contrast, analysis with geometric morphometrics revealed distinctive body shapes for pure coastal cutthroat trout and steelhead, and the combination of classical traits and geometric morphology was mostly accurate in distinguishing them. However, first generation hybrids and backcrosses overlapped completely with parental types, highlighting challenges in hybrid identification using phenotypic traits.

Implications of Large-Effect Loci for Conservation: A Review and Case Study with Pacific Salmon.

The increasing feasibility of assembling large genomic datasets for non-model species presents both opportunities and challenges for applied conservation and management. A popular theme in recent studies is the search for large-effect loci that explain substantial portions of phenotypic variance for a key trait(s). If such loci can be linked to adaptations, 2 important questions arise: 1) Should information from these loci be used to reconfigure conservation units (CUs), even if this conflicts with overall patterns of genetic differentiation? 2) How should this information be used in viability assessments of populations and larger CUs? In this review, we address these questions in the context of recent studies of Chinook salmon and steelhead (anadromous form of rainbow trout) that show strong associations between adult migration timing and specific alleles in one small genomic region. Based on the polygenic paradigm (most traits are controlled by many genes of small effect) and genetic data available at the time showing that early-migrating populations are most closely related to nearby late-migrating populations, adult migration differences in Pacific salmon and steelhead were considered to reflect diversity within CUs rather than separate CUs. Recent data, however, suggest that specific alleles are required for early migration, and that these alleles are lost in populations where conditions do not support early-migrating phenotypes. Contrasting determinations under the US Endangered Species Act and the State of California's equivalent legislation illustrate the complexities of incorporating genomics data into CU configuration decisions. Regardless how CUs are defined, viability assessments should consider that 1) early-migrating phenotypes experience disproportionate risks across large geographic areas, so it becomes important to identify early-migrating populations that can serve as reliable sources for these valuable genetic resources; and 2) genetic architecture, especially the existence of large-effect loci, can affect evolutionary potential and adaptability.

Linking climate niches across seasons to assess population vulnerability in a migratory bird.

Global loss of biodiversity has placed new urgency on the need to understand factors regulating species response to rapid environmental change. While specialists are often less resilient to rapid environmental change than generalists, species-level analyses may obscure the extent of specialization when locally adapted populations vary in climate tolerances. Until recently, quantification of the degree of climate specialization in migratory birds below the species level was hindered by a lack of genomic and tracking information, but recent technological advances have helped to overcome these barriers. Here we take a genome-wide genetic approach to mapping population-specific migratory routes and quantifying niche breadth within genetically distinct populations of a migratory bird, the willow flycatcher (Empidonax traillii), which exhibits variation in the severity of population declines across its breeding range. While our sample size is restricted to the number of genetically distinct populations within the species, our results support the idea that locally adapted populations of the willow flycatcher with narrow climatic niches across seasons are already federally listed as endangered or in steep decline, while populations with broader climatic niches have remained stable in recent decades. Overall, this work highlights the value of quantifying niche breadth within genetically distinct groups across time and space when attempting to understand the factors that facilitate or constrain the response of locally adapted populations to rapid environmental change.

Community oncologists' perceptions and utilization of large-panel genomic tumor testing.

PURPOSE: Large-panel genomic tumor testing (GTT) is an emerging technology with great promise but uncertain clinical value. Previous research has documented variability in academic oncologists' perceptions and use of GTT, but little is known about community oncologists' perceptions of GTT and how perceptions relate to clinicians' intentions to use GTT. METHODS: Community oncology physicians (N = 58) participating in a statewide initiative aimed at improving access to large-panel GTT completed surveys assessing their confidence in using GTT, attitudes regarding the value of GTT, perceptions of barriers to GTT implementation, and future intentions to use GTTs. Descriptive and multivariable regression analyses were conducted to characterize these perceptions and to explore the relationships between them. RESULTS: There was substantial variability in clinicians' perceptions of GTT. Clinicians generally had moderate confidence in their ability to use GTT, but lower confidence in patients' ability to understand test results and access targeted treatment. Clinicians had positive attitudes regarding the value of GTT. Clinicians' future intentions to use GTT were associated with greater confidence in using GTT and greater perceived barriers to implementing GTT, but not with attitudes about the value of GTT. CONCLUSIONS: Community oncologists' perceptions of large-panel genomic tumor testing are variable, and their future intentions to use GTT are associated with both their confidence in and perceived barriers to its use, but not with their attitudes towards GTT. More research is needed to understand other factors that determine how oncologists perceive and use GTT in clinical practice.

Polygenic Basis and the Role of Genome Duplication in Adaptation to Similar Selective Environments.

Genetic changes underlying adaptation vary greatly in terms of complexity and, within the same species, genetic responses to similar selective pressures may or may not be the same. We examine both complex (supergene) and simple (SNP) genetic variants occurring in populations of rainbow trout (Oncorhynchus mykiss) independently isolated from ocean access and compared them to each other and to an anadromous below-barrier population representing their ancestral source to search for signatures of both parallel and nonparallel adaptation. All landlocked populations displayed an increased frequency of a large inversion on chromosome Omy05, while 3 of the 4 populations exhibited elevated frequencies of another inversion located on chromosome Omy20. In addition, we identified numerous regions outside these 2 inversions that also show significant shifts in allele frequencies consistent with adaptive evolution. However, there was little concordance among above-barrier populations in these specific genomic regions under selection. In part, the lack of concordance appears to arise from ancestral autopolyploidy in rainbow trout that provides duplicate genomic regions of similar functional composition for selection to act upon. Thus, while selection acting on landlocked populations universally favors the resident ecotype, outside of the major chromosomal inversions, the resulting genetic changes are largely distinct among populations. Our results indicate that selection on standing genetic variation is likely the primary mode of rapid adaptation, and that both supergene complexes and individual loci contribute to adaptive evolution, further highlighting the diversity of adaptive genomic variation involved in complex phenotypic evolution.

The influence of uncertainty and uncertainty tolerance on attitudes and self-efficacy about genomic tumor testing.

Novel medical technologies, like large-panel genomic tumor testing (GTT), offer great promise but also substantial uncertainty regarding their clinical value and appropriate use. The goal of this study was to understand how clinicians' perceived uncertainty about GTT, and uncertainty tolerance (UT), a construct that describes trait-level differences in individuals' responses to uncertainty, influence attitudes and self-efficacy regarding GTT. Community-based oncologists participating in a study of large-panel GTT completed surveys assessing their perceptions of uncertainty about GTT, and their attitudes and self-efficacy regarding GTT. Multivariable regression analyses examined the relationship between oncologists' perceived uncertainty of GTT and their GTT-related attitudes and self-efficacy, and the potential moderating effect of individual differences in UT. Fifty-seven oncologists completed surveys. Greater perceived uncertainty about GTT was associated with more negative attitudes towards it. This association was moderated by UT, such that lower UT was associated with a stronger negative relationship between perceived uncertainty and attitudes. That is, oncologists who perceive GTT as uncertain, tended to have more negative attitudes, particularly if they were low in the trait of uncertainty tolerance. More research is warranted to understand how uncertainty and uncertainty tolerance influence clinicians' responses to GTT and other novel medical interventions.

Dispersal of a nearshore marine fish connects marine reserves and adjacent fished areas along an open coast.

Marine species with pelagic larvae typically exhibit little population structure, suggesting long-distance dispersal and high gene flow. Directly quantifying dispersal of marine fishes is challenging but important, particularly for the design of marine protected areas (MPAs). Here, we studied kelp rockfish (Sebastes atrovirens) sampled along ~25 km of coastline in a boundary current-dominated ecosystem and used genetic parentage analysis to identify dispersal events and characterize them, because the distance between sedentary parents and their settled offspring is the lifetime dispersal distance. Large sample sizes and intensive sampling are critical for increasing the likelihood of detecting parent-offspring matches in such systems and we sampled more than 6,000 kelp rockfish and analysed them with a powerful set of 96 microhaplotype markers. We identified eight parent-offspring pairs with high confidence, including two juvenile fish that were born inside MPAs and dispersed to areas outside MPAs, and four fish born in MPAs that dispersed to nearby MPAs. Additionally, we identified 25 full-sibling pairs, which occurred throughout the sampling area and included all possible combinations of inferred dispersal trajectories. Intriguingly, these included two pairs of young-of-the-year siblings with one member each sampled in consecutive years. These sibling pairs suggest monogamy, either intentional or accidental, which has not been previously demonstrated in rockfishes. This study provides the first direct observation of larval dispersal events in a current-dominated ecosystem and direct evidence that larvae produced within MPAs are exported both to neighbouring MPAs and to proximate areas where harvest is allowed.

Vegetarians' and omnivores' affective and physiological responses to images of food.

Many vegetarians report that meat is unpleasant, but little else is known about their affective responses to meat and non-meat foods. Here we explored affective responses to food images in vegetarians and omnivores and tested the hypothesis that vegetarians have global differences in affective processing (e.g., increased disgust sensitivity). We presented pictures of different food items and recorded participants' affective experience while we recorded peripheral physiology. We found that vegetarians' self-reported experience of meat meal images was less pleasant than omnivores', but that other food images were equally pleasant across the two groups. Moreover, vegetarians and omnivores had strikingly similar physiological responses to all food images - including meat meals. We interpret these results from a psychological constructionist perspective, which posits that individuals conceptualize changes in their bodily states in ways that match their beliefs, such that increased sympathetic nervous system activity may be conceptualized as an experience of excitement about a delicious meat meal for omnivores but as an experience of displeasure for a vegetarian who believes meat is cruel, wasteful, impure, or unhealthy. This interpretation is consistent with emerging neuroscience evidence that the brain constructs experience by predicting and making meaning of internal sensations based on past experience and knowledge.

Ecological genomics predicts climate vulnerability in an endangered southwestern songbird.

Few regions have been more severely impacted by climate change in the USA than the Desert Southwest. Here, we use ecological genomics to assess the potential for adaptation to rising global temperatures in a widespread songbird, the willow flycatcher (Empidonax traillii), and find the endangered desert southwestern subspecies (E. t. extimus) most vulnerable to future climate change. Highly significant correlations between present abundance and estimates of genomic vulnerability - the mismatch between current and predicted future genotype-environment relationships - indicate small, fragmented populations of the southwestern willow flycatcher will have to adapt most to keep pace with climate change. Links between climate-associated genotypes and genes important to thermal tolerance in birds provide a potential mechanism for adaptation to temperature extremes. Our results demonstrate that the incorporation of genotype-environment relationships into landscape-scale models of climate vulnerability can facilitate more precise predictions of climate impacts and help guide conservation in threatened and endangered groups.

Purging putative siblings from population genetic data sets: a cautionary view.

Interest has surged recently in removing siblings from population genetic data sets before conducting downstream analyses. However, even if the pedigree is inferred correctly, this has the potential to do more harm than good. We used computer simulations and empirical samples of coho salmon to evaluate strategies for adjusting samples to account for family structure. We compared performance in full samples and sibling-reduced samples of estimators of allele frequency (P^), population differentiation (F^ST) and effective population size (N^e). RESULTS: (i) unless simulated samples included large family groups together with a component of unrelated individuals, removing siblings generally reduced precision of P^ and F^ST; (ii) N^e based on the linkage disequilibrium method was largely unbiased using full random samples but became increasingly upwardly biased under aggressive purging of siblings. Under nonrandom sampling (some families over-represented), N^e using full samples was downwardly biased; removing just the right 'Goldilocks' fraction of siblings could produce an unbiased estimate, but this sweet spot varied widely among scenarios; (iii) weighting individuals based on the inferred pedigree (to produce a best linear unbiased estimator, BLUE) maximized precision of P^ when the inferred pedigree was correct but performed poorly when the pedigree was wrong; (iv) a variant of sibling removal that leaves intact small sibling groups appears to be more robust to errors in inferences about family structure. Our results illustrate the complex challenges posed by presence of family structure, suggest that no single optimal solution exists and argue for caution in adjusting population genetic data sets for the presence of putative siblings without fully understanding the consequences.

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens).

The little greenbul, a common rainforest passerine from sub-Saharan Africa, has been the subject of long-term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest-savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise FST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.

Bayesian pedigree inference with small numbers of single nucleotide polymorphisms via a factor-graph representation.

We develop a computational framework for addressing pedigree inference problems using small numbers (80-400) of single nucleotide polymorphisms (SNPs). Our approach relaxes the assumptions, which are commonly made, that sampling is complete with respect to the pedigree and that there is no genotyping error. It relies on representing the inferred pedigree as a factor graph and invoking the Sum-Product algorithm to compute and store quantities that allow the joint probability of the data to be rapidly computed under a large class of rearrangements of the pedigree structure. This allows efficient MCMC sampling over the space of pedigrees, and, hence, Bayesian inference of pedigree structure. In this paper we restrict ourselves to inference of pedigrees without loops using SNPs assumed to be unlinked. We present the methodology in general for multigenerational inference, and we illustrate the method by applying it to the inference of full sibling groups in a large sample (n=1157) of Chinook salmon typed at 95 SNPs. The results show that our method provides a better point estimate and estimate of uncertainty than the currently best-available maximum-likelihood sibling reconstruction method. Extensions of this work to more complex scenarios are briefly discussed.

Comment on 'Bayesian parentage analysis with systematic accountability of genotyping error, missing data and false matching'.

We show the software SOLOMON is improved by using the likelihood ratio instead of an ad hoc statistic. CODE: github.com/eriqande/solidmon/releases/tag/v0.1-bioinformatics

Next-generation sequencing for molecular ecology: a caveat regarding pooled samples.

We develop a model based on the Dirichlet-compound multinomial distribution (CMD) and Ewens sampling formula to predict the fraction of SNP loci that will appear fixed for alternate alleles between two pooled samples drawn from the same underlying population. We apply this model to next-generation sequencing (NGS) data from Baltic Sea herring recently published by (Corander et al., 2013, Molecular Ecology, 2931-2940), and show that there are many more fixed loci than expected in the absence of genetic structure. However, we show through coalescent simulations that the degree of population structure required to explain the fraction of alternatively fixed SNPs is extraordinarily high and that the surplus of fixed loci is more likely a consequence of limited representation of individual gene copies in the pooled samples, than it is of population structure. Our analysis signals that the use of NGS on pooled samples to identify divergent SNPs warrants caution. With pooled samples, it is hard to diagnose when an NGS experiment has gone awry; especially when NGS data on pooled samples are of low read depth with a limited number of individuals, it may be worthwhile to temper claims of unexpected population differentiation from pooled samples, pending verification with more reliable methods or stricter adherence to recommended sampling designs for pooled sequencing e.g. Futschik & Schlötterer 2010, Genetics, 186, 207; Gautier et al., 2013a, Molecular Ecology, 3766-3779). Analysis of the data and diagnosis of problems is easier and more reliable (and can be less costly) with individually barcoded samples. Consequently, for some scenarios, individual barcoding may be preferable to pooling of samples.

A role for migration-linked genes and genomic islands in divergence of a songbird.

Next-generation sequencing has made it possible to begin asking questions about the process of divergence at the level of the genome. For example, recently, there has been a debate around the role of 'genomic islands of divergence' (i.e. blocks of outlier loci) in facilitating the process of speciation-with-gene-flow. The Swainson's thrush, Catharus ustulatus, is a migratory songbird with two genetically distinct subspecies that differ in a number of traits known to be involved in reproductive isolation in birds (plumage coloration, song and migratory behaviour), despite contemporary gene flow along a secondary contact zone. Here, we use RAD-PE sequencing to test emerging hypotheses about the process of divergence at the level of the genome and identify genes and gene regions involved in differentiation in this migratory songbird. Our analyses revealed distinct genomic islands on 15 of the 23 chromosomes and an accelerated rate of divergence on the Z chromosome, one of the avian sex chromosomes. Further, an analysis of loci linked to traits known to be involved in reproductive isolation in songbirds showed that genes linked to migration are significantly more differentiated than expected by chance, but that these genes lie primarily outside the genomic islands. Overall, our analysis supports the idea that genes linked to migration play an important role in divergence in migratory songbirds, but we find no compelling evidence that the observed genomic islands are facilitating adaptive divergence in migratory behaviour.

Mapping migration in a songbird using high-resolution genetic markers.

Neotropic migratory birds are declining across the Western Hemisphere, but conservation efforts have been hampered by the inability to assess where migrants are most limited-the breeding grounds, migratory stopover sites or wintering areas. A major challenge has been the lack of an efficient, reliable and broadly applicable method for measuring the strength of migratory connections between populations across the annual cycle. Here, we show how high-resolution genetic markers can be used to identify genetically distinct groups of a migratory bird, the Wilson's warbler (Cardellina pusilla), at fine enough spatial scales to facilitate assessing regional drivers of demographic trends. By screening 1626 samples using 96 highly divergent single nucleotide polymorphisms selected from a large pool of candidates (~450 000), we identify novel region-specific migratory routes and timetables of migration along the Pacific Flyway. Our results illustrate that high-resolution genetic markers are more reliable, precise and amenable to high throughput screening than previously described intrinsic marking techniques, making them broadly applicable to large-scale monitoring and conservation of migratory organisms.

Genetic parentage reveals the (un)natural history of Central Valley hatchery steelhead.

Populations composed of individuals descended from multiple distinct genetic lineages often feature significant differences in phenotypic frequencies. We considered hatchery production of steelhead, the migratory anadromous form of the salmonid species Oncorhynchus mykiss, and investigated how differences among genetic lineages and environmental variation impacted life history traits. We genotyped 23,670 steelhead returning to the four California Central Valley hatcheries over 9 years from 2011 to 2019, confidently assigning parentage to 13,576 individuals to determine age and date of spawning and rates of iteroparity and repeat spawning within each year. We found steelhead from different genetic lineages showed significant differences in adult life history traits despite inhabiting similar environments. Differences between coastal and Central Valley steelhead lineages contributed to significant differences in age at return, timing of spawning, and rates of iteroparity among programs. In addition, adaptive genomic variation associated with life history development in this species varied among hatchery programs and was associated with the age of steelhead spawners only in the coastal lineage population. Environmental variation likely contributed to variations in phenotypic patterns observed over time, as our study period spanned both a marine heatwave and a serious drought in California. Our results highlight evidence of a strong genetic component underlying known phenotypic differences in life history traits between two steelhead lineages.

Genome-matched treatments and patient outcomes in the Maine Cancer Genomics Initiative (MCGI).

Genomic tumor testing (GTT) is an emerging technology aimed at identifying variants in tumors that can be targeted with genomically matched drugs. Due to limited resources, rural patients receiving care in community oncology settings may be less likely to benefit from GTT. We analyzed GTT results and observational clinical outcomes data from patients enrolled in the Maine Cancer Genomics Initiative (MCGI), which provided access to GTTs; clinician educational resources; and genomic tumor boards in community practices in a predominantly rural state. 1603 adult cancer patients completed enrollment; 1258 had at least one potentially actionable variant identified. 206 (16.4%) patients received a total of 240 genome matched treatments, of those treatments, 64% were FDA-approved in the tumor type, 27% FDA-approved in a different tumor type and 9% were given on a clinical trial. Using Inverse Probability of Treatment Weighting to adjust for baseline characteristics, a Cox proportional hazards model demonstrated that patients who received genome matched treatment were 31% less likely to die within 1 year compared to those who did not receive genome matched treatment (HR: 0.69; 95% CI: 0.52-0.90; p-value: 0.006). Overall, GTT through this initiative resulted in levels of genome matched treatment that were similar to other initiatives, however, clinical trials represented a smaller share of treatments than previously reported, and "off-label" treatments represented a greater share. Although this was an observational study, we found evidence for a potential 1-year survival benefit for patients who received genome matched treatments. These findings suggest that when disseminated and implemented with a supportive infrastructure, GTT may benefit cancer patients in rural community oncology settings, with further work remaining on providing genome-matched clinical trials.