An empirical verification of population assignment methods by marking and parentage data: hatchery and wild steelhead (Oncorhynchus mykiss) in Forks Creek, Washington, USA

Assignment tests are increasingly applied in ecology and conservation, although empirical comparisons of methods are still rare or are restricted to few of the available approaches. Furthermore, the performance of assignment tests in cases with low population differentiation, violations of Hardy-Weinberg equilibrium and unbalanced sampling designs has not been verified. The release of adult hatchery steelhead to spawn in Forks Creek in 1996 and 1997 provided an opportunity to compare the power of different assignment methods to distinguish their offspring from those of sympatric wild steelhead. We compared standard assignment methods requiring baseline samples (frequency, distance and Bayesian) and clustering approaches with and without baseline information, using six freely available computer programs. Assignments were verified by parentage data obtained for a subset of returning offspring. All methods provided similar assignment success, despite low differentiation between wild and hatchery fish (F(ST) = 0.02). Bayesian approaches with baseline data performed best, whereas the results of clustering methods were variable and depended on the samples included in the analysis and the availability of baseline information. Removal of a locus with null alleles and equalizing sample sizes had little effect on assignments. Our results demonstrate the robustness of most assignment tests to low differentiation and violations of assumptions, as well as their utility for ecological studies that require correct classification of different groups.     

Comparative Assessment of Genetic Variability in the Populations of Endemic and Endangered Yellow Catfish, Horabagrus brachysoma (Teleostei: Horabagridae), Based on Allozyme, RAPD, and Microsatellite Markers

The comparative assessment of genetic diversity using allozymes, random amplified polymorphic DNA (RAPD), and microsatellite markers was conducted in endemic and endangered yellow catfish (Horabagrus brachysoma) sampled from three locations in Western Ghats river systems of India. Among the three markers, microsatellites show more polymorphism, having 100% polymorphic loci, whereas allozymes show the least (56%). In RAPD, 60.5% of fragments were polymorphic. Observed heterozygosity and F(ST) values were very high in microsatellites, compared with the other markers. Microsatellite and RAPD markers reported a higher degree of genetic differentiation than allozymes among the populations depicted by pairwise F(ST)/G(ST), AMOVA, Nei's genetic distance, and UPGMA dendrogram. The three classes of markers demonstrated striking genetic differentiation between pairs of H. brachysoma populations. The data emphasize the need for fishery management, conservation, and rehabilitation of this species.     

Toward the development of highly homozygous diploid potato lines using the self-compatibility controlling Sli gene.

Cultivated diploid potatoes (2n = 2x = 24) are self-incompatible, but can be altered to become self-compatible using the Sli gene. Previously, a diploid clone 97H32-6 was selfed up to S3 using the Sli gene. To explore the usefulness of the Sli gene for the production of highly homozygous diploid potatoes, 2 S4 families from the above 97H32-6 derived S3 lines (inbred series A) and 3 S5 families by continuous selfings from a different F1 (= S0) plant (inbred series B) were developed. The level of heterozygosity and the location of heterozygous loci on the genetic map were investigated using RFLP and AFLP markers. The average heterozygosity levels of the originally heterozygous loci decreased from 100% in S0 to 10.7% in S4 and 8.6% in S5 (inbred series A and B, respectively). The average rate of reduction in heterozygosity per generation (38.4% and 38.5% for inbred series A and B, respectively) was lower than the theoretically expected rate (50%). However, none of the loci or chromosome sections was exclusively heterozygous in the advanced self-progeny. Thus, highly homozygous and seed-propagated diploid potatoes could be obtained by repeated selfing using the Sli gene.     

Selection of high heterozygosity popcorn varieties in Brazil based on SSR markers

We analyzed genetic structure and diversity among eight populations of popcorn, using SSR loci as genetic markers. Our objectives were to select SSR loci that could be used to estimate genetic diversity within popcorn populations, and to analyze the genetic structure of promising populations with high levels of heterozygosity that could be used in breeding programs. Fifty-seven alleles (3.7 alleles per locus) were detected; the highest effective number of alleles (4.21) and the highest gene diversity (0.763) were found for the Umc2226 locus. A very high level of population differentiation was found (F(ST) = 0.3664), with F(ST) for each locus ranging from 0.1029 (Umc1664) to 0.6010 (Umc2350). This analysis allowed us to identify SSR loci with high levels of heterozygosity and heterozygous varieties, which could be selected for production of inbred lines and for developing new cultivars.     

