Estimating genetic drift and effective population size from temporal shifts in dominant gene marker frequencies

Measurement of temporal change in allele frequencies represents an indirect method for estimating the genetically effective size of populations. When allele frequencies are estimated for gene markers that display dominant gene expression, such as, e.g. random amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP) markers, the estimates can be seriously biased. We quantify bias for previous allele frequency estimators and present a new expression that is generally less biased and provides a more precise assessment of temporal allele frequency change. We further develop an estimator for effective population size that is appropriate when dealing with dominant gene markers. Comparison with estimates based on codominantly expressed genes, such as allozymes or microsatellites, indicates that about twice as many loci or sampled individuals are required when using dominant markers to achieve the same precision.     

Development and Characterization of RAPD and Microsatellite Markers for Genetic Variation Analysis in the Critically Endangered Yellow Catfish <Emphasis Type="Italic">Horabagrus nigricollaris</Emphasis> (Teleostei: Horabagridae)

Random-amplified polymorphic DNA (RAPD) and microsatellite markers were developed and used for the analysis of genetic variability in the critically endangered yellow catfish Horabagrus nigricollaris, sampled from the Chalakkudy River, Kerala, India. Eight RAPD and five microsatellite markers were detected to genotype the species. In RAPD, the 73 fragments were 20.55% polymorphic, whereas 4 polymorphic loci (80%) were obtained in microsatellites. In microsatellites, the number of alleles across the 5 loci was 1-5, and the range of heterozygosity was 0.25-0.5. The mean observed number of alleles was 2.4, and the effective number was 1.775 per locus. The average heterozygosity across all investigated samples was 0.29, indicating a significant deficiency of heterozygotes in this species. RAPD and microsatellite methods report a low degree of gene diversity and lack of genetic heterogeneity in the population of H. nigricollaris, emphasizing the need for fishery management, conservation, and rehabilitation of this species.     

Distinguishing adaptive from nonadaptive genetic differentiation: comparison of Q(ST) and F(ST) at two spatial scales

Genetic differentiation in 20 hierarchically sampled populations of wild barley was analyzed with quantitative traits, allozymes and Random Amplified Polymorphic DNAs (RAPDs), and compared for three marker types at two hierarchical levels. Regional subdivision for both molecular markers was much lower than for quantitative traits. For both allozymes and RAPDs, most loci exhibited minor or no regional differentiation, and the relatively high overall estimates of the latter were due to several loci with exceptionally high regional differentiation. The allozyme- and RAPD-specific patterns of differentiation were concordant in general with one another, but not with quantitative trait differentiation. Divergent selection on quantitative traits inferred from very high regional Q(ST) was in full agreement with our previous results obtained from a test of local adaptation and multilevel selection analysis. In contrast, most variation in allozyme and RAPD variation was neutral, although several allozyme loci and RAPD markers were exceptional in their levels of regional differentiation. However, it is not possible to answer the question whether these exceptional loci are directly involved in the response to selection pressure or merely linked to the selected loci. The fact that Q(ST) and F(ST) did not differ at the population scale, that is, within regions, but differed at the regional scale, for which local adaptation has been previously shown, implies that comparison of the level of subdivision in quantitative traits, as compared with molecular markers, is indicative of adaptive population differentiation only when sampling is carried out at the appropriate scale.     

Combined Assessment of Genetic Variability in Populations of Brown Trout (Salmo trutta L.) Based on Allozymes, Microsatellites, and RAPD Markers

Genetic variability within and among four Spanish natural populations of Salmo trutta L. was evaluated on the basis of 25 enzyme loci, 3 microsatellite loci, and 9 randomly amplified polymorphic DNAs (RAPDs). A total of 21 allelic markers were found, 12 of which were reported by microsatellites, whereas enzyme and RAPD accounted only for 6 and 3, respectively. Genetic variation within samples was significantly higher for microsatellites and RAPD than for enzyme loci. Although all methods reported a high degree of allelic heterogeneity among samples, also revealing a high degree of gene diversity, genetic relationships depicted by UPGMA dendrograms closely agreed for all kinds of data. Microsatellite loci appeared to be the most feasible technique when searching for specific alleles for a population or an area, owing to the higher number of allelic variants found. Received July 1, 1998; accepted January 14, 1999     

