Analysis of genetic variation in the Belgian Blue cattle breed using DNA sequence polymorphism at the growth hormone, low density lipoprotein receptor, alpha-subunit of glycoprotein hormones and thyroglobulin loci

New DNA sequence polymorphisms were identified at four bovine autosomal loci: growth hormone, low density lipoprotein receptor, alpha-subunit of glycoprotein hormones and thyroglobulin. Assuming independent assortment between these polymorphisms, the probabilities to be heterozygous at these four loci are 0.48, 0.36, 0.10 and 0.77 respectively, within the Belgian Blue Cattle breed (BBCB). Nucleotide diversity was estimated, showing that animals from the BBCB are heterozygous for 1/1450 nucleotides, a value significantly lower than the 1/500 value found in man. Moreover, we have estimated that the mutation rate at the cytosines of CG dinucleotides is about 10 times higher than that for other nucleotides.     

Genetic Analysis of the Fungus, Bremia Lactucae, Using Restriction Fragment Length Polymorphisms

Restriction fragment length polymorphisms (RFLPs) were developed as genetic markers for Bremia lactucae, the biotrophic Oomycete fungus which causes lettuce downy mildew. By using 55 genomic and cDNA probes, a total of 61 RFLP loci were identified among three heterothallic isolates of B. lactucae. Of these 61 RFLP loci, 53 were heterozygous in at least one of the three strains and thus were informative for linkage analysis in at least one of two F1 crosses that were performed. Analysis of the cosegregation of these 53 RFLPs, eight avirulence loci and the mating type locus allowed the construction of a preliminary genetic linkage map consisting of 13 small linkage groups. Based on the extent of linkage detected among probes, the genome of B. lactucae can be estimated to be approximately 2000 cM. Linkage was detected between a RFLP locus and an avirulence gene, providing a potential starting point for chromosome walking to clone an avirulence gene. The high frequency of DNA polymorphism in naturally occurring isolates and the proper Mendelian segregation of loci detected by low copy number probes indicates that it will be possible to construct a detailed genetic map of B. lactucae using RFLPs as markers. The method of analysis employed here should be applicable to many other outbreeding, heterozygous species for which defined inbred lines are not available.     

Association of the widespread A149P hereditary fructose intolerance mutation with newly identified sequence polymorphisms in the aldolase B gene.

Hereditary fructose intolerance (HFI) is a potentially fatal autosomal recessive disease resulting from the catalytic deficiency of fructose 1-phosphate aldolase (aldolase B) in fructose-metabolizing tissues. The A149P mutation in exon 5 of the aldolase B gene, located on chromosome 9q21.3-q22.2, is widespread and the most common HFI mutation, accounting for 57% of HFI chromosomes. The possible origin of this mutation was studied by linkage to polymorphisms within the aldolase B gene. DNA fragments of the aldolase B gene containing the polymorphic marker loci from HFI patients homozygous for the A149P allele were amplified by PCR. Absolute linkage to a common PvuII RFLP allele was observed in 10 A149P homozygotes. In a more informative study, highly heterozygous polymorphisms were detected by direct sequence determination of a PCR-amplified aldolase B gene fragment. Two two-allele, single-base-pair polymorphisms, themselves in absolute linkage disequilibrium, in intron 8 (C at nucleotide 84 and A at nucleotide 105, or T at 84 and G at 105) of the aldolase B gene were identified. Mendelian segregation of these polymorphisms was confirmed in three families. Allele-specific oligonucleotide (ASO) hybridizations with probes for both sequence polymorphisms showed that 47% of 32 unrelated individuals were heterozygous at these loci; the calculated PIC value was .37. Finally, ASO hybridizations of PCR-amplified DNA from 15 HFI patients homozygous for the A149P allele with probes for these sequence polymorphisms revealed absolute linkage disequilibrium between the A149P mutation and the 84T/105G allele. These results are consistent with a single origin of the A149P allele and subsequent spread by genetic drift.     

Nuclear DNA polymorphisms in a wild population of yellow baboons (Papio hamadryas cynocephalus) from Mikumi National Park,Tanzania

