Genetic diversity of the black mangrove (Avicennia germinans L.) in Colombia

This study analyzed the genetic diversity and patterns of genetic structure in Colombian populations of Avicennia germinans L. using microsatellite loci. A lower genetic diversity was found on both the Caribbean (Ho = 0.439) and the Pacific coasts (Ho = 0.277) than reported for the same species in other locations of Central American Pacific, suggesting the deterioration of genetic diversity. All the populations showed high inbreeding coefficients (0.131–0.462) indicating heterozygotes deficience. The genetic structure between the Colombian coasts separated by Central American Isthmus was high (F_RT = 0.39) and the analyses of the genetic patterns of A. germinans revealed a clear differentiation of populations and no-recent gene flow evidence between coasts. Genetic structure was found within each coast (F_ST = 0.10 for the Caribbean coast and F_ST = 0.22 for the Pacific coast). The genetic patterns along the two coasts appear to reflect a forcing by local geomorphology and marine currents. Both coasts constitute a different Evolutionary Significant Unit, so we suggest for future transplantations plans that propagules or saplings of the populations of the Caribbean coast should not be mixed with those of the Pacific Colombian coast. Besides, we suggest that reforestation efforts should carefully distinguish propagules sources within each coast.     

Biochemical assessment of the taxonomic diversity of the operculate land snail, Cyclophorus fulguratus (Gastropoda: Cyclophoridae), from Thailand

Allozyme variation was examined in 223 samples of the operculate land snail Cyclophorus fulguratus from 13 localities across three regions of Thailand. Using horizontal starch gel electrophoresis, 13 allozyme presumed loci (12 polymorphic) were screened. Heterozygosity was moderate in C. fulguratus (H_exp = 0.008–0.127) with a high genetic heterogeneity among samples (F_st = 0.734). Populations showed a greater genetic differentiation in central Thailand (F_st = 0.380) than in northeastern Thailand (F_st = 0.108), suggesting frequent gene flow among populations in northeastern Thailand. C. fulguratus exhibits a strong pattern of isolation by distance over the entire tested species range in Thailand and may potentially have been involved in an extensive local fragmentation. Results of the distance analysis revealed that large genetic divergence has occurred among the central, northeastern and eastern Thailand groups [D = 0.361–0.701], strongly suggesting populations from these three geographical regions may actually represent or else be evolving into separate species.     

Population genetic structure of striped mullet, Mugil cephalus, along the coast of China, inferred by AFLP fingerprinting

The striped mullet, Mugil cephalus, is an economically important species for both aquaculture and commercial fisheries in China. In this study, the amplified fragment length polymorphism (AFLP) technique was employed to analyze population genetic diversity and genetic distance between different populations with the aim of elucidating the population structure and gene flow of M. cephalus along the coast of China. A total of 230 fragments with 150–500 bp were identified from 118 individuals by five AFLP primer combinations. The polymorphic loci within populations varied from 46.52% to 64.78%, with an average of 53.91%, and the average heterozygosity from 0.1829 to 0.2282, with an average of 0.2074. The UPGMA phenograms of 118 individuals were constructed based on Dice similarity coefficients and four clusters were recognized. AMOVA analysis revealed that 60.7% of genetic variations were within populations and 39.3% between populations. The estimated genetic distance (φ_ST) value over all polymorphic loci across the six populations was 0.393 (p < 0.0010), indicating a strong population structuring. The pairwise φ_ST value ranged from 0.1112 to 0.5358, with an average of 0.3693. The population pairwise gene flows (average N_m = 0.73) are low. In addition, the result of the Mantel test showed that there was a significant correlation between geographic and genetic distances (r = 0.5434, p = 0.0050). It was speculated that there exist at least four distinct geographic populational subdivisions of M. cephalus along the Chinese coast. This research has provided new molecular data which will aid our understanding of the genetic structure of this species.     

Evidence for adaptive phenotypic differentiation in Baltic Sea sticklebacks

