Estimates Of Effective Population Size Research Articles

Effect of sampling methods, effective population size and migration rate estimation in Glossina palpalis palpalis from Cameroon

Human and animal trypanosomiases are two major constraints to development in Africa. These diseases are mainly transmitted by tsetse flies in particular by Glossina palpalis palpalis in Western and Central Africa. To set up an effective vector control campaign, prior population genetics studies have proved useful. Previous studies on population genetics of G. p. palpalis using microsatellite loci showed high heterozygote deficits, as compared to Hardy–Weinberg expectations, mainly explained by the presence of null alleles and/or the mixing of individuals belonging to several reproductive units (Wahlund effect). In this study we implemented a system of trapping, consisting of a central trap and two to four satellite traps around the central one to evaluate a possible role of the Wahlund effect in tsetse flies from three Cameroon human and animal African trypanosomiases foci (Campo, Bipindi and Fontem). We also estimated effective population sizes and dispersal. No difference was observed between the values of allelic richness, genetic diversity and Wright’s FIS, in the samples from central and from satellite traps, suggesting an absence of Wahlund effect. Partitioning of the samples with Bayesian methods showed numerous clusters of 2–3 individuals as expected from a population at demographic equilibrium with two expected offspring per reproducing female. As previously shown, null alleles appeared as the most probable factor inducing these heterozygote deficits in these populations. Effective population sizes varied from 80 to 450 individuals while immigration rates were between 0.05 and 0.43, showing substantial genetic exchanges between different villages within a focus. These results suggest that the “suppression” with establishment of physical barriers may be the best strategy for a vector control campaign in this forest context.

Similarity thresholds used in DNA sequence assembly from short reads can reduce the comparability of population histories across species.

Comparing inferences among datasets generated using short read sequencing may provide insight into the concerted impacts of divergence, gene flow and selection across organisms, but comparisons are complicated by biases introduced during dataset assembly. Sequence similarity thresholds allow the de novo assembly of short reads into clusters of alleles representing different loci, but the resulting datasets are sensitive to both the similarity threshold used and to the variation naturally present in the organism under study. Thresholds that require high sequence similarity among reads for assembly (stringent thresholds) as well as highly variable species may result in datasets in which divergent alleles are lost or divided into separate loci (‘over-splitting’), whereas liberal thresholds increase the risk of paralogous loci being combined into a single locus (‘under-splitting’). Comparisons among datasets or species are therefore potentially biased if different similarity thresholds are applied or if the species differ in levels of within-lineage genetic variation. We examine the impact of a range of similarity thresholds on assembly of empirical short read datasets from populations of four different non-model bird lineages (species or species pairs) with different levels of genetic divergence. We find that, in all species, stringent similarity thresholds result in fewer alleles per locus than more liberal thresholds, which appears to be the result of high levels of over-splitting. The frequency of putative under-splitting, conversely, is low at all thresholds. Inferred genetic distances between individuals, gene tree depths, and estimates of the ancestral mutation-scaled effective population size (θ) differ depending upon the similarity threshold applied. Relative differences in inferences across species differ even when the same threshold is applied, but may be dramatically different when datasets assembled under different thresholds are compared. These differences not only complicate comparisons across species, but also preclude the application of standard mutation rates for parameter calibration. We suggest some best practices for assembling short read data to maximize comparability, such as using more liberal thresholds and examining the impact of different thresholds on each dataset.

Genome-wide estimates of coancestry, inbreeding and effective population size in the Spanish Holstein population.

Estimates of effective population size in the Holstein cattle breed have usually been low despite the large number of animals that constitute this breed. Effective population size is inversely related to the rates at which coancestry and inbreeding increase and these rates have been high as a consequence of intense and accurate selection. Traditionally, coancestry and inbreeding coefficients have been calculated from pedigree data. However, the development of genome-wide single nucleotide polymorphisms has increased the interest of calculating these coefficients from molecular data in order to improve their accuracy. In this study, genomic estimates of coancestry, inbreeding and effective population size were obtained in the Spanish Holstein population and then compared with pedigree-based estimates. A total of 11,135 animals genotyped with the Illumina BovineSNP50 BeadChip were available for the study. After applying filtering criteria, the final genomic dataset included 36,693 autosomal SNPs and 10,569 animals. Pedigree data from those genotyped animals included 31,203 animals. These individuals represented only the last five generations in order to homogenise the amount of pedigree information across animals. Genomic estimates of coancestry and inbreeding were obtained from identity by descent segments (coancestry) or runs of homozygosity (inbreeding). The results indicate that the percentage of variance of pedigree-based coancestry estimates explained by genomic coancestry estimates was higher than that for inbreeding. Estimates of effective population size obtained from genome-wide and pedigree information were consistent and ranged from about 66 to 79. These low values emphasize the need of controlling the rate of increase of coancestry and inbreeding in Holstein selection programmes.