Yellow baboons (Papio hamadryas cynocephalus) from Mikumi National Park, Tanzania were studied for polymorphisms in nuclear DNA. The study population consists of four social groups that inhabit overlapping home ranges and exchange males. As a result, these groups are considered to be members of a single interbreeding population. Human DNA clones were used as probes to screen five loci (AT3, REN, HEXB, VIM, and APOB) for restriction fragment length polymorphisms (RFLP). A total of 14 polymorphisms, at least one at each locus, was detected in a panel of 27 baboons tested using six restriction enzymes for each locus. Eleven of these RFLP systems have average heterozygosity values greater than 0.40. This initial screening demonstrates that human DNA clones can be used to detect significant numbers of informative DNA polymorphisms in single-copy nuclear genes of this species and suggests that the average proportion of nucleotides polymorphic across nuclear loci in this population may be between 1.0% and 1.3%. © 1993 Wiley-Liss, Inc.     

Fingerprinting trifoliate orange germ plasm accessions with isozymes, RFLPs, and inter-simple sequence repeat markers

 Trifoliate orange [Poncirus trifoliata (L.) Raf.] is frequently used as a parent in citrus rootstock breeding, but the origin and amount of genetic diversity in germ plasm collections are poorly understood. Most accessions are self-compatible, but produce a mixture of sexual and apomictic seedlings. Variation among 48 vegetatively propagated trifoliate orange accessions was assessed at seven isozyme loci, together with the restriction fragment length polymorphisms (RFLPs) detected by 38 probe-enzyme combinations and the inter-simple sequence repeat (ISSR) markers generated by 11 primers. Isozymes and RFLPs detected few polymorphisms among accessions, although genetic analysis has shown that the common phenotype is heterozygous for four isozyme and at least four RFLP loci. ISSR amplification generated multiple banding profiles with an average of 58 fragments/primer/accession. These fragments were repeatable across DNA samples extracted from different trees of the same accession or extracted at different times, and across separate PCR runs. Seventeen unique marker phenotypes were identified. The 48 trifoliate orange accessions were classified into four major groups based on polymorphic ISSR markers. All large-flowered accessions are in group 4, while small-flowered accessions are in group 3. Many ISSR markers segregated in progeny derived by open-pollination (probably mostly selfing) of a common accession, indicating that these ISSR markers are also heterozygous. Accessions having identical genotypes for a large number of heterozygous markers are unlikely to have diverged by recombination. Thus the limited divergence we detected among most accessions most likely originated by mutation. ‘Monoembryonic’ and ‘Simmons’ differed from other accessions only in the loss of specific markers, indicating that they originated as zygotic seedlings of individuals similar to the common genotype. Three accessions recently introduced from China have relatively different fingerprints with 3–14 unique ISSR markers, and probably represent a much more divergent germ plasm that may be a valuable breeding resource. Received: 8 August 1996 / Accepted: 21 March 1997     

Heterozygosity and fitness in a California population of the labyrinth spider Metepeira ventura (Araneae, Araneidae)

Abstract. The relationship between individual heterozygosity and characteristics likely to be associated with fitness was investigated in the labyrinth spider Metepeira ventura . Adult females and their egg sacs were collected at a coastal site in southern California, and three measures of bodily condition (carapace width, weight, residual index) and six measures of reproductive output (number of egg sacs, variation in egg number among sacs [coefficient of variation], total number of eggs, mean eggs/sac, mean eggs/sac divided by carapace width, mean eggs/sac divided by weight) were determined for each spider. The sample was polymorphic at three allozyme loci that were in Hardy–Weinberg equilibrium, and individual females were heterozygous at up to two of the three loci, forming three heterozygosity classes (0, 1, and 2). None of the bodily condition measures were significantly related to the number of heterozygous loci, while four of the reproductive output estimators (total number of eggs, mean eggs/sac, mean eggs/sac divided by carapace width, mean eggs/sac divided by weight) were significantly influenced by heterozygosity. In each significant case, values for class 2 females were less than those for class 0 and 1 females, whose values were usually more similar. Thus, while female bodily condition was comparable among classes, the most heterozygous females produced fewer total eggs and eggs per sac than their less heterozygous peers. The fact that females of M. ventura engage in a reproductive investment-number trade-off suggests that high-variability and low-variability females may be pursuing distinct reproductive strategies in the wild, with more heterozygous females being K -selected (smaller clutches, heavier eggs) and more homozygous females being r -selected (larger clutches, lighter eggs). Further investigation will be needed to assess more fully the fitness value of heterozygosity in M. ventura .     

Allele-specific PCR can improve the efficiency of experimental resolution of heterozygotes in resequencing studies