The evidence for adaptive phenotypic differentiation in mobile marine species remains scarce, partly due to the difficulty of obtaining quantitative genetic data to demonstrate the genetic basis of the observed phenotypic differentiation. Using a combination of phenotypic and molecular genetic approaches, we elucidated the relative roles of natural selection and genetic drift in explaining lateral plate number differentiation in threespine sticklebacks (Gasterosteus aculeatus) across the entire Baltic Sea basin (approximately 392 000 km²). We found that phenotypic differentiation (P_ST = 0.213) in plate number exceeded that in neutral markers (F_ST = 0.008), suggesting an adaptive basis for the observed differentiation. Because a close correspondence was found between plate phenotype and genotype at a quantitative trait loci (QTL; STN381) tightly linked to the gene (Ectodysplasin) underlying plate variation, the evidence for adaptive differentiation was confirmed by comparison of F_ST at the QTL (F_STQ = 0.089) with F_ST at neutral marker loci. Hence, the results provide a comprehensive demonstration of adaptive phenotypic differentiation in a high‐gene‐flow marine environment with direct, rather than inferred, verification for the genetic basis of this differentiation. In general, the results illustrate the utility of P_ST–F_ST–F_STQ comparisons in uncovering footprints of natural selection and evolution and add to the growing evidence for adaptive genetic differentiation in high‐gene‐flow marine environments, including that of the relatively young Baltic Sea.     

Genetic variation and hybridization in SwedishSchoenus (Cyperaceae)

Schoenus ferrugineus andS. nigricans have restricted distributions in Sweden and are almost exclusively confined to calcareous fen habitats. AtS. nigricans sites,S. ferrugineus is usually also present, and hybrids are frequently found. In this report, I used allozymes to estimate the amount of gene flow between the two species, and to compare the partitioning of genetic diversity in each of them. Thirteen loci were analysed at eight different enzyme systems. Seven loci were variable between or within the species. The two species had completely different alleles at two of the seven variable loci, whereas there was overlap at five loci. In all, 22 different alleles were found. Six of these alleles were confined toS. nigricans, and five alleles were confined toS. ferrugineus. Nei's genetic identity was 0.55.—InS. ferrugineus, three loci (23%) were polymorphic, and the average number of alleles per polymorphic locus was 2.0 (each polymorphic locus had two alleles). InS. nigricans, three loci (23%) were polymorphic, and the average number of alleles per polymorphic locus was 2.3.—The proportion of genetic diversity due to variation among sites (G _ST) was fairly similar in the two species, mean over loci = 0.12 inS. ferrugineus and 0.15 inS. nigricans. However, the proportion of genetic diversity due to variation among individuals within sites (G _IS) differed markedly between the two species, mean over loci = 0.54 inS. ferrugineus and 0.17 inS. nigricans. Accordingly, there was a much higher individual heterozygosity inS. nigricans than inS. ferrugineus. — Most hybrids were interpreted as F₁ hybrids. However, a small proportion, 0.5–1.6 %, were F_n hybrids or back-crosses.—On the Swedish mainland, all former occurrences ofS. nigricans are extinct, but viable hybrids are still present at a few sites in southernmost Sweden.     

Analysis of F ST outliers at allozyme loci in Pacific salmon: implications for natural selection

Natural selection has been invoked to explain the observed geographic distribution of allozyme allele frequencies for a number of teleost species. The effects of selection on allozyme loci in three species of Pacific salmon were tested. A simulation-based approach to estimate the null distribution of population differentiation (F _ST) and test for F _ST outliers was used. This approach showed that a majority of allozyme loci conform to neutral expectations predicted by the simulation model, with relatively few F _ST outliers found. No consistent F _ST outlier loci were found across species. Analysis of population sub-groups based on geography and genetic identity reduced the number of outlier loci for some species, indicating that large geographic groups may include genetically divergent populations and/or that there is geographic heterogeneity in selection pressure upon allozyme loci. Two outlier allozyme loci found in this analysis, lactate dehydrogenase-B and malic enzyme, have been shown to be influenced by selection in other teleost species. This approach is also useful in identifying allozyme loci (or other genetic markers) that meet assumptions for population genetic study.     

