Restricted gene flow in fragmented populations of a wind-pollinated tree

Fragmentation of natural populations can have negative effects at the genetic level, thus threatening their evolutionary potential. Many of the negative genetic impacts of population fragmentation can be ameliorated by gene flow and it has been suggested that in wind-pollinated tree species, high or even increased levels of gene flow are a feature of fragmented populations, although several studies have disputed this. We have used a combination of nuclear microsatellites and allele-specific PCR (AS-PCR) analysis of chloroplast single nucleotide polymorphisms (SNPs) to examine the levels and patterns of genetic diversity and population differentiation in fragmented populations of juniper (Juniperus communis) in Ireland and inform conservation programs for the species. Significant population differentiation was found for both chloroplast and nuclear markers, indicating restricted gene flow, particularly over larger geographic scales. For conservation purposes, the existence of genetically distinct clusters and geographically localised chloroplast haplotypes suggests that the concept of provenance should be taken into account when formulating augmentation or reintroduction strategies. Furthermore, the potential lack of seed dispersal and seedling establishment means that ex-situ approaches to seed and seedling management may have to be considered.     

Genetic effects of chronic habitat fragmentation on tree species: the case of Sorbus aucuparia in a deforested Scottish landscape

Sustainable forest restoration and management practices require a thorough understanding of the influence that habitat fragmentation has on the processes shaping genetic variation and its distribution in tree populations. We quantified genetic variation at isozyme markers and chloroplast DNA (cpDNA), analysed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) in severely fragmented populations of Sorbus aucuparia (Rosaceae) in a single catchment (Moffat) in southern Scotland. Remnants maintain surprisingly high levels of gene diversity (HE) for isozymes (HE = 0.195) and cpDNA markers (HE = 0.490). Estimates are very similar to those from non-fragmented populations in continental Europe, even though the latter were sampled over a much larger spatial scale. Overall, no genetic bottleneck or departures from random mating were detected in the Moffat fragments. However, genetic differentiation among remnants was detected for both types of marker (isozymes Theta n = 0.043, cpDNA Theta c = 0.131; G-test, P-value < 0.001). In this self-incompatible, insect-pollinated, bird-dispersed tree species, the estimated ratio of pollen flow to seed flow between fragments is close to 1 (r = 1.36). Reduced pollen-mediated gene flow is a likely consequence of habitat fragmentation, but effective seed dispersal by birds is probably helping to maintain high levels of genetic diversity within remnants and reduce genetic differentiation between them.     

Contact zone of two different chloroplast lineages and genetic guidelines for seed transfer in Quercus serrata and Quercus crispula

Within their natural distribution ranges, plant species exhibit a genetic structure that has been created by global climate change and natural selection over long periods. This genetic structure needs to be conserved for sustainable use of genetic resources. To conserve local forests with different genetic structures, genetic guidelines for seed and seedling transfer in individual species are necessary. Genetic guidelines have been published for 43 Japanese tree species using population genetic data; however, for practical use, more detailed genetic borders between important genetic lineages should be clarified to inform seed collection and planting. Thus, we investigated in detail the genetic borders between two important Japanese oak species, Quercus serrata and Quercus crispula, in the Chubu region of Japan using chloroplast and nuclear DNA markers, and we discuss the factors that influenced border creation using the results of species distribution modeling (SDM). The chloroplast DNA (cpDNA) haplotype was clearly different within the Chubu region of Japan but the difference in nuclear DNA between northern and southern haplotype populations was very small, both in Q. serrata and Q. crispula. The results of SDM showed that during the last glacial maximum (LGM) Q. serrata was distributed mostly along the coastline but Q. crispula was distributed not only along the coast but also in mountainous areas further inland. The cpDNA genetic borders of these two oak species are complex and seem to have been influenced by topography and their distribution during the LGM. We propose and discuss genetic guidelines for these two oak species based on the results of this study.     

Patterns of seed dispersal and pollen flow in Quercus garryana (Fagaceae) following post-glacial climatic changes