Resolution of the two haplotypes present in an individual that is heterozygous at a locus has been a difficult problem for nucleotide sequence-based population genetic studies. Here, we demonstrate a method in which allele-specific polymerase chain reaction (AS-PCR) and computational phasing are combined for relatively high-throughput, efficient resolution of phase in resequencing studies. Using data from multiple loci that were fully experimentally phased, we demonstrate that the popular computational tool PHASE can accurately phase heterozygous individuals with common SNPs (single nucleotide polymorphisms) and/or common haplotypes. However, we also demonstrate that experimental phasing with AS-PCR can efficiently supplement computational phasing, providing a rapid means to phase individuals with rare SNPs or haplotypes and with heterozygous insertion/deletion polymorphisms. By following simple stepwise procedures, AS-PCR can result in much more efficient and accurate experimental phasing of haplotypes than is possible with traditional methods such as cloning.     

Gene-environment interaction effects on the development of immune responses in the 1st year of life

Asthma is a common disease that results from both genetic and environmental risk factors. Children attending day care in the 1st year of life have lower risks for developing asthma, although the mechanism for this "day care" effect is largely unknown. We investigated the interactions between day care exposure in the 1st 6 mo of life and genotypes for 72 polymorphisms at 45 candidate loci and their effects on cytokine response profiles and on the development of atopic phenotypes in the 1st year of life in the Childhood Onset of Asthma (COAST) cohort of children. Six interactions (at four polymorphisms in three loci) with "day care" that had an effect on early-life immune phenotypes were significant at P<.001. The estimated false-discovery rate was 33%, indicating that an estimated four P values correspond to true associations. Moreover, the "day care" effect at some loci was accounted for by the increased number of viral infections among COAST children attending day care, whereas interactions at other loci were independent of the number of viral infections, indicating the presence of additional risk factors associated with day care environment. This study identified significant gene-environment interactions influencing the early patterning of the immune system and the subsequent development of asthma and highlights the importance of considering environmental risk factors in genetic analyses.     

Molecular analysis of the extent of asymmetry in two asymmetric somatic hybrids of tomato

Summary Two somatic hybrid plants generated from a single fusion event between Lycopersicon esculentum and irradiated L. pennellii protoplasts have been analyzed at the molecular level. Over 30 loci have been analyzed using isozymes and RFLPs. All loci tested on chromosomes 2–10 were heterozygous, while those loci on chromosome 12 were homozygous L. pennellii in both somatic hybrids. In one of the somatic hybrids, 2850, loci on chromosome 1 were also homozygous L. pennellii. The other somatic hybrid, 28F5, was heterozygous at all chromosome 1 loci tested, but exhibited altered stoichiometry of parental bands as compared to the sexual hybrid. Loci on chromosome 2 from both somatic hybrids have altered stoichiometry, with L. pennellii alleles being four times more abundant than expected. Both somatic hybrids contain the L. esculentum chloroplast genome, while only L. pennellii polymorphisms have been detected in the mitochondrial genome.     

Human lymphocyte polymorphisms detected by quantitative two-dimensional electrophoresis.

A survey of 186 soluble lymphocyte proteins for genetic polymorphism was carried out utilizing two-dimensional electrophoresis of 14C-labeled phytohemagglutinin (PHA)-stimulated human lymphocyte proteins. Nineteen of these proteins exhibited positional variation consistent with independent genetic polymorphism in a primary sample of 28 individuals. Each of these polymorphisms was characterized by quantitative gene-dosage dependence insofar as the heterozygous phenotype expressed approximately 50% of each allelic gene product as was seen in homozygotes. Patterns observed were also identical in monozygotic twins, replicate samples, and replicate gels. The three expected phenotypes (two homozygotes and a heterozygote) were observed in each of 10 of these polymorphisms while the remaining nine had one of the homozygous classes absent. The presence of the three phenotypes, the demonstration of gene-dosage dependence, and our own and previous pedigree analysis of certain of these polymorphisms supports the genetic basis of these variants. Based on this data, the frequency of polymorphic loci for man is: P = 19/186 = .102, and the average heterozygosity is .024. This estimate is approximately 1/3 to 1/2 the rate of polymorphism previously estimated for man in other studies using one-dimensional electrophoresis of isozyme loci. The newly described polymorphisms and others which should be detectable in larger protein surveys with two-dimensional electrophoresis hold promise as genetic markers of the human genome for use in gene mapping and pedigree analyses.     

Restriction fragment length polymorphism of ovine casein genes: close linkage between the αs1-,αs2-,ß- and k-casein loci