Inbreeding in Japanese quail estimated by pedigree and microsatellite analyses

Accurately estimating inbreeding is important because inbreeding reduces fitness and production traits in populations. We analyzed information from pedigrees and from microsatellite markers to estimate inbreeding in a line of Japanese quail derived from a randombred line (QO) and maintained for 17 generations by pedigreed matings of brothers to groups of sisters. Pedigree data were used to calculate the inbreeding coefficient (F(IT)), which is the level of inbreeding based on a reference ancestor. Data from analysis of 14 microsatellite markers in the inbred and QO lines were used to calculate the population differentiation (F(ST)) of the lines caused by inbreeding. The F(IT) was then calculated as F(IT) = F(IS) + (1 - F(IS)) x F(ST), where F(IS) is the level of inbreeding in the inbred line. Observed heterozygosity from analysis of the microsatellite markers of the QO and inbred lines was 0.43 and 0.21, respectively, and the number of alleles was 3.29 and 1.93, demonstrating a reduction of genetic diversity in the inbred line. The F(IT) of the inbred line calculated from the pedigree and microsatellite marker analyses was 0.69 +/- 0.07 and 0.57 +/- 0.33, respectively. These data suggest that pedigree analysis was more accurate than microsatellite marker analyses for estimating inbreeding in this line of Japanese quail.     

Population structure, genetic variation and morphological diversity in indigenous sheep of Ethiopia

We investigated genetic and morphological diversity and population structure of 14 traditional sheep populations originating from four ecological zones in Ethiopia (sub-alpine, wet highland, sub-humid lowland and arid lowland). All animals (n = 672) were genotyped for 17 microsatellite markers and scored for 12 morphological characters. The sheep were initially classified as fat-tailed (11 populations), thin-tailed (one population) and fat-rumped sheep (two populations). These classifications are thought to correspond to three consecutive introduction events of sheep from the Near-East into East Africa. For the 14 populations, allelic richness ranged from 5.87 to 7.51 and expected heterozygosity (H(E)) from 0.66 to 0.75. Genetic differentiations (F(ST) values) between all pairs of populations, except between sub-alpine populations, were significantly different from zero (P < 0.001). Cluster analysis of morphological characters and a dendrogram constructed from genetic distances were broadly consistent with the classification into fat-tailed, thin-tailed and fat-rumped sheep. Bayesian cluster analysis using microsatellite markers indicated that there has been further genetic differentiation after the initial introduction of sheep into Ethiopia. Investigation of factors associated with genetic variation showed that an isolation-by-distance model, independently of other factors, explained most of the observed genetic variation. We also obtained a strong indication of adaptive divergence in morphological characters, patterns of morphological variation being highly associated with ecology even when the effect of neutral genetic divergence (F(ST)) was parcelled out in partial Mantel tests. Using a combination of F(ST) values, Bayesian clustering analysis and morphological divergence, we propose a classification of Ethiopian sheep into six breed groups and nine breeds.     

Source population of dispersing rock-wallabies (Petrogale lateralis) identified by assignment tests on multilocus genotypic data