Genetic variation and population structure in Korean populations of sand dune speciesSalsola komarovi (Chenopodiaceae)

Salsola komarovi lljin is a herbaceous annual native to the sand dunes and beaches of Japan, northern China, Sakhalln and Korea. Starch-gel electrophoresis was conducted on leaves and stems collected from 300 plants in eight Korean populations. The mean number of alleles per locus (A _p=1.51), mean expected heterozygosity (He _p=0.116), and total genetic diversity (H _T=0.279) were comparable with those for species with similar life history and ecological traits. A general conformance of genotype frequencies to Hardy-Weinberg expectations (meanF _IS=−0.030) indicates thatS. komarovi is an outcrossing species. Slightly more than 20% of the genetic variation was found among populations (F _ST=0.204). In addition, significant differences in allele frequency were detected between populations at all 11 polymorphic loci (P<0.001). Nei's genetic identities range from 0.885 to 0.985 with a mean of 0.942. However, indirect estimates of the number of migrant per generation (0.97, calculated fromF _ST and 0.31, calculated from seven private alleles) indicate that the levels of gene flow is low among Korean populations. Although the species maintains a moderate level of genetic variation within populations, the small, isolated natural populations of the species have been severely destructed by human activities, particularly in summer season. If this is true, conservation efforts should be focused on those populations that currently maintain the most genetic diversity (e.g., populations of Cheju Island and coast of the southwestern Korean Peninsula).     

Diversity and genetic connectivity among populations of a threatened tree (<Emphasis Type="Italic">Dalbergia nigra</Emphasis>) in a recently fragmented landscape of the Brazilian Atlantic Forest

In this study we evaluated the influence of recent landscape fragmentation on the dynamics of remnant fragments from the Brazilian Atlantic Forest. This biome is one of the richest in the world and has been extensively deforested and fragmented. We sampled five populations of the threatened Dalbergia nigra, a tree endemic to the Brazilian Atlantic Forest, two located in a large reserve of continuous forest and three in fragments of different sizes and levels of disturbance. In order to assess historical changes, considering the longevity of the analyzed species, 119 adults and 116 saplings were genotyped for six microsatellite loci. Lower levels of genetic diversity were found in the most impacted fragments when compared to the most preserved population located inside the reserve, and there was significant genetic structure among the populations studied (pairwise F _ST = 0.031–0.152; pairwise D _EST = 0.039–0.301). However, genetic structure among saplings (F _ST = 0.056; D _EST = 0.231) was significantly lower than among adults (F _ST = 0.088; D _EST = 0.275). Estimates of contemporary gene flow based on assignment tests corroborated this result, suggesting that fragmentation led to an increase in gene flow. This connectivity among remnant fragments could mitigate the loss of genetic diversity through a metapopulation dynamic, but the high rate of habitat loss and the unknown long-term genetic effects add uncertainty. These results, taken together with the presence of private alleles in disturbed populations, highlight the importance of preserving the extant fragments.     

Comparison of genetic differentiation at marker loci and quantitative traits

The comparison of the degree of differentiation in neutral marker loci and genes coding quantitative traits with standardized and equivalent measures of genetic differentiation (F_ST and Q_ST, respectively) can provide insights into two important but seldom explored questions in evolutionary genetics: (i) what is the relative importance of random genetic drift and directional natural selection as causes of population differentiation in quantitative traits, and (ii) does the degree of divergence in neutral marker loci predict the degree of divergence in genes coding quantitative traits? Examination of data from 18 independent studies of plants and animals using both standard statistical and meta‐analytical methods revealed a number of interesting points. First, the degree of differentiation in quantitative traits (Q_ST) typically exceeds that observed in neutral marker genes (F_ST), suggesting a prominent role for natural selection in accounting for patterns of quantitative trait differentiation among contemporary populations. Second, the F_ST – Q_ST difference is more pronounced for allozyme markers and morphological traits, than for other kinds of molecular markers and life‐history traits. Third, very few studies reveal situations were Q_ST < F_ST, suggesting that selection pressures, and hence optimal phenotypes, in different populations of the same species are unlikely to be often similar. Fourth, there is a strong correlation between Q_ST and F_ST indices across the different studies for allozyme (r=0.81), microsatellite (r=0.87) and combined (r=0.75) marker data, suggesting that the degree of genetic differentiation in neutral marker loci is closely predictive of the degree of differentiation in loci coding quantitative traits. However, these interpretations are subject to a number of assumptions about the data and methods used to derive the estimates of population differentiation in the two sets of traits.     