Restriction fragment length polymorphism (RFLP) of ovine casein genes was investigated. Genomic DNA from 56 rams was digested with 10 restriction endonucleases and Southern blots probed with the four ovine casein cDNAs (alpha s1-, beta-, alpha s2- and kappa-Cn). Five enzymes, namely, BglI, PvuII, RsaI, TaqI and HindIII revealed nine different RFLPs. The inheritance of six of these polymorphisms was studied by segregation analysis of gametes in nine rams' families, and each of them could be related to the existence of alleles at the relevant casein locus. A close linkage between the four ovine casein genes was demonstrated since no recombination within the four pairs of loci examined, alpha s1-beta-Cn, alpha s1-kappa-Cn, beta-kappa-Cn and alpha s2-kappa-Cn, was observed in the progeny of double heterozygous rams. The casein genes are thus clustered in the ovine species as in the case of other mammals.     

Patterns of Gene Variation in Central and Marginal Populations of DROSOPHILA ROBUSTA

The central and marginal populations of D. robusta differ greatly in the level of inversion polymorphism; the marginal populations are monomorphic or nearly so and the central populations are highly polymorphic. This paper presents the frequencies of alleles at forty gene loci in various populations of D. robusta, studied by electrophoresis of proteins and enzymes. Population samples were obtained from eight widely separated populations of D. robusta which included the central, the extreme marginal and the intervening populations between the center and the margins. We find that the proportion of polymorphic loci and average heterozygosity per individual is slightly higher in the marginal populations than the central populations. In D. robusta on an average, 39% of the loci are polymorphic and the average proportion of loci heterozygous per individual is 11%. A breakdown of loci in three categories, viz, hydrolytic enzymes and some other enzymes, larval proteins and glycolytic and Kreb's cycle enzymes, shows that in all populations the level of polymorphism is highest in the hydrolytic enzymes, intermediate in larval proteins and least in the glycolytic and Kreb's cycle enzymes. On the average, the proportion of loci heterozygous per individual for three groups of loci is: hydrolytic enzymes and others (.164), larval proteins (.115) and glycolytic and Kreb's cycle enzymes (.037). We also observe that in all populations the level of polymorphism on the X chromosome is far less than the expected 38%; in salivary gland cells the euchromatic length of the X chromosome is 38% of the entire genome. Lower levels of polymorphism for the X chromosome loci are explained due to low probability of balanced polymorphisms for the X-linked loci since the conditions for establishment of balanced polymorphism for X-linked loci are more restrictive than for the autosomal loci.-The polymorphic loci can be grouped according to pattern of allele frequencies in different populations as follows: (1) The allele frequencies are similar in all populations at the XDH, Pep-1 and Hex-1 loci. (2) The alleles at the Est-1, Est-2, Amy loci and the AP-4(1.0) and the LAP-1(.90) alleles show north south clinal change in frequency. (3) There is north south and east west differentiation at the Pt-5, Pt-8 and Pt-9 loci and the allele AP-4(.81). (4) Polymorphism at loci such as Fum, B.Ox, Hex-8, Pep-2 and Pep-3 are restricted to only one or two of the populations. (5) Allele frequencies at the MDH and ODH loci fluctuate between populations. (6) Allele frequencies at many polymorphic loci such as Est-1, Est-2, LAP-1, AP-4, Pt-5, Pt-8, Pt-9, Pt-16, MDH, Fum change clinally within a gene arrangement. The pattern of gene variation in D. robusta is very complex and cannot be easily explained due to migration of neutral alleles between once-isolated populations or to semi-isolation of neutral alleles. The observations of the pattern of allele variation in different populations, high levels of polymorphism in the marginal populations which have small population size and low levels of polymorphism of the X chromosome loci all support the argument in favor of balancing selection as the main mechanism for the maintenance of these polymorphisms. Environmental factors must play a role in the maintenance of a great deal of these polymorphisms, since we observe clinal allele frequency changes even within a given inversion type.     

Genome-wide linkage analysis of inherited hydrocephalus in the H-Tx rat

