The effect of single nucleotide polymorphism identification strategies on estimates of linkage disequilibrium

At present there is tremendous interest in characterizing the magnitude and distribution of linkage disequilibrium (LD) throughout the human genome, which will provide the necessary foundation for genome-wide LD analyses and facilitate detailed evolutionary studies. To this end, a human high-density single-nucleotide polymorphism (SNP) marker map has been constructed. Many of the SNPs on this map, however, were identified by sampling a small number of chromosomes from a single population, and inferences drawn from studies using such SNPs may be influenced by ascertainment bias (AB). Through extensive simulations, we have found that AB is a potentially significant problem in estimating and comparing LD within and between populations. Specifically, the magnitude of AB is a function of the SNP discovery strategy, number of chromosomes used for SNP discovery, population genetic characteristics of the particular genomic region considered, amount of gene flow between populations, and demographic history of the populations. We demonstrate that a balanced SNP discovery strategy (where equal numbers of chromosomes are sampled from multiple subpopulations) is the optimal study design for generating broadly applicable SNP resources. Finally, we validate our theoretical predictions by comparing our results to publicly available data from ten genes sequenced in 24 African American and 23 European American individuals.     

Perl and Python codes to estimate Hardy–Weinberg‐ and linkage‐disequilibrium‐related parameters

We have developed two pieces of Perl and Python codes to calculate the natural log of lambda, a quantity that expresses deviations from Hardy–Weinberg equilibrium, and R, a weighted correlation coefficient between alleles of locus pairs as a quantity that expresses linkage disequilibrium. These quantities can be used in an estimation approach as tools to facilitate the characterization and the comparison of populations.     

High levels of linkage disequilibria between serologically defined class I bovine lymphocyte antigens (BOLA-A) and class II DQB restriction fragment length polymorphism (RFLP) in Norwegian cows

Serum defined BoLA-A antigens, together with BoLA-DQB RFLP patterns, were determined in 87 almost unrelated Norwegian cattle. Statistical analysis revealed strong linkage disequilibria between these loci at the population level. A total of 13 haplotypes were found to be present at frequencies significantly greater than those predicted on the basis of their component gene frequencies. Among these, the subgroups 1A and 1B of the DQ1 haplotype were found to be closely associated with the class I antigens A11 and w16, respectively. The association between A11 and DQ1A is of particular interest, as two independent studies, one employing class I serology, and the other RFLP analysis of the class II locus DQ, have previously indicated that A11 and DQ1A confer relative susceptibility to mastitis.     

Antigen and haplotype frequencies at three human leucocyte antigen loci (HLA-A, -B, -C) in the Pawaia of Papua New Guinea

The genetic profile of the Pawaia, a seminomadic, linguistic isolate from the highlands fringe of Papua New Guinea, is described in terms of antigen and haplotype frequencies at three class I human leucocyte antigen loci (HLA-A, -B, and -C). The Pawaia, like other Papua New Guinea populations, exhibit restricted polymorphisms at all three loci studied, both in the number of alleles segregating and in the level of average heterozygosity. An extremely high frequency (52.9%) of HLA-B27, the antigen implicated in the pathogenesis of seronegative arthropathies, was found. A novel HLA-C locus specificity, CNG, resulting probably from a gene duplication event, was also observed in significant numbers. Although the gene frequency comparisons suggest their strong affinities with the highlanders, the Pawaia haplotypes reveal significant admixture from other neighbouring groups as well. The usefulness of HLA haplotypes in tracing the movements of human populations in the New Guinea area is discussed.     

Forest dynamics,differential mortality and variable recruitment probabilities

Abstract. Both spatial and temporal variability in recruitment probabilities can lead to coexistence in gap-phase regenerating forests which would otherwise tend to be dominated by fewer species. Using modified Markov models, the potential roles were examined of temporal variability and differential mortality rates among species in the dynamics of a forest for which spatial variability has been rejected as a strong factor leading to coexistence. Differential longevity modifies results obtained from a simple Markov model: it exerts a strong influence on the equilibrium species composition, on the rate of community change and on the time a community requires to reach equilibrium. Simulations with varying transition probabilities mimicked a changing climate, producing four main results: 1. Unless the duration of climate states is very long or very short, forest composition is in a continual state of disequilibrium. 2. Species vary in their response times to changing climate. 3. The mean abundance of each species under a varying climate scenario is different from that expected from the mean climate state. 4. The rare, long-lived species was favored by climatic fluctuations at the expense of more common shorter lived species. Differential mortality rates provide an equilibrium-based mechanism for coexistence, and temporally fluctuating recruitment probabilities a non-equilibrium mechanism. Composition could be maintained by differential longevity among species and climatic fluctuations allowing periodic recruitment of the less common species.     

A hierarchical analysis of population genetic structure in Rhizobium leguminosarum bv. trifolii

Little is known about the population processes that shape the genetic diversity in natural populations of rhizobia. A sample of 912 Rhizobium leguminosarum biovar trifolii isolates were collected from naturalized red clover populations ( Trifolium pratense ) and analyzed for 15 allozyme loci to determine the levels and distribution of genetic diversity. Hierarchical analyses compared different sampling levels, geographical separation, and temporal separation. Total genetic diversity across all isolates was H = 0.426, with 57.6% of the total diversity found among isolates obtained from individual red clover plants. Relatively low genetic differentiation among populations and high differentiation among plants within populations was observed; this suggests that gene flow and founder effect act differently at geographical and local scales. Significant differences were observed in (i) allele frequencies among populations and among plants within populations, and (ii) the frequency distribution of the most widespread and the most abundant strains. When multilocus linkage disequilibrium was calculated, significant levels of disequilibrium were observed in the total sample and in three of the eight populations.     

Molecular diversity and multilocus organization of the parental lines used in the International Rice Molecular Breeding Program

One hundred and ninety three parental lines obtained from 26 countries for an international rice molecular breeding program were evaluated using 101 well-distributed simple sequence repeat (SSR) markers. An overall genetic diversity of 0.68 and an average of 6.3 alleles per locus were revealed, indicating a high level of genetic variation in these lines. Cluster analysis of the 193 accessions showed three major groups and nine subgroups. Group I corresponded to the classical indica subspecies, whereas groups II and III belong to the japonica subspecies. Indica and japonica differentiation accounted for only 6.5% of the total variation in the entire sample and 93.5% was due to within-subspecies diversity. Differentiation among eco-geographic regions accounted for 24% of the diversity within the subspecies. Larger amounts of the eco-geographical differentiation were resolved within japonica than within indica. The largest indica-japonica differentiation based on the single locus level was detected by markers on chromosomes 9 and 12, while the smallest differentiation was detected by markers on chromosomes 4 and 8. Furthermore, genetic differences at the single-locus and two-locus levels, as well as components due to allelic and gametic differentiation, were revealed between indica and japonica and among the main geographic regions. The multilocus analysis in genetic diversity showed a higher proportion of variation caused by predominant non-random associations of different loci within and among the classified subspecies and geographic subdivisions. The results suggest that selection for eco-geographical adaptation on multilocus associations was largely responsible for the maintenance of the extensive variation in the primary gene pool of rice.Communicated by Q. ZhangS.B. Yu and W.J. Xu contributed equally to the work