In recent years, Genomic Selection (GS) has accelerated genetic gain in dairy cattle breeds worldwide. Despite the evident genetic progress, several dairy populations have also encountered challenges ...such as heightened inbreeding rates and reduced effective population sizes. The challenge has been to find a balance between achieving substantial genetic gain while managing genetic diversity within the population, thereby mitigating the negative effects of inbreeding depression. This study aims to elucidate the impact of GS on pedigree and genomic rates of inbreeding (ΔF) and coancestry (ΔC) in Nordic Jersey (NJ) and Holstein (NH) cattle populations. Furthermore, key genetic metrics including the generation interval (L), effective population size (Ne), and future effective population size (FNe) were assessed between 2 time periods, before and after GS, and across distinct animal cohorts in both breeds: females, bulls, and approved semen-producing bulls (AI-sires).
Analysis of ΔF and ΔC revealed distinct trends across the studied periods and animal groups. Notably, there was a consistent increase in yearly ΔF for most animal groups in both breeds. An exception was observed in NH AI-sires, which demonstrated a slight decrease in yearly ΔF. Moreover, NJ displayed minimal changes in yearly ΔC between the periods, whereas NH exhibited elevated ΔC values across all animal groups. Particularly striking was the substantial increase in yearly ΔC within the NH female population, surging from 0.02% to 0.39% between the periods.
Implementation of GS resulted in a reduction of the generation interval across all animal cohorts in both NJ and NH breeds. However, the extent of reduction was more pronounced in males compared with females. This reduction in generation interval influenced generational changes in ΔF and ΔC. Bulls and AI-sires of both breeds exhibited reduced generational ΔF between periods, in contrast to females that demonstrated an opposing pattern.
Between the periods, NJ maintained a relatively stable Ne, 29.4 before and 30.3 after GS, while NH experienced a notable decline from 54.3 to 42.8. Female groups in both breeds displayed a negative Ne trend, while males demonstrated either neutral or positive Ne developments. Regarding FNe, NJ exhibited positive FNe development with an increase from 40.7 to 57.2. The opposite was observed in NH, where FNe decreased from 198.8 to 42.7.
In summary, it was evident that the genomic methods could detect differences between the populations and changes in ΔF and ΔC more efficiently than pedigree methods. GS implementation yielded positive outcomes within the NJ population regarding the rate of coancestry but the opposite was observed with NH. Moreover, analysis of ΔC data hints at the potential to decrease future ΔF through informed mating strategies. Conversely, NH faces more pressing concerns, even though ΔF remains comparatively modest in contrast to what has been observed in other Holstein populations. These findings underscore the necessity of genomic control of inbreeding and coancestry with strategic changes in the Nordic breeding schemes for dairy to ensure long-term sustainability in the forthcoming years.
•Genetic factors remain inadequately addressed in conservation management.•Effective population size (Ne)=50 does not prevent inbreeding depression.•Ne⩾100 is required to limit inbreeding depression ...to 10% over 5 generations.•Ne=500 is too low for retaining evolutionary potential; Ne⩾1000 is required.•IUCN Red List Criterion C thresholds for population size require doubling.
Conservation managers typically need to make prompt decisions based on limited information and resources. Consequently, generalisations have essential roles in guiding interventions. Here, we (i) critique information on some widely accepted generalisations and variables affecting them, (ii) assess how adequately genetic factors are currently incorporated into population viability analysis (PVA) models used to estimate minimum viable population sizes, and (iii) relate the above to population size thresholds of the IUCN Red List criteria for threatened species that were derived from genetic considerations. Evidence accumulated since 1980 shows that genetically effective population size (Ne)=50 is inadequate for preventing inbreeding depression over five generations in the wild, with Ne⩾100 being required to limit loss in total fitness to ⩽10%. Further, even Ne=500 is too low for retaining evolutionary potential for fitness in perpetuity; a better approximation is Ne⩾1000. Extrapolation from census population size (N) to Ne depends on knowing the ratio of Ne/N, yet this information is unavailable for most wild populations. Ratio averages (∼0.1–0.2) from meta-analyses are sufficient, provided adjustments are made for dissimilar life histories. Most PVA-based risk assessments ignore or inadequately model genetic factors. PVA should routinely include realistic inbreeding depression, and genetic impacts on evolutionary potential should be incorporated where appropriate. Genetic generalisations used in conservation, the treatment of genetics in PVAs, and sections of the IUCN Red List criteria derived from genetic considerations, all require revision to be more effective conservation tools.
DNA-based species delimitation may be compromised by limited sampling effort and species rarity, including "singleton" representatives of species, which hampers estimates of intra-versus interspecies ...evolutionary processes. In a case study of southern African chafers (beetles in the family Scarabaeidae), many species and subclades were poorly represented and 48.5% of species were singletons. Using cox1 sequences from >500 specimens and ~100 species, the Generalized Mixed Yule Coalescent (GMYC) analysis as well as various other approaches for DNA-based species delimitation (Automatic Barcode Gap Discovery (ABGD), Poisson tree processes (PTP), Species Identifier, Statistical Parsimony), frequently produced poor results if analyzing a narrow target group only, but the performance improved when several subclades were combined. Hence, low sampling may be compensated for by "clade addition" of lineages outside of the focal group. Similar findings were obtained in reanalysis of published data sets of taxonomically poorly known species assemblages of insects from Madagascar. The low performance of undersampled trees is not due to high proportions of singletons per se, as shown in simulations (with 13%, 40% and 52% singletons). However, the GMYC method was highly sensitive to variable effective population size (Ne), which was exacerbated by variable species abundances in the simulations. Hence, low sampling success and rarity of species affect the power of the GMYC method only if they reflect great differences in Ne among species. Potential negative effects of skewed species abundances and prevalence of singletons are ultimately an issue about the variation in Ne and the degree to which this is correlated with the census population size and sampling success. Clade addition beyond a limited study group can overcome poor sampling for the GMYC method in particular under variable Ne. This effect was less pronounced for methods of species delimitation not based on coalescent models.
Few complete human genomes from the European Early Upper Palaeolithic (EUP) have been sequenced. Using novel sampling and DNA extraction approaches, we sequenced the genome of a woman from “Peştera ...Muierii,” Romania who lived ∼34,000 years ago to 13.5× coverage. The genome shows similarities to modern-day Europeans, but she is not a direct ancestor. Although her cranium exhibits both modern human and Neanderthal features, the genome shows similar levels of Neanderthal admixture (∼3.1%) to most EUP humans but only half compared to the ∼40,000-year-old Peştera Oase 1. All EUP European hunter-gatherers display high genetic diversity, demonstrating that the severe loss of diversity occurred during and after the Last Glacial Maximum (LGM) rather than just during the out-of-Africa migration. The prevalence of genetic diseases is expected to increase with low diversity; however, pathogenic variant load was relatively constant from EUP to modern times, despite post-LGM hunter-gatherers having the lowest diversity ever observed among Europeans.
•Peştera Muierii woman is related to Europeans, but she is not a direct ancestor•Reduced diversity in Europe caused by Last Glaciation, not out-of-Africa bottleneck•Genetic load appears indifferent across 40,000 years of European history•New DNA extraction approach recovers up to 33 times more DNA from ancient remains
Svensson et al. sequence the complete genome of a woman from “Peştera Muierii,” Romania, who lived 34,000 years ago. Her genome is similar to modern-day Europeans, but she is not a direct ancestor. Her genome shows high levels of diversity, revealing that much loss of diversity in non-Africans occurred after she lived rather than before her time.
Effective population size is a fundamental parameter in population genetics, evolutionary biology, and conservation biology, yet its estimation can be fraught with difficulties. Several methods to ...estimate Ne from genetic data have been developed that take advantage of various approaches for inferring Ne. The ability of these methods to accurately estimate Ne, however, has not been comprehensively examined. In this study, we employ seven of the most cited methods for estimating Ne from genetic data (Colony2, CoNe, Estim, MLNe, ONeSAMP, TMVP, and NeEstimator including LDNe) across simulated datasets with populations experiencing migration or no migration. The simulated population demographies are an isolated population with no immigration, an island model metapopulation with a sink population receiving immigrants, and an isolation by distance stepping stone model of populations. We find considerable variance in performance of these methods, both within and across demographic scenarios, with some methods performing very poorly. The most accurate estimates of Ne can be obtained by using LDNe, MLNe, or TMVP; however each of these approaches is outperformed by another in differing demographic scenario. Knowledge of the approximate demography of population as well as the availability of temporal data largely improves Ne estimates.
Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have ...previously been hindered by a lack of means and genetic data to model changes in effective population size (Nₑ) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in Nₑ over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black‐and‐white flycatchers. Nₑ curves of Atlas, collared, pied and semicollared flycatcher converged 1–2 million years ago (Ma) at an Nₑ of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate Nₑ trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100–200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in Nₑ differed among species, and for pied flycatcher, the temporal dynamics of Nₑ differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per‐site filter of ≥10 reads and no more than 25% of missing data.
Frequency-dependent selection and demographic fluctuations play important roles in evolutionary and ecological processes. Under frequency-dependent selection, the average fitness of the population ...may increase or decrease based on interactions between individuals within the population. This should be reflected in fluctuations of the population size even in constant environments. Here, we propose a stochastic model that naturally combines these two evolutionary ingredients by assuming frequency-dependent competition between different types in an individual-based model. In contrast to previous game theoretic models, the carrying capacity of the population, and thus the population size, is determined by pairwise competition of individuals mediated by evolutionary games and demographic stochasticity. In the limit of infinite population size, the averaged stochastic dynamics is captured by deterministic competitive LotkaâVolterra equations. In small populations, demographic stochasticity may instead lead to the extinction of the entire population. Because the population size is driven by fitness in evolutionary games, a population of cooperators is less prone to go extinct than a population of defectors, whereas in the usual systems of fixed size the population would thrive regardless of its average payoff.
This contribution breaks with the tradition to restrict stochastic evolutionary game dynamics to populations of constant size and introduces a theoretical framework to investigate relevant and natural changes arising in populations that vary in size according to fitnessâa feature common to many real biological systems. Explicitly including ecological variation can result in significant effects on the stochastic evolutionary trajectories while providing a transparent link to the established, deterministic LotkaâVolterra systems.
•Not all genomic regions and lengths of hepatitis E virus (HEV) are equally suited for use in genetic analyses.•Generally, longer nucleic acid stretches yield better accuracy in predicting the viral ...subtype.•Analysis of molecular and epidemiological data of HEV in the Netherlands indicates the ability of the 493nt fragment to detect short term changes in the viral population size.•In the absence of whole genome sequences the 493bp fragment can be used for molecular typing and population structure characterization of HEV.
The aim of this study was to investigate to what extent fragments of the HEV genome could be used for accurate diagnostics and inference of viral population-scale processes. For this, we selected all the published whole genome sequences from the NCBI GenBank and trimmed them to various fragment lengths (ORF1,2,3, ORF1, ORF2, ORF3, 493nt in ORF2 and 148 nt in ORF2). Each of the fragment lengths was used to infer the richness and diversity of the viral sequence types, typing accuracy, and potential use in phylodynamics. The results obtained from the different fragments were compared. We observed that, generally, the longer the nucleic acid fragment used in typing, the better the accuracy in predicting the viral subtype. However, the dominant HEV subtypes circulating in Europe were relatively well classified even by the 493nt fragment, with false negative rates as low as 8 in 1000 typed sequences. Most fragments also give comparable results in analyses of population size, albeit with shorter fragments showing a broader 95% highest posterior density interval and less obvious increase of the viral effective population size. The reconstructed phylogenies of a heterochronous subset indicated a good concordance between all the fragments, with the major clades following similar branching patterns. Furthermore, we have used the HEV sequence data from the Netherlands available in the HEVnet database as a case study for reconstruction of population size changes in the past decades. This data showed that molecular and epidemiological results are concordant and point to an increase in the viral effective population size underlying the observed increase in incidence of acute HEV infection cases. In the absence of whole genome sequencing data, the 493bp fragment can be used for analyzing HEV strains currently circulating in Europe, as it is informative for describing short term population-scale processes.
Genetic diversity in livestock populations is a significant contributor to the sustainability of animal production. Also, genetic diversity allows animal production to become more responsive to ...environmental changes and market demands. The loss of genetic diversity can result in a plateau in production and may also result in loss of fitness or viability in animal production. In this study, we investigated the rate of inbreeding (ΔF), rate of coancestry (Δf), and effective population size (Ne) as important quantitative indicators of genetic diversity and evaluated the effect of the recent implementation of genomic selection on the loss of genetic diversity in North American Holstein and Jersey dairy cattle. To estimate the rate of inbreeding and coancestry, inbreeding and coancestry coefficients were calculated using the traditional pedigree method and genomic methods estimated from segment- and marker-based approaches. Furthermore, we estimated Ne from the rate of inbreeding and coancestry and extent of linkage disequilibrium. A total of 205,755 and 89,238 pedigreed and genotyped animals born between 1990 and 2018 inclusively were available for Holsteins and Jerseys, respectively. The estimated average pedigree inbreeding coefficients were 7.74 and 7.20% for Holsteins and Jerseys, respectively. The corresponding values for the segment and marker-by-marker genomic inbreeding coefficients were 13.61, 15.64, and 31.40% for Holsteins and 21.16, 22.54, and 42.62% for Jerseys, respectively. The average coancestry coefficients were 8.33 and 15.84% for Holsteins and 9.23 and 23.46% for Jerseys with pedigree and genomic measures, respectively. Generation interval for the whole 29-yr time period averaged approximately 5 yr for all selection pathways combined. The ΔF per generation based on pedigree, segment, and marker-by-marker genomic measures for the entire 29-yr period was estimated to be 0.75, 1.10, 1.16, and 1.02% for Holstein animals and 0.67, 0.62, 0.63, and 0.59% for Jersey animals, respectively. The Δf was estimated to be 0.98 and 0.98% for Holsteins and 0.73 and 0.78% for Jerseys with pedigree and genomic measures, respectively. These ΔF and Δf translated to an Ne that ranged from 43 to 66 animals for Holsteins and 64 to 85 animals for Jerseys. In addition, the Ne based on linkage disequilibrium was 58 and 120 for Holsteins and Jerseys, respectively. The 10-yr period that involved the application of genomic selection resulted in an increased ΔF per generation with ranges from 1.19 to 2.06% for pedigree and genomic measures in Holsteins. Given the rate at which inbreeding is increasing after the implementation of genomic selection, there is a need to implement measures and means for controlling the rate of inbreeding per year, which will help to manage and maintain farm animal genetic resources.