Inherited hydrocephalus in humans has received very little attention, most probably because known occurrences are sporadic and systematic investigation is difficult. The H-Tx rat, one of a number of rodent strains with inherited hydrocephalus, has a complex inheritance with more than one postulated susceptibility gene and 40% penetrance. The aim of this study was to perform a genome-wide scan on backcross progeny derived from H-Tx and Fisher F344 rats, to identify genomic regions associated with hydrocephalus. Penetrance of hydrocephalus in (H-Tx × F344) F₁× H-Tx was 12.3%. All severely hydrocephalic progeny (n = 185) and a subset of normal progeny (n = 128) were screened with 110 simple sequence length polymorphisms (SSLPs) with 83% coverage of the genome. A significant susceptibility locus was found on chromosome (Chr) 11 (LOD = 3.1). Three loci with suggestive linkage were found on Chr 17 (LOD = 2.4), on Chr 9 (LOD = 1.94), and on Chr 19 (LOD = 1.91). For the loci on Chr 11 and 19, hydrocephalus was associated with the heterozygous genotype, while the other two were recessive. Although none of the four loci was essential for the hydrocephalic phenotype, the additive effects of two, three, or four loci increased the penetrance in a linear fashion. Altogether these four loci accounted for 13.5% of the total variance. It is concluded that hydrocephalus in the H-Tx rat is associated with two, possibly four genetic loci, but that there may be additional undefined genetic and environmental influences. Received: 26 April 2000 / Accepted: 24 August 2000     

Genetic variation of the bovine thyroglobulin gene studied at the DNA level

The bovine thyroglobulin gene has been analysed for variation using restriction endonucleases. Six independent restriction fragment length polymorphisms have been identified. One of these results most probably from a 2.5-kb deletion, the others being compatible with point mutations. We determined that an individual taken at random within the Belgian White and Blue breed is, on average, heterozygous for one out of 1700 nucleotides within the thyroglobulin gene.     

Genetic Variation in Natural Populations of Five Drosophila Species and the Hypothesis of the Selective Neutrality of Protein Polymorphisms

We have studied genetic variation at 30-32 loci coding for enzymes in natural populations of five species of Drosophila. The average proportion of heterozygous loci per individual is 17.7 +/- 0.4%. The average proportion of polymorphic loci per population is 69.2 +/- 2.6% or 49.8 +/- 2.2%, depending on what criterion of polymorphism is used. The following generalizations are advanced: (1) The amount of genetic polymorphism varies considerably from locus to locus. (2) At a given locus, populations of the same species are very similar in the amount and pattern of genetic variation. (3) However, at some loci large differences sometimes occur between local populations of the same species. (4) The amount of variation at a given locus is approximately the same in all five species. (5) When different species are compared, the pattern of the variation is either essentially identical or totally different at a majority of loci. We have tested the hypothesis that protein polymorphisms are selectively neutral by examining four predictions derived from the hypothesis. Our results are at variance with every one of the predictions. We have measured the amount of genetic differentiation, D, between taxa of various degrees of evolutionary divergence. The average value of D is 0.033 for local populations, 0.228 for subspecies, 0.226 for semispecies, 0.538 for sibling species, and 1.214 for morphologically distinguishable species. Our results indicate that a substantial degree of genetic differentiation (22.8 allelic substitutions for every 100 loci) occurs between allopatric populations that have diverged to the point where they might become different species if they were to become sympatric. However, very little additional genetic change is required for the development of complete reproductive isolation. After the speciation process is completed, species continue to diverge genetically from each other.     

Extensive spatial genetic structure revealed by AFLP but not SSR molecular markers in the wind-pollinated tree, Fagus sylvatica

Studies of fine-scale spatial genetic structure (SGS) in wind-pollinated trees have shown that SGS is generally weak and extends over relatively short distances (less than 30-40 m) from individual trees. However, recent simulations have shown that detection of SGS is heavily dependent on both the choice of molecular markers and the strategy used to sample the studied population. Published studies may not always have used sufficient markers and/or individuals for the accurate estimation of SGS. To assess the extent of SGS within a population of the wind-pollinated tree Fagus sylvatica, we genotyped 200 trees at six microsatellite or simple sequence repeat (SSR) loci and 250 amplified fragment length polymorphisms (AFLP) and conducted spatial analyses of pairwise kinship coefficients. We re-sampled our data set over individuals and over loci to determine the effect of reducing the sample size and number of loci used for SGS estimation. We found that SGS estimated from AFLP markers extended nearly four times further than has been estimated before using other molecular markers in this species, indicating a persistent effect of restricted gene flow at small spatial scales. However, our SSR-based estimate was in agreement with other published studies. Spatial genetic structure in F. sylvatica and similar wind-pollinated trees may therefore be substantially larger than has been estimated previously. Although 100-150 AFLP loci and 150-200 individuals appear sufficient for adequately estimating SGS in our analysis, 150-200 individuals and six SSR loci may still be too few to provide a good estimation of SGS in this species.     

Enzyme polymorphisms in the European oyster, Ostrea edulis L.

