Effective Population Size,Extended Linkage Disequilibrium and Signatures of Selection in the Rare Dog Breed Lundehund

The Lundehund is an old dog breed with remarkable anatomical features including polydactyly in all four limbs and extraordinary flexibility of the spine. We genotyped 28 Lundehund using the canine Illumina high density beadchip to estimate the effective population size (N_e) and inbreeding coefficients as well as to identify potential regions of positive selection. The decay of linkage disequilibrium was slow with r² = 0.95 in 50 kb distance. The last 7-200 generations ago, Ne was at 10-13. An increase of N_e was noted in the very recent generations with a peak value of 19 for N_e at generation 4. The FROH estimated for 50-, 65- and 358-SNP windows were 0.87, 087 and 0.81, respectively. The most likely estimates for F_ROH after removing identical-by-state segments due to linkage disequilibria were at 0.80-0.81. The extreme loss of heterozygosity has been accumulated through continued inbreeding over 200 generations within a probably closed population with a small effective population size. The mean inbreeding coefficient based on pedigree data for the last 11 generations (F_Ped = 0.10) was strongly biased downwards due to the unknown coancestry of the founders in this pedigree data. The long-range haplotype test identified regions with genes involved in processes of immunity, olfaction, woundhealing and neuronal development as potential targets of selection. The genes QSOX2, BMPR1B and PRRX2 as well as MYOM1 are candidates for selection on the Lundehund characteristics small body size, increased number of digits per paw and extraordinary mobility, respectively.     

Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy

Increased inbreeding is an inevitable consequence of selection in livestock populations. The analysis of high‐density single nucleotide polymorphisms (SNPs) facilitates the identification of long and uninterrupted runs of homozygosity (ROH) that can be used to identify chromosomal regions that are identical by descent. In this work, the distribution of ROH of different lengths in five Italian cattle breeds is described. A total of 4095 bulls from five cattle breeds (2093 Italian Holstein, 749 Italian Brown, 364 Piedmontese, 410 Marchigiana and 479 Italian Simmental) were genotyped at 54K SNP loci. ROH were identified and used to estimate molecular inbreeding coefficients (F_ROH), which were compared with inbreeding coefficients estimated from pedigree information (F_PED) and using the genomic relationship matrix (F_GRM). The average number of ROH per animal ranged from 54 ± 7.2 in Piedmontese to 94.6 ± 11.6 in Italian Brown. The highest number of short ROH (related to ancient consanguinity) was found in Piedmontese, followed by Simmental. The Italian Brown and Holstein had a higher proportion of longer ROH distributed across the whole genome, revealing recent inbreeding. The F_PED were moderately correlated with F_ROH > 1 Mb (0.662, 0.700 and 0.669 in Italian Brown, Italian Holstein and Italian Simmental respectively) but poorly correlated with F_GRM (0.134, 0.128 and 0.448 for Italian Brown, Italian Holstein and Italian Simmental respectively). The inclusion of ROH > 8 Mb in the inbreeding calculation improved the correlation of F_ROH with F_PED and F_GRM. ROH are a direct measure of autozygosity at the DNA level and can overcome approximations and errors resulting from incomplete pedigree data. In populations with high linkage disequilibrium (LD) and recent inbreeding (e.g. Italian Holstein and Italian Brown), a medium‐density marker panel, such as the one used here, may provide a good estimate of inbreeding. However, in populations with low LD and ancient inbreeding, marker density would have to be increased to identify short ROH that are identical by descent more precisely.     

The estimates of effective population size based on linkage disequilibrium are virtually unaffected by natural selection

