Definition and properties of disequilibria within nuclear-mitochondrial-chloroplast and other nuclear-dicytoplasmic systems

Genetic Incompatibilities Between Mitochondria and Nuclear Genes: Effect on Gene Flow and Speciation

The process of speciation is, according to the biological species concept, the reduction in gene flow between genetically diverging populations. Most of the previous theoretical studies analyzed the effect of nuclear genetic incompatibilities on gene flow. There is, however, an increasing number of empirical examples suggesting that cytoplasmically inherited genetic elements play an important role in speciation. Here, we present a theoretical analysis of mitochondrial driven speciation, in which genetic incompatibilities occur between mitochondrial haplotypes and nuclear alleles. Four population genetic models with mainland-island structure were analyzed that differ with respect to the type of incompatibility and the underlying genetics. Gene flow reduction was measured on selectively neutral alleles of an unlinked locus and quantified by the effective migration rate. Analytical formulae for the different scenarios were derived using the fitness graph method. For the models with haploid genetics, we found that mito-nuclear incompatibilities (MtNI) are as strong as nuclear-nuclear incompatibilities (NNI) in reducing gene flow at the unlinked locus, but only if males and females migrate in equal number. For models with diploid genetics, we found that MtNI reduce gene flow stronger than NNI when incompatibilities are recessive, but weaker when they are dominant. For both haploid and diploid MtNI, we found that gene flow reduction is stronger if females are the migrating sex, but weaker than NNI when males are the migrating sex. These results encourage further examination on the role of mitochondria on genetic divergence and speciation and point toward specific factors (e.g., migrating sex) that could be the focus of an empirical test.

Cryptic genetic diversity in the mottled rabbitfish Siganus fuscescens with mitochondrial introgression at a contact zone in the South China Sea

The taxonomy of the mottled rabbitfish Siganus fuscescens species complex has long been challenging. In this study, we analyzed microsatellite genotypes, mitochondrial lineages, and morphometric data from 373 S. fuscescens individuals sampled from the northern Philippines and Hong Kong (South China Sea, Philippine Sea and Sulu Sea basins), to examine putative species boundaries in samples comprising three co-occurring mitochondrial lineages previously reported to characterize S. fuscescens (Clade A and Clade B) or S. canaliculatus (Clade C). We report the existence of two cryptic species within S. fuscescens in the northeast region of the South China Sea and northern Philippine Sea, supported by genetic and morphological differences. Individual-based assignment methods recovered concordant groupings of individuals into two nuclear genotype clusters (Cluster 1, Cluster 2) with (1) limited gene flow, if any, between them (FST = 0.241; P < 0.001); (2) low frequency of later-generation hybrids; (3) significant association with mitochondrial Clade A and Clade B, respectively; and (4) subtle yet significant body shape differences as inferred from geometric morphometric analysis. The divergence between mitochondrial Clade C and the two other clades was not matched by genetic differences at microsatellite marker loci. The occurrence of discordant mitonuclear combinations (20.5% of the total number of individuals) is thought to result from mitochondrial introgression, consistent with a scenario of demographic, and presumably spatial, post-Pleistocene expansion of populations from northern regions into a secondary contact zone in the South China Sea. Mitonuclear discordance due to introgression obscures phylogenetic relationships for recently-diverged lineages, and cautions against the use of mitochondrial markers alone for species identification within the mottled rabbitfish species complex in the South China Sea region.

Spatial genetic structure in wild cardoon, the ancestor of cultivated globe artichoke: Limited gene flow, fragmentation and population history

Nuclear and chloroplast markers and phenotypic characters were integrated to analyse the population genetic structure of wild cardoon, Cynara cardunculus var. sylvestris, the ancestor of cultivated globe artichoke, Cynara cardunculus var. scolymus on the island of Sardinia, Italy. The spatial scale ranged from a few metres to ∼200km. Wild cardoon appears to be genetically fragmented, with significant genetic divergence at various scales, indicating that gene flow is insufficient to counterbalance the effects of genetic drift or founder effects. Divergence between populations was higher for chloroplast (40%) than for nuclear markers (15%), suggesting that gene flow via seed was lower than via pollen. Two main genetic groups were detected; these correlated with differences in flowering time, capitula size, glossiness, and anthocyanin pigmentation. A complex population structure of wild cardoon emerged over small spatial scales, likely resulting from the interplay between gene dispersal, colonisation history and selective forces. Indeed, Sardinia appears to be a 'hybrid zone' of different gene pools. The island has unique diverse germplasm that has originated from hybridisation among different gene pools. The sampling of seeds from a few plants but from many sites is suggested as the best strategy to harvest the genetic diversity of wild cardoon.