RAPD polymorphism of wild emmer wheat populations, Triticum dicoccoides, in Israel

 Genetic diversity in random amplified polymorphic DNAs (RAPDs) was studied in 110 genotypes of the tetraploid wild progenitor of wheat, Triticum dicoccoides, from 11 populations sampled in Israel and Turkey. Our results show high level of diversity of RAPD markers in wild wheat populations in Israel. The ten primers used in this study amplified 59 scorable RAPD loci of which 48 (81.4%) were polymorphic and 11 monomorphic. RAPD analysis was found to be highly effective in distinguishing genotypes of T. dicoccoides originating from diverse ecogeographical sites in Israel and Turkey, with 95.5% of the 100 genotypes correctly classified into sites of origin by discriminant analysis based on RAPD genotyping. However, interpopulation genetic distances showed no association with geographic distance between the population sites of origin, negating a simple isolation by distance model. Spatial autocorrelation of RAPD frequencies suggests that migration is not influential. Our present RAPD results are non-random and in agreement with the previously obtained allozyme patterns, although the genetic diversity values obtained with RAPDs are much higher than the allozyme values. Significant correlates of RAPD markers with various climatic and soil factors suggest that, as in the case of allozymes, natural selection causes adaptive RAPD ecogeographical differentiation. The results obtained suggest that RAPD markers are useful for the estimation of genetic diversity in wild material of T. dicoccoides and the identification of suitable parents for the development of mapping populations for the tagging of agronomically important traits derived from T. dicoccoides. Received: 13 July 1998 / Accepted: 13 August 1998     

Estimating pairwise relatedness from dominant genetic markers

Knowledge of the genetic relatedness between a pair of individuals is important in many research areas of quantitative genetics, conservation genetics, evolution and ecology. Many estimators have been developed to estimate such pairwise relatedness (r) using codominant markers, such as microsatellites and enzymes. In contrast, only two estimators are proposed to use dominant markers, such as random amplified polymorphic DNAs (RAPDs) and amplified fragment length polymorphisms (AFLPs), in relatedness inference. They are both biased estimators, and their statistical properties and robustness to the sampling errors in allele frequency have not been investigated. In this short paper, I propose two new pairwise relatedness estimators for dominant markers, and compare them in precision, accuracy and robustness to sampling with the two previous estimators using simulations. It was found that the new estimator based on the least squares approach is unbiased when allele frequencies are known or estimated from a sample without correcting for sampling effects. It has, however, a low precision and as a result, an intermediate overall performance among the four estimators in terms of the mean squared deviation (MSD) of estimates from actual values of r. The new estimator based on a similarity index is slightly biased but has generally the lowest MSD among the four estimators compared, regardless of the number of loci, type of actual relationships, allele frequencies known or estimated from samples. Simulations also show that the confidence intervals estimated by bootstrapping are appropriate for different estimators provided that the number of loci used in the estimation is not small.     

Use of a subset of doubled-haploid lines for RAPD interval mapping in barley.

Molecular markers have been used in barley to locate genes and quantitative trait loci. Only a few RAPD markers have been located on barley marker maps. The objectives of this study were (i) to place RAPD markers in specific intervals on the barley linkage map developed from the cross Steptoe (S) x Morex (M), (ii) to examine the distribution of RAPD markers, and (iii) to compare markers amplified by Taq DNA polymerase with those amplified by the Stoffel fragment of Taq DNA polymerase. Screening of DNA from S and M with 362 decamer primers identified 85 that amplified 127 reliable RAPDs. A subset of 15 doubled-haploid (DH) lines from the 150 DH line mapping population was used to place these RAPD markers in intervals on the SM map. This subset can be used for rapid placement of any new markers on the SM linkage map. Most of the RAPD markers were dominant but four codominant RAPDs were identified. The RAPDs were not evenly distributed, with many clustered around the centromeric region of each chromosome. Two of these clusters were located in intervals larger than 15 cM. Testing of 38 to 42 additional DH lines provided more precise placement of eight of the markers in these clusters. Reliable RAPDs were detected with 44% of the primers tested with the Stoffel fragment, but with only 17% of the primers tested with Taq DNA polymerase. These RAPDs provide additional markers for use in barley improvement.     

The genetic structure of endangered populations in the Cranberry Fritillary,Boloria aquilonaris (Lepidoptera,Nymphalidae): RAPDs vs allozymes