Aim We examined the genetic structure of Quercus garryana to infer post‐glacial patterns of seed dispersal and pollen flow to test the hypotheses that (1) peripheral populations are genetically distinct from core populations and from one another; (2) genetic diversity declines towards the poleward edge of the species’ range; and (3) genetic diversity in the chloroplast genome, a direct measure of seed dispersal patterns, declines more sharply with increasing latitude than diversity in the nuclear genome. We address our findings in the context of known historical oak distribution from pollen core data derived from previously published research. Location Oak–savanna ecosystems from southern Oregon, USA (core populations/non‐glaciated range) northward to Vancouver Island, British Columbia, Canada (peripheral populations/glaciated range). Methods We genotyped 378 trees from 22 sites with five chloroplast and seven nuclear microsatellite loci. For both sets of markers, we estimated genetic diversity and differentiation using an analysis of molecular variance and generated Mantel correlograms to detect genetic and geographical distance correlations. For the nuclear markers, we also used a Bayesian approach to infer population substructure. Results There was a large degree of population differentiation revealed by six chloroplast haplotypes, with little (≤ 3) or no haplotype diversity within sites. Peripheral island locations shared the same, maternally inherited chloroplast haplotype, whereas locations in mainland Washington had greater haplotype diversity. In contrast, genetic diversity of the nuclear markers was high at all locations sampled. Populations clustered into two groups and were significantly positively correlated over large spatial scales (≤ 200 km), although allele richness decreased significantly with latitude. Population substructure was observed between core and peripheral populations because rare alleles were absent in peripheral localities and common allele frequencies differed. Main conclusions The observed pattern of chloroplast haplotype loss at the northern periphery suggests restricted seed dispersal events from mainland sites to peripheral islands. This pattern was unexpected, however, as refugial oak populations remained near the current post‐glacial range even during the Last Glacial Maximum. Using nuclear markers, we found high within‐population diversity and population differentiation only over large spatial scales, suggesting that pollen flow is relatively high among populations.     

Chloroplast DNA variation and population structure in the widespread forest tree, Eucalyptus grandis

Recognition of genetic structure of populations and the ability to identify vulnerable populations is useful for the formation of conservation management strategies for plants. Eucalyptus grandis is a tall forest tree that has a major area of occurrence in subtropical eastern Australia, with smaller populations located in the east coast tropics. Many widespread forest species exhibit population differentiation that corresponds to geographic regions. However, Eucalyptus grandis appears to be an exception based on isozyme and morphological data. This is intriguing given a large discontinuity between northern populations and those in the southern part of the species range. In this study, the distribution of a maternally inherited chloroplast locus was examined because it was more likely to reveal genetic structure due to the slower evolution of the chloroplast genome and limited dispersal of seed in eucalypts. As expected, the G _ST for chloroplast DNA was higher than that for nuclear DNA but indicated low population differentiation for a forest tree species. Phylogeographic analysis indicated that the 15 populations grouped into three broad geographical regions; however, overall population structure was weak suggesting that the large geographical disjunction in the distribution of E. grandis may be relatively recent. A paradigm for conservation management of E. grandis based on chloroplast DNA haplotype distribution would take into account the low differentiation among populations.     

Testing species delimitation in sympatric species complexes: the case of an African tropical tree, Carapa spp. (Meliaceae)

Plant species delimitation within tropical ecosystems is often difficult because of the lack of diagnostic morphological characters that are clearly visible. The development of an integrated approach, which utilizes several different types of markers (both morphological and molecular), would be extremely useful in this context. Here we have addressed species delimitation of sympatric tropical tree species that belong to Carapa spp. (Meliaceae) in Central Africa. We adopted a population genetics approach, sampling numerous individuals from three locations where sympatric Carapa species are known to exist. Comparisons between morphological markers (the presence or absence of characters, leaf-shape traits) and molecular markers (chloroplast sequences, ribosomal internal transcribed spacer region (ITS) sequences, and nuclear microsatellites) demonstrated the following: (i) a strong correlation between morphological and nuclear markers; (ii) despite substantial polymorphism, the inability of chloroplast DNA to discriminate between species, suggesting that cytoplasmic markers represent ineffective DNA barcodes; (iii) lineage sorting effects when using ITS sequences; and (iv) a complex evolutionary history within the genus Carapa, which includes frequent inter-specific gene flow. Our results support the use of a population genetics approach, based on ultra-polymorphic markers, to address species delimitation within complex taxonomic groups.     

Late Pleistocene speciation of three closely related tree peonies endemic to the Qinling–Daba Mountains,a major glacial refugium in Central China