Testing reticulate versus coalescent origins of Erica lusitanica using a species phylogeny of the northern heathers (Ericeae, Ericaceae)

Whilst most of the immense species richness of heathers (Calluna, Daboecia and Erica: Ericeae; Ericaceae) is endemic to Africa, particularly the Cape Floristic Region, the oldest lineages are found in the Northern Hemisphere. We present phylogenetic hypotheses for the major clades of Ericeae represented by multiple accessions of all northern Erica species and placeholder taxa for the large nested African/Madagascan clade. We identified consistent, strongly supported conflict between gene trees inferred from ITS and chloroplast DNA sequences with regard to the position of Erica lusitanica. We used coalescent simulations to test whether this conflict could be explained by coalescent stochasticity, as opposed to reticulation (e.g. hybridisation), given estimates of clade ages, generation time and effective population sizes (Ne). A standard approach, comparing overall differences between real and simulated trees, could not clearly reject coalescence. However, additional simulations showed that at the (higher) Ne necessary to explain conflict in E. lusitanica, further topological conflict would also be expected. Ancient hybridisation between ancestors of northern species is therefore a plausible scenario to explain the origin of E. lusitanica, and its morphological similarities to E. arborea. Assuming either process influences the results of species tree and further evolutionary inference. The coalescence scenario is equivocal with regard the standing hypothesis of stepping stone dispersal of Erica from Europe into Africa; whereas reticulate evolution in E. lusitanica would imply that the colonisation of Tropical East Africa by E. arborea instead occurred independently of dispersals within the rest of the African/Madagascan clade.

Defining population structure and genetic signatures of decline in the giant gartersnake (Thamnophis gigas): implications for conserving threatened species within highly altered landscapes

Anthropogenic habitat fragmentation can disrupt the ability of species to disperse across landscapes, which can alter the levels and distribution of genetic diversity within populations and negatively impact long-term viability. The giant gartersnake (Thamnophis gigas) is a state and federally threatened species that historically occurred in the wetland habitats of California’s Great Central Valley. Despite the loss of 93 % of historic wetlands throughout the Central Valley, giant gartersnakes continue to persist in relatively small, isolated patches of highly modified agricultural wetlands. Gathering information regarding genetic diversity and effective population size represents an essential component for conservation management programs aimed at this species. Previous mitochondrial sequence studies have revealed historical patterns of differentiation, yet little is known about contemporary population structure and diversity. On the basis of 15 microsatellite loci, we estimate population structure and compare indices of genetic diversity among populations spanning seven drainage basins within the Central Valley. We sought to understand how habitat loss may have affected genetic differentiation, genetic diversity and effective population size, and what these patterns suggest in terms of management and restoration actions. We recovered five genetic clusters that were consistent with regional drainage basins, although three northern basins within the Sacramento Valley formed a single genetic cluster. Our results show that northern drainage basin populations have higher connectivity than among central and southern basins populations, and that greater differentiation exists among the more geographically isolated populations in the central and southern portion of the species’ range. Genetic diversity measures among basins were significantly different, and were generally lower in southern basin populations. Levels of inbreeding and evidence of population bottlenecks were detected in about half the populations we sampled, and effective population size estimates were well below recommended minimum thresholds to avoid inbreeding. Efforts focused on maintaining and enhancing existing wetlands to facilitate dispersal between basins and increase local effective population sizes may be critical for these otherwise isolated populations.

Extremely low microsatellite diversity but distinct population structure in a long-lived threatened species, the Australian lungfish Neoceratodus forsteri (Dipnoi).

