Diving into broad-scale and high-resolution population genomics to decipher drivers of structure and climatic vulnerability in a marine invertebrate.

Species with widespread distributions play a crucial role in our understanding of climate change impacts on population structure. In marine species, population structure is often governed by both high connectivity potential and selection across strong environmental gradients. Despite the complexity of factors influencing marine populations, studying species with broad distribution can provide valuable insights into the relative importance of these factors and the consequences of climate-induced alterations across environmental gradients. We used the northern shrimp Pandalus borealis and its wide latitudinal distribution to identify current drivers of population structure and predict the species' vulnerability to climate change. A total of 1514 individuals sampled across 24° latitude were genotyped at high geographic (54 stations) and genetic (14,331 SNPs) resolutions to assess genetic variation and environmental correlations. Four populations were identified in addition to finer substructure associated with local adaptation. Geographic patterns of neutral population structure reflected predominant oceanographic currents, while a significant proportion of the genetic variation was associated with gradients in salinity and temperature. Adaptive landscapes generated using climate projections suggest a larger genomic offset in the southern extent of the P. borealis range, where shrimp had the largest adaptive standing genetic variation. Our genomic results combined with recent observations point to further deterioration in southern regions and an impending vulnerable status in the regions at higher latitudes for P. borealis. They also provide rare insights into the drivers of population structure and climatic vulnerability of a widespread meroplanktonic species, which is crucial to understanding future challenges associated with invertebrates essential to ecosystem functioning.

Reduction in size-at-maturity in unprecedentedly strong cohorts of redfish (Sebastes mentella and S. fasciatus) in the Gulf of St. Lawrence and Laurentian Channel.

Life-history traits, such as size-at-maturity, are key parameters to model population dynamics used to inform fisheries management. Fishery-induced evolution, density-dependent effects, and global warming have been shown to affect size- and age-at-maturity, and resulting spawning stock biomass (SSB) in a wide range of commercial fish stocks. Marked changes in redfish biomass and environmental conditions in the Gulf of St. Lawrence and Laurentian Channel over the past decade called for a review and update of size-at-maturity for commercially important deepwater redfish Sebastes mentella and Acadian redfish Sebastes fasciatus stocks. Following a 25-year moratorium, local redfish biomass has recently reached unprecedented levels, co-occurring with an overall warming of bottom water temperatures. Our objectives were (1) to perform a histological assessment of redfish reproduction stages, including the validation and fine-tuning of a robust visual chart to facilitate monitoring of size-at-maturity and SSB in a transforming environment, and (2) to evaluate changes in size-at-maturity in unprecedentedly strong cohorts of redfish, and consequences for stock status assessment and fisheries management. Each specimen was genetically identified to species, and gonad reproduction stages were determined by histology and macroscopic appearances. The present study enabled a robust visual chart for continued and cost-effective monitoring of redfish reproduction stages to be refined and validated, and has shown a large decrease in redfish length when 50% of the individuals are considered mature that led to an increase in estimates of SSB during the 2011-2021 period for S. mentella and S. fasciatus. These changes modified the perception of stock status, thus having significant implications for fisheries management. Given that fishery-induced evolution and community structure changes along with global warming are affecting numerous stocks worldwide, the present study outlines a major and global challenge for scientists and resources managers. As shown by our results, the monitoring and frequent updates of life-history traits in transforming environments are needed to provide reliable science advice for sustainable fisheries.

Contrasting new and available reference genomes to highlight uncertainties in assemblies and areas for future improvement: an example with monodontid species.

BACKGROUND: Reference genomes provide a foundational framework for evolutionary investigations, ecological analysis, and conservation science, yet uncertainties in the assembly of reference genomes are difficult to assess, and by extension rarely quantified. Reference genomes for monodontid cetaceans span a wide spectrum of data types and analytical approaches, providing the context to derive broader insights related to discrepancies and regions of uncertainty in reference genome assembly. We generated three beluga (Delphinapterus leucas) and one narwhal (Monodon monoceros) reference genomes and contrasted these with published chromosomal scale assemblies for each species to quantify discrepancies associated with genome assemblies. RESULTS: The new reference genomes achieved chromosomal scale assembly using a combination of PacBio long reads, Illumina short reads, and Hi-C scaffolding data. For beluga, we identified discrepancies in the order and orientation of contigs in 2.2-3.7% of the total genome depending on the pairwise comparison of references. In addition, unsupported higher order scaffolding was identified in published reference genomes. In contrast, we estimated 8.2% of the compared narwhal genomes featured discrepancies, with inversions being notably abundant (5.3%). Discrepancies were linked to repetitive elements in both species. CONCLUSIONS: We provide several new reference genomes for beluga (Delphinapterus leucas), while highlighting potential avenues for improvements. In particular, additional layers of data providing information on ultra-long genomic distances are needed to resolve persistent errors in reference genome construction. The comparative analyses of monodontid reference genomes suggested that the three new reference genomes for beluga are more accurate compared to the currently published reference genome, but that the new narwhal genome is less accurate than one published. We also present a conceptual summary for improving the accuracy of reference genomes with relevance to end-user needs and how they relate to levels of assembly quality and uncertainty.