The genetic population structure of the Cranberry Fritillary Boloria aquilonaris was studied using both RAPDs (random amplified polymorphic DNA) and allozymes. In Belgium, B. aquilonaris has a naturally fragmented distribution that has been accentuated due to human activity during the last century. The genetic population structure of this butterfly was analysed at the regional (several Ardenne uplands) and at the landscape level (several populations within an Ardenne upland). Both population genetic markers confirmed results from a previous CMR study at the landscape scale. At the regional scale however, important incongruences were observed between RAPDs and allozymes. The average gene diversity for the RAPD data was twice that of the allozyme data. The degree of population subdivision was also much greater for RAPDs than for allozymes. The UPGMA clusters produced by each of these markers differed significantly. We believe that, given the higher rate of mutation of RAPDs and the greater number of loci assayed by this method, RAPDs reveal a more accurate and recent population genetic structure than allozymes.     

Inheritance of random amplified polymorphic DNA (RAPD-PCR) markers and their use in population genetic studies of Calomys musculinus (Rodentia,Muridae), the reservoir of Argentine hemorrhagic fever

In order to assess the reliability of RAPD markers in the estimation of the genetic structure of natural populations of the murid rodent Calomys musculinus (reservoir of Junin virus, ethiological agent of Argentine Hemorrhagic Fever), we have analyzed the heritability of RAPD bands in 10 parents and their offspring (33 individuals). Fourteen out of a total of 119 bands obtained were absent in the parental patterns, but consistently amplified in offspring from some families. These bands can be eliminated from analyses. Overall degree of band sharing between individuals, including non-parental bands, correctly grouped members of a family in the same cluster in a UPGMA tree, with a high bootstrap percentage. Results support the usefulness of RAPDs as hereditable markers. One hundred polymorphic RAPD loci were identified in three natural populations of C. musculinus. Mean expected heterozygosity in three natural populations ranged from 0.206 to 0.220. Allele frequency based and phenotype based measures of genetic differentiation among natural populations of C. musculinus gave similar results (Weir and Cockerham's theta = 0.133; Excoffier et al.'s phi = 0.127). These values were considerably higher than those found previously using allozymes as genetic markers, and are compatible with moderate to low levels of gene flow among populations.     

ABC inference of multi‐population divergence with admixture from unphased population genomic data

Rapidly developing sequencing technologies and declining costs have made it possible to collect genome‐scale data from population‐level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three‐population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation‐based cross‐validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.     

Random amplified polymorphic DNA (RAPD) and quantitative trait analyses across a major phylogeographical break in the Mediterranean ragwort Senecio gallicus Vill. (Asteraceae)

Random amplified polymorphic DNA (RAPD) and quantitative trait variation of the widespread and ephemeral Senecio gallicus were surveyed in 11 populations sampled from the Iberian Peninsula and southern France. The aim of the study was to compare population relationships and levels of geographical differentiation with chloroplast (cp) DNA and allozyme variation assessed previously in the same populations. Employing multivariate statistics, a moderate level of intraspecific differentiation was observed among populations from Iberian coastal and inland regions for both RAPDs and quantitative traits. However, RAPDs provided greater resolution in identifying additional population structure within the hypothesized, Pleistocene refugial source area of the species in coastal Iberia. A major part of the geographical subdivision in RAPD and quantitative traits was concordant with the coastal vs. inland divergence as previously inferred from cpDNA haplotype frequencies, but strongly contrasted with the geographical uniformity of the species for allozymes. This concordance across various nuclear and cytoplasmic markers (RAPDs/quantitative traits, cpDNA) suggests that geographical uniformity for allozymes is more attributable to low rates of evolution and/or small genome sampling rather than high rates of pollen dispersal, slow rates of nuclear lineage sorting, or indirect balancing selection. The present study underscores the value of using additional classes of nuclear markers for narrowing the numbers of competing causal hypotheses about intraspecific cpDNA-allozyme discrepancies and their underlying evolutionary processes.     

Intra-specific genetic diversity in wild olives (Olea europaea ssp cuspidata) in Hormozgan Province, Iran