The effective population size (N_e) is a key parameter to quantify the magnitude of genetic drift and inbreeding, with important implications in human evolution. The increasing availability of high-density genetic markers allows the estimation of historical changes in N_e across time using measures of genome diversity or linkage disequilibrium between markers. Directional selection is expected to reduce diversity and N_e, and this reduction is modulated by the heterogeneity of the genome in terms of recombination rate. Here we investigate by computer simulations the consequences of selection (both positive and negative) and recombination rate heterogeneity in the estimation of historical N_e. We also investigate the relationship between diversity parameters and N_e across the different regions of the genome using human marker data. We show that the estimates of historical N_e obtained from linkage disequilibrium between markers (N_eLD) are virtually unaffected by selection. In contrast, those estimates obtained by coalescence mutation-recombination-based methods can be strongly affected by it, which could have important consequences for the estimation of human demography. The simulation results are supported by the analysis of human data. The estimates of N_eLD obtained for particular genomic regions do not correlate, or they do it very weakly, with recombination rate, nucleotide diversity, proportion of polymorphic sites, background selection statistic, minor allele frequency of SNPs, loss of function and missense variants and gene density. This suggests that N_eLD measures mainly reflect demographic changes in population size across generations.     

Variable rates of random genetic drift in protected populations of English yew: implications for gene pool conservation

Protecting populations in their natural habitat allows for the maintenance of naturally evolved adaptations and ecological relationships. However, the conservation of genetic resources often requires complementary practices like gene banks, translocations or reintroductions. In order to minimize inbreeding depression and maximize the adaptive potential of future populations, populations chosen for ex situ conservation should be selected according to criteria that will result in a reduction of global coancestry in the population. Generally, large populations should reveal lower coancestry and higher genetic variation than small populations. If detailed knowledge about coancestry is lacking, census population number (N _c ) can be used as a proxy for required characteristics. However, a simple measure of N _c may be misleading in particular cases as genetic processes rely on effective population size (N _e ) rather than N _c and these two measures may differ substantially due to demographic processes. We used an example of English yew to address whether N _c can be a good predictor of genetic parameters when used in conservation programs. Using microsatellite markers, we estimated allelic richness, inbreeding and coancestry coefficients of six relatively large yew populations in Poland. Each population was characterized by N _e using the linkage disequilibrium method. Our results showed that populations of English yew were subject to substantial divergence and genetic drift, with both being inversely proportional to the effective subpopulation size (N _e ). Additionally, allelic richness appeared proportional to N _e but not to N _c . However, the N _e /N ratio differed greatly among populations, which was possibly due to different population histories. From the results we concluded that choosing source populations based only on their census size can be fairly misleading. Implications for conservation are briefly discussed.     

Comparative Analysis of Genome Diversity in Bullmastiff Dogs

Management and preservation of genomic diversity in dog breeds is a major objective for maintaining health. The present study was undertaken to characterise genomic diversity in Bullmastiff dogs using both genealogical and molecular analysis. Genealogical analysis of diversity was conducted using a database consisting of 16,378 Bullmastiff pedigrees from year 1980 to 2013. Additionally, a total of 188 Bullmastiff dogs were genotyped using the 170,000 SNP Illumina CanineHD Beadchip. Genealogical parameters revealed a mean inbreeding coefficient of 0.047; 142 total founders (f); an effective number of founders (f_e) of 79; an effective number of ancestors (f_a) of 62; and an effective population size of the reference population of 41. Genetic diversity and the degree of genome-wide homogeneity within the breed were also investigated using molecular data. Multiple-locus heterozygosity (MLH) was equal to 0.206; runs of homozygosity (ROH) as proportion of the genome, averaged 16.44%; effective population size was 29.1, with an average inbreeding coefficient of 0.035, all estimated using SNP Data. Fine-scale population structure was analysed using NETVIEW, a population analysis pipeline. Visualisation of the high definition network captured relationships among individuals within and between subpopulations. Effects of unequal founder use, and ancestral inbreeding and selection, were evident. While current levels of Bullmastiff heterozygosity, inbreeding and homozygosity are not unusual, a relatively small effective population size indicates that a breeding strategy to reduce the inbreeding rate may be beneficial.     

Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds

Single nucleotide polymorphism (SNP) genotyping tools, which can analyse thousands of SNPs covering the whole genome, have opened new opportunities to estimate the inbreeding level of animals directly using genome information. One of the most commonly used genomic inbreeding measures considers the proportion of the autosomal genome covered by runs of homozygosity (ROH), which are defined as continuous and uninterrupted chromosome portions showing homozygosity at all loci. In this study, we analysed the distribution of ROH in three commercial pig breeds (Italian Large White, n = 1968; Italian Duroc, n = 573; and Italian Landrace, n = 46) and four autochthonous breeds (Apulo-Calabrese, n = 90; Casertana, n = 90; Cinta Senese, n = 38; and Nero Siciliano, n = 48) raised in Italy, using SNP data generated from Illumina PorcineSNP60 BeadChip. We calculated ROH-based inbreeding coefficients (F_ROH) using ROH of different minimum length (1, 2, 4, 8, 16 Mbp) and compared them with several other genomic inbreeding coefficients (including the difference between observed and expected number of homozygous genotypes (F_HOM)) and correlated all these genomic-based measures with the pedigree inbreeding coefficient (F_PED) calculated for the pigs of some of these breeds. Autochthonous breeds had larger mean size of ROH than all three commercial breeds. F_HOM was highly correlated (0.671 to 0.985) with F_ROH measures in all breeds. Apulo-Calabrese and Casertana had the highest F_ROH values considering all ROH minimum lengths (ranging from 0.273 to 0.189 and from 0.226 to 0.152, moving from ROH of minimum size of 1 Mbp (F_ROH1) to 16 Mbp (F_ROH16)), whereas the lowest F_ROH values were for Nero Siciliano (from 0.072 to 0.051) and Italian Large White (from 0.117 to 0.042). F_ROH decreased as the minimum length of ROH increased for all breeds. Italian Duroc had the highest correlations between all F_ROH measures and F_PED (from 0.514 to 0.523) and between F_HOM and F_PED (0.485). Among all analysed breeds, Cinta Senese had the lowest correlation between F_ROH and F_PED. This might be due to the imperfect measure of F_PED, which, mainly in local breeds raised in extensive production systems, cannot consider a higher level of pedigree errors and a potential higher relatedness of the founder population. It appeared that ROH better captured inbreeding information in the analysed breeds and could complement pedigree-based inbreeding coefficients for the management of these genetic resources.     

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome‐wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.     

Genome‐wide linkage disequilibrium and past effective population size in three Korean cattle breeds

The routine collection and use of genomic data are useful for effectively managing breeding programs for endangered populations. Linkage disequilibrium (LD) using high‐density DNA markers has been widely used to determine population structures and predict the genomic regions that are associated with economic traits in beef cattle. The extent of LD also provides information about historical events, including past effective population size (N_e), and it allows inferences on the genetic diversity of breeds. The objective of this study was to estimate the LD and N_e in three Korean cattle breeds that are genetically similar but have different coat colors (Brown, Brindle and Jeju Black Hanwoo). Brindle and Jeju Black are endangered breeds with small populations, whereas Brown Hanwoo is the main breeding population in Korea. DNA samples from these cattle breeds were genotyped using the Illumina BovineSNP50 Bead Chip. We examined 13 cattle breeds, including European taurines, African taurines and indicines, and hybrids to compare their LD values. Brown Hanwoo consistently had the lowest mean LD compared to Jeju Black, Brindle and the other 13 cattle breeds (0.13, 0.19, 0.21 and 0.15–0.22 respectively). The high LD values of Brindle and Jeju Black contributed to small N_e values (53 and 60 respectively), which were distinct from that of Brown Hanwoo (531) for 11 generations ago. The differences in LD and N_e for each breed reflect the breeding strategy applied. The N_e for these endangered cattle breeds remain low; thus, effort is needed to bring them back to a sustainable tract.     

Effective Population Size,Genetic Variation,and Their Relevance for Conservation: The Bighorn Sheep in Tiburon Island and Comparisons with Managed Artiodactyls