Spatial distribution and speciation of mercury in a recovering deepwater redfish (Sebastes mentella) population from St. Lawrence Estuary and Gulf, Canada.

Mercury (Hg) pollution poses a significant threat to the environment, particularly in the form of methylmercury (MeHg). However, little is known about the distribution and influencing factors of Hg in deep-sea (>200m) fish, which is crucial for assessing potential health risks to fish and humans. In Canada, the deepwater redfish (Sebastes mentella) has been designated as an endangered species. After a 25-year fishing moratorium, the redfish population in the St. Lawrence Estuary and Gulf is recovering, and resuming of commercial fishing and human consumption are expected. This study aimed to investigate the distribution of MeHg and total Hg (THg) in the muscle of redfish, as well as the factors influencing its distribution, and to assess the potential human health risks associated with redfish consumption. The redfish samples (n = 123) were collected by Fisheries and Oceans Canada in 2019. The concentrations of THg and MeHg in redfish muscle were determined to be 93.3 ± 183 ng/g (mean ± SD, wet weight) and 78.2 ± 149 ng/g, respectively. Large redfish (>30 cm) accumulated 20 to 30 times more Hg than small redfish (17-30 cm). Small redfish from the Estuary-Western Gulf had higher levels of MeHg and THg than those from the Laurentian Channel and the Northeast Gulf, but the Hg availability to redfish among the three areas were similar. Significant predictors of MeHg concentrations in redfish muscle were determined to be fish length, muscle moisture, δ15N, and N%. MeHg consumption by the general population with an average fish consumption rate is not anticipated to have adverse effects. This study establishes a baseline for future Hg monitoring in the deep water environments in this region. Further research is required to elucidate the cause-effect relationships between various environmental/biological parameters and Hg accumulation in deep-sea biota.

Population genomics and history of speciation reveal fishery management gaps in two related redfish species (Sebastes mentella and Sebastes fasciatus).

Understanding the processes shaping population structure and reproductive isolation of marine organisms can improve their management and conservation. Using genomic markers combined with estimation of individual ancestries, assignment tests, spatial ecology, and demographic modeling, we (i) characterized the contemporary population structure, (ii) assessed the influence of space, fishing depth, and sampling years on contemporary distribution, and (iii) reconstructed the speciation history of two cryptic redfish species, Sebastes mentella and S. fasciatus. We genotyped 860 individuals in the Northwest Atlantic Ocean using 24,603 filtered single nucleotide polymorphisms (SNPs). Our results confirmed the clear genetic distinctiveness of the two species and identified three ecotypes within S. mentella and five populations in S. fasciatus. Multivariate analyses highlighted the influence of spatial distribution and depth on the overall genomic variation, while demographic modeling revealed that secondary contact models best explained inter- and intragenomic divergence. These species, ecotypes, and populations can be considered as a rare and wide continuum of genomic divergence in the marine environment. This acquired knowledge pertaining to the evolutionary processes driving population divergence and reproductive isolation will help optimizing the assessment of demographic units and possibly to refine fishery management units.

Connecting tree-ring phenotypes, genetic associations and transcriptomics to decipher the genomic architecture of drought adaptation in a widespread conifer.

As boreal forests face significant threats from climate change, understanding evolutionary trajectories of coniferous species has become fundamental to adapting management and conservation to a drying climate. We examined the genomic architecture underlying adaptive variation related to drought tolerance in 43 populations of a widespread boreal conifer, white spruce (Picea glauca [Moench] Voss), by combining genotype-environment associations, genotype-phenotype associations, and transcriptomics. Adaptive genetic variation was identified by correlating allele frequencies for 6,153 single nucleotide polymorphisms from 2,606 candidate genes with temperature, precipitation and aridity gradients, and testing for significant associations between genotypes and 11 dendrometric and drought-related traits (i.e., anatomical, growth response and climate-sensitivity traits) using a polygenic model. We identified a set of 285 genes significantly associated with a climatic factor or a phenotypic trait, including 110 that were differentially expressed in response to drought under greenhouse-controlled conditions. The interlinked phenotype-genotype-environment network revealed eight high-confidence genes involved in white spruce adaptation to drought, of which four were drought-responsive in the expression analysis. Our findings represent a significant step toward the characterization of the genomic basis of drought tolerance and adaptation to climate in conifers, which is essential to enable the establishment of resilient forests in view of new climate conditions.