Biogeography of "Cyprinella lutrensis": intensive genetic sampling from the Pecos River 'melting pot' reveals a dynamic history and phylogenetic complexity

Thorough sampling is necessary to delineate lineage diversity for polytypic "species" such as Cyprinella lutrensis. We conducted extensive mtDNA sampling (cytochrome b and ND4) from the Pecos River, Rio Grande, and South Canadian River, New Mexico. Our study emphasized the Pecos River due to its complex geological history and potential to harbor multiple lineages. We used geometric-morphometric, morphometric, and meristic analyses to test for phenotypic divergence and combined nucDNA with mtDNA to test for cytonuclear disequilibrium and combined our sequences with published data to conduct a phylogenetic re-assessment of the entire C. lutrensis clade. We detected five co-occurring mtDNA lineages in the Pecos River, but no evidence for cytonuclear disequilibrium or phenotypic divergence. Recognized species were interspersed amongst divergent lineages of "C. lutrensis". Allopatric divergence among drainages isolated in the Late Miocene and Pliocene apparently produced several recognized species and major divisions within "C. lutrensis". Pleistocene re-expansion and subsequent re-fragmentation of a centralized lineage founded younger, divergent lineages throughout the Rio Grande basin and Edwards Plateau. There is also evidence of recent introductions to the Rio Grande, Pecos and South Canadian Rivers. Nonetheless, deeply divergent lineages have coexisted since the Pleistocene.

Association between Chloroplast and Mitochondrial DNA sequences in Chinese Prunus genotypes (Prunus persica, Prunus domestica, and Prunus avium)

BACKGROUND

The nuclear DNA is conventionally used to assess the diversity and relatedness among different species, but variations at the DNA genome level has also been used to study the relationship among different organisms. In most species, mitochondrial and chloroplast genomes are inherited maternally; therefore it is anticipated that organelle DNA remains completely associated. Many research studies were conducted simultaneously on organelle genome. The objectives of this study was to analyze the genetic relationship between chloroplast and mitochondrial DNA in three Chinese Prunus genotypes viz., Prunus persica, Prunus domestica, and Prunus avium.

RESULTS

We investigated the genetic diversity of Prunus genotypes using simple sequence repeat (SSR) markers relevant to the chloroplast and mitochondria. Most of the genotypes were genetically similar as revealed by phylogenetic analysis. The Y2 Wu Xing (Cherry) and L2 Hong Xin Li (Plum) genotypes have a high similarity index (0.89), followed by Zi Ye Li (0.85), whereas; L1 Tai Yang Li (plum) has the lowest genetic similarity (0.35). In case of cpSSR, Hong Tao (Peach) and L1 Tai Yang Li (Plum) genotypes demonstrated similarity index of 0.85 and Huang Tao has the lowest similarity index of 0.50. The mtSSR nucleotide sequence analysis revealed that each genotype has similar amplicon length (509 bp) except M5Y1 i.e., 505 bp with CCB256 primer; while in case of NAD6 primer, all genotypes showed different sizes. The MEHO (Peach), MEY1 (Cherry), MEL2 (Plum) and MEL1 (Plum) have 586 bps; while MEY2 (Cherry), MEZI (Plum) and MEHU (Peach) have 585, 584 and 566 bp, respectively. The CCB256 primer showed highly conserved sequences and minute single polymorphic nucleotides with no deletion or mutation. The cpSSR (ARCP511) microsatellites showed the harmonious amplicon length. The CZI (Plum), CHO (Peach) and CL1 (Plum) showed 182 bp; whileCHU (Peach), CY2 (Cherry), CL2 (Plum) and CY1 (Cherry) showed 181 bp amplicon lengths.

CONCLUSIONS

These results demonstrated high conservation in chloroplast and mitochondrial genome among Prunus species during the evolutionary process. These findings are valuable to study the organelle DNA diversity in different species and genotypes of Prunus to provide in depth insight in to the mitochondrial and chloroplast genomes.