The Australian lungfish is a unique living representative of an ancient dipnoan lineage, listed as ‘vulnerable’ to extinction under Australia’s Environment Protection and Biodiversity Conservation Act 1999. Historical accounts indicate this species occurred naturally in two adjacent river systems in Australia, the Burnett and Mary. Current day populations in other rivers are thought to have arisen by translocation from these source populations. Early genetic work detected very little variation and so had limited power to answer questions relevant for management including how genetic variation is partitioned within and among sub-populations. In this study, we use newly developed microsatellite markers to examine samples from the Burnett and Mary Rivers, as well as from two populations thought to be of translocated origin, Brisbane and North Pine. We test whether there is significant genetic structure among and within river drainages; assign putatively translocated populations to potential source populations; and estimate effective population sizes. Eleven polymorphic microsatellite loci genotyped in 218 individuals gave an average within-population heterozygosity of 0.39 which is low relative to other threatened taxa and for freshwater fishes in general. Based on F ST values (average over loci = 0.11) and STRUCTURE analyses, we identify three distinct populations in the natural range, one in the Burnett and two distinct populations in the Mary. These analyses also support the hypothesis that the Mary River is the likely source of translocated populations in the Brisbane and North Pine rivers, which agrees with historical published records of a translocation event giving rise to these populations. We were unable to obtain bounded estimates of effective population size, as we have too few genotype combinations, although point estimates were low, ranging from 29 - 129. We recommend that, in order to preserve any local adaptation in the three distinct populations that they be managed separately.

Population differentiation determined from putative neutral and divergent adaptive genetic markers in Eulachon (Thaleichthys pacificus, Osmeridae), an anadromous Pacific smelt.

Twelve eulachon (Thaleichthys pacificus, Osmeridae) populations ranging from Cook Inlet, Alaska and along the west coast of North America to the Columbia River were examined by restriction-site-associated DNA (RAD) sequencing to elucidate patterns of neutral and adaptive variation in this high geneflow species. A total of 4104 single-nucleotide polymorphisms (SNPs) were discovered across the genome, with 193 putatively adaptive SNPs as determined by F(ST) outlier tests. Estimates of population structure in eulachon with the putatively adaptive SNPs were similar, but provided greater resolution of stocks compared with a putatively neutral panel of 3911 SNPs or previous estimates with 14 microsatellites. A cline of increasing measures of genetic diversity from south to north was found in the adaptive panel, but not in the neutral markers (SNPs or microsatellites). This may indicate divergent selective pressures in differing freshwater and marine environments between regional eulachon populations and that these adaptive diversity patterns not seen with neutral markers could be a consideration when determining genetic boundaries for conservation purposes. Estimates of effective population size (N(e)) were similar with the neutral SNP panel and microsatellites and may be utilized to monitor population status for eulachon where census sizes are difficult to obtain. Greater differentiation with the panel of putatively adaptive SNPs provided higher individual assignment accuracy compared to the neutral panel or microsatellites for stock identification purposes. This study presents the first SNPs that have been developed for eulachon, and analyses with these markers highlighted the importance of integrating genome-wide neutral and adaptive genetic variation for the applications of conservation and management.

Genetic Diversity of White Sharks, Carcharodon carcharias, in the Northwest Atlantic and Southern Africa.

The white shark, Carcharodon carcharias, is both one of the largest apex predators in the world and among the most heavily protected marine fish. Population genetic diversity is in part shaped by recent demographic history and can thus provide information complementary to more traditional population assessments, which are difficult to obtain for white sharks and have at times been controversial. Here, we use the mitochondrial control region and 14 nuclear-encoded microsatellite loci to assess white shark genetic diversity in 2 regions: the Northwest Atlantic (NWA, N = 35) and southern Africa (SA, N = 131). We find that these 2 regions harbor genetically distinct white shark populations (Φ ST = 0.10, P < 0.00001; microsatellite F ST = 0.1057, P < 0.021). M-ratios were low and indicative of a genetic bottleneck in the NWA (M-ratio = 0.71, P < 0.004) but not SA (M-ratio = 0.85, P = 0.39). This is consistent with other evidence showing a steep population decline occurring in the mid to late 20th century in the NWA, whereas the SA population appears to have been relatively stable. Estimates of effective population size ranged from 22.6 to 66.3 (NWA) and 188 to 1998.3 (SA) and evidence of inbreeding was found (primarily in NWA). Overall, our findings indicate that white population dynamics within NWA and SA are determined more by intrinsic reproduction than immigration and there is genetic evidence of a population decline in the NWA, further justifying the strong domestic protective measures that have been taken for this species in this region. Our study also highlights how assessment of genetic diversity can complement other sources of information to better understand the status of threatened marine fish populations.

