Linkage Disequilibrium Between Allozymes in Natural Populations of Lodgepole Pine

Multilocus structure in thePinus contorta – Pinus banksianacomplex

Natural hybridization between lodgepole pine ( Pinus contorta Dougl. ex Loud.) and jack pine ( Pinus banksiana Lamb.) in western Canada provides a distinct opportunity to infer evolutionary and demographic determinants of nonrandom associations of the alleles at different loci occurring in this Pinus contorta – Pinus banksiana complex (PCBC). Here, we investigated multilocus associations among and within 40 PCBC populations sampled from central and northwestern Alberta, using 39 random amplified polymorphic DNA (RAPD) markers. For each of the 40 PCBC populations, we examined distributions of 741 pairs of linkage (gametic) disequilibria (LD) between the 39 RAPDs and measured the “average” multilocus associations using the summary statistics that allow for packaging of individual LD in each population. We then partitioned the variance of LD in the total population to assess the causes of multilocus population structure. The results showed that (i) LD were more prominent in hybrid populations than in parental populations; (ii) multilocus Wahlund effect was a much more important determinant of population structure than its single-locus counterpart, particularly in hybrid populations; and (iii) considerable mutlilocus associations across the populations within each taxa group was due to the presence of different multilocus haplotypes in different populations. Such results are best explained by the fact that PCBC populations are geographically and ecologically marginal and are produced at the balance between mixing of two distinct gene pools creating new recombinants and selection in favor of parental gametes, but against the recombinants.

Highly variable SSR markers in Douglas-fir: Mendelian inheritance and map locations

Twenty-two highly variable SSR markers were developed in Douglas-fir [ Pseudotsuga menziesii (Mirb.) Franco] from five SSR-enriched genomic libraries. Fifteen PCR primer pairs amplified a single codominant locus, while seven primer pairs occasionally amplified two loci. The Mendelian inheritance of all 22 SSRs was confirmed via segregation analyses in several Douglas-fir families. The mean observed heterozygosity and the mean number of alleles per locus were 0.855 (SE=0.020) and 23 (SE=1.6), respectively. Twenty markers were used in genetic linkage analysis and mapped to ten known linkage groups. Because of their high polymorphism and unambiguous phenotypes, 15 single-locus markers were selected as the most suitable for DNA fingerprinting and parentage analysis. Only three SSRs were sufficient to achieve an average probability of exclusion from paternity of 0.998 in a Douglas-fir seed orchard block consisting of 59 parents.

Complex patterns of mating revealed in a Eucalyptus regnans seed orchard using allozyme markers and the neighbourhood model

The neighbourhood model apportions offspring of individual mother plants to self-fertilization, outcrossing to males within a circumscribed area around the mother plant (the neighbourhood), and outcrossing to males outside the neighbourhood. Formerly the model was applied only to haploid pollen gametes in the offspring of conifers, but is extended so that it can be used with genotypic data from diploid offspring of both angiosperms and gymnosperms. In addition, it is shown that the mating parameters can be estimated without independent estimates of allele frequencies in the pollen pools outside the neighbourhood; thus the model might be applied effectively to natural populations exposed to unknown external pollen sources. Parameters of the neighbourhood mating model were estimated for a 10-year-old seed orchard population of the insect-pollinated tree, Eucalyptus regnans, in southeast Australia, which contained a mixture of two geographical provenances (Victoria and Tasmania). The mating patterns revealed were complex. Crosses between trees of the same provenance occurred three times more often than crosses between trees of different provenances. Levels of self-fertilization and patterns of mating within neighbourhoods were influenced by provenance origin, crop fecundity and orchard position (central vs. edge) of mother trees. Gene dispersal, however, was extensive, with approximately 50% of effective pollen gametes coming from males more than 40 m away from mother trees (average distance between neighbouring trees was 7.4 m). Thus, insect pollinators are efficient promoters of cross-fertilization in this orchard, with the result that the effective number of males mating with each female is large.

Joint linkage and linkage disequilibrium mapping in natural populations