Determining the factors promoting speciation is a major task in ecological and evolutionary research and can be aided by phylogeographic analysis. The Qinling–Daba Mountains (QDM) located in central China form an important geographic barrier between southern subtropical and northern temperate regions, and exhibit complex topography, climatic, and ecological diversity. Surprisingly, few phylogeographic analyses and studies of plant speciation in this region have been conducted. To address this issue, we investigated the genetic divergence and evolutionary histories of three closely related tree peony species (Paeonia qiui, P. jishanensis, and P. rockii) endemic to the QDM. Forty populations of the three tree peony species were genotyped using 22 nuclear simple sequence repeat markers (nSSRs) and three chloroplast DNA sequences to assess genetic structure and phylogenetic relationships, supplemented by morphological characterization and ecological niche modeling (ENM). Morphological and molecular genetic analyses showed the three species to be clearly differentiated from each other. In addition, coalescent analyses using DIYABC conducted on nSSR variation indicated that the species diverged from each other in the late Pleistocene, while ecological niche modeling (ENM) suggested they occupied a larger area during the Last Glacial Maximum (LGM) than at present. The combined genetic evidence from nuclear and chloroplast DNA and the results of ENM indicate that each species persisted through the late Pleistocene in multiple refugia in the Qinling, Daba, and Taihang Mountains with divergence favored by restricted gene flow caused by geographic isolation, ecological divergence, and limited pollen and seed dispersal. Our study contributes to a growing understanding of the origin and population structure of tree peonies and provides insights into the high level of plant endemism present in the Qinling–Daba Mountains of Central China.     

Genetic diversity and gene flow in a Caribbean tree <Emphasis Type="Italic">Pterocarpus officinalis</Emphasis> Jacq.: a study based on chloroplast and nuclear microsatellites

We analysed the molecular diversity of Pterocarpus officinalis, a tree species distributed in Caribbean islands, South and Central America to quantify the genetic variation within island, to assess the pattern of differentiation and infer levels of gene flow; with the overall goal of defining a strategy of conservation. Two hundred two individuals of 9 populations were analysed using three chloroplast and six nuclear microsatellite markers. The observed heterozygosity varied markedly among the populations for nuclear (H _Onuc = 0.20–0.50) and chloroplast microsatellites (H _cp = 0.22–0.68). The continental population from French Guyana showed a higher value of H _Onuc than island populations, and the differences were significant in some cases. The fixation index F _IS ranged from −0.043 to 0.368; a significant heterozygote deficit was detected in 7 populations. The heterozygosity excess method suggested that two populations in Guadeloupe have undergone a recent bottleneck. Global and pairwise F _ST were high for both nuclear (F _STnuc = 0.29) and chloroplast microsatellites (F _STcp = 0.58). The neighbour-joining tree based on both markers, presented a differentiation pattern that can be explained by the seed dispersal by flotation and marine stream. The comparison of Bayesian approach and the method based on allelic frequency demonstrate a very limited number of migrants between populations.     

The impact of pollination syndrome and habitat on gene flow: a comparative study of two Streptocarpus (Gesneriaceae) species

Gene flow through pollen and seed dispersal is important in terms of population differentiation and eventually speciation. Seed and pollen flow are affected in turn by habitats and pollen vectors. We examined the effect of different pollinators and habitats on gene flow by comparing two species of Streptocarpus, using microsatellite and chloroplast RFLP markers. Populations of the forest-dwelling S. primulifolius were highly differentiated according to nuclear microsatellite data and had mutually exclusive chloroplast haplotypes. This result is congruent with infrequent seed dispersal and limited between-population foraging by the long-tongued fly pollinator Stenobasipteron wiedemanni. In contrast, populations of S. dunnii growing in exposed crags had lower levels of population differentiation according to both nuclear and chloroplast data, congruent with a hypothesis of more effective between population seed dispersal and greater pollen-mediated gene flow due to the sunbird pollinator Nectarinia famosa. The population genetic behavior of these species is reflected in their taxonomy and phylogenetic position; S. primulifolius belongs to a taxonomically complex clade in which recent speciation is evident, while the clade containing S. dunnii is characterized by taxonomically well-defined species on longer phylogenetic branches. Our study shows that pollinator movements and seed dispersal patterns are a major determinant of the evolutionary trajectories of these species.     