The amount of genetic diversity in a finite biological population mostly depends on the interactions among evolutionary forces and the effective population size (N _e) as well as the time since population establishment. Because the N _e estimation helps to explore population demographic history, and allows one to predict the behavior of genetic diversity through time, N _e is a key parameter for the genetic management of small and isolated populations. Here, we explored an N _e-based approach using a bighorn sheep population on Tiburon Island, Mexico (TI) as a model. We estimated the current (N _crnt) and ancestral stable (N _stbl) inbreeding effective population sizes as well as summary statistics to assess genetic diversity and the demographic scenarios that could explain such diversity. Then, we evaluated the feasibility of using TI as a source population for reintroduction programs. We also included data from other bighorn sheep and artiodactyl populations in the analysis to compare their inbreeding effective size estimates. The TI population showed high levels of genetic diversity with respect to other managed populations. However, our analysis suggested that TI has been under a genetic bottleneck, indicating that using individuals from this population as the only source for reintroduction could lead to a severe genetic diversity reduction. Analyses of the published data did not show a strict correlation between H _E and N _crnt estimates. Moreover, we detected that ancient anthropogenic and climatic pressures affected all studied populations. We conclude that the estimation of N _crnt and N _stbl are informative genetic diversity estimators and should be used in addition to summary statistics for conservation and population management planning.     

Inbreeding and extinction: Effects of rate of inbreeding

Deleterious alleles may be removed (purged) bynatural selection in populations undergoinginbreeding. However, there is controversyregarding the effectiveness of selection inreducing the risk of extinction due toinbreeding, especially in relation to the rateof inbreeding. We evaluated the effect of therate of inbreeding on reducing extinction risk,in populations of Drosophila melanogastermaintained using full-sib mating (160replicates), or at effective population sizes(N _e) of 10 (80) or 20 (80).Extinction rates in the populations maintainedusing full-sib mating occurred at lower levelsof inbreeding than in the larger populations,whereas the two larger populations did notdiffer significantly from each other.Inbreeding coefficients at 50% extinction were0.62, 0.79 and 0.77 for the full-sib (N _e = 2.6), N _e = 10 and N _e = 20 treatments, respectively. Populations of N _e = 20 that remained extant after 60 generations, showed inbreeding depression, with the mean fitness of these populations being only 45% of the outbredcontrols. There was considerable variationamong the 31 inbred populations in fitness, butnone of the N _e = 20 populations hadfitness that was higher than the outbredcontrol. We conclude that purging may slow therate of extinction slightly, but it cannot berelied on to eliminate the deleterious effectsof inbreeding.     

INCREASED PROBABILITY OF EXTINCTION DUE TO DECREASED GENETIC EFFECTIVE POPULATION SIZE: EXPERIMENTAL POPULATIONS OF CLARKIA PULCHELLA

We established replicated experimental populations of the annual plant Clarkia pulchella to evaluate the existence of a causal relationship between loss of genetic variation and population survival probability. Two treatments differing in the relatedness of the founders, and thus in the genetic effective population size (N_e), were maintained as isolated populations in a natural environment. After three generations, the low N_e treatment had significantly lower germination and survival rates than did the high N_e treatment. These lower germination and survival rates led to decreased mean fitness in the low N_e populations: estimated mean fitness in the low N_e populations was only 21% of the estimated mean fitness in the high N_e populations. This inbreeding depression led to a reduction in population survival: at the conclusion of the experiment, 75% of the high N_e populations were still extant, whereas only 31% of the low N_e populations had survived. Decreased genetic effective population size, which leads to both inbreeding and the loss of alleles by genetic drift, increased the probability of population extinction over that expected from demographic and environmental stochasticity alone. This demonstrates that the genetic effective population size can strongly affect the probability of population persistence.     

The Effects of Demography and Long-Term Selection on the Accuracy of Genomic Prediction with Sequence Data

