Heredity

The European flat oyster (Ostrea edulis), like numerous other oyster species, has been exploited for millennia and cultivated and translocated for centuries. Following a severe population decline, and in the context of ongoing conservation and restoration programs, genetic considerations must now be addressed to avoid mistakes. The objective of our study was to complement population genetic studies conducted at various scales along European coasts. Our sampling primarily targeted the French Atlantic, English Channel, and Mediterranean coasts, aiming to provide a fine-scale genetic characterization of populations in these regions. By integrating SNP array and low-coverage sequencing datasets, we obtained a comprehensive overview of the population genetic structure of Ostrea edulis across western Europe. Most previously identified clusters in Western Europe were confirmed. In France, populations assigned to these clusters exhibited notable within-patch homogeneity. However, two key findings emerged: (1) an extensive overlap zone between the Atlantic and western Mediterranean clusters, spanning at least from southern Portugal to southern France, and (2) the detection of a novel, clearly distinct cryptic cluster east of the English Channel, whose geographic range remains to be better delineated. These insights are critical for informing management decisions, particularly as restoration and conservation plans are currently being implemented across the species range.

Ectomycorrhizal (EM) fungi are vital to forest ecosystems, supporting tree growth and survival. However, their inclusion in conservation policy and action remains limited and little is known about the status of their genetic diversity, which is essential for their long-term survival and adaptation. The Global Biodiversity Framework adopted a genetic indicator based on the effective population size, Ne, to monitor genetic diversity in all species. To date, it is still uncertain how Ne, a key parameter, can be reliably assessed in species with complex life history traits. Ectomycorrhizal fungi are a highly diverse group of taxa displaying haplodiplontic life cycles with partially clonal reproduction. Here, we review the literature to understand how these life history traits might affect Ne and its estimation in six species of EM fungi. We estimated Ne in 19 populations using eight genetic and genomic datasets from selected studies. We compared Ne estimates using Linkage Disequilibrium (LD) and Sibship Frequency (SF) methods. We tested how Ne estimates change due to partial clonality and genetic structure gradients and whether the number of genetic markers influence the precision of the estimates. We show a systematic bias in Ne estimations when large clones are present and when populations are not correctly delimited. We found both methods are not robust to these factors, which makes them unreliable for conservation assessment purposes in EM fungi. This study provides new perspectives for further research into the links between life history traits and the effective population size of ectomycorrhizal fungi.

Genomic offset (GO) is increasingly used to predict population maladaptation risk under climate change, with larger offsets assumed to indicate greater vulnerability. Despite rapid adoption in conservation planning, it remains unclear how sensitive GO estimates are to key methodological choices, including SNP set composition, genotype-environment association (GEA) methods, and the specific GO metric used. Empirical validation against observed population dynamics also remains limited. Here, we evaluate the methodological robustness and predictive performance of GO using multidecadal demographic monitoring data from Lycaeides butterflies, a system with short generation times and high fecundity that may facilitate rapid adaptive responses. GO estimates were broadly consistent across SNP sets, regardless of composition or size, with climate-associated and randomly selected SNPs yielding largely concordant values. Consistency across GEA methods was moderate and depended on the SNP set used. In contrast, GO metrics differed substantially in the magnitude of maladaptation estimated, suggesting they capture distinct biological signals and should not be treated as interchangeable. Crucially, GO was a poor predictor of observed population trends, regardless of SNP set composition, GO metric, or GEA method, both at sites used to fit GEA models and when extrapolated to independent demographic sites. These findings suggest that, while GO provides a valuable conceptual framework for assessing potential maladaptation, its quantitative estimates and predictive power are sensitive to methodological choices and species-specific biological context. We therefore urge careful alignment of GO metric assumptions with conservation objectives, along with rigorous empirical validation, before GO estimates are used to inform management decisions.