Genetic variability and the effect of habitat fragmentation in spotted suslik Spermophilus suslicus populations from two different regions

Endangered species worldwide exist in remnant populations, often within fragmented landscapes. Although assessment of genetic diversity in fragmented habitats is very important for conservation purposes, it is usually impossible to evaluate the amount of diversity that has actually been lost. Here, we compared population structure and levels of genetic diversity within populations of spotted suslik Spermophilus suslicus, inhabiting two different parts of the species range characterized by different levels of habitat connectivity. We used microsatellites to analyze 10 critically endangered populations located at the western part of the range, where suslik habitat have been severely devastated due to agriculture industrialization. Their genetic composition was compared with four populations from the eastern part of the range where the species still occupies habitat with reasonable levels of connectivity. In the western region, we detected extreme population structure (F _ST = 0.20) and levels of genetic diversity (Allelic richness ranged from 1.45 to 3.07) characteristic for highly endangered populations. Alternatively, in the eastern region we found significantly higher allelic richness (from 5.09 to 5.81) and insignificant population structure (F _ST = 0.03). As we identified a strong correlation between genetic and geographic distance and a lack of private alleles in the western region, we conclude that extreme population structure and lower genetic diversity is due to recent habitat loss. Results from this study provide guidelines for conservation and management of this highly endangered species.     

Genetic structure of refugial populations of the temperate plant Shortia rotundifolia (Diapensiaceae) on a subtropical island

Continental island systems harbour relict biota and populations that might have migrated during glacial periods due to the formation of landbridges. Here we analysed the genetic structure of relict populations of the temperate plant Shortia rotundifolia on the subtropical island of Iriomotejima, Japan. This plant, which inhabits riparian environments, is designated “near threatened”. Only five extant populations have been found on the island. Our analyses of 10 nuclear microsatellite loci detected genetic diversity of H _E = 0.488 and H _O = 0.358 for all populations of S. rotundifolia on the island. A high inbreeding coefficient for all populations together (F _IS = 0.316) and each population separately (F _IS = 0.258–0.497) might be attributable to crossing among closely related descendants within a population, an idea that is supported by the relatedness coefficient. These results and an examination of the populations’ demographic histories suggest that the extant populations on Iriomotejima have not experienced a recent population bottleneck. The five extant populations were genetically differentiated (F _ST = 0.283; P < 0.001), suggesting low seed dispersal by gravity and/or low pollen flow via pollinators in the riparian environment. In addition, population differentiation was not related to genetic distance, implying that at one time, ancestral populations might have been distributed over a wider area of the island. However, population fragmentation and range contraction might have occurred at random during the postglacial period.     

Asynchronous spawning in sympatric populations of a hard coral reveals cryptic species and ancient genetic lineages

Genetic subdivision within a species is a vital component of the evolution of biodiversity. In some species of Acropora corals in Western Australia, conspecific individuals spawn in two seasons 6 months apart, which has the potential to impede gene flow and result in genetic divergence. Genetic comparison of sympatric spring and autumn spawners of Acropora samoensis was conducted to assess the level of reproductive isolation and genetic divergence between the spawning groups based on multiple loci (13 microsatellite loci, the mitochondrial control region and two nuclear introns). Bayesian clustering and principal coordinate analysis of the microsatellite loci showed a high level of genetic differentiation between the spawning groups (F’_ST = 0.30; P < 0.001), as did the sequence data from PaxC and Calmodulin (Φ_ST = 0.97 and 0.31, respectively). At the PaxC locus, the autumn and spring spawners were associated with two divergent lineages that were separated by an evolutionary distance of 1.7% and statistical tests indicate divergent selection in PaxC, suggesting this gene may play a role in coral spawning. This study indicates that the autumn and spring spawners represent two cryptic species, and highlights the importance of asynchronous spawning as a mechanism influencing speciation in corals.     

