Dataset Information

Estimating contemporary effective population size in non-model species using linkage disequilibrium across thousands of loci.

ABSTRACT: Contemporary effective population size (Ne) can be estimated using linkage disequilibrium (LD) observed across pairs of loci presumed to be selectively neutral and unlinked. This method has been commonly applied to data sets containing 10-100 loci to inform conservation and study population demography. Performance of these Ne estimates could be improved by incorporating data from thousands of loci. However, these thousands of loci exist on a limited number of chromosomes, ensuring that some fraction will be physically linked. Linked loci have elevated LD due to limited recombination, which if not accounted for can cause Ne estimates to be downwardly biased. Here, we present results from coalescent and forward simulations designed to evaluate the bias of LD-based Ne estimates ([Ncirc ]e). Contrary to common perceptions, increasing the number of loci does not increase the magnitude of linkage. Although we show it is possible to identify some pairs of loci that produce unusually large r(2) values, simply removing large r(2) values is not a reliable way to eliminate bias. Fortunately, the magnitude of bias in [Ncirc ]e is strongly and negatively correlated with the process of recombination, including the number of chromosomes and their length, and this relationship provides a general way to adjust for bias. Additionally, we show that with thousands of loci, precision of [Ncirc ]e is much lower than expected based on the assumption that each pair of loci provides completely independent information.

SUBMITTER: Waples RK

PROVIDER: S-EPMC5026758 | biostudies-literature | 2016 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Estimating contemporary effective population size in non-model species using linkage disequilibrium across thousands of loci.

Waples R K RK Larson W A WA Waples R S RS

Heredity 20160824 4

Contemporary effective population size (Ne) can be estimated using linkage disequilibrium (LD) observed across pairs of loci presumed to be selectively neutral and unlinked. This method has been commonly applied to data sets containing 10-100 loci to inform conservation and study population demography. Performance of these Ne estimates could be improved by incorporating data from thousands of loci. However, these thousands of loci exist on a limited number of chromosomes, ensuring that some frac ...[more]

PMID: 27553452

Similar Datasets

Project description: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 (Ne) and inbreeding coefficients as well as to identify potential regions of positive selection. The decay of linkage disequilibrium was slow with r2=0.95 in 50 kb distance. The last 7-200 generations ago, Ne was at 10-13. An increase of Ne was noted in the very recent generations with a peak value of 19 for Ne 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 FROH 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 (FPed=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. The objectives of the present study were to genotype 28 Lundehund using the canine Illumina high density beadchip (Illumina, San Diego, CA, USA) to estimate the effective population size (Ne) from data on linkage disequilibria (LD). We identified runs of homozygosity (ROH) as regions with a local loss of genetic variation. Inbreeding coefficients were calculated on the basis of pedigree data (FPed) and genotype information (FROH, FIS) to compare the results of these different methods for this population. We applied the long-range haplotype test to search for potential selective sweeps.

Project description:BackgroundThe use of molecular genetic technologies for broodstock management and selective breeding of aquaculture species is becoming increasingly more common with the continued development of genome tools and reagents. Several laboratories have produced genetic maps for rainbow trout to aid in the identification of loci affecting phenotypes of interest. These maps have resulted in the identification of many quantitative/qualitative trait loci affecting phenotypic variation in traits associated with albinism, disease resistance, temperature tolerance, sex determination, embryonic development rate, spawning date, condition factor and growth. Unfortunately, the elucidation of the precise allelic variation and/or genes underlying phenotypic diversity has yet to be achieved in this species having low marker densities and lacking a whole genome reference sequence. Experimental designs which integrate segregation analyses with linkage disequilibrium (LD) approaches facilitate the discovery of genes affecting important traits. To date the extent of LD has been characterized for humans and several agriculturally important livestock species but not for rainbow trout.ResultsWe observed that the level of LD between syntenic loci decayed rapidly at distances greater than 2 cM which is similar to observations of LD in other agriculturally important species including cattle, sheep, pigs and chickens. However, in some cases significant LD was also observed up to 50 cM. Our estimate of effective population size based on genome wide estimates of LD for the NCCCWA broodstock population was 145, indicating that this population will respond well to high selection intensity. However, the range of effective population size based on individual chromosomes was 75.51 - 203.35, possibly indicating that suites of genes on each chromosome are disproportionately under selection pressures.ConclusionsOur results indicate that large numbers of markers, more than are currently available for this species, will be required to enable the use of genome-wide integrated mapping approaches aimed at identifying genes of interest in rainbow trout.

Dataset Information

Estimating contemporary effective population size in non-model species using linkage disequilibrium across thousands of loci.

Publications

Estimating contemporary effective population size in non-model species using linkage disequilibrium across thousands of loci.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets