Population differentiation in an annual legume: genetic architecture

Abstract. The presence or absence of epistasis, or gene interaction, is explicitly assumed in many evolutionary models. Although many empirical studies have documented a role of epistasis in population divergence under laboratory conditions, there have been very few attempts at quantifying epistasis in the native environment where natural selection is expected to act. In addition, we have little understanding of the frequency with which epistasis contributes to the evolution of natural populations. In this study we used a quantitative genetic design to quantify the contribution of epistasis to population divergence for fitness components of a native annual legume, Chamaecrista fasciculata . The design incorporated the contrast of performance of F₂ and F₃ segregating progeny of 18 interpopulation crosses with the F₁ and their parents. Crosses were conducted between populations from 100 m to 2000 km apart. All generations were grown for two seasons in the natural environment of one of the parents. The F₁ often outperformed the parents. This F₁ heterosis reveals population structure and suggests that drift is a major contributor to population differentiation. The F₂ generation demonstrated that combining genes from different populations can sometimes have unexpected positive effects. However, the F₃ performance indicated that combining genes from different populations decreased vigor beyond that due to the expected loss of heterozygosity. Combined with previous data, our results suggest that both selection and drift contribute to population differentiation that is based on epistatic genetic divergence. Because only the F₃ consistently expressed hybrid breakdown, we conclude that the epistasis documented in our study reflects interactions among linked loci.     

A hierarchical analysis of minisatellite DNA diversity in Gambel oak (Quercus gambelii Nutt.; Fagaceae)

Gambel oak is a short shrub to medium-sized tree of southwest North America for which the potential to grow in large clonal stands has been proposed. Here we use three different multilocus VNTR DNA probes (synthesized via PCR) to demonstrate that clones growing together can be identified, but individual clones exceeding 50 m in diameter do not commonly occur in the study population of Quercus gambelii . Further, using a hierarchical sampling scheme (three circular transects with diameters of 9, 1100, and 28 000 m), we estimate that for nonclonal individuals: (i) mean number of bands analysed per individual = 22.92; (ii) mean similarity (band sharing) between individuals = 0.322; (iii) mean probability that two randomly chosen individuals share all bands = 4.93 × 10^–11; and (iv) mean estimated heterozygosity = 0.796. F _ST calculated for the two nonclonal levels of the sampling hierarchy was 0.023, indicating that little genetic differentiation exists between them. These results support previous findings that, due to life-history traits of the genus and this species (out-crossing, wind-pollinated, animal-dispersed, long-lived woody perennials), gene flow is high and genetic subdivision of populations is low in oaks. At this study site, the clonal nature demonstrated for Gambel oak appears to have little detectable effect on these population genetic characteristics.     

Comparing the genetic population structure of two species of arctic lemmings: more local differentiation in Lemmus trimucronatus than in Dicrostonyx groenlandicus

Collared and brown lemmings ( Dicrostonyx groenlandicus and Lemmus trimucronatus ) are two largely sympatric and ecologically comparable species of arctic microtine rodents, differing however in some respects which allow us to hypothesise differences in the genetic structure of their populations. Collared lemmings are particularly well adapted to life at high latitude, they occasionally emerge to the surface of the snow and may disperse over larger distances than brown lemmings – possibly even over snow and ice. This should result in more local differentiation among populations of brown lemmings than among populations of collared lemmings. We compared the genetic population structure between the two lemming species in a fragmented landscape with small islands in the central Canadian Arctic using four microsatellite loci and partial mitochondrial control region sequences. Both types of genetic markers showed higher differentiation ( F _ST values) among local populations for brown lemmings than for collared lemmings. We discuss to what extent the observed genetic differences may be explained by differences in dispersal rates in addition to differences in average effective population size.     

Population structure of Moroccan water frogs: genetic cohesion despite a fragmented distribution