Association between chloroplast DNA and mitochondrial DNA haplotypes in Prunus spinosa L. (Rosaceae) populations across Europe

Chloroplast DNA (cpDNA) and mitochondrial DNA (mtDNA) were studied in 24 populations of Prunus spinosa sampled across Europe. The cpDNA and mtDNA fragments were amplified using universal primers and subsequently digested with restriction enzymes to obtain the polymorphisms. Combinations of all the polymorphisms resulted in 33 cpDNA haplotypes and two mtDNA haplotypes. Strict association between the cpDNA haplotypes and the mtDNA haplotypes was detected in most cases, indicating conjoint inheritance of the two genomes. The most frequent and abundant cpDNA haplotype (C20; frequency, 51 %) is always associated with the more frequent and abundant mtDNA haplotype (M1; frequency, 84 %). All but two of the cpDNA haplotypes associated with the less frequent mtDNA haplotype (M2) are private haplotypes. These private haplotypes are phylogenetically related but geographically unrelated. They form a separate cluster on the minimum-length spanning tree.

Comparing relative rates of pollen and seed gene flow in the island model using nuclear and organelle measures of population structure

We describe a method for comparing nuclear and organelle population differentiation (F(ST)) in seed plants to test the hypothesis that pollen and seed gene flow rates are equal. Wright's infinite island model is used, with arbitrary levels of self-fertilization and biparental organelle inheritance. The comparison can also be applied to gene flow in animals. Since effective population sizes are smaller for organelle genomes than for nuclear genomes and organelles are often uniparentally inherited, organelle F(ST) is expected to be higher at equilibrium than nuclear F(ST) even if pollen and seed gene flow rates are equal. To reject the null hypothesis of equal seed and pollen gene flow rates, nuclear and organelle F(ST)'s must differ significantly from their expected values under this hypothesis. Finite island model simulations indicate that infinite island model expectations are not greatly biased by finite numbers of populations (>/=100 subpopulations). The power to distinguish dissimilar rates of pollen and seed gene flow depends on confidence intervals for fixation index estimates, which shrink as more subpopulations and loci are sampled. Using data from the tropical tree Corythophora alta, we rejected the null hypothesis that seed and pollen gene flow rates are equal but cannot reject the alternative hypothesis that pollen gene flow is 200 times greater than seed gene flow.

Mitochondrial DNA diversity, population structure, and gender association in the gynodioecious plant Silene vulgaris

A highly variable mitochondrial DNA (mtDNA) restriction fragment length polymorphism (RFLP) locus is used to assess the population structure of mitochondrial genomes in the gynodioecious plant Silene vulgaris at two spatial scales. Thirteen mtDNA haplotypes were identified within 250 individuals from 18 populations in a 20-km diameter region of western Virginia. The population structure of these mtDNA haplotypes was estimated as thetaST = 0.574 (+/- 0.066 SE) and, surprisingly, genetic differentiation among populations was negatively correlated with geographic distance (Mantel r = -0.246, P < 0.002). Additionally, mtDNA haplotypes were spatially clumped at the scale of meters within one population. Gender in S. vulgaris is determined by an interaction between autosomal male fertility restorers and cytoplasmic male sterility (CMS) factors, and seed fitness is affected by an interaction between gender and population sex ratio; thus, selection acting on gender could influence the distribution of mtDNA RFLP haplotypes. The sex ratio (females:hermaphrodites) varied among mtDNA haplotypes across the entire metapopulation, possibly because the haplotypes were in linkage disequilibrium with different CMS factors. The gender associated with some of the most common haplotypes varied among populations, suggesting that there is also population structure in male fertility restorer genes. In comparison with reports of mtDNA variation from other published studies, we found that S. vulgaris exhibits a large number of mtDNA haplotypes relative to that observed in other species.

Phylogeography of Kandelia candel in East Asiatic mangroves based on nucleotide variation of chloroplast and mitochondrial DNAs