The climatic oscillations of the Quaternary have likely affected the demographic history of many species, and PSMC (Pairwise Sequentially Markovian Coalescent) has been widely used to investigate these histories. However, it is increasingly acknowledged that PSMC trajectories are difficult to interpret. First, they are influenced by connectivity changes, even without population size changes. Second, most PSMC curves exhibit a few humps when tens of cycles occurred during the Pleistocene. Finally, responses to ancient habitat change have been shown to be species-specific. To address these issues, we simulated structured populations where connectivity (or population size and connectivity) varied according to successive interglacial and glacial periods during the last 2.6 million years. We computed the IICR (Inverse Instantaneous Coalescence Rate), the function that PSMC estimates, and ran PSMC. We further varied the generation length and assumed that some species were positively or negatively affected by glacials. We found that the IICR carries information regarding the demographic oscillations, but that PSMC fails to recover it for times older than 300 ky. For the last 200 ky, PSMC was often able to reproduce qualitatively the demographic oscillations. We also tested SNIF (Structured Non-stationary Inferential Framework), which produced good results using the IICR curve as an input but not when using the PSMC curve. Altogether, our study suggests that the humps older than 300 ky in PSMC histories are unlikely to represent trends of population size or connectivity. However, improving the estimation of the IICR could potentially help reconstruct some of these past demographic changes.

Japanese chum salmon supported by one of the world largest hatchery programs have experienced severe declines in marine survival and egg size. To investigate the underlying mechanisms, we analyzed a 21-year time series (1999-2019) of reproductive traits of age-4 chum salmon from 13 rivers together with climate and salmon abundance data using a bootstrap-supported Bayesian network. Here, we assumed that environmental variables can affect the chum salmon populations, but not vice versa, and that there could be maternal effect on reproductive traits, but not the other way around. These constraints enabled us to infer the causal links that shaped the biogeography of North Pacific chum salmon. Global warming caused a decline in Japanese chum salmon abundance, resulting in the increase of the competing Russian chum, which in turn decreased the female body size, fecundity, and egg size of Japanese chum. These findings suggest that climate-driven warming may have exposed genetic effects of hatchery practices, contributing to fitness decline in Japanese chum salmon and the ecological reorganization of chum salmon populations in the North Pacific.

Infection can substantially reduce host fitness and influence population dynamics, yet it is often difficult to detect and quantify in wild animal populations. Molecular tools offer a valuable means of identifying cryptic infection in natural systems. Using whole-genome bisulfite sequencing, we examined whether infection with the parasitic leech Ozobranchus margoi is associated with DNA methylation variation in loggerhead sea turtles (Caretta caretta), while also assessing the potential value of this variation as a biomarker of parasite infection. In nesting females, we identified infection-associated differentially methylated CpG sites associated with genes implicated in immune signalling and cellular regulation. Offspring of infected females also showed infection-associated methylation patterns, despite not being directly exposed to the parasite themselves. Differential methylation analyses identified genes involved in immunity, neurodevelopment and metabolic activity, with limited overlap in associated genes and no overlap in differentially methylated sites between generations. Maternal and offspring genome-wide methylation levels showed a non-linear association that differed subtly with maternal infection status, indicating that infection modifies intergenerational methylation associations. Finally, methylation profiles showed strong discriminatory power for maternal infection status in both maternal and hatchling samples using machine learning models, supporting their potential as candidate biomarkers of cryptic infection. Together, these results show that parasite infection is associated with distinct, generation-specific DNA methylation signatures, and highlight the potential value of epigenetic data for monitoring cryptic infection states in conservation-relevant systems.

Animals fight for resources to obtain fitness benefits; most contests are intrasexual, and males tend to fight more than females. Although the genetic basis of male aggression is well studied, we know little about the genetic variation of female aggression. Female aggression varies with reproductive status and is potentially influenced not only by her genotype, but also by the genotype of her mate. Here we measured both male and female aggression in a set of Drosophila melanogaster inbred lines by competing each line against a standard competitor. Aggression varied among lines for both sexes, but male and female aggression were not correlated. Female aggression for many lines increased with mating, as expected, but not all lines changed aggression. However, when females were mated to males of different lines, male genotype did not affect the post-mating change in aggression, suggesting that ejaculate-mediated effects do not vary across these lines. The aggression level of the standard opponent was positively correlated with that of focal individuals indicating that individuals modulate their behavior according to the genotype of their opponent.