A new strategy for studying the genome structure and organization of natural populations is proposed on the basis of a combined analysis of linkage and linkage disequilibrium using known polymorphic markers. This strategy exploits a random sample drawn from a panmictic natural population and the open-pollinated progeny of the sample. It is established on the principle of gene transmission from the parental to progeny generation during which the linkage between different markers is broken down due to meiotic recombination. The strategy has power to simultaneously capture the information about the linkage of the markers (as measured by recombination fraction) and the degree of their linkage disequilibrium created at a historic time. Simulation studies indicate that the statistical method implemented by the Fisher-scoring algorithm can provide accurate and precise estimates for the allele frequencies, recombination fractions, and linkage disequilibria between different markers. The strategy has great implications for constructing a dense linkage disequilibrium map that can facilitate the identification and positional cloning of the genes underlying both simple and complex traits.

Simultaneous estimation of pollen contamination and pollen fertilities of individual trees in conifer seed orchards using multilocus genetic data

Seed orchards of forest trees are established to produce genetically-improved seed for reforestation. Genetic efficiency requires seed orchards to be (1) reproductively isolated form surrounding trees, (2) that there be similar fertilities among all orchard trees, and (3) minimum inbreeding. Each aspect of seed orchard reproduction can be simultaneously estimated using the observed frequency of each multilocus gametic type contributed through fertilizing pollen and the expected multilocus gametic segregation frequencies of orchard tree and of the contaminating population. These genetic estimates are directly relevant to the genetics of the tree breeding program. The flexibility of sampling seed - the basic data for these techniques - allows great scope for hypothesis testing, including tests of the accuracy of predictions of biophysical models of pollen movement. A simple example and a white spruce seed orchard case study are presented to illustrate the estimation technique and to investigate its sensitivity.

A genetic linkage map of Picea abies Karst., based on RAPD markers, as a tool in population genetics

Norway spruce (Picea abies Karst.) is a most important species among European forest trees for both economical and ecological reasons. However, this species has suffered from a lack of information on the genetic side due to the scarcity of linkage data. In this study we have used a population of 72 megagametophytes from a single tree in a natural Italian stand to produce a genetic linkage map by means of RAPD markers. Ninety-six random decamers used as primers yielded 185 polymorphic loci showing Mendelian inheritance. Analysis of the segregation by multipoint analysis allowed us to define 17 major linkage groups covering a total distance of 3584 cM, with an average spacing between markers of 22 cM. Possible uses of a genetic linkage map with respect to population ecology and genetics are discussed.

Multilocus structure in Pinus contorta Dougl

We studied isozyme variation at 21 loci in 66 populations from three subspecies of Pinus contorta Dougl.; 35 in spp. latifolia, 20 in spp contorta and 11 in spp. murrayana. The objectives were to assess gametic disequilibria and multilocus structure. There was considerable differentiation of allele frequencies at 19 polymorphic loci across the 66 populations and within the subspecies. Allele frequencies at many loci correlated with geographic variables. Genetic variability varied considerably among populations within subspecies but the subspecies means were similar. The mean number of polymorphic loci and the mean heterozygosity over 19 polymorphic loci were, respectively, 13 and 0.194 in latifolia, 12 and 0.196 in murrayana, and 12 and 0.180 in contorta. The mean heterozygosity correlated with longitude and altitude across the 66 populations and with latitude in latifolia. Gametic disequilibria were evident in 40 populations; 29 in latifolia, eight in murrayana and three in contorta. Gametic disequilibria correlated with latitude across the 66 populations and with longitude in latifolia. The single-locus F ST averaged 0.0339 in latifolia, 0.0567 in murrayana, and 0.0764 in contorta. The multilocus F STM was 0.1227 in latifolia, 0.2926 in murrayana, and 0.3328 in contorta. Multilocus Wahlund and founder effects, migration patterns, and natural selection, probably played significant roles in generating and maintaining the multilocus genetic structure in P. contorta in general and the subspecies latifolia in particular.

Associations between nuclear loci and chloroplast DNA genotypes in wild barley

Associations among alleles at nine nuclear loci and three chloroplast DNA (cpDNA) genotypes were assessed in a sample of 247 accessions of the wild barley, Hordeum vulgare ssp. spontaneum. Alleles at two of the nine nuclear loci are marked by length variations in the intergenic spacer region of ribosomal DNA (rDNA), and those of the other seven loci are well characterized allozymes. The three chloroplast DNA (cpDNA) genotypes are marked by restriction fragment length polymorphisms resulting from three polymorphic restriction sites detected by Southern blot hybridization. The analyses were performed by dividing the nine nuclear loci into a series of two-locus subsets and constructing log-linear models to characterize associations between the subsets of two nuclear loci and the cpDNA genotypes. Statistically significant associations were detected between six of the nine nuclear loci and the cpDNA genotypes, either individually as pairwise correlations, or through interaction with another nuclear locus to form three-variate complexes. Although the sample size of the present study was inadequate for statistical evaluation of higher order interactions, the results suggest the existence of interactions in which more than two nuclear loci are involved in associations with cpDNA genotypes. The observed cytonuclear associations appear to result from interplay among a number of evolutionary forces including a mating system of predominant selfing, differentiation among gene pools of local populations, and adaptation of barley genotypes to specific environmental conditions.

