Population Subdivision and Migration Assessment of Mangalica Pig Breeds Based on Pedigree Analysis.

In conserving the genetic diversity of domestic animal breeds, strategies that emphasise between-breed diversity may not be optimal, as they neglect within-breed variation. The aim of the present study was to assess the extent of population subdivision in three Mangalica pig breeds and the contribution of migration to their substructure. Wright's FST coefficient was calculated based on genealogical data with breeding animals born between 1981 and 2023, with three colour variants (Blonde, Swallow-Belly and Red). These Wright's FST coefficients were analysed using multidimensional scaling to reveal the population substructure. The average FST coefficient was 0.04 for the Blonde breed and 0.047 for the Swallow-Belly and Red Mangalica breeds, while these parameters were lower in the active herds at 0.03 and 0.04, respectively. The migration of individuals between herds was 61.63% for the Blonde breed and 75.53% and 63.64% for the Swallow-Belly and Red Magalica breeds, respectively. No population substructure was observed in any of the Mangalica breeds, which can be explained by the extensive migration between herds.

The fitness consequences of genetic divergence between polymorphic gene arrangements.

Inversions restrict recombination when heterozygous with standard arrangements, but often have few noticeable phenotypic effects. Nevertheless, there are several examples of inversions that can be maintained polymorphic by strong selection under laboratory conditions. A long-standing model for the source of such selection is divergence between arrangements with respect to recessive or partially recessive deleterious mutations, resulting in a selective advantage to heterokaryotypic individuals over homokaryotypes. This paper uses a combination of analytical and numerical methods to investigate this model, for the simple case of an autosomal inversion with multiple independent nucleotide sites subject to mildly deleterious mutations. A complete lack of recombination in heterokaryotypes is assumed, as well as constancy of the frequency of the inversion over space and time. It is shown that a significantly higher mutational load will develop for the less frequent arrangement. A selective advantage to heterokaryotypes is only expected when the two alternative arrangements are nearly equal in frequency, so that their mutational loads are very similar in size. The effects of some Drosophila pseudoobscura polymorphic inversions on fitness traits seem to be too large to be explained by this process, although it may contribute to some of the observed effects. Several population genomic statistics can provide evidence for signatures of a reduced efficacy of selection associated with the rarer of two arrangements, but there is currently little published data that are relevant to the theoretical predictions.

Genomic Insights Into the Evolution and Demographic History of the SARS-CoV-2 Omicron Variant: Population Genomics Approach.