The ability to confidently identify or exclude a population as the source of an individual has numerous powerful applications in molecular ecology. Several alternative assignment methods have recently been developed and are yet to be fully evaluated with empirical data. In this study we tested the efficacy of different assignment methods by using a translocated rock-wallaby (Petrogale lateralis) population, of known provenance. Specimens from the translocated population (n = 43), its known source population (n = 30) and four other nearby populations (n = 19-32) were genotyped for 11 polymorphic microsatellite loci. The results identified Bayesian clustering, frequency and Bayesian methods as the most consistent and accurate, correctly assigning 93-100% of individuals up to a significance threshold of P = 0.01. Performance was variable among the distance-based methods, with the Cavalli-Sforza and Edwards chord distance performing best, whereas Goldstein et al.'s (deltamu)2 consistently performed poorly. Using Bayesian clustering, frequency and Bayesian methods we then attempted to determine the source of rock-wallabies which have recently recolonized an outcrop (Gardners) 8 km from the nearest rock-wallaby population. Results indicate that the population at Gardners originated via a recent dispersal event from the eastern end of Mt. Caroline. This is only the second published record of dispersal by rock-wallabies between habitat patches and is the longest movement recorded to date. Molecular techniques and methods of analysis are now available to allow detailed studies of dispersal in rock-wallabies and should also be possible for many other taxa.     

Microsatellite Analysis of Homozygosity Progression of Heterozygous Genotypes Segregating in the Rice Subspecies Cross Pei��ai64s/Nipponbare

Progression to homozygosity of heterozygous genotypes was studied in a cross of the rice subspecies Pei'ai64s and Nipponbare, using a set of 157 polymorphic microsatellite (SSR) markers. The segregation of heterozygous genotypes ranged from 49.13% in the F(2) population to 4.52% in the F(6) population (progression value 11.15%). The heterozygous genotypes were widely distributed in 180 F(2) plants, 330 F(6) lines, and 157 SSR markers. Homozygosity progression showed a wide distribution in plants and SSR markers but not in chromosomes. The segregation of heterozygous genotypes was not significant between populations but varied greatly in F(2) plants, F(6) lines, and SSR markers. The correlation between the progression to homozygosity and the heterozygosity of SSR markers was significant at the chromosome level. The segregation of heterozygous genotypes in plants, SSR markers, and chromosomes was not completely in accordance with Mendel's law. This information will help rice geneticists and breeders to understand heterozygous genotype segregation at the DNA level and to screen special markers for breeding.     

Characterization of 15 microsatellite loci in the pulmonate snail Pseudosuccinea columella (Mollusca, Gastropoda)

We characterized 15 new variable microsatellites in the freshwater snail Pseudosuccinea (Lymnaea) columella, as well as conditions for multiplexing and simultaneously genotyping sets of loci. Two to six alleles were detected per locus over the six populations studied. Gene diversity ranged from 0.000 to 0.498, but essentially no heterozygous individuals were observed. This resulted in extremely high F(IS) estimates, and therefore high selfing rates. The F(ST) estimates ranged from 0.18 to 1 among populations, but was generally high. These markers will constitute efficient tools for investigating the population structure of this invasive species. Cross-species amplification was on the whole unsuccessful.     

Sex-biased dispersal in a migratory bat: a characterization using sex-specific demographic parameters

We studied the noctule bat (Nyctalus noctula), in which the mitochondrial F(ST) is about 10 times that revealed by nuclear markers, to address two questions. We first verified whether random dispersal of one sex is compatible with highly contrasted mitochondrial and nuclear population structures. Using computer simulations, we then assessed the power of multilocus population differentiation tests when the expected population structure departs only slightly from panmixia. Using an island model with sex-specific demographic parameters, we found that random male dispersal is consistent with the population structure observed in the noctule. However, other parameter combinations are also compatible with the data. We computed the minimum sex bias in dispersal (at least 69% of the dispersing individuals are males), a result that would not be available if we had used more classical population genetic models. The power of multilocus population differentiation tests was unexpectedly high, the tests being significant in almost 100% of the replicates, although the observed population structure infered from nuclear markers was extremely low (F(ST) = 0.6%).     

Effects of microsatellite null alleles on assignment testing