Estimation of the variance effective population size in age structured populations

The variance effective population size for age structured populations is generally hard to estimate and the temporal method often gives biased estimates. Here, we give an explicit expression for a correction factor which, combined with estimates from the temporal method, yield approximately unbiased estimates. The calculation of the correction factor requires knowledge of the age specific offspring distribution and survival probabilities as well as possible correlation between survival and reproductive success. In order to relax these requirements, we show that only first order moments of these distributions need to be known if the time between samples is large, or individuals from all age classes which reproduce are sampled. A very explicit approximate expression for the asymptotic coefficient of standard deviation of the estimator is derived, and it can be used to construct confidence intervals and optimal ways of weighting information from different markers. The asymptotic coefficient of standard deviation can also be used to design studies and we show that in order to maximize the precision for a given sample size, individuals from older age classes should be sampled since their expected variance of allele frequency change is higher and easier to estimate. However, for populations with fluctuating age class sizes, the accuracy of the method is reduced when samples are taken from older age classes with high demographic variation. We also present a method for simultaneous estimation of the variance effective and census population size.

Using parentage analysis to estimate rates of straying and homing in Chinook salmon (Oncorhynchus tshawytscha).

We used parentage analysis based on microsatellite genotypes to measure rates of homing and straying of Chinook salmon (Oncorhynchus tshawytscha) among five major spawning tributaries within the Wenatchee River, Washington. On the basis of analysis of 2248 natural-origin and 11594 hatchery-origin fish, we estimated that the rate of homing to natal tributaries by natural-origin fish ranged from 0% to 99% depending on the tributary. Hatchery-origin fish released in one of the five tributaries homed to that tributary at a far lower rate than the natural-origin fish (71% compared to 96%). For hatchery-released fish, stray rates based on parentage analysis were consistent with rates estimated using physical tag recoveries. Stray rates among major spawning tributaries were generally higher than stray rates of tagged fish to areas outside of the Wenatchee River watershed. Within the Wenatchee watershed, rates of straying by natural-origin fish were significantly affected by spawning tributary and by parental origin: progeny of naturally spawning hatchery-produced fish strayed at significantly higher rates than progeny whose parents were themselves of natural origin. Notably, none of the 170 offspring that were products of mating by two natural-origin fish strayed from their natal tributary. Indirect estimates of gene flow based on FST statistics were correlated with but higher than the estimates from the parentage data. Tributary-specific estimates of effective population size were also correlated with the number of spawners in each tributary.

A test of the central-marginal hypothesis using population genetics and ecological niche modelling in an endemic salamander (Ambystoma barbouri).

The central-marginal hypothesis (CMH) predicts that population size, genetic diversity and genetic connectivity are highest at the core and decrease near the edges of species' geographic distributions. We provide a test of the CMH using three replicated core-to-edge transects that encompass nearly the entire geographic range of the endemic streamside salamander (Ambystoma barbouri). We confirmed that the mapped core of the distribution was the most suitable habitat using ecological niche modelling (ENM) and via genetic estimates of effective population sizes. As predicted by the CMH, we found statistical support for decreased genetic diversity, effective population size and genetic connectivity from core to edge in western and northern transects, yet not along a southern transect. Based on our niche model, habitat suitability is lower towards the southern range edge, presumably leading to conflicting core-to-edge genetic patterns. These results suggest that multiple processes may influence a species' distribution based on the heterogeneity of habitat across a species' range and that replicated sampling may be needed to accurately test the CMH. Our work also emphasizes the importance of identifying the geographic range core with methods other than using the Euclidean centre on a map, which may help to explain discrepancies among other empirical tests of the CMH. Assessing core-to-edge population genetic patterns across an entire species' range accompanied with ENM can inform our general understanding of the mechanisms leading to species' geographic range limits.