Phylogeography of two East Asian species in Croomia (Stemonaceae) inferred from chloroplast DNA and ISSR fingerprinting variation

The genus Croomia (Stemonaceae) comprises three herbaceous perennial species that are distributed in temperate-deciduous forests in Southeastern North America (C. pauciflora) and East Asia (C. japonica, C. heterosepala). The two Asian species have abutting ranges in South Japan, but C. japonica also occurs disjunctively on the adjacent Asiatic mainland in East China. In our phylogenetic analysis of Croomia, based on chloroplast (cp) DNA sequence variation of the trnL-F region, and rooted with Stemona spp., the two Asian species are identified as sister that likely diverged in the Mid-to-Late Pleistocene (0.84–0.13 mya), whereas the divergence of C. pauciflora dates back to the Late Plio-/Pleistocene (<2.6 mya). Phylogeographical analysis of the two East Asian species detected seven cpDNA (trnL-F) haplotypes across 16 populations surveyed, and all of those were fixed for a particular cpDNA haplotype (H_E = 0.0, G_ST = 1). A survey of inter-simple sequence repeats (ISSRs) markers also detected remarkably low levels of within-population diversity (C. japonica: H_E = 0.085; C. heterosepala: H_E = 0.125), and high levels of inter-population differentiation (C. japonica: Φ_ST = 0.736; C. heterosepala: Φ_ST = 0.550), at least partly due to pronounced regional genetic substructure within both species. Non-overlapping distributions of cpDNA haplotypes and strong genetic (cpDNA/ISSR) differentiation among populations and/or regions accord with findings of a nested clade analysis, which inferred allopatric fragmentation as the major process influencing the spatial haplotype distribution of the two species. Based on mismatch distribution analysis and neutrality tests, we do not find evidence of population expansion in both species. Overall, we conclude that components of temperate-deciduous forest types in South Japan and East China are particularly sensitive to range fragmentation, isolation, and enhanced (incipient) species formation through climate-induced expansions of other forest types over glacial and interglacial periods of the (Late) Quaternary.     

Identifying chromosomal selection‐sweep regions in facial eczema selection‐line animals using an ovine 50K‐SNP array

Facial eczema (FE) is a hepato‐mycotoxicosis found mainly in New Zealand sheep and cattle. When genetics was found to be a factor in FE susceptibility, resistant and susceptible selection lines of Romney sheep were established to enable further investigations of this disease trait. Using the Illumina OvineSNP50 BeadChip, we conducted a selection‐sweep experiment on these FE genetic lines. Two analytical methods were used to detect selection signals, namely the Peddrift test (Dodds & McEwan, 1997) and fixation index F_ST (Weir & Hill, 2002). Of 50 975 single nucleotide polymorphism (SNP) markers tested, there were three that showed highly significant allele frequency differences between the resistant and susceptible animals (Peddrift nominal P < 0.000001). These SNP loci are located on chromosomes OAR1, OAR11 and OAR12 that coincide precisely with the three highest genomic F_ST peaks. In addition, there are nine less significant Peddrift SNPs (nominal P ≤ 0.000009) on OAR6 (n = 2), OAR9 (n = 2), OAR12, OAR19 (n = 2), OAR24 and OAR26. In smoothed F_ST (five‐SNP moving average) plots, the five most prominent peaks are on OAR1, OAR6, OAR7, OAR13 and OAR19. Although these smoothed F_ST peaks do not coincide with the three most significant Peddrift SNP loci, two (on OAR6 and OAR19) overlap with the set of less significant Peddrift SNPs above. Of these 12 Peddrift SNPs and five smoothed F_ST regions, none is close to the FE candidate genes catalase and ABCG2; however, two on OAR1 and one on OAR13 fall within suggestive quantitative trait locus regions identified in a previous genome screen experiment. The present studies indicated that there are at least eight genomic regions that underwent a selection sweep in the FE lines.     