Demographic factors are intrinsically crucial to evaluate species' extinction risk. However, measuring them remains difficult and time-consuming and the use of genomic summary statistics has been advocated to assess the conservation status of a species. In the present study, we estimated (i) the census number (Nc), (ii) effective population size (Ne) over three different time periods, recent, historical and ancient, (iii) neutral genetic diversity ({pi}4), and (iv) a measure of the efficacy of purifying selection ({pi}0/{pi}4) for 101 plant species using population genomic sequencing data. Twenty-one species are from the Plant Species with Extremely Small Populations (PSESP) program of SW China. Threatened species exhibited significantly lower Ne, Nc, {pi}4, and weaker purifying selection, but had a higher Ne/Nc ratio than non-threatened ones. Nc was the main determinant in identifying conservation status, and contemporary neutral genetic diversity was predominantly influenced by historical Ne. In the absence of demographic information, genetic parameters are a good proxy of conservation status, likely because currently threatened species also had a low historical population size. In summary, our findings suggest that direct estimates of Nc are more useful than {pi}4, although the latter remains a valuable conservation indicator. Hence, efforts such as the PSESP should be extended.

The potential for environmental change to compound the detrimental effects of inbreeding depression in small and isolated populations is a significant concern in conservation biology. Previous evidence syntheses suggested that environmental stress exacerbates inbreeding depression, but were based on limited data. Here, we comprehensively test the relationship between inbreeding depression and environmental stress in natural populations using Bayesian mixed-effects meta-analysis on a large, high-quality data set of 2127 inbreeding depression effect sizes from animals and plants. Our results show that inbreeding depression is significantly higher in benign than in stressful environments. Analyses of both inbreeding depression and stress-induced changes in genetic load supported a unimodal (humped) relationship between the costs of inbreeding and stress intensity, with a peak at intermediate stress. At the highest levels of stress there was, on average, a significantly greater inbreeding load in benign than in stressful environments. We suggest that the lower cost of inbreeding associated with extreme stress results from constraints on the expression of inbreeding depression as fitness and phenotypes decline towards zero. Our findings help to resolve long-standing uncertainty around how inbreeding and environmental change interact, revealing that inbreeding responses vary non-linearly with environmental stress intensity, but showing that stress does not generally amplify inbreeding depression. As such, they will inform both the management of populations of conservation concern and predictions of species responses to global environmental change.

Many biological characteristics arise by interactions between more than one biological organism or unit. Fertilization success in sexually reproducing species represents such an extended phenotype where both mates are required to be fertile for a successful outcome. Consequently, predictive models should account for the joint nature of reproductive performance while offering interpretable estimates for individual mate contributions. Recent advances in genomics and machine learning (ML) provide standardized, high-dimensional genetic information on one hand and computational tools capable of modeling complex biological systems on the other. Here, we construct and evaluate two-tower (TT) machine learning architectures for genomic prediction of binary reproductive outcomes and recovery of sex-specific fertility liabilities. Simulated datasets, generated under a range of genetic architectures, were utilized to compare multilayer perceptron (TT-MLP), convolutional neural network (TT-CNN), and L1-regularized linear (TT-LASSO) two-tower models. Simulation scenarios varied sex-specific heritabilities, genetic correlations, infertility prevalence, mating structure, and sex-specific infertility rates. Models were evaluated with regard to their ability to predict reproductive success at pair level and also recover true underlying genetic values for male and female fertility. Prediction accuracy increased with the underlying heritable component as expected, while sex-specific tower-scores successfully recovered latent fertility liabilities despite models being trained only on observed joint outcomes. TT-LASSO achieved the highest overall classification performance, whereas TT-MLP provided more balanced and consistent recovery of sex-specific genetic values across scenarios. An additional simulation, incorporating genotype-dependent mate compatibility demonstrated advantages of fully-connected neural networks for capturing non-additive interactions. These results indicate that two-tower frameworks provide a powerful approach for modeling reproductive traits, enabling simultaneous prediction of aggregate reproductive outcomes and sex-specific fertility liabilities from genotypic information.