Microsatellite null alleles are found to a varying degree across all taxa. They are problematic as they may inflate measures of genetic differentiation and create false homozygotes. Although there are several methods for correcting allele frequencies for null alleles and enable estimations of F(ST), much less is known about how null alleles affect assignment testing. Data presented here, based on simulations, show that the percentage of correctly assigned individuals in model-based clustering and Bayesian assignment methods were slightly, though significantly, reduced in the presence of null alleles (frequency range from 0.000 to 0.913). The bias in assignment tests caused by null alleles lead to a slight reduction in the power to correctly assigned individuals (0.2 and 1.0 percent units for STRUCTURE- and 2.4 percent units for GENECLASS-based assignment tests). Further, the presence of null alleles caused a small, however, significant overestimation of F(ST). Consequently, microsatellite loci affected by null alleles would probably not alter the overall outcome of assignment testing and could therefore be included in these types of studies. Nevertheless, loci prone to null alleles should be used with caution as they lower the power of assignment tests and alter the accuracy of F(ST), and loci less prone to null alleles should always be preferred.     

High interpopulation genetic differentiation and unidirectional linear migration patterns in Myricaria laxiflora (Tamaricaceae), an endemic riparian plant in the Three Gorges Valley of the Yangtze River

Myricaria laxiflora is restricted to the riverbanks of the Yangtze River valley and will be completely lost owing to the construction of the Three Gorges Dam. Genetic diversity and structure of nine natural and one ex situ populations were investigated using amplified fragment length polymorphisms (AFLPs). A moderate level of gene diversity was found in natural populations, while the ex situ population had the highest. The F statistics calculated by different approaches consistently revealed a high genetic differentiation among natural populations, contributing >45% of the total gene diversity. The Bayesian-based analysis differentiated nine independent populations in accordance with the sites sampled. Estimates of gene flow by F(ST) and coalescent-based simulation analysis indicated a restricted recurrent gene exchange among populations (Nm = 0.290-0.401), whereas genetic distance-based clustering and coalescent-based assignment analyses revealed significant genetic isolation among populations. The migration pattern in M. laxiflora is best explained by a classical metapopulation model, but with a unique unidirectional direction underlined by hydrochoric force that drove dispersal of seeds and propagules from upstream toward downstream populations. Previous efforts in preserving genomic integrity in ex situ conservation were evaluated, and the results provide valuable information to formulate conservation guidelines for successfully reintroducing M. laxiflora to the wild.     

An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure

There is great interest in using amplified fragment length polymorphism (AFLP) markers because they are inexpensive and easy to produce. It is, therefore, possible to generate a large number of markers that have a wide coverage of species genomes. Several statistical methods have been proposed to study the genetic structure using AFLPs but they assume Hardy-Weinberg equilibrium and do not estimate the inbreeding coefficient, F(IS). A Bayesian method has been proposed by Holsinger and colleagues that relaxes these simplifying assumptions but we have identified two sources of bias that can influence estimates based on these markers: (i) the use of a uniform prior on ancestral allele frequencies and (ii) the ascertainment bias of AFLP markers. We present a new Bayesian method that avoids these biases by using an implementation based on the approximate Bayesian computation (ABC) algorithm. This new method estimates population-specific F(IS) and F(ST) values and offers users the possibility of taking into account the criteria for selecting the markers that are used in the analyses. The software is available at our web site (http://www-leca.ujf-grenoble.fr/logiciels.htm). Finally, we provide advice on how to avoid the effects of ascertainment bias.     

Number of SNPS loci needed to detect population structure

The study of the association of polymorphic genetic markers with common diseases is one of the most powerful tools in modern genetics. Interest in single nucleotide polymorphisms (SNPs) has steadily grown over the last decade. SNPs are currently the most developed markers in the human genome because they have a number of advantages over other marker types. One of the critical problems responsible for 'spurious' association findings in case-control studies is population stratification. There are many statistical approaches developed for detecting population heterogeneity. However the power to detect population structure by known methods is highly dependent on the number of loci utilised. We performed an analysis of SNPs data available in the public domain from The Single Nucleotide Consortia Ltd. (TSCL). Three populations, Afro-American, Asian and Caucasian, were compared. Estimation of the minimum number of SNPs loci necessary for detection of the population structure was performed. Two clustering approaches, distance-based and model-based, were compared. The model-based approach was superior when compared with the distance-based method. We found more than 65 random SNPs loci are required for identifying distinct geographically separated populations. Increasing the number of markers to over 100 raises the probability of correct assignment of a particular individual to an origin group to over 90%, even with conventional clustering methods.     

Evaluation of genetic diversity of Fusarium head blight resistance in European winter wheat

Genetic diversity in relation to Fusarium head blight (FHB) resistance was investigated among 295 European winter wheat cultivars and advanced breeding lines using 47 wheat SSR markers. Twelve additional wheat lines with known FHB resistance were included as reference material. At least one SSR marker per chromosome arm, including SSR markers reported in the literature with putative associations with QTLs for FHB resistance, were assayed to give an even distribution of SSR markers across the wheat genome. A total of 404 SSR alleles were detected. The number of alleles per locus ranged from 2 to 21, with an average of 8.6 alleles. The polymorphism information content of the SSR markers ranged from 0.13 (Xwmc483) to 0.87 (Xwmc607), with an average of 0.54. Cluster analysis was performed by both genetic distance-based and model-based methods. In general, the dendrogram based on unweighted pair-group method with arithmetic averages showed similar groupings to the model-based analysis. Seven clusters were identified by the model-based method, which did not strictly correspond to geographical origin. The FHB resistance level of the wheat lines was evaluated in field trials conducted over multiple years or locations by assessing the following traits: % FHB severity, % FHB incidence, % diseased kernels, in spray inoculation trials, and % FHB spread and % wilted tips, in point inoculation trials. Association analysis between SSR markers and the FHB disease traits detected markers significantly associated with FHB resistance, including some that have not been previously reported. The percentage of variance explained by each individual marker was, however, rather low. Haplotype analysis revealed that the FHB-resistant European wheat lines do not contain the 3BS locus derived from Sumai 3. The information generated in this study will assist in the selection of parental lines in order to increase the efficiency of breeding efforts for FHB resistance.     

Introgression of the avian naked neck gene assisted by DNA fingerprints

Theoretical predictions suggest that DNA markers can be useful tools for genomic selection in gene introgression programmes. An experiment was carried out to evaluate the efficiency of using multi-locus DNA markers in an introgression programme designed to transfer the naked neck gene from a donor to a recipient chicken line. The donor line was a commercial egg layer chicken stock heterozygous at the naked neck locus (Na/na⁺), while the recipients were from a Cornish broiler line. These two lines differ markedly in their average body weight, a quantitative trait that can also represent the comprehensive differences between the genomes of the two lines involved. Three groups of naked neck BC₁ individuals were selected according to the following criteria: (i) low band-sharing with their donor grandsires evaluated by multi-locus DNA markers, (ii) high body weight at six weeks of age, and (iii) selection at random as a control group. Birds from each of these groups were mated at random to individuals from the heavier Cornish line to produce three groups of BC₂ individuals whose body weights were recorded weekly from three to seven weeks of age. Results indicated that BC₂ birds obtained from BC₁ parents selected for band-sharing levels and those selected for body weight, performed equally well at 4–7 weeks of age; both were 3.1–3.9% heavier than birds from the randomly selected group. The additional genome recovery of the heavier broiler line, obtained by DNA markers, was found to be in agreement with theoretically predicted values.     

Molecular characterization and genetic structure of the Nero Siciliano pig breed

Nero Siciliano is an autochthonous pig breed that is reared mainly in semi-extensive systems in northeastern Sicily. Despite its economic importance and well-appreciated meat products, this breed is currently endangered. Consequently, an analysis of intra-breed variability is a fundamental step in preserving this genetic resource and its breeding system. In this work, we used 25 microsatellite markers to examine the genetic composition of 147 unrelated Nero Siciliano pigs. The total number of alleles detected (249, 9.96 per locus) and the expected heterozygosity (0.708) indicated that this breed had a high level of genetic variability. Bayesian cluster analysis showed that the most likely number of groups into which the sample could be partitioned was nine. Based on the proportion of each individuals genome derived from ancestry, pigs with at least 70% of their genome belonging to one cluster were assigned to that cluster. The cluster size ranged from 7 to 17 (n = 108). Genetic variability in this sub-population was slightly lower than in the whole sample, genetic differentiation among clusters was moderate (F(ST) 0.125) and the F(IS) value was 0.011. NeighborNet and correspondence analysis revealed two clusters as the most divergent. Molecular coancestry analysis confirmed the good within-breed variability and highlighted the clusters that retained the highest genetic diversity.     

A population genetics study of Anopheles darlingi (Diptera: Culicidae) from Colombia based on random amplified polymorphic DNA-polymerase chain reaction and amplified fragment lenght polymorphism markers

The genetic variation and population structure of three populations of Anopheles darlingi from Colombia were studied using random amplified polymorphic markers (RAPDs) and amplified fragment length polymorphism markers (AFLPs). Six RAPD primers produced 46 polymorphic fragments, while two AFLP primer combinations produced 197 polymorphic fragments from 71 DNA samples. Both of the evaluated genetic markers showed the presence of gene flow, suggesting that Colombian An. darlingi populations are in panmixia. Average genetic diversity, estimated from observed heterozygosity, was 0.374 (RAPD) and 0.309 (AFLP). RAPD and AFLP markers showed little evidence of geographic separation between eastern and western populations; however, the F ST values showed high gene flow between the two western populations (RAPD: F ST = 0.029; Nm: 8.5; AFLP: F ST = 0.051; Nm: 4.7). According to molecular variance analysis (AMOVA), the genetic distance between populations was significant (RAPD:phiST = 0.084; AFLP:phiST = 0.229, P < 0.001). The F ST distances and AMOVAs using AFLP loci support the differentiation of the Guyana biogeographic province population from those of the Chocó-Magdalena. In this last region, Chocó and Córdoba populations showed the highest genetic flow.     

Nei's to Bayes': comparing computational methods and genetic markers to estimate patterns of genetic variation in Tolpis (Asteraceae)

Accurate determination of patterns of genetic variation provides a powerful inferential tool for studies of evolution and conservation. For more than 30 years, enzyme electrophoresis was the preferred method for elucidating these patterns. As a result, evolutionary geneticists have acquired considerable understanding of the relationship between patterns of allozyme variation and aspects of evolutionary process. Myriad molecular markers and statistical analyses have since emerged, enabling improved estimates of patterns of genetic diversity. With these advances, there is a need to evaluate results obtained with different markers and analytical methods. We present a comparative study of gene statistic estimates (F(ST), G(ST), F(IS), H(S), and H(T)) calculated from an intersimple sequence repeat (ISSR) and an allozyme data set derived from the same populations using both standard and Bayesian statistical approaches. Significant differences were found between estimates, owing to the effects of marker and analysis type. Most notably, F(ST) estimates for codominant data differ between Bayesian and standard approaches. Levels of statistical significance are greatly affected by methodology and, in some cases, are not associated with similar levels of biological significance. Our results suggest that caution should be used in equating or comparing results obtained using different markers and/or methods of analysis.     

Phylogeography of North African Atlas cedar (Cedrus atlantica, Pinaceae): Combined molecular and fossil data reveal a complex Quaternary history

Northwest Africa is a major hotspot of plant biodiversity, but very little is known about the Quaternary range dynamics of plant species in this region. Here we investigate the range-wide population structure and phylogeography of Atlas cedar (Cedrus atlantica), an emblematic forest tree endemic to Morocco and Algeria. We genotyped 261 individuals from 11 populations using AFLP markers. Data were analyzed using both conventional F(ST)-based techniques and Bayesian clustering. Overall population differentiation was high (F(ST) = 0.25). Two major groups of populations were identified, one distributed through the Rif and Middle Atlas mountains in Morocco and the other through the Algerian Tell Atlas and Aurès mountains as well as the Middle Atlas. Combined molecular and fossil data indicate that C. atlantica survived the Last Glacial Maximum in at least three disjunct refugia along the coast of the Mediterranean Sea, whereas the Middle Atlas, today the core of the species range, has been colonized relatively recently (<10000 yr BP). The colonization history of individual populations has left clear imprints in their present-day diversity, which may vary greatly even between nearby stands. Our study illustrates how integrating different data sources and analytical approaches can help elucidate complex range dynamics that would otherwise remain undeciphered.