Vivipary with precocious seedlings in mangrove plants was thought to be a hindrance to long-range dispersal. To examine the extent of seedling dispersal across oceans, we investigated the phylogeny and genetic structure among East Asiatic populations of Kandelia candel based on organelle DNAs. In total, three, 28 and seven haplotypes of the chloroplast DNA (cpDNA) atpB-rbcL spacer, cpDNA trnL-trnF spacer, and mitochondrial DNA (mtDNA) internal transcribed spacer (ITS) were identified, respectively, from 202 individuals. Three data sets suggested consistent phylogenies recovering two differentiated lineages corresponding to geographical regions, i.e. northern South-China-Sea + East-China-Sea region and southern South-China-Sea region (Sarawak). Phylogenetically, the Sarawak population was closely related to the Ranong population of western Peninsula Malaysia instead of other South-China-Sea populations, indicating its possible origin from the Indian Ocean Rim. No geographical subdivision was detected within the northern geographical region. An analysis of molecular variance (AMOVA) revealed low levels of genetic differentiation between and within mainland and island populations (phiCT = 0.015, phiSC = 0.037), indicating conspicuous long-distance seedling dispersal across oceans. Significant linkage disequilibrium excluded the possibility of recurrent homoplasious mutations as the major force causing phylogenetic discrepancy between mtDNA and the trnL-trnF spacer within the northern region. Instead, relative ages of alleles contributed to non-random chlorotype-mitotype associations and tree inconsistency. Widespread distribution and random associations (chi2 = 0.822, P = 0.189) of eight hypothetical ancestral cytotypes indicated the panmixis of populations of the northern geographical region as a whole. In contrast, rare and recently evolved alleles were restricted to marginal populations, revealing some preferential directional migration.

Organelle DNA phylogeography of Cycas taitungensis, a relict species in Taiwan

The phylogegraphic pattern of Cycas taitungensis, an endemic species with two remaining populations in Taiwan, was investigated based on genetic variability and phylogeny of the atpB-rbcL noncoding spacer of chloroplast DNA (cpDNA) and the ribosomal DNA (rDNA) internal transcribed spacer (ITS) of mitochondrial DNA (mtDNA). High levels of genetic variation at both organelle loci, due to frequent intramolecular recombination, and low levels of genetic differentiation were detected in the relict gymnosperm. The apportionment of genetic variation within and between populations agreed with a migrant-pool model, which describes a migratory pattern with colonists recruited from a random sample of earlier existing populations. Phylogenies obtained from cpDNA and mtDNA were discordant according to neighbour-joining analyses. In total four chlorotypes (clades I-IV) and five mitotypes (clades A-E) were identified based on minimum spanning networks of each locus. Significant linkage disequilibrium in mitotype-chlorotype associations excluded the possibility of the recurrent homoplasious mutations as the major force causing phylogenetic inconsistency. The most abundant chlorotype I was associated with all mitotypes and the most abundant mitotype C with all chlorotypes; no combinations of rare mitotypes with rare chlorotypes were found. According to nested clade analyses, such nonrandom associations may be ascribed to relative ages among alleles associated with the geological history through which cycads evolved. Nested in networks as interior nodes coupled with wide geographical distribution, the most dominant cytotypes of CI and EI may represent ancestral haplotypes of C. taitungensis with a possible long existence prior to the Pleistocene glacial maximum. In contrast, rare chlorotypes and mitotypes with restricted and patchy distribution may have relatively recent origins. Newly evolved genetic elements of mtDNA, with a low frequency, were likely to be associated with the dominant chlorotype, and vice versa, resulting in the nonrandom mitotype-chlorotype associations. Paraphyly of CI and EI cytotypes, leading to the low level of genetic differentiation between cycad populations, indicated a short period for isolation, which allowed low possibilities of the attainment of coalescence at polymorphic ancestral alleles.

Cytonuclear disequilibrium and genetic drift in a natural population of ponderosa pine

We measured the cytonuclear disequilibrium between 11 nuclear allozyme loci and both mitochondrial and chloroplast DNA haplotypes in a natural population of ponderosa pine (Pinus ponderosa, Laws). Three allozyme loci showed significant associations with mtDNA variation, while two other loci showed significant association with cpDNA. However, the absolute number of individuals involved in any of the associations was small, such that in none of the nuclear-organellar combinations was the difference between observed and expected numbers >11 individuals. Patterns of association were not consistent across loci or organellar genomes, suggesting that they are not the result of mating patterns, which would act uniformly on all loci. This pattern of disequilibria is consistent with the action of genetic drift and with existing knowledge of the structure of this population and thus does not imply the action of other evolutionary processes. The overall magnitude (normalized disequilibrium) of associations was greater for maternally inherited mtDNA than for paternally inherited cpDNA, though this difference was neither large nor significant. Such significant disequilibria involving the paternally inherited organelle indicate that not only are there a limited number of seed parents, but the effective number of pollen parents is also limited.