This study analyses the allozymic variation of 20 presumptive loci in eight populations of Rana saharica from Morocco. Populations were collected from the very different climatic zones of this country: the Rif area, the Atlas mountains and the desert. Moroccan water-frog populations are genetically well differentiated from the geographically closed Algerian populations. Thus, to check if such a differentiation process is taking place within Moroccan water frogs, we attempted to analyse the genetic structure and patterns of gene flow of Moroccan populations, by means of estimates of Hardy-Weinberg equilibrium, F-statistics and indirect measures of gene flow. F_st(0.250) and F_is(0.254) values were similar, which means that both intra and interpopulation differentiation contribute equally to the amount of genetic divergence revealed. F_is values indicated some degree of structure within ponds, which is possibly related to the homing behaviour of some amphibians. On the other hand, F_st and genetic distances between populations were not very high. Despite the low levels of gene flow estimated, together with the homing behaviour revealed and the spatially discontinuous distribution, it was found that genetic differentiation among populations was not as high as expected. The likelihood of genetic homogeneity being the consequence of continuous population extinction and recolonization events is discussed.     

Responses to elevated temperature and CO2 in the perennial grass Agrostis curtisii in relation to population origin

1. Evolutionary responses to climate change will depend on the presence of heritable variation within species populations for traits that increase fitness under the changing conditions. Patterns of ecotypic differentiation in relation to latitude in some species suggest that such variation exists in relation to temperature responses. Response to elevated CO₂, whether heritable or not, is not expected to be related to latitudinal or climatic differences within temperate regions.
2. To test these ideas, seeds were collected from 10 populations of the outbreeding perennial grass Agrostis curtisii across its range in Europe from south Wales to Portugal. Plants were grown under ambient and elevated temperature and CO₂ conditions, in a factorial design, in solardomes; two half sibs from each population were planted in separate pots in each of the two replicate domes with each combination of treatments. One half sib was harvested at the end of the first summer, the second at the end of the second summer.
3. Survival was uniformly high and flowering uniformly low across treatments and populations.
4. Responses to temperature and CO₂ treatments varied over time for almost all populations. Treatment effects were not significant on plants harvested in year 1, although there was a trend towards higher shoot biomass under the elevated temperature and CO₂ treatment. In year 2 shoot biomass was significantly higher under the elevated temperature treatment across all populations and there was a strong trend towards decreased biomass under elevated CO₂.
5. There were no significant correlations of plant response to either CO₂ or temperature with climate at origin.
6. These results warn of the dangers of extrapolating evolutionary plant responses to CO₂ from short-term experiments.     

Genetic diversity and multilocus structure in sexual Otiorhynchus alpicola populations (Goleoptera: Gurculionidae)