High connectivity among argali sheep from Afghanistan and adjacent countries: Inferences from neutral and candidate gene microsatellites

We quantified population connectivity and genetic variation in the Marco Polo subspecies of argali mountain sheep (Ovis ammon polii) by genotyping 9 neutral and 8 candidate gene microsatellite loci in 172 individuals noninvasively sampled across five study areas in Afghanistan, China, and Tajikistan. Heterozygosity and allelic richness were generally high (mean H = 0.67, mean A = 6.1), but were significantly lower in the China study area (H = 0.61, P < 0.001; A = 4.9, P < 0.01). One marker in an immune system gene (TCRG4) showed an excess of rare alleles compared to neutral expectations. Another immune system gene (GLYCAM-1) showed excessive differentiation (high F _ST) between study areas. Estimates of genetic differentiation were similar (F _ST = 0.035 vs. 0.033) with and without the two loci deviating from neutrality, suggesting that selection is not a primary driver of overall molecular variation, and that candidate gene loci can be used for connectivity monitoring, as long as selection tests are conducted to avoid biased gene flow estimates. Adequate protection of argali and maintenance of inter-population connectivity will require monitoring and international cooperation because argali exhibit high gene flow across international borders.     

Twelve polymorphic expressed sequence tags‐derived markers for Plasmopara halstedii,the causal agent of sunflower downy mildew

Twelve expressed sequence tags‐derived markers were isolated from Plasmopara halstedii (Oomycetes), the causal agent of sunflower downy mildew. A total of 25 single nucleotide polymorphisms and five indels were detected by single‐strand conformation polymorphism analysis and developed for high‐throughput genotyping of 32 isolates. There was a high level of genetic diversity (H_E = 0.484). Observed heterozygosity ranged from 0 to 0.143 indicating that P. halstedii is probably a selfing species. These markers were also useful in detecting significant genetic variations among French populations (F_ST = 0.193) and between French and Russian populations (F_ST = 0.23). Cross‐amplification tests on three closely related species indicated that no loci amplified in other Oomycete species.     

Low levels of quantitative and molecular genetic differentiation among natural populations of Medicago ciliaris Kroch. (Fabaceae) of different Tunisian eco-geographical origin

In this paper, we analyze the genetic variability in four Tunisian natural populations of Medicago ciliaris using 19 quantitative traits and six polymorphic microsatellite loci. We investigated the amplification transferability of 30 microsatellites developed in the model legume M. truncatula to M. ciliaris. Results revealed that about 56.66% of analyzed markers are valuable genetic markers for M. ciliaris. The most genetic diversity at quantitative traits and microsatellite loci was found to occur within populations (>80%). Low differentiations among populations at quantitative traits Q _ST  = 0.146 and molecular markers F _ST  = 0.18 were found. The majority of measured traits exhibited no significant difference in the level of Q _ST and F _ST . Furthermore, significant correlations established between these traits and eco-geographical factors suggested that natural selection should be invoked to explain the level of phenotypic divergence among populations rather than drift. There was no significant correlation between population differentiation at quantitative traits and molecular markers. Significant spatial genetic structure consistent with models of isolation by distance was detected within all studied populations. The site-of-origin environmental factors explain about 9.07% of total phenotypic genetic variation among populations. The eco-geographical factors that influence more the variation of measured traits among populations are the soil texture and altitude. Nevertheless, there were no consistent pattern of associations between gene diversity (He) and environmental factors.     

Comparative phylogeography of bagrid catfishes in Taiwan