Linkage disequilibrium and phylogenetic congruence between chloroplast and mitochondrial haplotypes in Silene vulgaris

Both the chloroplast and mitochondrial genomes are used extensively in studies of plant population genetics and systematics. In the majority of angiosperms, the chloroplast DNA (cpDNA) and mitochondrial DNA (mtDNA) are each primarily transmitted maternally, but rare biparental transmission is possible. The extent to which the cpDNA and mtDNA are in linkage disequilibrium is argued to be dependent on the fidelity of co-transmission and the population structure. This study reports complete linkage disequilibrium between cpDNA and mtDNA haplotypes in 86 individuals from 17 populations of Silene vulgaris, a gynodioecious plant species. Phylogenetic analysis of cpDNA and mtDNA haplotypes within 14 individuals supports a hypothesis that the evolutionary histories of the chloroplasts and mitochondria are congruent within S. vulgaris, as might be expected if this association persists for long periods. This provides the first documentation of the evolutionary consequences of long-term associations between chloroplast and mitochondrial genomes within a species. Factors that contribute to the phylogenetic and linkage associations, as well as the potential for intergenomic hitchhiking resulting from selection on genes in one organellar genome are discussed.

The cytonuclear effects of facultative apomixis. II. Definitions and dynamics of disequilibria in tetraploid populations

We develop a cytonuclear framework for tetraploid populations in which a diallelic nuclear marker exhibits tetrasomic inheritance. This system requires two separate parameterizations, with six cytonuclear disequilibria (nonrandom associations) in tetraploid individuals and four in their diploid gametes. Double reduction during meiosis adds further complexity by causing gametic output to vary with the distance of the nuclear locus from the centromere. We derive and analyze dynamical solutions for the disequilibria under generalized mixed mating, with any combination of apomixis, selfing, and outcrossing, with and without double reduction. As in comparable diploid systems, all disequilibria ultimately decay to zero, unless nuclear and cytoplasmic alleles are nonrandomly associated and outcrossing is absent, in which case permanent associations result. Selfing and apomixis retard the decay of disequilibria (or approach to equilibrium), and often to the same extent. In contrast, double reduction can accelerate the loss of tetraploid cytonuclear associations, but only negligibly in hybrid zones, and this loss is never faster than in diploids. Only in the absence of allelic associations or outcrossing is the asymptotic approach to equilibrium differentially affected by apomixis and selfing or slower under tetrasomic than disomic inheritance. To facilitate empirical applications, we also examine tetraploid hybrid zone dynamics and offer practical guidelines for experimental design and data analysis, showing how the consequences of the mating system alone provide a valuable baseline for drawing evolutionary inferences from the observed patterns of cytonuclear associations.

The cytonuclear effects of facultative apomixis. I. Disequilibrium dynamics in diploid populations

We comprehensively analyze the cytonuclear effects of generalized mixed mating, including all combinations of selfing, outcrossing, and apomixis, the asexual production of seeds. After first deriving the time-dependent solutions for nonrandom associations (disequilibria) between a diallelic cytoplasmic marker and the alleles and genotypes at a diploid nuclear locus, we delimit all possible dynamical behaviors and the conditions under which each occurs. As in standard mixed mating systems, all disequilibria ultimately decay to zero except when outcrossing is absent, in which case permanent disequilibria result if the allelic association is initially nonzero. When at least some outcrossing is present, any initial allelic association decays at a constant geometric rate, whereas genotypic disequilibria may first increase in magnitude or change sign. Although selfing and apomixis tend to retard the decay of disequilibria (or approach to equilibrium) and often to the same extent, apomixis can have a stronger effect under some conditions. We also determine the dynamics of cytonuclear disequilibria in specific examples that may be of particular interest for empirical studies of hybrid zones. The results suggest several practical guidelines for experimental design and data analysis and show how the cytonuclear disequilibrium dynamics under mating system alone furnish a quantitative baseline for null hypotheses against which to test for the presence of other evolutionary forces.

The effects of pollen and seed migration on nuclear-dicytoplasmic systems. I. Nonrandom associations and equilibrium structure with both maternal and paternal cytoplasmic inheritance