The Wangs garden lizard, Calotes wangi, is a widely distributed agamid species in Southern China and Northern Vietnam and exhibits pronounced colour variation and rapid body colour change. Despite increasing interest in the genomic basis of colour variation, chromosome-level genomic resources remain limited in agamid lizards. Here, we generated a chromosome-level reference genome of C. wangi using PacBio HiFi sequencing and Hi-C scaffolding. The final genome assembly was approximately 1.66 Gb in size and comprised 6 macrochromosomes and 11 microchromosomes, with a contig N50 of 110.09 Mb and 98.9% complete BUSCO genes. A total of 20,442 protein-coding genes were annotated. Comparative genomic analyses identified 297 significantly expanded gene families, with enriched functions associated with steroid metabolism, chromatin regulation, and epigenetic processes. This high-quality genome assembly provides an important genomic resource for future studies of colour variation, phenotypic plasticity, and evolutionary diversification in agamid lizards.

Understanding spatial and temporal connectivity among individuals with different migration strategies is essential for migratory ecology and effective conservation, yet it often requires integrating multiple data sources. In Northeast Atlantic European sea bass (Dicentrarchus labrax), electronic tagging has revealed partial migration, with both resident and long-distance migrants showing fidelity to summer feeding and winter spawning areas. However, the role of regional spawning-site philopatry in shaping migration patterns and stock connectivity remains unclear. Here, we combine reconstructed migration trajectories with genome-wide analyses of gene flow and recent relatedness in 708 individuals sampled from 10 French Atlantic locations. We identify a seasonally shifting genetic discontinuity between the Bay of Biscay (BOB) and Northern (NS) stocks, located off western Brittany during winter spawning and displaced northeastward into the central English Channel during summer feeding. Despite seasonal mixing in the English Channel, an association between individual genetic composition and spawning-site selection supports regional spawning-site philopatry. Analyses of long genomic segments shared identical-by-descent reveal substantially greater connectivity within stocks than between stocks, indicating that philopatry constrains effective gene flow despite seasonal mixing. Reanalysis of independent genomic data further shows that sea bass from the northern Atlantic range predominantly belong to the Northern stock. Together, these results show how seasonal movements reshape spatial genetic structure while maintaining demographic subdivision, with direct implications for fisheries management.

Long-term studies of isolated animal populations have greatly improved the understanding of various evolutionary processes. However, potentially elevated inbreeding in those compared to wild populations is a common concern. Conventionally, inbreeding has been investigated using reconstructed pedigrees, but nowadays it can be done directly at the genomic level. Here, we utilize genomic data from an intensively studied isolated rhesus macaque (Macaca mulatta) population on the small island Cayo Santiago (Puerto Rico), which was founded in 1938 with wild animals from India. We quantified inbreeding levels by inferring runs of homozygosity (ROH), i.e., long identical haplotypes inherited from both parents. We identified ROH in 97 ~5x-coverage genomes from Cayo Santiago and, for comparison, in 79 rhesus macaque genomes from five wild populations from China. Notably, this conventionally considered low-coverage data proved sufficient to infer ROHs >4 centimorgans long after imputing the genomes using a reference panel. Our results revealed that the ROH-derived effective population size on Cayo Santiago, 420 individuals, falls within the ranges we inferred in wild populations. Moreover, a general scarcity of individuals with long ROH in both the Cayo and wild populations indicates very few cases of close-kin breeding, suggesting that mechanisms to avoid close-kin breeding operate in rhesus macaques, both in wild and isolated populations. Taken together, our results suggest that Cayo Santiago remains a representative study population.

We present a chromosome-level genome assembly of a female great grey owl (Strix nebulosa lapponica). The assembly comprises two pseudo-haplotypes of 1554 Mb and 1242 Mb, with 83.2% and 91.4% scaffolded into 40 autosomal chromosomes, in addition to the W and Z sex chromosomes both placed in hap1. Assembly completeness is high (BUSCO 99.2% and 94.8%), with 18,493 and 17,279 annotated protein-coding genes for hap1 and hap2, respectively. This genome establishes a reference for investigating genetic variation and chromosome evolution in great grey owls. Compared with the previous S. nebulosa assembly, this assembly includes both sex chromosomes, separates regions that were previously collapsed, and resolves 82 chromosomes total. While larger chromosomes show broadly conserved synteny across owl assemblies, the recovery of additional conserved microchromosome-associated genes suggests that ONT reads improved resolution of the smallest chromosomes relative to HiFi-based assemblies.