The use of dense SNPs to predict the genetic value of an individual for a complex trait is often referred to as “genomic selection” in livestock and crops, but is also relevant to human genetics to predict, for example, complex genetic disease risk. The accuracy of prediction depends on the strength of linkage disequilibrium (LD) between SNPs and causal mutations. If sequence data were used instead of dense SNPs, accuracy should increase because causal mutations are present, but demographic history and long-term negative selection also influence accuracy. We therefore evaluated genomic prediction, using simulated sequence in two contrasting populations: one reducing from an ancestrally large effective population size (N_e) to a small one, with high LD common in domestic livestock, while the second had a large constant-sized N_e with low LD similar to that in some human or outbred plant populations. There were two scenarios in each population; causal variants were either neutral or under long-term negative selection. For large N_e, sequence data led to a 22% increase in accuracy relative to ∼600K SNP chip data with a Bayesian analysis and a more modest advantage with a BLUP analysis. This advantage increased when causal variants were influenced by negative selection, and accuracy persisted when 10 generations separated reference and validation populations. However, in the reducing N_e population, there was little advantage for sequence even with negative selection. This study demonstrates the joint influence of demography and selection on accuracy of prediction and improves our understanding of how best to exploit sequence for genomic prediction.     

Effects of population characteristics and structure on estimates of effective population size in a house sparrow metapopulation

Effective population size (N_e) is a key parameter to understand evolutionary processes and the viability of endangered populations as it determines the rate of genetic drift and inbreeding. Low N_e can lead to inbreeding depression and reduced population adaptability. In this study, we estimated contemporary N_e using genetic estimators (LDNE, ONeSAMP, MLNE and CoNe) as well as a demographic estimator in a natural insular house sparrow metapopulation. We investigated whether population characteristics (population size, sex ratio, immigration rate, variance in population size and population growth rate) explained variation within and among populations in the ratio of effective to census population size (N_e/N_c). In general, N_e/N_c ratios increased with immigration rates. Genetic N_e was much larger than demographic N_e, probably due to a greater effect of immigration on genetic than demographic processes in local populations. Moreover, although estimates of genetic N_e seemed to track N_c quite well, the genetic N_e‐estimates were often larger than N_c within populations. Estimates of genetic N_e for the metapopulation were however within the expected range (<N_c). Our results suggest that in fragmented populations, even low levels of gene flow may have important consequences for the interpretation of genetic estimates of N_e. Consequently, further studies are needed to understand how N_e estimated in local populations or the total metapopulation relates to actual rates of genetic drift and inbreeding.     

The inbreeding effective number and the effective number of alleles in a population that varies in size