Life-stage-dependent supergene haplotype frequencies and metapopulation neutral genetic patterns of Atlantic cod, Gadus morhua, from Canada's Northern cod stock region and adjacent areas.

Among highly migratory fish species, nursery areas occupied by juveniles often differ from adult habitats. To better understand the spatial dynamics of Canada's Northern cod stock, juveniles caught off the east coast of Newfoundland and Labrador were compared to adults from the same region as well as individuals from other areas in Atlantic Canada using double-digest restriction site-associated DNA sequencing-derived single nucleotide polymorphisms. A reduced proportion of homozygotes with a chromosomal inversion located in linkage group 1 (LG1) was detected between juvenile and adult samples in the Northern cod stock region, potentially indicating age-dependent habitat use or ontogenetic selection for attributes associated with the many genes located in LG1. No selectively neutral genetic differences were found between samples from the Northern cod stock; nevertheless, significant differences were found between some of these samples and cod collected from St. Pierre Bank, Bay of Fundy, Browns Bank and the southern Scotian Shelf. Clustering analysis of variants at neutral loci provided evidence for three major genetic units: (a) the Newfoundland Atlantic Coast, (b) eastern and southern Gulf of St. Lawrence and Burgeo Bank and (c) the Bay of Fundy, Browns Bank and southern Scotian Shelf. Both adaptive and neutral population structure within the Northern cod stock should be considered by managers to promote demographic rebuilding of the stock.

Hydroxyacetophenone defenses in white spruce against spruce budworm.

We review a recently discovered white spruce (Picea glauca) chemical defense against spruce budworm (Choristoneura fumiferana) involving hydroxyacetophenones. These defense metabolites detected in the foliage accumulate variably as the aglycons, piceol and pungenol, or the corresponding glucosides, picein and pungenin. We summarize current knowledge of the genetic, genomic, molecular, and biochemical underpinnings of this defense and its effects on C. fumiferana. We present an update with new results on the ontogenic variation and the phenological window of this defense, including analysis of transcript responses in P. glauca to C. fumiferana herbivory. We also discuss this chemical defense from an evolutionary and a breeding context.

The large repertoire of conifer NLR resistance genes includes drought responsive and highly diversified RNLs.

The NLRs or NBS-LRRs (nucleotide-binding, leucine-rich-repeat) form the largest resistance gene family in plants, with lineage-specific contingents of TNL, CNL and RNL subfamilies and a central role in resilience to stress. The origin, evolution and distribution of NLR sequences has been unclear owing in part to the variable size and diversity of the RNL subfamily and a lack of data in Gymnosperms. We developed, searched and annotated transcriptomes assemblies of seven conifers and identified a resource of 3816 expressed NLR sequences. Our analyses encompassed sequences data spanning the major groups of land plants and determinations of NLR transcripts levels in response to drought in white spruce. We showed that conifers have among the most diverse and numerous RNLs in tested land plants. We report an evolutionary swap in the formation of RNLs, which emerged from the fusion of an RPW8 domain to a NB-ARC domain of CNL. We uncovered a quantitative relationship between RNLs and TNLs across all land plants investigated, with an average ratio of 1:10. The conifer RNL repertoire harbours four distinct groups, with two that differ from Angiosperms, one of which contained several upregulated sequences in response to drought while the majority of responsive NLRs are downregulated.

Association genetics of acetophenone defence against spruce budworm in mature white spruce.

BACKGROUND: Outbreaks of spruce budworm (SBW, Choristoneura fumiferana Clem.) cause major recurrent damage in boreal conifers such as white spruce (Picea glauca [Moench] Voss) and large losses of forest biomass in North America. Although defensive phenolic compounds have recently been linked to chemical resistance against SBW, their genetic basis remains poorly understood in forest trees, especially in conifers. Here, we used diverse association genetics approaches to discover genes and their variants that may control the accumulation of acetophenones, and dissect the genetic architecture of these defence compounds against SBW in white spruce mature trees. RESULTS: Out of 4747 single nucleotide polymorphisms (SNPs) from 2312 genes genotyped in a population of 211 unrelated individuals, genetic association analyses identified 35 SNPs in 33 different genes that were significantly associated with the defence traits by using single-locus, multi-locus and multi-trait approaches. The multi-locus approach was particularly effective at detecting SNP-trait associations that explained a large fraction of the phenotypic variance (from 20 to 43%). Significant genes were regulatory including the NAC transcription factor, or they were involved in carbohydrate metabolism, falling into the binding, catalytic or transporter activity functional classes. Most of them were highly expressed in foliage. Weak positive phenotypic correlations were observed between defence and growth traits, indicating little or no evidence of defence-growth trade-offs. CONCLUSIONS: This study provides new insights on the genetic architecture of tree defence traits, contributing to our understanding of the physiology of resistance mechanisms to biotic factors and providing a basis for the genetic improvement of the constitutive defence of white spruce against SBW.