Otiorhynchus alpicola Boh. (Coleoptera: Curculionidae) is distributed on high mountains of central, southern and south-eastern Europe. On the Balkan peninsula, this species is patchily distributed on mountain peaks at heights over 1800 m. To examine the organization of isoenzyme variation of the ten sexual (diploid) populations, starch gel electrophoresis of enzymes was used. The average heterozygosity over 14 gene loci was 0.11. For diploid O. alpicola, F-statistics were used to assess population heterogeneity and substructuring. The data indicate that, compared with other insects, this high-altitude weevil species is genetically very differentiated (the average F _ST= 0.35). In addition, gene-flow among populations was extremely low (the estimates of N_em from Wright's F _ST and Slatkin private-allele methods were 0.47 and 0.83, respectively). Pairwise linkage disequilibrium (LD) parameters were estimated from zygotic frequencies using Burrows' method. The average rate of significant LD in analysed populations of O. alpicola was about 27%. The pattern of LD over the diploid populations indicates that stochastic factors might be a primary cause of the observed multilocus associations.     

Long-term stability and effective population size in North Sea and Baltic Sea cod (Gadus morhua)

DNA from archived otoliths was used to explore the temporal stability of the genetic composition of two cod populations, the Moray Firth (North Sea) sampled in 1965 and 2002, and the Bornholm Basin (Baltic Sea) sampled in 1928 and 1997. We found no significant changes in the allele frequencies for the Moray Firth population, while subtle but significant genetic changes over time were detected for the Bornholm Basin population. Estimates of the effective population size ( N _e ) generally exceeded 500 for both populations when employing a number of varieties of the temporal genetic method. However, confidence intervals were very wide and N _e 's most likely range in the thousands. There was no apparent loss of genetic variability and no evidence of a genetic bottleneck for either of the populations. Calculations of the expected levels of genetic variability under different scenarios of N _e showed that the number of alleles commonly reported at microsatellite loci in Atlantic cod is best explained by N _e 's exceeding thousand. Recent fishery-induced bottlenecks can, however, not be ruled out as an explanation for the apparent discrepancy between high levels of variability and recently reported estimates of N _e  << 1000. From life history traits and estimates of survival rates in the wild, we evaluate the compatibility of the species' biology and extremely low N _e / N ratios. Our data suggest that very small N _e 's are not likely to be of general concern for cod populations and, accordingly, most populations do not face any severe threat of losing evolutionary potential due to genetic drift.     

Heterozygosity and components of fitness in a strain of rainbow trout

We examined the relationship between multilocus heterozygosity and several components of fitness in the Arlee strain of rainbow trout (Salmo gairdneri) . There is a significant negative association between multilocus heterozygosity and standardized rate of oxygen consumption (VO₂) in 168 fish produced from a pooled mating between 55 females and 60 males. Fish which are more heterozygous also tend to be longer, heavier and have a significantly greater condition factor than fish which are less heterozygous. However, no relationship is evident between VO₂, weight, length, condition factor and multilocus heterozygosity within two full-sib families. These results indicate that the larger size of more heterozygous fish may partially be attributed to the increased metabolic efficiency in heterozygotes compared to homozygotes and that these differences may only be detectable at the population level. More heterozygous females of the same age also have larger eggs. Energy conversion may be more efficient in more heterozygous females resulting in the production of larger eggs. The prediction that fish with faster developmental rates are metabolically more efficient and, therefore, have lower VO₂ was also tested. Earlier hatching fish from five families do not have significantly lower VO₂ than their later hatching full-sibs.     

Evidence for coadaptation in geographic populations of Drosophila bipectinata

In order to test whether there is genetic coadaptation in geographic populations of Drosophila bipectinata with respect to body size, reciprocal crosses were made among five strains derived from ecogeographically different localities in India. Wing length was used as an index of body size, and was measured in all the five strains, and their crosses in F and F₂ generations. The statistical analysis of the data show that there is significant interpopulation variation in body size and in all the crosses, there is an increase in body size in F₁ generation when compared with mid-parent value. Further, there is a decrease in body size in F₂ generation as compared to F₁ in most of the crosses with increased variability. These results provide evidence for genetic coadaptation in geographic populations of D. bipectinata .     

Mechanistic basis of differences in water-use efficiency between a CAM and a C3 species of Peperomia (Piperaceae)

Under ecologically realistic environmental conditions, the water-use efficiency (WUE) of Peperomia scandens , a CAM plant, was higher than that of the C₃ congener P. obtusifolia . This difference has been attributed to differences in stomatal activity between C₃ and CAM plants, coupled with differences in the evaporative demand of the atmosphere during which the stomata are open. This explanation has apparently not, however, been experimentally tested. Thus, WUEs were compared in these species in two experiments in which the atmospheric evaporative demand was identical (or nearly so) during the period of stomatal opening (i.e. during the night for the CAM plant and during the day for the C₃ species). In both experiments, the WUE of the CAM species was higher than that of the C₃ species. These results suggest that factors other than differences in atmospheric environmental conditions must also be responsible for the observed differences in WUE. Because CO₂ uptake rates of the CAM species were substantially lower than those of the C₃ species, the lower WUE in the CAM species resulted primarily from lower transpiration rates. Lower rates of water loss in P. scandens , relative to rates in P. obtusifolia , were ascribed, in part, to lower stomatal densities. Thus, leaf morphological differences, in addition to differences in atmospheric evaporative demand, help to explain the high WUE typically measured in CAM plants.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.     

Population genetic structure of wild rice Oryza glumaepatula distributed in the Amazon flood area influenced by its life-history traits

Wild plant species develop their own way of living to adapt to the specific environment of their habitats. Their life-history traits strongly affect the genetic structure of the population. The wild species Oryza glumaepatula Steud. growing in the Amazon basin seems to have characteristic life-history traits suited for the flood condition. At the vegetative growth stage, the culms frequently break at internodes. With no roots anchoring on the ground, plant bodies floating in the water move downriver by water current and wind. To examine the association between the life-history traits and genetic population structure of Amazonian O. glumaepatula , we analysed allozyme variability at 29 loci of 16 enzymes using 37 populations from five regions. Allozymes were not so variable (total gene diversity H _E = 0.044) compared with Asian wild rice, O. rufipogon Griff. The bottleneck effect and rare opportunity of interspecies gene flow may prevent the development of allozyme variability. Population genotypes tended to be differentiated among geographically isolated regions. Observed heterozygosities were much lower than expected heterozygosities, or gene diversity ( H _O = 0.003 for whole population) and F _IS over polymorphic loci was 0.931, indicating that O. glumaepatula has developed an inbreeding system. But, the intrapopulation gene diversity ( H _S) was higher than interpopulation gene diversity ( D _ST), as generally observed in outbreeding populations. The migration ability of O. glumaepatula makes long-distance seed dispersal possible. This might have led to frequent gene flow among populations.     

Genetic variability in Irish populations of the invasive zebra mussel, Dreissena polymorpha: discordant estimates of population differentiation from allozymes and microsatellites

1. The recent arrival and explosive spread of the zebra mussel, Dreissena polymorpha (Pallas), in Ireland provided a rare opportunity to study the population genetics of an invasive species.
2. Eight polymorphic allozyme loci ( ACO-1, ACO-2 , EST-D, GPI, IDH-2, MDH, OPDH and PGM ) were used to investigate genetic diversity and population structure in five Irish populations, and the results were compared with those from a previous microsatellite study on the same samples.
3. The mean number of alleles per locus (2.7 ± 0.1) was similar to the mean for the same loci in European populations, suggesting that Irish founder populations were large and/or multiple colonization events took place after foundation. A deficiency of heterozygotes was observed in all populations, but was uneven across loci.
4. Pairwise comparisons, using Fisher's exact tests and F _ST values, revealed significant genetic differentiation among populations. The overall multilocus F _ST estimate was 0.118 ± 0.045, which contrasted with an estimate of 0.015 ± 0.007 from five microsatellite loci on the same samples in a previous study.
5. Assuming that microsatellites can be used as a neutral baseline, the discordant results from allozymes and microsatellites suggest that selection may be acting on some allozyme loci, specifically ACO-1, ACO-2 , IDH-2 and MDH, which contributed most to the significant differentiation between samples.     

Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

Genetic structure of Lake Magadi tilapia populations

The control region of the mitochondrial DNA haplotype frequencies were significantly different among the two separate lagoon populations of Oreochromis alcalicus grahami in Lake Magadi and of O. a. alcalicus from lake Natron, and DNA fingerprint similarity indices were significantly higher for intra-population comparisons of the two Magadi lagoon populations and the Lake Natron population than the inter-population similarity indices among these populations. A modified F _st measure indicated population sub-division and the phylogeographic partitioning of the VNTR fragments observed were unique to specific populations further indicating substantial genetic differentiation. The lagoon populations within Lake Magadi demonstrated the same degree of genetic differentiation as either of these populations did to the outgroup (the Lake Natron population). There appears to be limited gene flow between Lake Magadi tilapia populations and this population structure has important implications for protecting locally adapted populations within this unique ecosystem.     

Gene flow and effective population size in Lepanthes (Orchidaceae): a case for genetic drift

Genetic drift can play an important role in population differentiation, particularly when effective population sizes are small and gene flow is limited. Such conditions are suspected to be common in the species-rich Orchidaceae. We investigated the likelihood of genetic drift in natural populations of three endemic species of Lepanthes (Orchidaceae) from Puerto Rico. We estimated effective population size, Ne, using three ecologically based methods. Two of the three estimates were based on variance in reproductive potential and the third was based on coalescence time. All estimates of Ne were usually <40% of the standing population size, resulting in values of <20 individuals per population. Based on starch gel electrophoresis of isozymes, Nm estimates suggest restricted gene flow among populations in the range of one or less successful migrant per generation. Genetic differentiation among populations is expected under such conditions from random genetic drift. Indeed we observed high genetic differentiation among populations (L. rubripetala, F_ST> G_ST, 6; θ.248, 0.266, 0.293; L. rupestris, 0.148, 0.169, 0.138; L. eltoroensis, 0.251, 0.219, 0.218, respectively). Genetic drift is likely to be important for population differentiation in Lepanthes as a result of small effective population sizes and restricted gene flow.     

Genetic diversity and genetic structure of populations of Ranunculus japonicus Thunb. (Ranunculaceae)

Abstract In order to clarify the genetic diversity and population structure of Ranunculus japonicus , allozymic analysis was conducted on 60 populations in southwestern Japan. Considerable genetic variati ons were detected among the populations of R. japonicus . The genetic diversities within species ( H _es = 0.215) and within populations ( H _ep = 0.172) were slightly higher than those of other perennial herbs with widespread distribution and outcrossing plants. Significantly higher values of fixation index were detected in some populations, which might have arisen from restricted mating partners. The majority of genetic variation (approx. 80%) resided within a population and a moderate level of genetic differentiation ( G _ST = 0.203) was observed among populations. The F _ST value (0.203) suggests the existence of a substantial population structure in this species. The highly significant correlation between geographic distance and F _ST values indicates that isolation by distance has played an important role in the construction of the genetic structure of this species.     

Crossability variation in four South American Hordeum species

Andersson, 2. 1995. Crossability variation in four South American Hordeum species— Nord. J. Bot. IS: 355–364. Copenhagen. ISSN 0107–055X.
Intra- and interpopulation crosses within the four South American Hordeum species, H. cordobense, H. chilense, H. muticum and H. stenostachys were performed. Seeds were obtained in all intrapopulation crosses within all four species. Hordeum cordobense showed little variation in crossability of interpopulation crosses and the F₁-families had rather high fertility. A wide range in crossability between populations was observed in the three species H. chilense, H. muticum and H. stenostachys . Some combinations failed to set seed while others had high seed set. The fertility of the F₁-_- families varied from almost sterile to fully fertile. All parental plants, intra and interpopulation hybrids in the four species had very high meiotic pairing frequency with 12–14 chiasmata per cell. Hordeum cordobense is a homogeneous species while the other three species are heterogeneous and genetically diverse.     

Genetic structure within and among regional populations of the Eurasian badger (Meles meles) from Denmark and the Netherlands

The Eurasian badger Meles meles has a wide distribution area ranging from Japan to Ireland. In western Europe badger habitats are severely disturbed by anthropogenic factors, leading to fragmentation into subpopulations and formation of a metapopulation substructuring of once continuous panmictic populations. We have examined the genetic structure of Dutch and Danish badger populations on a relatively small scale (within countries) and a larger scale (between countries). The levels of genetic variation of populations were moderate and did not differ significantly among populations (overall H _O=0.30, overall H _E=0.34). Considerable genetic differentiation between the Dutch and Danish populations was found (overall F _ST=0.32, mean pairwise Dutch–Danish F _ST=0.42), indicating a large-scale substructuring of these western European badger populations. Further analysis showed that the Danish badger population can be substructured into three clusters [ P _{( k =3)}=0.99], but the Dutch populations cluster into one more or less panmictic population [ P _{( k =1)}=0.99] with little or no substructuring. The presence of migration barriers, such as roads, together with the peninsular geography of Denmark, may have led to this structuring of badger populations. In contrast, measures that improve migration and connection to other populations from neighboring countries may have prevented substructuring of the Dutch badger population.     

Population genetic structure in Myrtus communis L. in a chronically fragmented landscape in the Mediterranean: can gene flow counteract habitat perturbation?

Ancient managed landscapes provide ideal opportunities to assess the consequences of habitat fragmentation on the patterns of genetic diversity and gene flow in long-lived plant species. Using amplified fragment length polymorphism (AFLP) and allozyme markers, we quantified seed-mediated gene flow and population genetic diversity and structure in 14 populations of Myrtus communis (myrtle), a common endozoochorous shrub species of forest patches in lowland agricultural Mediterranean areas. Overall, allozyme diversity for myrtle was low (P₉₅   =   25%; A   =   1.411; H_e = 0.085) compared to other known populations, and a significant portion of populations (57%) had lower levels of allelic diversity and/or heterozygosity than expected at random, as shown by simulated resampling of the whole diversity of the landscape. We found significant correlations between allozyme variability and population size and patch isolation, but no significant inbreeding in any population. Genetic differentiation among populations for both allozyme and AFLP markers was significant (Φ_ST = 0.144 and Φ_ST = 0.142, respectively) but an isolation-by-distance pattern was not detected. Assignment tests on AFLP data indicated a high immigration rate in the populations ( ca. 20–22%), likely through effective seed dispersal across the landscape by birds and mammals. Our results suggest that genetic isolation is not the automatic outcome of habitat destruction since substantial levels of seed-mediated gene flow are currently detectable. However, even moderate rates of gene flow seem insufficient in this long-lived species to counteract the genetic erosion and differentiation imposed by chronic habitat destruction.