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.

Organic UV absorbents in the deepwater redfish (Sebastes mentella) from the St. Lawrence Estuary and Gulf: Distribution and human health risk assessment

UV absorbents (UVAs), such as organic UV filters (UVFs) and benzotriazole UV stabilizers (BZT-UVs), are used in a wide range of consumer and industrial products and they are contaminants of emerging concern in the environment. However, their occurrence and fate in the deep-sea environments are inadequately understood. This study investigated the occurrence and distribution of five UVFs and ten BZT-UVs in the muscle (n = 127) of 2019-collected deepwater redfish (Sebastes mentella) from the St. Lawrence Estuary and Gulf (SLEG) (Canada) to better understand the accumulation of these contaminants in deep-sea fish. Small redfish (<30 cm) tended to have higher concentrations of UVAs in the muscle than that of larger specimens (>30 cm). The UVF 2-hydroxy-4-methoxybenzophenone (BP3) was the most frequently detected (present in 34 % of all samples) target UVA, with concentrations as high as 413 ng/g (dry weight). According to the δ15N and δ13C data, pelagic-eating redfish, and individuals with lower trophic levels had higher lipid content and accumulated more BP3 in their muscles. Four BZT-UVs were detected in redfish muscle, but the detection frequency was lower than 30 %. The estimated hazard quotient for these contaminants was <2.3 × 10-2 for general Canadian populations, indicating that they are unlikely to pose health risks to humans through redfish consumption. Factors influencing UVAs bioaccumulation in redfish, as well as the effects UVAs may have on deep-sea species, should be researched further.

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.

The lack of genetic variation underlying thermal transcriptomic plasticity suggests limited adaptability of the Northern shrimp, Pandalus borealis

Introduction

Genetic variation underlies the populations’ potential to adapt to and persist in a changing environment, while phenotypic plasticity can play a key role in buffering the negative impacts of such change at the individual level.

Methods

We investigated the role of genetic variation in the thermal response of the northern shrimp Pandalus borealis, an ectotherm species distributed in the Arctic and North Atlantic Oceans. More specifically, we estimated the proportion transcriptomic responses explained by genetic variance of female shrimp from three origins after 30 days of exposure to three temperature treatments.

Results

We characterized the P. borealis transcriptome (170,377 transcripts, of which 27.48% were functionally annotated) and then detected a total of 1,607 and 907 differentially expressed transcripts between temperatures and origins, respectively. Shrimp from different origins displayed high but similar level of transcriptomic plasticity in response to elevated temperatures. Differences in transcript expression among origins were not correlated to population genetic differentiation or diversity but to environmental conditions at origin during sampling.

Discussion

The lack of genetic variation explaining thermal plasticity suggests limited adaptability in this species’ response to future environmental changes. These results together with higher mortality observed at the highest temperature indicate that the thermal niche of P. borealis will likely be restricted to higher latitudes in the future. This prediction concurs with current decreases in abundance observed at the southern edge of this species geographical distribution, as it is for other cold-adapted crustaceans.

Maximizing the reliability and the number of species assignments in metabarcoding studies using a curated regional library and a public repository

Biodiversity assessments relying on DNA have increased rapidly over the last decade. However, the reliability of taxonomic assignments in metabarcoding studies is variable and affected by the reference databases and the assignment methods used. Species level assignments are usually considered as reliable using regional libraries but unreliable using public repositories. In this study, we aimed to test this assumption for metazoan species detected in the Gulf of St. Lawrence in the Northwest Atlantic. We first created a regional library (GSL-rl) by data mining COI barcode sequences from BOLD, and included a reliability ranking system for species assignments. We then estimated 1) the accuracy and precision of the public repository NCBI-nt for species assignments using sequences from the regional library and 2) compared the detection and reliability of species assignments of a metabarcoding dataset using either NCBI-nt or the regional library and popular assignment methods. With NCBI-nt and sequences from the regional library, the BLAST-LCA (least common ancestor) method was the most precise method for species assignments, but the accuracy was higher with the BLAST-TopHit method (&gt;80% over all taxa, between 70% and 90% amongst taxonomic groups). With the metabarcoding dataset, the reliability of species assignments was greater using GSL-rl compared to NCBI-nt. However, we also observed that the total number of reliable species assignments could be maximized using both GSL-rl and NCBI-nt with different optimized assignment methods. The use of a two-step approach for species assignments, i.e., using a regional library and a public repository, could improve the reliability and the number of detected species in metabarcoding studies.

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 (Piceaglauca [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.

see also the Perspective by Lars Opgenoorth and Christian Rellstab

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.

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.

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.