Extracts of the digestive diverticula of more than 300 individuals from four geographically separated populations of the European oyster, Ostrea edulis L. have been examined by electrophoresis for esterase and phosphoglucose isomerase polymorphisms. Three regions of esterase variability have been detected in the electrophoretograms. It is inferred that the variants in two of these regions are governed by two co-dominant alleles at each of two loci. Two phenotypes of phosphoglucose isomerase, which is inferred to be a dimeric enzyme encoded by two alleles at a single locus, have been observed. Allele frequencies and phenotype distributions are compared in the four populations. It is concluded that the populations differ genetically at the loci investigated, the magnitude of the differences being generally greatest between different 'physiological races'.     

Population genetics of short tandem repeat (STR) loci

To investigate the population genetics of short tandem repeat (STR) polymorphisms in human populations, we have studied the allele frequency distributions of four STR loci (HUMTH01, HUMVWA31, HUMF13A1 and HUMFES) in 16 different population surveys which can be categorised within three broadly defined ethnic groups: Caucasian, Asian (Indian subcontinent), and African (Afro-Caribbean and US black). We have observed that allele frequency distributions of populations within ethnic groups are similar; consequently, genetic distances are an order of magnitude lower than between ethnic groups. Inbreeding coefficients (F-statistics) and calculations of the number of mean heterozygous loci per individual, along with estimates of variance, did not suggest that the populations were substructured. This included a study of an immigrant Asian population known to comprise at least three different sub-groups. Finally, an indication of the discriminating power is given by calculation of likelihood ratios (LR) of each individual tested across all four loci. Approximately 70% of Caucasians give an LR of greater than 10,000; the test is even more discriminating in Afro-Caribbeans-approximately 90% of tests are greater than 10,000.Editor's commentsThe authors present data generated by the move from VNTR to STR loci for human identification. The data they present for samples within major racial groupings address some of the concerns about population substructuring discussed by Balding and Nichols in this volume.     

The strength of selection on ultraconserved elements in the human genome

Ultraconserved elements are stretches of consecutive nucleotides that are perfectly conserved in multiple mammalian genomes. Although these sequences are identical in the reference human, mouse, and rat genomes, we identified numerous polymorphisms within these regions in the human population. To determine whether polymorphisms in ultraconserved elements affect fitness, we genotyped unrelated human DNA samples at loci within these sequences. For all single-nucleotide polymorphisms tested in ultraconserved regions, individuals homozygous for derived alleles (alleles that differ from the rodent reference genomes) were present, viable, and healthy. The distribution of allele frequencies in these samples argues against strong, ongoing selection as the force maintaining the conservation of these sequences. We then used two methods to determine the minimum level of selection required to generate these sequences. Despite the lack of fixed differences in these sequences between humans and rodents, the average level of selection on ultraconserved elements is less than that on essential genes. The strength of selection associated with ultraconserved elements suggests that mutations in these regions may have subtle phenotypic consequences that are not easily detected in the laboratory.     

Computation vs. cloning: evaluation of two methods for haplotype determination

Nuclear sequence data, often from multiple loci, are increasingly being employed in analyses of population structure and history, yet there has been relatively little evaluation of methods for accurately and efficiently separating the alleles or haplotypes in heterozygous individuals. We compared the performance of a computational method of haplotype reconstruction and standard cloning methods using a highly variable intron (ornithine decarboxylase, intron 6) in three closely related species of dabbling ducks (genus Anas). Cloned sequences from 32 individuals were compared to results obtained from phase 2.1.1 . phase correctly identified haplotypes in 28 of 30 heterozygous individuals when the underlying model assumed no recombination. Haplotypes of the remaining two individuals were also inferred correctly except for unique polymorphisms, the phase of which was appropriately indicated as uncertain (phase probability = 0.5). For a larger set of 232 individuals, results were essentially identical regardless of the recombination model used and haplotypes for all 30 of the tested heterozygotes were correctly inferred, with the exception of uncertain phase for unique polymorphisms in one individual. In contrast, initial sequences of one clone per sample yielded accurate haplotype determination in only 26 of 30 individuals; polymerase chain reaction (PCR)/cloning errors resulting from misincorporation of individual nucleotides could be recognized and avoided by comparison to direct sequences, but errors due to PCR recombination resulted in incorrect haplotype reconstruction in four individuals. The accuracy of haplotypes reconstructed by phase , even when dealing with a relatively small number of samples and numerous variable sites, suggests broad utility of computational approaches for reducing the cost and improving the efficiency of data collection from nuclear sequence loci.