Genetic control and evolutionary potential of a constitutive resistance mechanism against the spruce budworm (Choristoneura fumiferana) in white spruce (Picea glauca).

Insect herbivory may drive evolution by selecting for trees with heritable resistance against defoliation. The spruce budworm (Choristoneura fumiferana, SBW) is a highly damaging forest insect pest that can affect population structure of white spruce (Picea glauca) in North America. Resistance against SBW was recently described in white spruce and was linked to three constitutive resistance biomarkers: the phenolic compounds piceol and pungenol, and expression of a beta-glucosidase encoding gene (Pgßglu-1). We investigated the phenotypic variability and heritability of these resistance biomarkers and of picein, the precursor of piceol, in the foliage of 874 trees belonging to 33 full-sib families and 71 clonal lines under evaluation in seven field locations in Eastern Canada. We aimed to (i) determine their genetic control, (ii) estimate the genetic and phenotypic correlations among defense biomarkers, and (iii) determine whether their constitutive levels are associated with detrimental trade-offs on growth. Quantitative genetics analyses indicated that all four traits are moderately to highly heritable. The full-sib and clonal analyses showed that additive and non-additive genetic effects play major and minor roles, respectively. Positive genetic and phenotypic correlations between resistance biomarkers and primary growth indicated that there is no trade-off between total height and height increment and resistance traits, contradicting the GDBH (Growth Differentiation Balance Hypothesis). Our findings about the predominant additive genetic basis of the resistance biomarkers show that adaptive evolution of white spruce natural populations to resist to SBW is possible and that potentially important gains could also be expected from artificial selection.

Evolution of the biosynthesis of two hydroxyacetophenones in plants.

Acetophenones are phenolic metabolites of plant species. A metabolic route for the biosynthesis and release of 2 defence-related hydroxyacetophenones in white spruce (Picea glauca) was recently proposed to involve 3 phases: (a) biosynthesis of the acetophenone aglycons catalysed by a currently unknown set of enzymes, (b) formation and accumulation of the corresponding glycosides catalysed by a glucosyltransferase, and (c) release of the aglycons catalysed by a glucosylhydrolase (PgßGLU-1). We tested if this biosynthetic model is conserved across Pinaceae and land plant species. We assayed and surveyed the literature and sequence databases for possible patterns of the presence of the acetophenone aglycons piceol and pungenol and their glucosides, as well as sequences and expression of Pgßglu-1 orthologues. In the Pinaceae, the 3 phases of the biosynthetic model are present and differences in expression of Pgßglu-1 gene orthologues explain some of the interspecific variation in hydroxyacetophenones. The phylogenetic signal in the metabolite phenotypes was low across species of 6 plant divisions. Putative orthologues of PgßGLU-1 do not form a monophyletic group in species producing hydroxyacetophenones. The biosynthetic model for acetophenones appears to be conserved across Pinaceae, whereas convergent evolution has led to the production of acetophenone glucosides across land plants.

Insect herbivory (Choristoneura fumiferana, Tortricidea) underlies tree population structure (Picea glauca, Pinaceae).

Variation in insect herbivory can lead to population structure in plant hosts as indicated by defence traits. In annual herbaceous, defence traits may vary between geographic areas but evidence of such patterns is lacking for long-lived species. This may result from the variety of selection pressures from herbivores, long distance gene flow, genome properties, and lack of research. We investigated the antagonistic interaction between white spruce (Picea glauca) and spruce budworm (SBW, Choristoneura fumiferana) the most devastating forest insect of eastern North America in common garden experiments. White spruces that are able to resist SBW attack were reported to accumulate the acetophenones piceol and pungenol constitutively in their foliage. We show that levels of these acetophenones and transcripts of the gene responsible for their release is highly heritable and that their accumulation is synchronized with the most devastating stage of SBW. Piceol and pungenol concentrations negatively correlate with rate of development in female SBW and follow a non-random geographic variation pattern that is partially explained by historical damage from SBW and temperature. Our results show that accumulation of acetophenones is an efficient resistance mechanism against SBW in white spruce and that insects can affect population structure of a long-lived plant.