Conserving the small milkwort, <Emphasis Type="Italic">Comesperma polygaloides</Emphasis>, a vulnerable subshrub in a fragmented landscape

The conservation of remnant grassland vegetation on the Victorian volcanic plain (VVP) is crucial for the persistence of local biodiversity. Recent habitat loss has restricted the grassland to only a small percentage of its former range. Along with grassland habitats, species that occur on the VVP are in decline and many are legally protected. Comesperma polygaloides is a grassland species of the VVP that also occurs outside of the region in woodland habitats. We use 12 neutral microsatellite loci and two chloroplast regions to understand genotypic patterns of C. polygaloides in southeastern Australia. We found separate genetic clusters but they do not follow geographic boundaries. There are fewer alleles (2.96) and effective alleles (2.01) than expected from 12 microsatellite markers compared to other species. Even with the low number of alleles per locus there was a moderate level of genetic diversity detected (I = 0.69; H_o = 0.43; H_e = 0.40). Populations of the VVP could not be differentiated from populations elsewhere using neutral markers or chloroplast analyses. The genetic structure discovered was not consistent with the level of fragmentation observed. There may be several reasons for the observed lack of genetic structure: the species is more common than perceived, plants are long-lived and can reproduce clonally, and the bioregion is relatively young, geologically. Results indicate that restoration projects and long-term viability of C. polygaloides will be improved by composite seed sourcing, alleviating the risk of insufficient genetic diversity posed by an over-emphasis on local provenancing.     

Increased inbreeding and strong kinship structure in Taxus baccata estimated from both AFLP and SSR data

Habitat fragmentation can have severe genetic consequences for trees, such as increased inbreeding and decreased effective population size. In effect, local populations suffer from reduction of genetic variation, and thus loss of adaptive capacity, which consequently increases their risk of extinction. In Europe, Taxus baccata is among a number of tree species experiencing strong habitat fragmentation. However, there is little empirical data on the population genetic consequences of fragmentation for this species. This study aimed to characterize local genetic structure in two natural remnants of English yew in Poland based on both amplified fragment length polymorphism (AFLP) and microsatellite (SSR) markers. We introduced a Bayesian approach that estimates the average inbreeding coefficient using AFLP (dominant) markers. Results showed that, in spite of high dispersal potential (bird-mediated seed dispersal and wind-mediated pollen dispersal), English yew populations show strong kinship structure, with a spatial extent of 50–100 m, depending on the population. The estimated inbreeding levels ranged from 0.016 to 0.063, depending on the population and marker used. Several patterns were evident: (1) AFLP markers showed stronger kinship structure than SSRs; (2) AFLP markers provided higher inbreeding estimates than SSRs; and (3) kinship structure and inbreeding were more pronounced in denser populations regardless of the marker used. Our results suggest that, because both kinship structure and (bi-parental) inbreeding exist in populations of English yew, gene dispersal can be fairly limited in this species. Furthermore, at a local scale, gene dispersal intensity can be more limited in a dense population.     

Characterization of chloroplast microsatellite loci from whole chloroplast genome of Camellia taliensis and their utilization for evaluating genetic diversity of Camellia reticulata (Theaceae)

This study characterized chloroplast microsatellite markers for Camellia reticulata, a famous ornamental and edible economic tree species only distributed in Southwestern China. Thirty-two chloroplast microsatellite markers were originally annotated in the whole chloroplast genome of Camellia taliensis, ten polymorphic microsatellite markers were tested in 96 individuals from four natural populations of C. reticulata. Alleles numbered from two to four, and average value of Shannon's Information index, Nei's gene diversity, total genetic diversity, genetic diversity within populations, gene differentiation coefficient and gene flow index were 0.408, 0.225, 0.217, 0.102, 0.530 and 0.443, respectively. Our results showed high genetic differentiation and limited gene flow among the studied populations, which may be explained by recent fragmentation of the remnant populations due to human destruction and disturbance of natural habitats of the species. The identified cpSSR markers will be useful for the future studies on population genetics, conservation biology and phylogeography of C. reticulata.     

Analysis of microsatellite loci in tree of heaven (<Emphasis Type="Italic">Ailanthus altissima</Emphasis> (Mill.) Swingle) using SSR-GBS