A phylogeographic analysis of two bagrid catfishes in Taiwan was conducted using sequence data from a portion of the mitochondrial DNA (mtDNA) control region. For Pseudobagrus brevianalis, which is most probably endemic to Taiwan, a total of eight haplotypes were detected in 189 specimens from nine river systems covering its entire distribution range, from northern to central western areas of the island. Obvious genetic differentiation was observed among its populations (average F _ST = 0.753); in particular, the northernmost Tamsui River population was fixed for a single endemic haplotype. Nested clade phylogeographic analysis (NCPA) indicated that the dispersal center of mtDNA was the area around the Touchien River and Holong River, north to the Miaoli Plateau, in northwestern Taiwan, suggesting both northward and southward dispersal in this species. There was no evidence for the validity of P. taiwanensis, the nominal species described from Taiwan, morphologically similar to P. brevianalis. We confirmed that P. adiposalis was distributed discontinuously in three river systems; analysis of 42 specimens from the rivers indicated a total of four haplotypes and population differentiation (average F _ST = 0.876). Fixation into a largely differentiated haplotype in the northernmost Tamsui River population was also found in this species, but different processes for this phylogeographic pattern were implied for the two species. Comparison with P. ussuriensis, a widespread continental species morphologically similar to P. adiposalis, suggested the possibility that P. adiposalis is a group of local populations of P. ussuriensis in Taiwan. Two migration routes vs. random fixation scenarios for the population structure of P. adiposalis are discussed with information on other fishes and the geological history of the island.     

Diversity and genetic structure of the endangered cycad Dioon sonorense (Zamiaceae) from Sonora, Mexico: Evolutionary and conservation implications

The diversity and genetic structure of the cycad Dioon sonorense (De Luca, Sabato & Vázq. Torres) Chemnick, T.J. Greg. & S. Salas-Mor. were evaluated from the four known populations in tropical dry forests of Sonora, Mexico. The estimates of genetic diversity in the 19 loci analysed were: mean number of alleles per locus (A = 2.0), percentage of polymorphic loci (P = 81.6) and expected heterozygosity (H_E = 0.314). The subpopulations were found to be in Hardy–Weinberg equilibrium at the local level (F_is), even though we found a global excess of homozygotes (F_it = 0.13). Genetic variation attributed to differences among populations was 15%. Our results indicate that the geographical isolation caused by the historical effects of the Pleistocene among populations has probably generated clinal variation of the allelic frequencies at two loci, in relation to their latitudinal distribution. We compare our results with the long-lived southern hemisphere trees Fitzroya and Nothofagus in relationship to the Pleistocene glaciations. We also infer that the populations have become fragmented recently due to increasing pressure of habitat conversion, disturbance and illegal extraction of plants and seeds for the international horticultural trade. We therefore recommend three aspects: (i) propagation of this species for reintroduction, (ii) declare the localities of its northern distribution as a protected area, and (iii) an amendment of the current International Union for Conservation of Nature Red List categorisation for this species to ‘Critically Endangered’.     

Comparison of quantitative and molecular variation in agroforestry populations of the shea tree (Vitellaria paradoxa C.F. Gaertn) in Mali

In this study we investigated the within- and between-population genetic variation using microsatellite markers and quantitative traits of the shea tree, Vitellaria paradoxa, an important agroforestry tree species of the Sudano–Sahelian region in Africa. Eleven populations were sampled across Mali and in northern Côte d’Ivoire. Leaf size and form and growth traits were measured in a progeny test at the nursery stage. Eight microsatellites were used to assess neutral genetic variation. Low levels of heterozygosity were recorded (1.6–3.0 alleles/locus; H_E = 0.25–0.42) and the fixation index (F_IS = −0.227–0.186) was not significantly different from zero suggesting that Hardy–Weinberg equilibrium is encountered in all populations sampled. Quantitative traits exhibited a strong genetic variation between populations and between families within populations. The degree of population differentiation of the quantitative traits (Q_ST = 0.055–0.283, Q_STmean = 0.189) strongly exceeds that in eight microsatellite loci (F_ST = −0.011–0.142, F_STmean = 0.047). Global and pairwise F_ST values were very low and not significantly different from zero suggesting agroforestry practices are amplifying gene flow (Nm = 5.07). The population means for quantitative traits and the rainfall variable were not correlated, showing variation was not linked with this climatic cline. It is suggested that this marked differentiation for quantitative traits, independent of environmental clines and despite a high gene flow, is a result of local adaptation and human selection of shea trees. This process has induced high linkage disequilibrium between underlying loci of polygenic characters.