We determine the nuclear-dicytoplasmic effects of unidirectional gene flow via pollen and seeds upon a mixed-mating plant population, focusing on nuclear-mitochondrial-chloroplast systems where mitochondria are inherited maternally and chloroplasts paternally, as in many conifers. After first delineating the general effects of admixture (via seeds or individuals) on the nonrandom associations in such systems, we derive the full dicytonuclear equilibrium structure, including when disequilibria may be indicators of gene flow. Substantial levels of permanent two- and three-locus disequilibria can be generated in adults by (i) nonzero disequilibria in the migrant pools or (ii) intermigrant admixture effects via different chloroplast frequencies in migrant pollen and seeds. Additionally, three-locus disequilibria can be generated by higher-order intermigrant effects such as different chloroplast frequencies in migrant pollen and seeds coupled with nuclear-mitochondrial disequilibria in migrant seeds, or different nuclear frequencies in migrant pollen and seeds coupled with mitochondrial-chloroplast disequilibria in migrant seeds. Further insight is provided by considering special cases with seed or pollen migration alone, complete random mating or selfing, or migrant pollen and seeds lacking disequilibria or intermigrant admixture effects. The results complete the theoretical foundation for a new method for estimating pollen and seed migration using joint cytonuclear or dicytonuclear data.

Symbiont survival and host-symbiont disequilibria under differential vertical transmission

Interspecific genetic interactions in host-symbiont systems raise intriguing coevolutionary questions and may influence the effectiveness of public health and management policies. Here we present an analytical and numerical investigation of the effects of host genetic heterogeneity in the rate of vertical transmission of a symbiont. We consider the baseline case with a monomorphic symbiont and a single diallelic locus in its diploid host, where vertical transmission is the sole force. Our analysis introduces interspecific disequilibria to quantify nonrandom associations between host genotypes and alleles and symbiont presence/absence. The transient and equilibrium behavior is examined in simulations with randomly generated initial conditions and transmission parameters. Compared to the case where vertical transmission rates are uniform across host genotypes, differential transmission (i) increases average symbiont survival from 50% to almost 60%, (ii) dramatically reduces the minimum average transmission rate for symbiont survival from 0.5 to 0.008, and (iii) readily creates permanent host-symbiont disequilibria de novo, whereas uniform transmission can neither create nor maintain such associations. On average, heterozygotes are slightly more likely to carry and maintain the symbiont in the population and are more randomly associated with the symbiont. Results show that simple evolutionary forces can create substantial nonrandom associations between two species.

The exact test for cytonuclear disequilibria

We extend the analysis of the statistical properties of cytonuclear disequilibria in two major ways. First, we develop the asymptotic sampling theory for the nonrandom associations between the alleles at a haploid cytoplasmic locus and the alleles and genotypes at a diploid nuclear locus, when there are an arbitrary number of alleles at each marker. This includes the derivation of the maximum likelihood estimators and their sampling variances for each disequilibrium measure, together with simple tests of the null hypothesis of no disequilibrium. In addition to these new asymptotic tests, we provide the first implementation of Fisher's exact test for the genotypic cytonuclear disequilibria and some approximations of the exact test. We also outline an exact test for allelic cytonuclear disequilibria in multiallelic systems. An exact test should be used for data sets when either the marginal frequencies are extreme or the sample size is small. The utility of this new sampling theory is illustrated through applications to recent nuclear-mtDNA and nuclear-cpDNA data sets. The results also apply to population surveys of nuclear loci in conjunction with markers in cytoplasmically inherited microorganisms.

Sampling theory for cytonuclear disequilibria

We examine the statistical properties of cytonuclear disequilibria within a system including one diploid nuclear locus and one haploid cytoplasmic locus, each with two alleles. The results provide practical guidelines for the design and interpretation of cytonuclear surveys seeking to utilize the novel evolutionary information recorded in the observed pattern of cytonuclear associations. Important applications include population studies of nuclear allozymes in conjunction with genes from mitochondria, chloroplasts, or cytoplasmically inherited microorganisms. Our attention focuses on the allelic and genotypic disequilibria, which respectively measure the nonrandom associations between the cytotypes and the nuclear alleles and genotypes. We first derive the maximum likelihood estimators and their approximate large sample variances for each disequilibrium measure. These are each in turn used to set up an asymptotic test of the null hypothesis of no disequilibrium. We then calculate the minimum sample sizes required to detect the disequilibria under specified alternate hypotheses. The work also incorporates the deviation from Hardy-Weinberg equilibrium at the nuclear locus, which can significantly affect the results. The practical utility of this new sampling theory is illustrated through applications to two nuclear-mitochondrial data sets.