Wild olive (O. europaea ssp cuspidata) plants grow in various regions of Iran and are expected to have considerable genetic diversity due to adaptation to the various environmental conditions. We examined the genetic diversity of four populations of wild olive growing in Hormozgan Province located in southern Iran by using 30 RAPDs and 10 ISSR markers. The mean value of polymorphism for RAPD loci was 73.71%, while the value for ISSR loci was 81.74%. The Keshar population had the highest value of intra-population polymorphism for both RAPD and ISSR loci (66.86 and 62.71%, respectively), while the Tudar population had the lowest values (20.35 and 28.81%, respectively). Similarly, the highest and lowest number of effective alleles, Shannon index and Nei's genetic diversity were also found for these two populations. The highest value of H(pop)/H(sp) within population genetic diversity for RAPD and ISSR loci was found for the Keshar population (H(pop) = 0.85 and H(sp) = 0.90). OPA04-750, OPA13-650 and OPA02-350 RAPD bands were specific for Tudar, Bondon and Keshar populations, respectively, while no specific ISSR bands were observed. Analysis of molecular variance as well as the pairwise F(ST) test showed significant differences for RAPD and ISSR markers among the populations. The NJ and UPGMA trees also separated the wild olive populations from each other, indicating their genetic distinctness. UPGMA clustering of the four wild olive populations placed the Tudar population far from the other populations; Keshar and Bokhoon population samples revealed more similarity and were grouped together. We conclude that there is high genetic diversity among O. europaea ssp cuspidata populations located in southern Iran. We also found RAPD and ISSR markers to be useful molecular tools to discriminate and evaluate genetic variations in wild olive trees.     

Optimal sampling strategy for estimation of spatial genetic structure in tree populations

Fine-scale spatial genetic structure (SGS) in natural tree populations is largely a result of restricted pollen and seed dispersal. Understanding the link between limitations to dispersal in gene vectors and SGS is of key interest to biologists and the availability of highly variable molecular markers has facilitated fine-scale analysis of populations. However, estimation of SGS may depend strongly on the type of genetic marker and sampling strategy (of both loci and individuals). To explore sampling limits, we created a model population with simulated distributions of dominant and codominant alleles, resulting from natural regeneration with restricted gene flow. SGS estimates from subsamples (simulating collection and analysis with amplified fragment length polymorphism (AFLP) and microsatellite markers) were correlated with the 'real' estimate (from the full model population). For both marker types, sampling ranges were evident, with lower limits below which estimation was poorly correlated and upper limits above which sampling became inefficient. Lower limits (correlation of 0.9) were 100 individuals, 10 loci for microsatellites and 150 individuals, 100 loci for AFLPs. Upper limits were 200 individuals, five loci for microsatellites and 200 individuals, 100 loci for AFLPs. The limits indicated by simulation were compared with data sets from real species. Instances where sampling effort had been either insufficient or inefficient were identified. The model results should form practical boundaries for studies aiming to detect SGS. However, greater sample sizes will be required in cases where SGS is weaker than for our simulated population, for example, in species with effective pollen/seed dispersal mechanisms.     

Identification of RAPD markers linked to genetic factors controlling the milling energy requirement of barley

Doubled haploid (DH) populations of barley have been used in combination with PCR-based polymorphic-assay procedures to identify molecular markers linked to genes controlling the milling energy requirement of the grain. Milling energy (ME) is a quantitative trait and locating individual quantitative trait loci (QTLs) involved the construction of bulks by combining DNA from DH families representing the extreme members of the distribution for ME. In addition, the individuals had alternative alleles at theRrn2 locus that has previously been shown to be linked to an ME QTL. The DNA bulks were screened with Randomly Amplified Polymorphic DNA (RAPD) markers and polymorphic amplification products tested for linkage to genes influencing the expression of ME in a DH population. Several markers were identified which are linked to a QTL controlling ME and the recombination fraction determined by maximum likelihood procedures. The results indicate that DHs in combination with RAPDs and bulked segregant analysis provide an efficient method for locating QTLs in barely. Furthermore, this approach is applicable to mapping other QTLs in a range of organisms from which DH or recombinant inbred lines can be extracted.     

A maximum-likelihood estimation of pairwise relatedness for autopolyploids

Relatedness between individuals is central to ecological genetics. Multiple methods are available to quantify relatedness from molecular data, including method-of-moment and maximum-likelihood estimators. We describe a maximum-likelihood estimator for autopolyploids, and quantify its statistical performance under a range of biologically relevant conditions. The statistical performances of five additional polyploid estimators of relatedness were also quantified under identical conditions. When comparing truncated estimators, the maximum-likelihood estimator exhibited lower root mean square error under some conditions and was more biased for non-relatives, especially when the number of alleles per loci was low. However, even under these conditions, this bias was reduced to be statistically insignificant with more robust genetic sampling. We also considered ambiguity in polyploid heterozygote genotyping and developed a weighting methodology for candidate genotypes. The statistical performances of three polyploid estimators under both ideal and actual conditions (including inbreeding and double reduction) were compared. The software package POLYRELATEDNESS is available to perform this estimation and supports a maximum ploidy of eight.     

分子水平的遗传多样性及其测量方法