Microsatellite markers are still the marker of choice for many research questions in the field of forest genetics. However, the number of available markers is often low for species that have not been studied intensively like the tree of heaven (Ailanthus altissima). During the last decade, next-generation sequencing (NGS) has offered advanced techniques for efficiently identifying microsatellite markers and accurately genotyping samples. Here, we identify new microsatellite markers for the tree of heaven by applying an NGS-based method using the Illumina MiSeq platform. NGS technology was proved to be an effective method for fast and cost-efficient identification of microsatellite markers by implementing a genotyping-by-sequencing approach based on Illumina amplicon sequencing (SSR-GBS). We screened three populations from Eastern Austria for genetic variation at 19 newly identified microsatellite loci. We tested two different genotyping approaches: (1) considering only allele lengths (forming a so-called “allele length dataset”), (2) taking also single nucleotide polymorphisms (SNPs) within the amplified fragments into account (forming a so-called “SNP dataset”). The results revealed higher values for all genetic diversity parameters, as well as a better resolution of genetic assignment, when the latter approach was followed. Thus, by taking advantage of sequence information which is provided by SSR-GBS, one may achieve considerable gains in performance using the same marker set. The developed markers provide a cost-efficient tool for genotyping populations of tree of heaven and the approach presented here promises to be of high value for medium throughput genotyping applications in non-model forest tree species. We will use this method to widen the perspectives for further population genetic investigations of the tree of heaven.     

Molecular,quantitative and abiotic variables for the delineation of evolutionary significant units: case of sandalwood (<Emphasis Type="Italic">Santalum austrocaledonicum</Emphasis> Vieillard) in New Caledonia

Various approaches have been developed to define conservation units for plant and animal species. In this study we combined nuclear microsatellites (from a previous published study) and chloroplast microsatellites (assessed in the present study), leaf and seed morphology traits and abiotic variables (climate and soil) to define evolutionary significant units (ESU) of Santalum austrocaledonicum, a tree species growing in New Caledonia. Results for chloroplast microsatellites showed that the total population heterozygosity was␣high, (H _cp = 0.84) but varied between islands. Differentiation was strong in the total population (F _stcp = 0.66) but also within the main island Grande Terre (F _stcp = 0.73) and within Iles Loyauté (F _stcp = 0.52), highlighting a limited gene flow between populations. These results confirmed those obtained with nuclear microsatellites. The cluster analysis on molecular markers discriminated two main groups constituted by the populations of Grande Terre and the populations of Iles Loyauté. A principal component analysis of leaf and seed morphology traits singled out the populations of Iles Loyauté and the western populations of Grande Terre. Quantitative genetic analyses showed that the variation between populations was under genetic control (broad sense heritability close to 80%). A high correlation between rainfall and morphological traits suggested an impact of climate on this variation. The integration of these results allows to define two ESUs, one corresponding to Grande Terre and Ile des Pins and the other the Iles Loyauté archipelago. This study stresses the need to restore some populations of Grande Terre that are currently threatened by their small size.     

Multiscale spatial genetic structure within and between populations of wild cherry trees in nuclear genotypes and chloroplast haplotypes

Spatial genetic structure (SGS) of plants mainly depends on the effective population size and gene dispersal. Maternally inherited loci are expected to have higher genetic differentiation between populations and more intensive SGS within populations than biparentally inherited loci because of smaller effective population sizes and fewer opportunities of gene dispersal in the maternally inherited loci. We investigated biparentally inherited nuclear genotypes and maternally inherited chloroplast haplotypes of microsatellites in 17 tree populations of three wild cherry species under different conditions of tree distribution and seed dispersal. As expected, interpopulation genetic differentiation was 6–9 times higher in chloroplast haplotypes than in nuclear genotypes. This difference indicated that pollen flow 4–7 times exceeded seed flow between populations. However, no difference between nuclear and chloroplast loci was detected in within‐population SGS intensity due to their substantial variation among the populations. The SGS intensity tended to increase as trees became more aggregated, suggesting that tree aggregation biased pollen and seed dispersal distances toward shorter. The loss of effective seed dispersers, Asian black bears, did not affect the SGS intensity probably because of mitigation of the bear loss by other vertebrate dispersers and too few tree generations after the bear loss to alter SGS. The findings suggest that SGS is more variable in smaller spatial scales due to various ecological factors in local populations.     