Biological invasions on islands provide a natural framework to study how dispersal and connectivity influence evolutionary and ecological processes. The avian nest parasitic fly, Philornis downsi - first recorded in Darwin's finch nests in 1997 - causes high mortality in endemic land birds, yet its inter-island and sex-specific patterns of dispersal and genetic structure remain poorly understood. We use low-coverage whole genome sequencing to investigate genome-wide patterns of genetic diversity, directional migration and effective population size in P. downsi across five major Galapagos Islands and its native range in mainland Ecuador. We find evidence for a genetic bottleneck in the Galapagos, isolation by distance, and evidence that the island closest to the Ecuadorian mainland, San Cristobal, is genetically divergent from the other four islands sampled, despite retaining the highest genetic diversity. No evidence was found for sex-biased dispersal; however, sex-biased genetic structure was detected using only markers from inferred autosomal scaffolds. We found asymmetric gene flow with higher migration rates from San Cristobal westward to the other islands, matching the direction of both southeast trade winds and major cargo shipping routes. Our results suggest both natural and human-mediated colonisation of P. downsi from the mainland through San Cristobal to the other islands, followed by high inter-island dispersal among closely situated sink islands. Our findings are critical for prioritising islands for control strategies that will reduce P. downsi impacts on vulnerable endemic birds and underscore the value of understanding directional migration patterns for managing invasive species in metapopulations.

Plant regeneration is a dynamic developmental process that spans from cell dedifferentiation to organ reconstruction in response to inductive cues, such as wounding stress and hormonal signals. Although this capacity varies widely both between and within species, a comprehensive understanding of the genetic and epigenetic bases of this variation remains incomplete. To address this issue, we revisited published datasets on natural variation in in vitro regeneration capacity in Arabidopsis thaliana. Using quantitative genetic approaches, including meta-analyses of genome-wide association studies (GWAS) and multi-locus models, we dissected the genetic architecture underlying regeneration traits. Our results showed that shoot regeneration capacity is primarily explained by allelic variation in the cis-regulatory region of WUSCHEL (WUS), a key regulator of shoot meristem formation. Notably, these polymorphisms are also associated with epigenetic variants of the DNA transposon ATDNA2T9C, which is located within the regulatory region. Furthermore, allelic variation in ARABIDOPSIS RESPONSE REGULATOR 2 (ARR2), a positive regulator of cytokinin signaling, is associated with callus formation and greening traits and may promote shoot formation through genetic interactions with WUS alleles. Although in vitro regeneration is controlled by complex, multilayered gene regulatory networks, our results suggest that, in A. thaliana, natural variation in regeneration capacity is largely shaped by a small number of major-effect modifiers together with epigenetic variation and genetic interactions, despite the substantial heterogeneity observed among natural populations.

The global study of freshwater clams in the genus Corbicula is frequently confounded by invasive androgenetic lineages that experience mitochondrial DNA capture and clonal propagation. In contrast, the endemic Corbicula of Sulawesi's ancient lakes reproduce sexually, offering a uniquely reliable system for mitochondrial population genetics. This study provides the first population-genetic framework for two endemic species, Corbicula possoensis (Lake Poso) and C. linduensis (Lake Lindu), using the cytochrome c oxidase subunit I (COI) marker. We analysed 90 newly generated COI sequences from C. possoensis (six stations) and C. linduensis (three stations), integrated with reference sequences from GenBank, to assess genetic diversity, population structure, and phylogeographic patterns. Hierarchical AMOVA revealed deep divergence between the two lakes ({Phi}_CT = 0.607), consistent with prolonged independent isolation rather than a single shared vicariance event, as the two species do not form a sister pair in the phylogeny. Within Lake Poso, C. possoensis exhibited exceptionally high genetic diversity (24 haplotypes; h = 0.876; {pi} = 0.016) and pronounced micro-geographic structuring into three phylogeographic zones (North: Tentena and Siuri; East: Tando Nceppo and Busogo Beach; Southwest: Bancea and Pendolo), each characterised by distinct haplogroups. Remarkably, the maximum divergence between zones (K2P = 2.33%) approached the interspecific distance between C. possoensis and C. linduensis (K2P = 2.42%), indicating that within-lake mitochondrial divergence has reached near-interspecific levels. Conversely, C. linduensis displayed near-panmixia and extreme genetic depauperation (3 haplotypes; h = 0.246; {pi} = 0.0004), indicating long-term demographic stasis within a restricted habitat. The deep phylogeographic zonation in C. possoensis suggests that its discrete populations should be treated as separate Management Units (MUs) in conservation planning to preserve locally adapted gene complexes, whereas the severely depauperate gene pool of C. linduensis renders it critically vulnerable to environmental disturbance and invasive species, warranting urgent IUCN Red List assessment. To validate these mitochondrial boundaries and inform future conservation strategies, multi-marker and genome-wide reassessments are strongly recommended.