Spatial autocorrelation analysis of the distribution of genotypes within populations of lodgepole pine

Spatial autocorrelation analyses of point samples within two populations of lodgepole pine (Pinus contorta ssp. latifolia) indicate that single-locus mature tree and pollen genotypes are distributed in a nearly random fashion for most of the allozyme loci assayed. This lack of structure in the distributions of most genotypes is consistent with outcrossing rates that are very nearly 1.0 and with estimates indicating that both pollen and seed are dispersed over long distances in lodgepole pine. However, spatial autocorrelation of genotypes for a few loci suggests that genotypes at these loci may be under natural selection.

Linkage disequilibrium in natural and experimental populations of Drosophila melanogaster

We have studied linkage disequilibrium in Drosophila melanogaster in two samples from a wild population and in four large laboratory populations derived from the wild samples. We have assayed four polymorphic enzyme loci, fairly closely linked in the third chromosome: Sod Est-6, Pgm, and Odh. The assay method used allows us to identify the allele associations separately in each of the two homologous chromosomes from each male sampled. We have detected significant linkage disequilibrium between two loci in 16.7% of the cases in the wild samples and in 27.8% of the cases in the experimental populations, considerably more than would be expected by chance alone. We have also found three-locus disequilibria in more instances than would be expected by chance. Some disequilibria present in the wild samples disappear in the experimental populations derived from them, but new ones appear over the generations. The effective population sizes required to generate the observed disequilibria by randomness range from 40 to more than 60,000 individuals in the natural population, depending on which locus pair is considered, and from 100 to more than 60,000 in the experimental populations. These population sizes are unrealistic; the fact that different locus-pairs yield disparate estimates within the same population argues against the likelihood that the disequilibria may have arisen as a consequence of population bottlenecks. Migration, or population mixing, cannot be excluded as the process generating the disequilibria in the wild samples, but can in the experimental populations. We conclude that linkage disequilibrium in these populations is most likely due to natural selection acting on the allozymes, or on loci very tightly linked to them.

Adaptation at specific loci. V. Metabolically adjacent enzyme loci may have very distinct experiences of selective pressures

The polymorphic phosphoglucomutase (PGM) and glucose-6-phosphate dehydrogenase (G6PD) loci have been studied in parallel to experimental work on the phosphoglucose isomerase (PGI) polymorphism in Colias butterflies. PGI, PGM and G6PD are also autosomal in Colias. PGM and G6PD are loosely linked (and represent the first identified autosomal linkage group in Colias); they assort independently from PGI. Recombination occurs in both sexes. Neither PGM nor G6PD shows large, consistent differences in flight capacity through the day among its genotypes, as PGI does. PGM shows some change of allele frequencies, and match to Hardy-Weinberg expectation, with air temperature in middle and latter parts of the season, but not early in the season. G6PD may show some heterozygote excess over Hardy-Weinberg expectation early in the day, but more testing is needed. No evidence for differential survivorship was seen at PGM or G6PD, in contrast to PGI. At the PGM and G6PD loci, male heterozygotes are advantaged in mating with females, but without the evidence of female choice which occurs for PGI. These effects are not correlated among the three loci. There is no assortative mating at G6PD (nor at PGI). There is minor positive assortative mating of PGM heterozygotes, but it is too weak to account for the PGM-genotype-specific male mating advantage. No trends of multilocus genotype frequencies involving PGI are seen. Certain PGM-G6PD two-locus genotypes are over-represented, and others under-represented, in wild adult samples, particularly among males and uniformly among successfully mating males. Our results emphasize that enzyme loci sharing a substrate need not have common experience of the existence or strength of natural selection, and suggest initial food-resource processing and allocation as a possible context for fitness-related effects of the PGM and G6PD polymorphisms.