Using long-term data for a reintroduced population to empirically estimate future consequences of inbreeding

Meeting Report on the Symposium "Evolutionary Applications" at the 3rd Joint Congress on Evolutionary Biology

The symposium "Evolutionary Applications" took place on June 28, 2024 in the virtual part of the 3rd Joint Congress on Evolutionary Biology. It was contributed to the conference by the European Society for Evolutionary Biology (ESEB). The symposium highlighted research on evolutionary biology applied to address questions and contemporary problems in medicine and public health, conservation biology, and food production and agriculture. Each of the six talks covered a different application and a different organism: domestication of cheese-making fungi, restoration of long-lived bird populations, evolution of herbicide resistance, coral reef conservation, gene drive systems targeting Malaria vectors, and antibiotic resistance evolution in bacteria. By including speakers who are active in a consortium or work in an NGO, the symposium also showed how to make the step from scientific findings to practical application. The symposium furthermore featured a range of scientific methods, ranging from genomic analyses and mathematical modeling to laboratory evolution and field experiments. Speakers from across 15 time zones highlighted the potential of virtual symposia to foster global collaboration in evolutionary biology.

What Can Genome Sequence Data Reveal About Population Viability?

Biologists have long sought to understand the impacts of deleterious genetic variation on fitness and population viability. However, our understanding of these effects in the wild is incomplete, in part due to the rarity of sufficient genetic and demographic data needed to measure their impact. The genomics revolution is promising a potential solution by predicting the effects of deleterious genetic variants (genetic load) bioinformatically from genome sequences alone bypassing the need for costly demographic data. After a historical perspective on the theoretical and empirical basis of our understanding of the dynamics and fitness effects of deleterious genetic variation, we evaluate the potential for these new genomic measures of genetic load to predict population viability. We argue that current genomic analyses alone cannot reliably predict the effects of deleterious genetic variation on population growth, because these depend on demographic, ecological and genetic parameters that need more than just genome sequence data to be measured. Thus, while purely genomic analyses of genetic load promise to improve our understanding of the composition of the genetic load, they are currently of little use for evaluating population viability. Demographic data and ecological context remain crucial to our understanding of the consequences of deleterious genetic variation for population fitness. However, when combined with such demographic and ecological data, genomic information can offer important insights into genetic variation and inbreeding that are crucial for conservation decision making.

Using experimental reintroductions to resolve the roles of habitat quality and metapopulation dynamics on patch occupancy in fragmented landscapes

Declines of species in fragmented landscapes can potentially be reversed either by restoring connectivity or restoring local habitat quality. Models fitted to snapshot occupancy data can be used to predict the effectiveness of these actions. However, such inferences can be misleading if the reliability of the habitat and landscape metrics used is unknown. The only way to unambiguously resolve the roles of habitat quality and metapopulation dynamics is to conduct experimental reintroductions to unoccupied patches so that habitat quality can be measured directly from data on vital rates. We, therefore, conducted a 15-year study that involved reintroducing a threatened New Zealand bird to unoccupied forest fragments to obtain reliable data on their habitat quality and reassess initial inferences made by modeling occupancy against habitat and landscape metrics. Although reproductive rates were similar among fragments, subtle differences in adult survival rates resulted in λ (finite rate of increase) estimations of <0.9 for 9 of the 12 fragments that were previously unoccupied. This was the case for only 1 of 14 naturally occupied fragments. This variation in λ largely explained the original occupancy pattern, reversing our original conclusion from occupancy modeling that this occupancy pattern was isolation driven and suggesting that it would be detrimental to increase connectivity without improving local habitat quality. These results illustrate that inferences from snapshot occupancy should be treated with caution and subjected to testing through experimental reintroductions in selected model systems.

Capturing the dynamics of small populations: A retrospective assessment using long-term data for an island reintroduction

The art of population modelling is to incorporate factors essential for capturing a population's dynamics while otherwise keeping the model as simple as possible. However, it is unclear how optimal model complexity should be assessed, and whether this optimal complexity has been affected by recent advances in modelling methodology. This issue is particularly relevant to small populations because they are subject to complex dynamics but inferences about those dynamics are often constrained by small sample sizes. We fitted Bayesian hierarchical models to long-term data on vital rates (survival and reproduction) for the toutouwai Petroica longipes population reintroduced to Tiritiri Matangi, a 220-ha New Zealand island, and quantified the performance of those models in terms of their likelihood of replicating the observed population dynamics. These dynamics consisted of overall growth from 33 (±0.3) to 160 (±6) birds from 1992-2018, including recoveries following five harvest events for further reintroductions to other sites. We initially included all factors found to affect vital rates, which included inbreeding, post-release effects (PRE), density-dependence, sex, age and random annual variation, then progressively removed these factors. We also compared performance of models where data analysis and simulations were done simultaneously to those produced with the traditional two-step approach, where vital rates are estimated first then fed into a separate simulation model. Parametric uncertainty and demographic stochasticity were incorporated in all projections. The essential factors for replicating the population's dynamics were density-dependence in juvenile survival and PRE, i.e. initial depression of survival and reproduction in translocated birds. Inclusion of other factors reduced the precision of projections, and therefore the likelihood of matching observed dynamics. However, this reduction was modest when the modelling was done in an integrated framework. In contrast, projections were much less precise when done with a two-step modelling approach, and the cost of additional parameters was much higher under the two-step approach. These results suggest that minimization of complexity may be less important than accounting for covariances in parameter estimates, which is facilitated by integrating data analysis and population projections using Bayesian methods.

Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi

Inbreeding can lead to a loss of heterozygosity in a population and when combined with genetic drift may reduce the adaptive potential of a species. However, there is uncertainty about whether resequencing data can provide accurate and consistent inbreeding estimates. Here, we performed an in-depth inbreeding analysis for hihi (Notiomystis cincta), an endemic and nationally vulnerable passerine bird of Aotearoa New Zealand. We first focused on subsampling variants from a reference genome male, and found that low-density data sets tend to miss runs of homozygosity (ROH) in some places and overestimate ROH length in others, resulting in contrasting homozygosity landscapes. Low-coverage resequencing and 50 K SNP array densities can yield comparable inbreeding results to high-coverage resequencing approaches, but the results for all data sets are highly dependent on the software settings employed. Second, we extended our analysis to 10 hihi where low-coverage whole genome resequencing, RAD-seq and SNP array genotypes are available. We inferred ROH and individual inbreeding to evaluate the relative effects of sequencing depth versus SNP density on estimating inbreeding coefficients and found that high rates of missingness downwardly bias both the number and length of ROH. In summary, when using genomic data to evaluate inbreeding, studies must consider that ROH estimates are heavily dependent on analysis parameters, data set density and individual sequencing depth.