遗传多样性水平是一项很重要的数据。目前从分子水平量化遗传多样性的方法以等位酶分析、ＲＦＬＰ分析和ＲＡＰＤ分析为代表。遗传变异在基因组中并非随机分布,所以取样方式对分析结果的影响不容忽视。本文叙述了遗传变异的产生和分布,并以此为基础比较了上述３种分析方法的理论与实践。     

Population genetic study in the Balearic endemic plant species Digitalis minor (Scrophulariaceae) using RAPD markers

Digitalis minor (Scrophulariaceae) is a cardenolide-producing plant endemic to the eastern Balearic Islands (Mallorca, Menorca, and Cabrera) that occurs in two morphologically distinct varieties: D. minor var. minor (pubescent) and D. minor var. palaui (glabrous). Levels and patterns of genetic diversity in 162 individuals from 17 D. minor populations across the entire geographic range were assessed using random amplified polymorphic DNA (RAPD) markers. Comigrating RAPD fragments tested were found to be homologous by Southern hybridization in both var. minor and var. palaui. To avoid bias in parameter estimation, analyses of population genetic structure were restricted to those RAPD bands that fulfilled the 3/N criterion (observed frequencies were less than 1 - [3/N] in each population) either among or within each island. Analyses of molecular variance (AMOVAs) with distances among individuals corrected for the dominant nature of RAPD (genotypic analysis) showed low values (1.57-17.55%) of between-population variability, indicating a relatively restricted population differentiation as expected for an outcrossing species such as D. minor. Nested AMOVAs demonstrated, however, a not significant partitioning of genetic diversity among Mallorca, Menorca, and Cabrera islands. Estimates of the Wright, Weir, and Cockerham and the Lynch and Milligan F(ST) from null allele frequencies corroborated AMOVA partitioning and provided evidence for population differentiation in D. minor. Our RAPD data did not show significant differences between pubescent and glabrous populations of D. minor, suggesting a failure to find a correlation between the RAPD loci and this morphological trait.     

Comparative evaluation of a new effective population size estimator based on approximate bayesian computation

We describe and evaluate a new estimator of the effective population size (N(e)), a critical parameter in evolutionary and conservation biology. This new "SummStat" N(e) estimator is based upon the use of summary statistics in an approximate Bayesian computation framework to infer N(e). Simulations of a Wright-Fisher population with known N(e) show that the SummStat estimator is useful across a realistic range of individuals and loci sampled, generations between samples, and N(e) values. We also address the paucity of information about the relative performance of N(e) estimators by comparing the SummStat estimator to two recently developed likelihood-based estimators and a traditional moment-based estimator. The SummStat estimator is the least biased of the four estimators compared. In 32 of 36 parameter combinations investigated using initial allele frequencies drawn from a Dirichlet distribution, it has the lowest bias. The relative mean square error (RMSE) of the SummStat estimator was generally intermediate to the others. All of the estimators had RMSE > 1 when small samples (n = 20, five loci) were collected a generation apart. In contrast, when samples were separated by three or more generations and N(e) < or = 50, the SummStat and likelihood-based estimators all had greatly reduced RMSE. Under the conditions simulated, SummStat confidence intervals were more conservative than the likelihood-based estimators and more likely to include true N(e). The greatest strength of the SummStat estimator is its flexible structure. This flexibility allows it to incorporate any potentially informative summary statistic from population genetic data.     

RFLP and RAPD mapping in flax (Linum usitatissimum)

A map of flax (Linum usitatissimum) using restriction fragment length polymorphisms (RFLPs) and random amplified polymorphic DNAs (RAPDs), and comprising 15 linkage groups containing 94 markers, has been developed covering about 1000 cM. The mapping populations were the F₂ populations from two crosses between diverse cultivars. From one cross, CI1303 and Stormont Cirrus, 20 RFLP and 520 RAPD markers were analyzed. Thirteen RFLP and 80 RAPD markers were on the 15 linkage groups, in addition to one sequence-tagged site (STS). Seven polymorphic RAPD markers were found to have unusual segregation patterns. RAPDs were expressed as dominant markers, but for these markers a prevalence of the progeny lacked a band rather than the expected one-fourth ratio. However, these exceptions may be related to the instability of the genome of Stormont Cirrus in which stable and heritable genomic changes can be induced by environmental factors. The current map could be used for the identification of markers linked to loci controlling the ability to generate heritable changes in response to environmental growth conditions, and to develop anchor loci with STSs for a more general application. Received: 20 March 1999 / Accepted: 16 December 1999