Assessment of recent bottlenecks and estimation of effective population size in the Ethiopian wild sorghum using simple sequence repeat allele diversity and mutation models

Since the immediate wild relatives ofSorghum bicolor(L.) Moench are indigenous to Ethiopia, studying their population biology is timely for undertaking conservation measures. A study was conducted to investigate the occurrence of population bottlenecks and to estimate the long-term effective population size (Ne) in wild relatives of sorghum. For this, 40 samples of wild sorghum were collected from two remotely located populations that were allopatric to the cultivated sorghum. The presence of bottlenecks was investigated using heterozygosity excess/deficiency, mode shift and allelic diversity based on nine polymorphic simple sequence repeat (SSR) loci. We also estimated theNeof the studied populations using two different methods employing SSR mutation models. The expected heterozygosity was found to be 0.41 and 0.71 and allelic richness was 3.0 and 4.9, in Awash and Gibe populations, respectively. Neither the heterozygosity excess nor the mode-shift methods detected signatures of bottlenecks in the studied populations. The effective size of the two wild sorghum populations studied also showed no risk of population reduction in these regions of Ethiopia. Therefore, these allopatric wild sorghum populations can survive by occupying patches by the roadsides and fences, areas within abandoned farm lands, forests, etc., which shows that their wild characteristics of adaptation have been adequate for them to survive from extinction despite extensive deforestation of their habitat for modern agriculture and frequent grazing by livestock. However, this does not guarantee the survival of these species for the future andex situconservation measures or policies could help maintain their diversity.

Neon: An R Package to Estimate Human Effective Population Size and Divergence Time from Patterns of Linkage Disequilibrium between SNPS

Objective: Estimating the effective population size (Ne) is crucial to understanding how populations evolved, expanded or shrunk. One possible approach is to compare DNA diversity, so as to obtain an average Ne over many past generations; however as the population sizes change over time, another possibility is to describe this change. Linkage Disequilibrium (LD) patterns contain information about these changes, and, whenever a large number of densely linked markers are available, can be used to monitor fluctuating population size through time. Here, we present a new R package, NeON that has been designed to explore population’s LD patterns to reconstruct two key parameters of human evolution: the effective population size and the divergence time between populations. Methods: NeON starts with binary or pairwise-LD PLINK files, and allows (a) to assign a genetic map position using HapMap (NCBI release 36 or 37) (b) to calculate the effective population size over time exploiting the relationship between Ne and the average squared correlation coefficient of LD (r2LD) within predefined recombination distance categories, and (c) to calculate the confidence interval about Ne based on the observed variation of the estimator across chromosomes; the outputs of the functions are both numerical and graphical. This package also offers the possibility to estimate the divergence time between populations given the Ne values calculated from the within-population LD data and a matrix of between-populations FST. These routines can be adapted to any species whenever genetic map positions are available. Results and Conclusion: The functions contained in the R package NeON provide reliable estimates of effective population sizes of human chromosomes from LD patterns of genome-wide SNPs data, as it is shown here for the populations contained in the CEPH panel. The NeON package enables to accommodate variable numbers of individuals, populations and genetic markers, allowing analyzing those using standard personal computers.

Genetic structure, diversity, and interisland dispersal in the endangered Mariana Common Moorhen (Gallinula chloropus guami)

ABSTRACT The Mariana Common Moorhen (Gallinula chloropus guami) is a highly endangered taxon, with fewer than 300 individuals estimated to occur in the wild. The subspecies is believed to have undergone population declines attributable to loss of wetland habitats on its native islands in the Mariana Islands. We analyzed mitochondrial DNA (mtDNA) sequences (control region and ND2 genes) and nuclear microsatellite loci in Mariana Common Moorhens from Guam and Saipan, the two most distal islands inhabited by the subspecies. Our analyses revealed similar nuclear genetic diversity and effective population size estimates on Saipan and Guam. Birds from Guam and Saipan were genetically differentiated (microsatellites: FST = 0.152; control region: FST = 0.736; ND2: FST = 0.390); however, assignment tests revealed the presence of first-generation dispersers from Guam onto Saipan (1 of 27 sampled birds) and from Saipan onto Guam (2 of 28 sampled birds), suggesting the capability for long-distance interpopulation mov...