The coast redwood (Sequoia sempervirens) is a long-lived, hexaploid conifer of high ecological, cultural, and economic value whose range has been greatly reduced by historical logging. Effective restoration and conservation depend on understanding patterns of genetic differentiation across the redwood range to delineate populations for management prioritization. Yet, past range-wide studies provided only a partial picture of population structure in coast redwood as they relied on a limited set of genetic markers or limited sampling, as sequencing was done on the same range-wide provenance collection. Here, we analyze 334,029 SNPs from a new range-wide set of 224 individuals using a dosage-based approach that accounts for polyploidy. Principal coordinates and neighbor-joining analyses reveal clear latitudinal genetic differentiation, with a distinct break south of San Francisco Bay. Outlier SNP analysis indicates new candidate loci involved in salinity tolerance, climate stress response, and nutrient uptake, suggesting potential local adaptation. These results point to the central role of geography in shaping genetic variation in coast redwood and give scientific basis for designing new conservation strategies and future experiments, including assisted migration, provenance trials, and restoration planning aimed at preserving the species into the future.

The unprecedented rate of climate warming threatens many species, and assessing their vulnerability to climate change represents a critical challenge in conservation biology. The Apollo butterfly, an emblematic mountain species, is expected to be impacted by climate change. Here, we analysed thousands of SNPs from 101 localities across Apollo French distribution. We identified 93 SNPs strongly associated with climate variation using five genotype-environment association analyses. We forecasted future climate maladaptation of French Apollo populations using four genomic offset methods and integrated these results with neutral and adaptive genetic diversity, genetic structure and adaptive climatic niches to infer their vulnerability to climate change. Jura and Alps populations exhibited the lowest risk of vulnerability to climate change, with low genomic offsets, high genetic diversity and connectivity, whereas Auvergne populations showed the highest genomic offsets and lowest neutral and adaptive genetic diversity. Only a reduced percentage (<1%) of the current distribution is predicted to face climatic conditions outside the current range, suggesting that adaptive variability required to adapt to future climates may already be present, and that assisted gene flow could represent an effective conservation strategy. Finally, we discuss some of the main challenges of genomic forecasts, particularly for declining non-model species.

Aedes aegypti is a major vector of arboviruses, and belongs to subfamily Culicinae, a diverse group of mosquitoes with homomorphic sex-determining chromosomes. Males are the heterogametic sex with a dominant male-determining locus (M locus). The M locus and its counterpart m locus are embedded in a region of suppressed recombination, with a large portion of this recombination desert showing significant molecular differentiation despite homomorphy. We developed a mathematical framework to examine M-linked genome editors that specifically target the m-chromosome during spermatogenesis, mimicking the naturally occurring sex ratio distorters (SRDs) in Culicinae that produce male-biased meiotic drives. Unlike previous models for species with heteromorphic sex chromosomes (e.g., X and Y), we incorporate features stemming from the homomorphic nature of the Ae. aegypti sex chromosomes such as varied linkage to the M locus, making the degree of super-Mendelian inheritance readily tunable. We evaluated in silico SRDs with a range of M-linkage and editing efficiencies and established the theoretical foundation for developing highly efficient SRDs that outperform several methods of population suppression. These SRDs can be tuned to mitigate impact on a neighboring population. The framework developed here is suitable for exploring SRD-mediated genetic biocontrol of pests with homomorphic sex chromosomes.