Paternally inherited chloroplast polymorphism in Pinus: estimation of diversity and population subdivision, and tests of disequilibrium with a maternally inherited mitochondrial polymorphism

We have surveyed a chloroplast DNA restriction fragment length polymorphism in 745 individuals, distributed rangewide in eight allopatric natural populations of jack pine (Pinus banksiana Lamb.) and eight allopatric natural populations of lodgepole pine (Pinus contorta Dougl.). The polymorphic region of the chloroplast genome is located near duplicated psbA genes. Fourteen length variants were found in the survey, and these variants distinguished the two species qualitatively. Variant diversities were high in both species (hes = 0.43 in jack pine; hes = 0.44 in lodgepole pine). Population subdivision was weak within and among lodgepole pine subspecies and in jack pine (i.e., theta values were less than 0.05). This weak subdivision is compatible with theoretical predictions for paternally inherited markers in wind-pollinated outcrossers, as well as for polymorphisms with high length mutation rates. If these populations are at a drift-migration equilibrium, the chloroplast DNA restriction fragment data and previous mitochondrial frequency data from the same individuals are consistent with gene flow that is differential through seeds and pollen. The new data have permitted the first empirical tests of disequilibrium between maternally and paternally inherited factors. As expected, these tests failed to detect convincing evidence of nonrandom association between chloroplast and mitochondrial variants.

Comparative effects of pollen and seed migration on the cytonuclear structure of plant populations. II. Paternal cytoplasmic inheritance

We continue our study of the effects of pollen and seed migration on the cytonuclear structure of mixed-mating plant populations by analyzing two deterministic continent-island models under the critical assumption of paternal cytoplasmic inheritance. The major results of this study that contrast with our previous conclusions based on maternal cytoplasmic inheritance are (i) pollen gene flow can significantly affect the cytonuclear structure of the island population, and in particular can help to generate cytonuclear disequilibria that greatly exceed the magnitude of those that would be produced by seed migration or mixed mating alone; (ii) with simultaneous pollen and seed migration, nonzero cytonuclear disequilibria will be maintained not only when there is disequilibrium in the immigrant pollen or seeds, but also through a variety of intermigrant admixture effects when the two pools of immigrants differ appropriately in their cytonuclear compositions; (iii) either immigrant pollen or immigrant seeds can generate disequilibria de novo in populations with initially random cytonuclear associations, but pollen migration alone generally produces lower levels of disequilibrium than does comparable seed migration, especially at high levels of self-fertilization when the overall fraction of immigrant pollen is low; (iv) the equilibrium state of the island population will be influenced by the rate of pollen gene flow whenever there is either allelic disequilibrium in the immigrant pollen or simultaneous seed migration coupled with different cytoplasmic or nuclear allele frequencies in immigrant pollen and seeds or nonzero allelic disequilibrium in either immigrant pool. The estimation of pollen migration should therefore be facilitated with paternal cytoplasmic inheritance relative to the case of maternal cytoplasmic inheritance. These basic conclusions hold whether the population is censused as seeds or as adults, but with simultaneous pollen and seed migration, the relationship between census time and the ability to detect nonrandom cytonuclear associations is complex. When migration is through pollen alone, however, the cytonuclear structure of the island population is independent of the life stage censused.

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.

Paternal leakage of mitochondrial DNA inPinus

We studied mitochondrial DNA restriction fragment length polymorphism in 11 parents and 125 seedlings of 23 controlled matings within and between jack pine (Pinus banksiana Lamb.) and lodgepole pine (P. contorta Dougl.). A potential mitochondrial distinction between these two conifers was evident in the parental samples. Only maternal mitochondrial restriction fragments were observed in a majority of the seedlings, which is consistent with results from angiosperms and other members of the genusPinus L. However, we detected exclusively paternal mitochondrial DNA in six of the seedlings. These unusual seedlings were not attributable to heteroplasmy or contamination of the experimental material, indicating that mitochondrial inheritance was not strictly maternal. Paternal mitochondrial leakage inPinus may permit novel insights into the transmission genetics and evolution of organellar polymorphisms.