Background: A thorough understanding of the patterns of genetic subdivision in a pathogen can provide crucial information that is necessary to prevent disease spread. For SARS-CoV-2, the availability of millions of genomes makes this task analytically challenging, and traditional methods for understanding genetic subdivision often fail. Objective: The aim of our study was to use population genomics methods to identify the subtle subdivisions and demographic history of the Omicron variant, in addition to those captured by the Pango lineage. Methods: We used a combination of an evolutionary network approach and multivariate statistical protocols to understand the subdivision and spread of the Omicron variant. We identified subdivisions within the BA.1 and BA.2 lineages and further identified the mutations associated with each cluster. We further characterized the overall genomic diversity of the Omicron variant and assessed the selection pressure for each of the genetic clusters identified. Results: We observed concordant results, using two different methods to understand genetic subdivision. The overall pattern of subdivision in the Omicron variant was in broad agreement with the Pango lineage definition. Further, 1 cluster of the BA.1 lineage and 3 clusters of the BA.2 lineage revealed statistically significant signatures of selection or demographic expansion (Tajima's D<-2), suggesting the role of microevolutionary processes in the spread of the virus. Conclusions: We provide an easy framework for assessing the genetic structure and demographic history of SARS-CoV-2, which can be particularly useful for understanding the local history of the virus. We identified important mutations that are advantageous to some lineages of Omicron and aid in the transmission of the virus. This is crucial information for policy makers, as preventive measures can be designed to mitigate further spread based on a holistic understanding of the variability of the virus and the evolutionary processes aiding its spread.

The effects of inversion polymorphisms on patterns of neutral genetic diversity.

The strong reduction in the frequency of recombination in heterozygotes for an inversion and a standard gene arrangement causes the arrangements to become partially isolated genetically, resulting in sequence divergence between them and changes in the levels of neutral variability at nucleotide sites within each arrangement class. Previous theoretical studies on the effects of inversions on neutral variability have assumed either that the population is panmictic or that it is divided into 2 populations subject to divergent selection. Here, the theory is extended to a model of an arbitrary number of demes connected by migration, using a finite island model with the inversion present at the same frequency in all demes. Recursion relations for mean pairwise coalescent times are used to obtain simple approximate expressions for diversity and divergence statistics for an inversion polymorphism at equilibrium under recombination and drift, and for the approach to equilibrium following the sweep of an inversion to a stable intermediate frequency. The effects of an inversion polymorphism on patterns of linkage disequilibrium are also examined. The reduction in effective recombination rate caused by population subdivision can have significant effects on these statistics. The theoretical results are discussed in relation to population genomic data on inversion polymorphisms, with an emphasis on Drosophila melanogaster. Methods are proposed for testing whether or not inversions are close to recombination-drift equilibrium, and for estimating the rate of recombinational exchange in heterozygotes for inversions; difficulties involved in estimating the ages of inversions are also discussed.

Evolutionary Genomics of a Subdivided Species.

The ways in which genetic variation is distributed within and among populations is a key determinant of the evolutionary features of a species. However, most comprehensive studies of these features have been restricted to studies of subdivision in settings known to have been driven by local adaptation, leaving our understanding of the natural dispersion of allelic variation less than ideal. Here, we present a geographic population-genomic analysis of 10 populations of the freshwater microcrustacean Daphnia pulex, an emerging model system in evolutionary genomics. These populations exhibit a pattern of moderate isolation-by-distance, with an average migration rate of 0.6 individuals per generation, and average effective population sizes of ∼650,000 individuals. Most populations contain numerous private alleles, and genomic scans highlight the presence of islands of excessively high population subdivision for more common alleles. A large fraction of such islands of population divergence likely reflect historical neutral changes, including rare stochastic migration and hybridization events. The data do point to local adaptive divergence, although the precise nature of the relevant variation is diffuse and cannot be associated with particular loci, despite the very large sample sizes involved in this study. In contrast, an analysis of between-species divergence highlights positive selection operating on a large set of genes with functions nearly nonoverlapping with those involved in local adaptation, in particular ribosome structure, mitochondrial bioenergetics, light reception and response, detoxification, and gene regulation. These results set the stage for using D. pulex as a model for understanding the relationship between molecular and cellular evolution in the context of natural environments.

ddRAD-Seq reveals evolutionary insights into population differentiation and the cryptic phylogeography of Hyporhamphus intermedius in Mainland China.

Species differentiation and local adaptation in heterogeneous environments have attracted much attention, although little is known about the mechanisms involved. Hyporhamphus intermedius is an anadromous, brackish-water halfbeak that is widely distributed in coastal areas and hyperdiverse freshwater systems in China, making it an interesting model for research on phylogeography and local adaptation. Here, 156 individuals were sampled at eight sites from heterogeneous aquatic habitats to examine environmental and genetic contributions to phenotypic divergence. Using double-digest restriction-site-associated DNA sequencing (ddRAD-Seq) in the specimens from the different watersheds, 5498 single nucleotide polymorphisms (SNPs) were found among populations, with obvious population differentiation. We find that present-day Mainland China populations are structured into distinct genetic clusters stretching from southern and northern ancestries, mirroring geography. Following a transplant event in Plateau Lakes, there were virtually no variations of genetic diversity occurred in two populations, despite the fact two main splits were unveiled in the demographic history. Additionally, dorsal, and anal fin traits varied widely between the southern group and the others, which highlighted previously unrecognized lineages. We then explore genotype-phenotype-environment associations and predict candidate loci. Subgroup ranges appeared to correspond to geographic regions with heterogeneous hydrological factors, indicating that these features are likely important drivers of diversification. Accordingly, we conclude that genetic and phenotypic polymorphism and a moderate amount of genetic differentiation occurred, which might be ascribed to population subdivision, and the impact of abiotic factors.

Population subdivision promoted by a sea-level-change-driven bottleneck: A glimpse from the evolutionary history of the mangrove plant Aegiceras corniculatum.

Historic climate changes drive geographical populations of coastal plants to contract and recover dynamically, even die out completely. Species suffering from such bottlenecks usually lose intraspecific genetic diversity, but how do these events influence population subdivision patterns of coastal plants? Here, we investigated this question in the typical coastal plant: mangrove species Aegiceras corniculatum. Inhabiting the intertidal zone of the tropical and subtropical coast of the Indo-West Pacific oceans, its populations are deemed to be greatly shaped by historic sea-level fluctuations. Using dual methods of Sanger and Illumina sequencing, we found that the 18 sampled populations were structured into two groups, namely, the "Indo-Malayan" group, comprising three subgroups (the northern South China Sea, Gulf of Bengal, and Bali), and the "Pan-Australasia" group, comprising the subgroups of the southern South China Sea and Australasia. Based on the approximate Bayesian computations and Stairway Plot, we inferred that the southern South China Sea subgroup, which penetrates the interior of the "Indo-Malayan" group, originated from the Australasia subgroup, accompanied by a severe bottleneck event, with a spot of gene flow from both the Australasia and "Indo-Malayan" groups. Geographical barriers such as the Sundaland underlie the genetic break between Indian and Pacific Oceans, but the discontinuity between southern and northern South China Sea was originated from genetic drift in the bottleneck event. Hence, we revealed a case evidencing that the bottleneck event promoted population subdivision. This conclusion may be applicable in other taxa beyond coastal plants.

Metapopulation structure modulates sexual antagonism.

Despite the far-reaching evolutionary implications of sexual conflict, the effects of metapopulation structure, when populations are subdivided into several demes connected to some degree by migration, on sexual conflict dynamics are unknown. Here, we used experimental evolution in an insect model system, the seed beetle Callosobruchus maculatus, to assess the independent and interacting effects of selection histories associated with mating system (monogamy vs. polygamy) and population subdivision on sexual conflict evolution. We confirm traditional predictions from sexual conflict theory by revealing increased resistance to male harm in females from populations with a history of intense sexual selection (polygamous populations) compared to females from populations with a history of relaxed sexual selection (monogamous populations). However, selection arising from metapopulation structure reversed the classic pattern of sexually antagonistic coevolution and led to reduced resistance in females from polygamous populations. These results underscore that population spatial structure moderates sexual selection and sexual conflict, and more broadly, that the evolution of sexual conflict is contingent on ecological context. The findings also have implications for population dynamics, conservation biology, and biological control.

Genetic nurturing, missing heritability, and causal analysis in genetic statistics.

Genetic nurturing, the effect of parents' genotypes on offspring phenotypes through parental phenotypic transmission, can be modeled in terms of gene-culture interactions. This paper first uses a simple one-locus, two-phenotype gene-culture cotransmission model to compute the effect of genetic nurturing in terms of regression of children's phenotypes on transmitted and nontransmitted alleles. With genetic nurturing, interpreting heritability and hence the meaning of "missing heritability" becomes problematic. Other factors, for example, population subdivision and assortative mating, generate similar signals to those of genetic nurturing, namely, correlation between parents' nontransmitted alleles and children's phenotypes. Corrections must be made for these to isolate the signal of genetic nurturing. Finally, a unified causal framework is constructed for genetic nurturing, population subdivision, and assortative mating. Causal and noncausal paths from transmitted and nontransmitted alleles to children's phenotypes are identified and investigated in the presence of genetic nurturing, population subdivision, and assortative mating. Using causal analysis, assumptions made in inferring direct and indirect effects are then clarified and evaluated in a broader causal context.

Is population subdivision different from speciation? From phylogeography to species delimitation.

Species-level diversity and the underlying mechanisms that lead to the formation of new species, that is, speciation, have often been confounded with intraspecific diversity and population subdivision. The delineation between intraspecific and interspecific divergence processes has received much less attention than species delimitation. The ramifications of confounding speciation and population subdivision are that the term speciation has been used to describe many different biological divergence processes, rendering the results, or inferences, between studies incomparable. Phylogeographic studies have advanced our understanding of how spatial variation in the pattern of biodiversity can begin, become structured, and persist through time. Studies of species delimitation have further provided statistical and model-based approaches to determine the phylogeographic entities that merit species status. However, without a proper understanding and delineation between the processes that generate and maintain intraspecific and interspecific diversity in a study system, the delimitation of species may still not be biologically and evolutionarily relevant. I argue that variation in the continuity of the divergence process among biological systems could be a key factor leading to the enduring contention in delineating divergence patterns, or species delimitation, meriting future comparative studies to help us better understand the nature of biological species.

Genetic analyses reveal population structure and recent decline in leopards (Panthera pardus fusca) across the Indian subcontinent.

BACKGROUND: Large carnivores maintain the stability and functioning of ecosystems. Currently, many carnivore species face declining population sizes due to natural and anthropogenic pressures. The leopard, Panthera pardus, is probably the most widely distributed and highly adaptable large felid globally, still persisting in most of its historic range. However, we lack subspecies-level data on country or regional scale on population trends, as ecological monitoring approaches are difficult to apply on such wide-ranging species. We used genetic data from leopards sampled across the Indian subcontinent to investigate population structure and patterns of demographic decline. METHODS: We collected faecal samples from the Terai-Arc landscape of northern India and identified 56 unique individuals using a panel of 13 microsatellite markers. We merged this data with already available 143 leopard individuals and assessed genetic structure at country scale. Subsequently, we investigated the demographic history of each identified subpopulations and compared genetic decline analyses with countrywide local extinction probabilities. RESULTS: Our genetic analyses revealed four distinct subpopulations corresponding to Western Ghats, Deccan Plateau-Semi Arid, Shivalik and Terai region of the north Indian landscape, each with high genetic variation. Coalescent simulations with microsatellite loci revealed a possibly human-induced 75-90% population decline between ∼120-200 years ago across India. Population-specific estimates of genetic decline are in concordance with ecological estimates of local extinction probabilities in these subpopulations obtained from occupancy modeling of the historic and current distribution of leopards in India. CONCLUSIONS: Our results confirm the population decline of a widely distributed, adaptable large carnivore. We re-iterate the relevance of indirect genetic methods for such species in conjunction with occupancy assessment and recommend that detailed, landscape-level ecological studies on leopard populations are critical to future conservation efforts. Our approaches and inference are relevant to other widely distributed, seemingly unaffected carnivores such as the leopard.

Performing Parentage Analysis in the Presence of Inbreeding and Null Alleles.

Parentage analysis is an important method that is used widely in zoological and ecological studies. Current mathematical models of parentage analyses usually assume that a population has a uniform genetic structure and that mating is panmictic. In a natural population, the geographic or social structure of a population, and/or nonrandom mating, usually leads to a genetic structure and results in genotypic frequencies deviating from those expected under the Hardy-Weinberg equilibrium (HWE). In addition, in the presence of null alleles, an observed genotype represents one of several possible true genotypes. The true father of a given offspring may thus be erroneously excluded in parentage analyses, or may have a low or negative LOD score. Here, we present a new mathematical model to estimate parentage that includes simultaneously the effects of inbreeding, null alleles, and negative amplification. The influences of these three factors on previous model are evaluated by Monte-Carlo simulations and empirical data, and the performance of our new model is compared under controlled conditions. We found that, for both simulated and empirical data, our new model outperformed other methods in many situations. We make available our methods in a new, free software package entitled parentage This can be downloaded via http://github.com/huangkang1987/parentage.

Volatility in the effective size of a freshwater gastropod population.

Despite the utility of gastropod models for the study of evolutionary processes of great generality and importance, their effective population size has rarely been estimated in the field. Here, we report allele frequency variance at three allozyme-encoding loci monitored over 7 years in a population of the invasive freshwater pulmonate snail Physa acuta (Draparnaud 1805), estimating effective population size with both single-sample and two-sample approaches. Estimated Ne declined from effectively infinite in 2009 to approximately 40-50 in 2012 and then rose back to infinity in 2015, corresponding to a striking fluctuation in the apparent census size of the population. Such volatility in Ne may reflect cryptic population subdivision.

Worldwide population structure of the wheat rust fungus Puccinia striiformis in the past.

Puccinia striiformis is a basidiomycete causing yellow rust on wheat. The availability of historic samples of this pathogen from the 'Stubbs collection' enabled us to investigate past population structure and temporal dynamics on a global scale. A set of 212 single genotype urediniospore isolates, representing samples collected from five continents between 1958 and 1991, were genotyped using 19 polymorphic microsatellite markers. The population genetic analyses revealed the existence of seven genetic groups in the past worldwide P. striiformis population. This genetic grouping generally corresponded with geographical sample origin except for the Middle East, where six of the seven genetic groups were represented. The presence of many genetic groups in the Middle Eastern population reflected a low differentiation from the populations in East Africa (FST=0.052) and in South Asia (FST=0.064). A high diversity and recombinant population structure was observed in China and South Asia, while a clonal population structure was observed in NW Europe, East Africa and the Mediterranean region. The high genetic diversity in the Himalayan region supported recent studies suggesting a putative center of diversity for P. striiformis in this area. Four of the 89 multilocus genotypes detected were resampled in different geographical regions suggesting long-distance migration in the past. Comparison of the past populations with more recent ones, represented by 309 isolates mainly collected between 2001 and 2009, revealed temporal divergence for all populations except for Northwest Europe. Overall, we observed a clear subdivision within the worldwide population structure of P. striiformis and migration in the past.

Population subdivision of the surf clam Mactra chinensis in the East China Sea: Changjiang River outflow is not the sole driver.

The northwestern Pacific, characterized by unique tectonic and hydrological settings, has greatly intrigued marine phylogeographers. However, current studies mostly focus on the influence of Pleistocene isolation of sea basins in population structure of species in the region, leaving the contribution of other factors (such as freshwater outflow and environmental gradients) largely unexploited. Here we shed light on the question by investigating phylogeography of the surf clam Mactra chinensis in the East China Sea (ECS). Genetic information was acquired from 501 specimens collected from its main distribution in the region, represented by mitochondrial cytochrome oxidase I (COI) and nine polymorphic microsatellite loci. A shallow and star-like phylogeny was revealed for all COI haplotypes, indicating the origin of populations from a single refugium. Although no divergent lineages existed, population subdivision was detected in both data sets. The most striking pattern was the significant differentiation between populations north and south of a biogeographic boundary-the Changjiang Estuary, suggesting a barrier effect of the freshwater outflow to gene flow. For the northern group, substructure was revealed by COI result as one southernmost population was significant different from other ones. Clear latitude gradations in allele frequencies were revealed by microsatellite analyses, likely influenced by environmental gradient factors such as temperature. Our results demonstrate that genetic subdivision can arise for populations within the ECS despite they have a single origin, and multiple mechanisms including Changjiang River outflow, environmental gradient factors and life-history traits may act in combination in the process.

The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations.

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and FST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.

Emergence of long-term balanced polymorphism under cyclic selection of spatially variable magnitude.

A fundamental question in evolutionary biology is what promotes genetic variation at nonneutral loci, a major precursor to adaptation in changing environments. In particular, balanced polymorphism under realistic evolutionary models of temporally varying environments in finite natural populations remains to be demonstrated. Here, we propose a novel mechanism of balancing selection under temporally varying fitnesses. Using forward-in-time computer simulations and mathematical analysis, we show that cyclic selection that spatially varies in magnitude, such as along an environmental gradient, can lead to elevated levels of nonneutral genetic polymorphism in finite populations. Balanced polymorphism is more likely with an increase in gene flow, magnitude and period of fitness oscillations, and spatial heterogeneity. This polymorphism-promoting effect is robust to small systematic fitness differences between competing alleles or to random environmental perturbation. Furthermore, we demonstrate analytically that protected polymorphism arises as spatially heterogeneous cyclic fitness oscillations generate a type of storage effect that leads to negative frequency dependent selection. Our findings imply that spatially variable cyclic environments can promote elevated levels of nonneutral genetic variation in natural populations.

Clines in quantitative traits: the role of migration patterns and selection scenarios.

The existence, uniqueness, and shape of clines in a quantitative trait under selection toward a spatially varying optimum is studied. The focus is on deterministic diploid two-locus n-deme models subject to various migration patterns and selection scenarios. Migration patterns may exhibit isolation by distance, as in the stepping-stone model, or random dispersal, as in the island model. The phenotypic optimum may change abruptly in a single environmental step, more gradually, or not at all. Symmetry assumptions are imposed on phenotypic optima and migration rates. We study clines in the mean, variance, and linkage disequilibrium (LD). Clines result from polymorphic equilibria. The possible equilibrium configurations are determined as functions of the migration rate. Whereas for weak migration, many polymorphic equilibria may be simultaneously stable, their number decreases with increasing migration rate. Also for intermediate migration rates polymorphic equilibria are in general not unique, however, for loci of equal effects the corresponding clines in the mean, variance, and LD are unique. For sufficiently strong migration, no polymorphism is maintained. Both migration pattern and selection scenario exert strong influence on the existence and shape of clines. The results for discrete demes are compared with those from models in which space varies continuously and dispersal is modeled by diffusion. Comparisons with previous studies, which investigated clines under neutrality or under linkage equilibrium, are performed. If there is no long-distance migration, the environment does not change abruptly, and linkage is not very tight, populations are almost everywhere close to linkage equilibrium.

Population structure and the rate of evolution.

The way population size, population structure (with migration), and spatially dependent selection (where there is no globally optimal allele), combine to affect the substitution rate is poorly understood. Here, we consider a two patch model where mutant alleles are beneficial in one patch and deleterious in the other patch. We assume that the spatial average of selection on mutant alleles is zero. We take each patch to maintain a finite number of N adults each generation, hence random genetic drift can independently occur in each patch. We show that the principal way the population size, N, when large, affects the substitution rate, R∞, is through its dependence on two composite parameters. These are the scaled migration rate M (∝ population size × migration rate), and the scaled selection intensity S (∝population size × beneficial effect of a mutant). Any relation between S and M that arises for ecological/evolutionary reasons can strongly influence the way the substitution rate, R∞, depends on the population size, N. In the simplest situation, both M and S are proportional to N, and this is shown to lead to R∞ increasing with N when S is not large. The behaviour, that R∞ increases with N, is not inevitable; a more complex relation between S and M can lead to the opposite or other behaviours. In particular, let us assume that dM/dN is positive, as would occur if the migration rate were constant, S is not large, and S depends on M (i.e., S=S(M)). We then find that if S(M) satisfies S(M)>((1+M)/1+2M)S(0) then the substitution rate, R∞, increases with N, but if S(M)<((1+M)/1+2M)S(0) then R∞ decreases with N.

Turnover and accumulation of genetic diversity across large time-scale cycles of isolation and connection of populations.

Major climatic and geological events but also population history (secondary contacts) have generated cycles of population isolation and connection of long and short periods. Recent empirical and theoretical studies suggest that fast evolutionary processes might be triggered by such events, as commonly illustrated in ecology by the adaptive radiation of cichlid fishes (isolation and reconnection of lakes and watersheds) and in epidemiology by the fast adaptation of the influenza virus (isolation and reconnection in hosts). We test whether cyclic population isolation and connection provide the raw material (standing genetic variation) for species evolution and diversification. Our analytical results demonstrate that population isolation and connection can provide, to populations, a high excess of genetic diversity compared with what is expected at equilibrium. This excess is either cyclic (high allele turnover) or cumulates with time depending on the duration of the isolation and the connection periods and the mutation rate. We show that diversification rates of animal clades are associated with specific periods of climatic cycles in the Quaternary. We finally discuss the importance of our results for macroevolutionary patterns and for the inference of population history from genomic data.