Conservation Genetics of Six Species of GenusDionda(Cyprinidae) in the Southwestern United States

We examined allelic variation at nuclear-encoded microsatellites and sequences of mitochondrial (mt)DNA in 10 geographic samples representing 6 nominal species of the cyprinid genus Dionda. Species of Dionda are found in springs and spring-fed headwaters in the southwestern United States and Mexico and are of particular interest to conservation and management, in part because of their limited distribution and habitat specificity, and in part as indicator species of habitat quality. All 10 samples examined appear to be discrete, demographically independent populations, with greater observed FST values between or among samples within species (0.123–0.280) than threshold values above which demographic independence is indicated. All 10 exhibited microsatellite and mtDNA variation comparable to or lower than that found in other cyprinids considered to be threatened or endangered; across microsatellites, average number of alleles across populations ranged from of 2.09 to 9.76, allelic richness from 2.24 to 8.45, and gene diversity from 0.0211 to 0.606; for mtDNA, the number of haplotypes across populations ranged from 1 to 14. Estimates of historical and present-day genetic demography indicated that all 10 populations have experienced order-of-magnitude declines in effective population size, with lower bounds of time intervals for the declines in 9 of the populations ranging from 6 to 65 years. Estimates of average long-term effective population size (536 in Dionda argentosa from San Felipe Creek to 2335 in D. texensis) and effective number of breeders (22 in D. flavipinnis from Fessenden Spring to 555 in D. diaboli from Devils River) also indicated recent declines in effective size, and at least 5 of the populations appear to have undergone recent, severe bottlenecks (mean Mc range 0.806–0.848, P value range 0.000–0.0350). The observation that all 10 populations are demographically independent indicates that local extirpations likely would not be replaced by new migrants and that loss of any of the populations would represent loss of a unique genetic entity. Conservation recommendations for each of the populations are briefly discussed.

DNA analysis of ancient dogs of the Americas: Identifying possible founding haplotypes and reconstructing population histories

As dogs have traveled with humans to every continent, they can potentially serve as an excellent proxy when studying human migration history. Past genetic studies into the origins of Native American dogs have used portions of the hypervariable region (HVR) of mitochondrial DNA (mtDNA) to indicate that prior to European contact the dogs of Native Americans originated in Eurasia. In this study, we summarize past DNA studies of both humans and dogs to discuss their population histories in the Americas. We then sequenced a portion of the mtDNA HVR of 42 pre-Columbian dogs from three sites located in Illinois, coastal British Columbia, and Colorado, and identify four novel dog mtDNA haplotypes. Next, we analyzed a dataset comprised of all available ancient dog sequences from the Americas to infer the pre-Columbian population history of dogs in the Americas. Interestingly, we found low levels of genetic diversity for some populations consistent with the possibility of deliberate breeding practices. Furthermore, we identified multiple putative founding haplotypes in addition to dog haplotypes that closely resemble those of wolves, suggesting admixture with North American wolves or perhaps a second domestication of canids in the Americas. Notably, initial effective population size estimates suggest at least 1000 female dogs likely existed in the Americas at the time of the first known canid burial, and that population size increased gradually over time before stabilizing roughly 1200 years before present.

High levels of genetic diversity through pollen flow of the coniferous Araucaria angustifolia: a landscape level study in Southern Brazil

In this study, genetic diversity, inbreeding, spatial genetic structure (SGS), and pollen dispersal are analyzed using ten microsatellite loci from two populations of the dioecious, wind-pollinated, coniferous tree Araucaria angustifolia in Southern Brazil. The study populations include an undisturbed 7 ha A. angustifolia cluster, where all adult trees were mapped and sampled, and an adjacent, long-abandoned, open agricultural area with an aggregated A. angustifolia population. Seeds were collected from 13 seed trees inside the forest cluster and from eight seed trees in the open forest. Our results showed that the adults present high levels of heterozygosity (H o = 0.91) and an absence of inbreeding. However, significant SGS was detected up to 90 m in the forest cluster suggesting that near-neighbor trees are related. The estimate of effective population size was lower than the total number of trees in the cluster (N e/N = 0.19), which can be explained by the presence of SGS in the stand. Substantial external pollen flow was detected in the forest cluster (26 %) and open forest (20 %), indicating that the reproductive population size is greater than the sampled populations, explaining the high genetic diversity in this population. Our results indicate that this site has potential for in and ex situ conservation due to high levels of genetic diversity and gene immigration resulting from pollen flow. Conservation strategies for A. angustifolia should focus not only on forest fragments but should also include the preservation of isolated trees throughout the landscape.