Consider a population that does not change in size. If it is assumed that there are an infinite number of possible neutral alleles at a locus and u is the probability that a particular gene mutates to some other gene in one generation, the effective number of alleles n_e is computed to be 4N_eu + 1, where N_e is the inbreeding effective population number. It is assumed in this paper that the number of individuals in a monoecious population, or the numbers of males and females in a dioecious population, are states in a finite irreducible Markov chain. In general it is impossible to obtain a single value of n_e. In some cases where the computation of n_e is possible, the results are as follows. When the population is monoecious, N_e is the reciprocal of the asymptotic average, over population sizes, of the probabilities that two gametes uniting to form an individual came from the same individual one generation earlier. In dioecious populations, N_e is the reciprocal of the long-run average of the probabilities that two homologous genes in separate individuals of one generation came from the same individual one generation earlier. Special cases are discussed.     

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F_ROH) in three local dairy cattle breeds. F_ROH values were compared with F estimated from the genomic relationship matrix (F_GRM), based on the difference between observed v. expected number of homozygous genotypes (F_HOM) and the genomic homozygosity of individual i (F_{MOL i}). The molecular coancestry coefficient (f_{MOL ij}) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F_ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F_HOM and F_ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F_ROH and other estimated F coefficients were generally lower ranged from 0.45 (F_{MOL i}−F_ROH) in Cinisara to 0.17 (F_GRM−F_ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.     

Estimation of effective population size in continuously distributed populations: there goes the neighborhood

Use of genetic methods to estimate effective population size (N_e) is rapidly increasing, but all approaches make simplifying assumptions unlikely to be met in real populations. In particular, all assume a single, unstructured population, and none has been evaluated for use with continuously distributed species. We simulated continuous populations with local mating structure, as envisioned by Wright''s concept of neighborhood size (NS), and evaluated performance of a single-sample estimator based on linkage disequilibrium (LD), which provides an estimate of the effective number of parents that produced the sample (N_b). Results illustrate the interacting effects of two phenomena, drift and mixture, that contribute to LD. Samples from areas equal to or smaller than a breeding window produced estimates close to the NS. As the sampling window increased in size to encompass multiple genetic neighborhoods, mixture LD from a two-locus Wahlund effect overwhelmed the reduction in drift LD from incorporating offspring from more parents. As a consequence, never approached the global N_e, even when the geographic scale of sampling was large. Results indicate that caution is needed in applying standard methods for estimating effective size to continuously distributed populations.     

Simple method for calculating confidence limits around inbreeding rate and effective population size estimates from pedigrees

An important concept in population genetics is effective population size (N_e), which describes the expected rate of loss of genetic variability from a population. One way to estimate N_e is using a pedigree. However, there are no methods for comparing the N_e estimated from a pedigree with that expected from life-history models. In the paper we show how N_e can be estimated from the change in inbreeding rate (f) estimated from a pedigree. The mean individual inbreeding rate in a population at a given time must be calculated from averaged values for males and females, where each age class is weighted by its reproductive value. We show an exact method for placing confidence intervals around f and N_e using a binomial distribution, and present a method for approximating this interval for large N_es using a Poisson distribution. These confidence intervals can be used to compare f and N_e from a pedigree to expected values from demographic models, and to compare N_es of two populations.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

Estimates of linkage disequilibrium and effective population sizes in Chinese Merino (Xinjiang type) sheep by genome-wide SNPs

Knowledge of linkage disequilibrium (LD) is important for effective genome-wide association studies and accurate genomic prediction. Chinese Merino (Xinjiang type) is well-known fine wool sheep breed. However, the extent of LD across the genome remains unexplored. In this study, we calculated autosomal LD based on genome-wide SNPs of 635 Chinese Merino (Xinjiang type) sheep by Illumina Ovine SNP50 BeadChip. A moderate level of LD (r ²?≥?0.25) across the whole genome was observed at short distances of 0–10 kb. Further, the ancestral effective population size (N _e ) was analyzed by extent of LD and found that N _e increased with the increase of generations and declined rapidly within the most recent 50 generations, which is consistent with the history of Chinese Merino sheep breeding, initiated in 1971. We also noted that even when the effective population size was estimated across different single chromosomes, N _e only ranged from 140.36 to 183.33 at five generations in the past, exhibiting a rapid decrease compared with that at ten generations in the past. These results indicated that the genetic diversity in Chinese Merino sheep recently decreased and proper protective measures should be taken to maintain the diversity. Our datasets provided essential genetic information to track molecular variations which potentially contribute to phenotypic variation in Chinese Merino sheep.     

Inbreeding depression and drift load in small populations at demographic disequilibrium

Inbreeding depression is a major driver of mating system evolution and has critical implications for population viability. Theoretical and empirical attention has been paid to predicting how inbreeding depression varies with population size. Lower inbreeding depression is predicted in small populations at equilibrium, primarily due to higher inbreeding rates facilitating purging and/or fixation of deleterious alleles (drift load), but predictions at demographic and genetic disequilibrium are less clear. In this study, we experimentally evaluate how lifetime inbreeding depression and drift load, estimated by heterosis, vary with census (N_c) and effective (estimated as genetic diversity, H_e) population size across six populations of the biennial Sabatia angularis as well as present novel models of inbreeding depression and heterosis under varying demographic scenarios at disequilibrium (fragmentation, bottlenecks, disturbances). Our experimental study reveals high average inbreeding depression and heterosis across populations. Across our small sample, heterosis declined with H_e, as predicted, whereas inbreeding depression did not vary with H_e and actually decreased with N_c. Our theoretical results demonstrate that inbreeding depression and heterosis levels can vary widely across populations at disequilibrium despite similar H_e and highlight that joint demographic and genetic dynamics are key to predicting patterns of genetic load in nonequilibrium systems.