Comparative effects of pollen and seed migration on the cytonuclear structure of plant populations. I. Maternal cytoplasmic inheritance

We explicitly solve and analyze a series of deterministic continent-island models to delimit the effects of pollen and seed migration on cytonuclear frequencies and disequilibria in random-mating, mixed-mating and self-fertilized populations. Given the critical assumption of maternal cytoplasmic inheritance, five major findings are (i) nonzero cytonuclear disequilibria will be maintained in the island population if and only if at least some migration occurs each generation through seeds with nonrandom cytonuclear associations; (ii) immigrant seeds with no cytonuclear disequilibria can strongly affect the genetic structure of the island population by generating significant and long-lasting transient associations; (iii) with all else being equal, substantially greater admixture disequilibria are generally found with higher rates of seed migration into, or higher levels of self-fertilization within, the island population (with the possible exception of the heterozygote disequilibrium); (iv) pollen migration can either enhance or reduce the cytonuclear disequilibria caused by seed migration, or that due to mixed-mating in the absence of seed migration, but the effect is usually small and appears primarily to make a noticeable difference in predominantly outcrossing populations; and (v) pollen migration alone cannot generate even transient disequilibria de novo in populations with completely random associations. This same basic behavior is exhibited as long as there is some random outcrossing in the island population. Self-fertilized populations represent a special case, however, in that they are necessarily closed to pollen migration, and nonzero disequilibria can be maintained even in the absence of seed migration. All of these general results hold whether the population is censused as adults or as seeds, but the ability to detect nonrandom cytonuclear associations can depend strongly on the life stage censused in populations with a significant level of random outcrossing. We suggest how these models might be used for the estimation of seed and pollen migration.

The effects of admixture and population subdivision on cytonuclear disequilibria

We examine the generation of cytonuclear disequilibria by admixture and continued gene flow. General formulas analogous to the nuclear case are first derived showing that the allelic and genotypic disequilibria from admixture or population subdivision equal their expected value across the contributing (sub) populations plus the covariance across these sources between the cytoplasmic gene frequency and the relevant nuclear frequency. A detailed study is then presented of the cytonuclear dynamics, in a random-mating population under two different migration scenarios. In both cases closed-form solutions are given for all variables as a function of the initial conditions and relevant migration parameters. The dynamics of the gene frequencies and allelic disequilibria, which dominate each system, are the same as those involving two unlinked nuclear loci, while the dynamics of the genotypic disequilibria and cytonuclear frequencies have no nuclear counterpart. The continent-island formulation focuses on a population receiving continued immigration from a large source of constant composition. A major discovery is that cytonuclear disequilibria can transiently build up on the "island" to levels far exceeding those found at equilibrium. In contrast, the admixture formulation focuses on the dynamics within two populations undergoing continued intermigration. Although in this case all cytonuclear associations must ultimately decay to zero, long-term transient disequilibria can develop which are many times their initial admixture values. For both migration scenarios it is shown that the time of population censusing relative to migration and reproduction dramatically affects both the amount and pattern of the nonrandom associations produced. The empirical relevance of these models is discussed in light of nuclear-mitochondrial data from a hybrid zone between European and North American eels and from a zone of racial admixture in humans.

Chloroplast DNA restriction fragment length polymorphism in Sequoia sempervirens D. Don Endl., Pseudotsuga menziesii (Mirb.) Franco, Calocedrus decurrens (Torr.), and Pinus taeda L

The extent and type of chloroplast DNA restriction fragment length polymorphism was determined among individual tree samples of coast redwood, Douglas fir, incense-cedar, and loblolly pine. A total of 107 trees was surveyed for three restriction enzymes (BamHI, EcoRI, HindIII) and six chloroplast DNA probes from petunia (P3, P4, P6, P8, P10, S8). The probes comprise 64% of the petunia chloroplast genome. Polymorphisms were detected in all species but loblolly pine. Coast redwood and incense-cedar had a small number of rare variants, whereas Douglas fir had one highly polymorphic region of insertions/deletions in sequences revealed by the P6 probe from petunia. The mutation hotspot is currently being studied by DNA sequence analysis.