Estimations of genomic linkage disequilibrium and effective population sizes in three sheep populations

As one of the most important measures in population genetics, linkage disequilibrium (LD) is the basis of effective population size estimation, genomic selection, genome-wide association study and quantitative trait locus mapping. Here, we characterized the pattern of LD in three sheep populations (Sunite, German Mutton Merino, and Dorper) and subsequently estimated persistence of phase among them using the Illumina Ovine SNP50k Chip. In addition, we calculated the corresponding effective population sizes using the metric r2 estimated. The results showed that the extent of LD was highest in Dorpers, followed by German Mutton Merinos and Sunites. Average r2 values were 0.13, 0.20, and 0.22 at adjacent marker distances of 57.28, 59.01, and 59.93kb for Sunite, German Mutton Merino, and Dorper sheep, respectively. The average r2 values at the same genomic distance in German Mutton Merino and Dorper sheep were similar and higher than that in Sunite sheep, and the rate of LD decay over distance was faster than the rates in the other two sheep breeds. The correlation of phase for different breed comparisons ranged from 0.784 to 0.825. The effective population sizes in Sunite, German Mutton Merino, and Dorper sheep 7 generations ago were approximately 207, 74, and 67, respectively. We also detected a bottleneck in the effective population size for Sunite sheep that was consistent with the timing for the origin of the domestication of Chinese sheep. Our findings provide a theoretical reference for future genome marker-assisted evaluation and selection in sheep.

Estrutura populacional da raça Girolando

O objetivo neste estudo foi avaliar a estrutura genética da população de bovinos da raça Girolando no Brasil. Analisou-se o arquivo de pedigree de 26.969 animais, composto de 3.031 machos e 23.938 fêmeas. O nível de conteúdo de informação do pedigree na geração atual foi 61%, mostrando ser de qualidade moderada. O coeficiente de endogamia médio e o coeficiente de relação médio da população Girolando foram 0,11 e 0,13%, respectivamente. O tamanho efetivo da população, considerando a geração completa traçada, foi 188, acima do nível crítico. Do total de 9.457 ancestrais que contribuíram para a população de referência, 457 explicaram 50% da variabilidade genética da população. O número efetivo de fundadores foi 551 e o de ancestrais 393. O intervalo médio de geração foi de 5,26 anos, sendo ligeiramente maior nas trilhas gaméticas mãe-filho e pai-filha. A partir dos coeficientes estimados, pode-se concluir que a endogamia nos rebanhos da raça Girolando foi de pequena magnitude e que as práticas de acasalamento foram adequadas durante o período avaliado. No entanto, é importante continuar com o monitoramento desses coeficientes a fim de prevenir perda de variabilidade genética

Lack of detectable genetic differentiation between den populations of the Prairie Rattlesnake (Crotalusviridis) in a fragmented landscape

Numerous studies have reported genetic fragmentation of species whose habitat has been modified by roads and other anthropogenic features, but it is still not clear how most species respond to roads and whether genetic effects can be detected over a limited number of generations. We used road-crossing models and population genetic analysis (based on microsatellite loci) to make inferences about functional connectivity between populations of the Prairie Rattlesnake (Crotalus viridis (Rafinesque, 1818)) on opposite sides of the Trans-Canada Highway near Medicine Hat (Alberta, Canada). The road-crossing model predicted a high probability of mortality while crossing the Trans-Canada Highway. However, model-based genetic clustering methods (STRUCTURE and BAPS) did not detect structure; a nonmodel-based clustering method (DAPC) found structure, but most groups consisted of individuals captured throughout the study area. Estimates of effective population size were immeasurably large and power to detect genic differentiation was diminished if the effective size exceeded 500; this reduction in power was intensified when the number of loci was reduced (from eight to five to account for null alleles). Our results corroborate accounts of long-distance migration by this species and indicate that genetic fragmentation may not be easily detectable over this spatial and temporal scale.