Limited structure and widespread diversity suggest potential buffers to genetic erosion in a threatened grassland shrub Pimelea spinescens (Thymelaeaceae)

Pimelea spinescens is a critically endangered species of the temperate grasslands of southeastern Australia. Two subspecies are recognised. Subspecies, P. spinescens subsp. spinescens, formerly common and widespread, is found in isolated remnants of previously extensive grasslands. The second subspecies, subsp. pubiflora, is thought to have been historically rare with only two geographically-isolated extant populations. The grassland communities exist now as fragmented remnants representing <1 % of their extent prior to European settlement in the early 1800s. Conservation management strategies for species in these critically endangered ecosystems rely on an understanding of genetic diversity and population structure to ensure long-term evolutionary potential. We used chloroplast DNA (cpDNA) and microsatellite markers to examine the population genetic structure of both subspecies. Analysis of cpDNA revealed 14 haplotypes with high divergence between the single haplotype found in subsp. pubiflora and most remaining haplotypes restricted to subsp. spinescens. Microsatellites also indicated high genetic differentiation between subspecies but little evidence of sub-structuring within either subspecies. Results suggest that seed dispersal has not been as limited as previously thought. In this fragmented habitat, a lack of genetic structure suggests buffers to genetic erosion, with current patterns reflecting genetic diversity prior to fragmentation. Plant longevity and the presence of seed banks may contribute to the maintenance of these patterns, resulting in a lag between fragmentation and genetic erosion. Whilst factors such as longevity and seed banks may be preserving historic genetic diversity, management is required to ensure the maintenance of this diversity into the future.     

Origins and genetic conservation of tropical trees in agroforestry systems: a case study from the Peruvian Amazon

Hundreds of native tree species are currently found in extensive agroforestry ecosystems in the Peruvian Amazon, forming an important reservoir of biodiversity. To further promote conservation, farmers are encouraged to supplement intra-specific genetic diversity in these populations with seed collected from local forests. For some tree species, however, this approach may be inappropriate, as stands of these taxa already found on-farm may not be of local origin. Despite this issue being of importance for conservation, little information is available on the history of cultivated trees in the region, a situation that we here rectify for the important fruit tree Inga edulis. Based on nuclear SSR and chloroplast DHPLC analyses of closely geographically matched natural and planted stands at five sites, it appears that cultivated material of I. edulis is primarily of non-local origin, indicating that conservation based on new wide-scale infusions from local wild stands into farms may be inappropriate in the region. Although nuclear and chloroplast diversity were both lower in planted stands, values were still relatively high (∼80 and 70% of natural stands, respectively), indicating that when farmers plant trees, good collection practice of seed from already cultivated I. edulis should be an effective means for ensuring long-term conservation on farms.     

Standardized genetic diversity‐life history correlates for improved genetic resource management of Neotropical trees

Aim

Life history traits and range size are key correlates of genetic diversity in trees. We used a standardized sampling protocol to explore how life history traits and range size relate to the magnitude, variance and structuring (both between‐ and within‐population) of genetic diversity in Neotropical tree species.

Location

The Neotropics

Methods

We present a meta‐analysis of new population genetic data generated for 23 Neotropical tree species (=2,966 trees, 86 populations) across a shared and broad geographic area. We compared established population genetic metrics across these species (e.g., genetic diversity, population structure, fine‐scale genetic structure), plus we estimated the rarely used variance in genetic diversity among populations. We used a multivariate, maximum likelihood, multimodel inference approach to explore the relative influence of life history traits and range size on patterns of neutral genetic diversity.

Results

We found that pioneer and narrow range species had lower levels but greater variance in genetic diversity—signs of founder effects and stronger genetic drift. Animal‐dispersed species had lower population differentiation, indicating extensive gene flow. Abiotically dispersed and pioneer species had stronger fine‐scale genetic structure, suggesting restricted seed dispersal and family cohort establishment.

Main conclusions

Our multivariable and multispecies approach allows ecologically relevant conclusions, since knowing whether one parameter has an effect, or one species shows a response in isolation, is dependent on the combination of traits expressed by a species. Our study demonstrates the influence of ecological processes on the distribution of genetic variation in tropical trees, and will help guide genetic resource management, and contribute to predicting the impacts of land use change.