Range‐wide analysis of genetic structure in a widespread,highly mobile species (Odocoileus hemionus) reveals the importance of historical biogeography

Highly mobile species that thrive in a wide range of habitats are expected to show little genetic differentiation across their range. A limited but growing number of studies have revealed that patterns of broad‐scale genetic differentiation can and do emerge in vagile, continuously distributed species. However, these patterns are complex and often shaped by both historical and ecological factors. Comprehensive surveys of genetic variation at a broad scale and at high resolution are useful for detecting cryptic spatial genetic structure and for investigating the relative roles of historical and ecological processes in structuring widespread, highly mobile species. In this study, we analysed 10 microsatellite loci from over 1900 samples collected across the full range of mule deer (Odocoileus hemionus), one of the most widely distributed and abundant of all large mammal species in North America. Through both individual‐ and population‐based analyses, we found evidence for three main genetic lineages, one corresponding to the ‘mule deer’ morphological type and two to the ‘black‐tailed deer’ type. Historical biogeographic events likely are the primary drivers of genetic divergence in this species; boundaries of the three lineages correspond well with predictions based on Pleistocene glacial cycles, and substructure within each lineage demonstrates island vicariance. However, across large geographic areas, including the entire mule deer lineage, we found that genetic variation fit an isolation‐by‐distance pattern rather than discrete clusters. A lack of genetic structure across wide geographic areas of the continental west indicates that ecological processes have not resulted in restrictions to gene flow sufficient for spatial genetic structure to emerge. Our results have important implications for our understanding of evolutionary mechanisms of divergence, as well as for taxonomy, conservation and management.     

Balancing genetic uniqueness and genetic variation in determining conservation and translocation strategies: a comprehensive case study of threatened dwarf galaxias,Galaxiella pusilla (Mack) (Pisces: Galaxiidae)

Genetic markers are widely used to define and manage populations of threatened species based on the notion that populations with unique lineages of mtDNA and well‐differentiated nuclear marker frequencies should be treated separately. However, a danger of this approach is that genetic uniqueness might be emphasized at the cost of genetic diversity, which is essential for adaptation and is potentially boosted by mixing geographically separate populations. Here, we re‐explore the issue of defining management units, focussing on a detailed study of Galaxiella pusilla, a small freshwater fish of national conservation significance in Australia. Using a combination of microsatellite and mitochondrial markers, 51 populations across the species range were surveyed for genetic structure and diversity. We found an inverse relationship between genetic differentiation and genetic diversity, highlighting a long‐term risk of deliberate isolation of G. pusilla populations based on protection of unique lineages. Instead, we adopt a method for identifying genetic management units that takes into consideration both uniqueness and genetic variation. This produced a management framework to guide future translocation and re‐introduction efforts for G. pusilla, which contrasted to the framework based on a more traditional approach that may overlook important genetic variation in populations.     

Habitat predictors of genetic diversity for two sympatric wetland‐breeding amphibian species

Population genetic diversity is widely accepted as important to the conservation and management of wildlife. However, habitat features may differentially affect evolutionary processes that facilitate population genetic diversity among sympatric species. We measured genetic diversity for two pond‐breeding amphibian species (Dwarf salamanders, Eurycea quadridigitata; and Southern Leopard frogs, Lithobates sphenocephalus) to understand how habitat characteristics and spatial scale affect genetic diversity across a landscape. Samples were collected from wetlands on a longleaf pine reserve in Georgia. We genotyped microsatellite loci for both species to assess population structures and determine which habitat features were most closely associated with observed heterozygosity and rarefied allelic richness. Both species exhibited significant population genetic structure; however, structure in Southern Leopard frogs was driven primarily by one outlier site. Dwarf salamander allelic richness was greater at sites with less surrounding road area within 0.5 km and more wetland area within 1.0 and 2.5 km, and heterozygosity was greater at sites with more wetland area within 0.5 km. In contrast, neither measure of Southern Leopard frog genetic diversity was associated with any habitat features at any scale we evaluated. Genetic diversity in the Dwarf salamander was strongly associated with land cover variables up to 2.5 km away from breeding wetlands, and/or results suggest that minimizing roads in wetland buffers may be beneficial to the maintenance of population genetic diversity. This study suggests that patterns of genetic differentiation and genetic diversity have associations with different habitat features across different spatial scales for two syntopic pond‐breeding amphibian species.     

A trait‐based approach to predict population genetic structure in bees

Understanding population genetic structure is key to developing predictions about species susceptibility to environmental change, such as habitat fragmentation and climate change. It has been theorized that life‐history traits may constrain some species in their dispersal and lead to greater signatures of population genetic structure. In this study, we use a quantitative comparative approach to assess if patterns of population genetic structure in bees are driven by three key species‐level life‐history traits: body size, sociality, and diet breadth. Specifically, we reviewed the current literature on bee population genetic structure, as measured by the differentiation indices Nei's G_ST, Hedrick's G′_ST, and Jost's D. We then used phylogenetic generalised linear models to estimate the correlation between the evolution of these traits and patterns of genetic differentiation. Our analyses revealed a negative and significant effect of body size on genetic structure, regardless of differentiation index utilized. For Hedrick's G′_ST and Jost's D, we also found a significant impact of sociality, where social species exhibited lower levels of differentiation than solitary species. We did not find an effect of diet specialization on population genetic structure. Overall, our results suggest that physical dispersal or other functions related to body size are among the most critical for mediating population structure for bees. We further highlight the importance of standardizing population genetic measures to more easily compare studies and to identify the most susceptible species to landscape and climatic changes.     

Interspecific interactions influence contrasting spatial genetic structures in two closely related damselfly species

Spatial genetic structure (SGS) is largely determined by colonization history, landscape and ecological characteristics of the species. Therefore, sympatric and ecologically similar species are expected to exhibit similar SGSs, potentially enabling prediction of the SGS of one species from that of another. On the other hand, due to interspecific interactions, ecologically similar species could have different SGSs. We explored the SGSs of the closely related Calopteryx splendens and Calopteryx virgo within Finland and related the genetic patterns to characteristics of the sampling localities. We observed different SGSs for the two species. Genetic differentiation even within short distances in C. splendens suggests genetic drift as an important driver. However, we also observed indication of previous gene flow (revealed by a negative relationship between genetic differentiation and increasing potential connectivity of the landscape). Interestingly, genetic diversity of C. splendens was negatively related to density of C. virgo, suggesting that interspecific interactions influence the SGS of C. splendens. In contrast, genetic differentiation between C. virgo subpopulations was low and only exhibited relationships with latitude, pointing to high gene flow, colonization history and range margin effects as the drivers of SGS. The different SGSs of the two ecologically similar species caution indirect inferences of SGS based on ecologically similar surrogate species.     

Range‐wide population genetic structure of the Caribbean sea fan coral,Gorgonia ventalina

The population structure of benthic marine organisms is of central relevance to the conservation and management of these often threatened species, as well as to the accurate understanding of their ecological and evolutionary dynamics. A growing body of evidence suggests that marine populations can be structured over short distances despite theoretically high dispersal potential. Yet the proposed mechanisms governing this structure vary, and existing empirical population genetic evidence is of insufficient taxonomic and geographic scope to allow for strong general inferences. Here, we describe the range‐wide population genetic structure of an ecologically important Caribbean octocoral, Gorgonia ventalina. Genetic differentiation was positively correlated with geographic distance and negatively correlated with oceanographically modelled dispersal probability throughout the range. Although we observed admixture across hundreds of kilometres, estimated dispersal was low, and populations were differentiated across distances <2 km. These results suggest that populations of G. ventalina may be evolutionarily coupled via gene flow but are largely demographically independent. Observed patterns of differentiation corroborate biogeographic breaks found in other taxa (e.g. an east/west divide near Puerto Rico), and also identify population divides not discussed in previous studies (e.g. the Yucatan Channel). High genotypic diversity and absence of clonemates indicate that sex is the primary reproductive mode for G. ventalina. A comparative analysis of the population structure of G. ventalina and its dinoflagellate symbiont, Symbiodinium, indicates that the dispersal of these symbiotic partners is not coupled, and symbiont transmission occurs horizontally.     

Sampling related individuals within ponds biases estimates of population structure in a pond‐breeding amphibian

Effective conservation and management of pond‐breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding‐site connectivity. Population‐level studies of pond‐breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within‐pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome‐wide SNPs generated using a double‐digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within‐pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.     

Recent phylogeographic structure in a widespread ‘weedy’ Neotropical tree species,Cordia alliodora (Boraginaceae)

Aim Although hundreds of tree species have broad geographic ranges in the Neotropics, little is known about how such widespread species attained disjunct distributions around mountain, ocean and xeric barriers. Here, we examine the phylogeographic structure of a widespread and economically important tree, Cordia alliodora, to: (1) test the roles of vicariance and dispersal in establishing major range disjunctions, (2) determine which geographic regions and/or habitats contain the highest levels of genetic diversity, and (3) infer the geographic origin of the species. Location Twenty‐five countries in Central and South America, and the West Indies. Methods Chloroplast simple sequence repeats (cpSSR; eight loci) were assayed in 67 populations (240 individuals) sampled from the full geographic range of C. alliodora. Chloroplast (trnH–psbA) and nuclear (internal transcribed spacer, ITS) DNA sequences were sampled from a geographically representative subset. Genetic structure was determined with samova , structure and haplotype networks. Analysis of molecular variance (AMOVA) and rarefaction analyses were used to compare regional haplotype diversity and differentiation. Results Although the ITS region was polymorphic it revealed limited phylogeographic structure, and trnH–psbA was monomorphic. However, structure analysis of cpSSR variation recovered three broad demes spanning Central America (Deme 1), the Greater Antilles and the Chocó (Deme 2), and the Lesser Antilles and cis‐Andean South America (Deme 3). samova showed two predominant demes (Deme 1 + 2 and Deme 3). The greatest haplotype diversity was detected east of the Andes, while significantly more genetic variation was partitioned among trans‐Andean populations. Populations experiencing high precipitation seasonality (dry ecotype) had greater levels of genetic variation. Main conclusions Cordia alliodora displayed weak cis‐ and trans‐Andean phylogeographic structure based on DNA sequence data, indicative of historical dispersal around this barrier and genetic exchange across its broad range. The cpSSR data revealed phylogeographic structure corresponding to three biogeographic zones. Patterns of genetic diversity are indicative of an origin in the seasonally dry habitats of South America. Therefore, C. alliodora fits the disperser hypothesis for widespread Neotropical species. Dispersal is evident in the West Indies and the northern Andean cordilleras. The dry ecotype harbours genetic variation that is likely to represent the source for the establishment of populations under future warmer and drier climatic scenarios.     

Next‐generation sequencing (NGS)‐based identification of induced mutations in a doubly mutagenized tomato (Solanum lycopersicum) population

The identification of mutations in targeted genes has been significantly simplified by the advent of TILLING (Targeting Induced Local Lesions In Genomes), speeding up the functional genomic analysis of animals and plants. Next‐generation sequencing (NGS) is gradually replacing classical TILLING for mutation detection, as it allows the analysis of a large number of amplicons in short durations. The NGS approach was used to identify mutations in a population of Solanum lycopersicum (tomato) that was doubly mutagenized by ethylmethane sulphonate (EMS). Twenty‐five genes belonging to carotenoids and folate metabolism were PCR‐amplified and screened to identify potentially beneficial alleles. To augment efficiency, the 600‐bp amplicons were directly sequenced in a non‐overlapping manner in Illumina MiSeq, obviating the need for a fragmentation step before library preparation. A comparison of the different pooling depths revealed that heterozygous mutations could be identified up to 128‐fold pooling. An evaluation of six different software programs (camba , crisp , gatk unified genotyper , lofreq , snver and vipr ) revealed that no software program was robust enough to predict mutations with high fidelity. Among these, crisp and camba predicted mutations with lower false discovery rates. The false positives were largely eliminated by considering only mutations commonly predicted by two different software programs. The screening of 23.47 Mb of tomato genome yielded 75 predicted mutations, 64 of which were confirmed by Sanger sequencing with an average mutation density of 1/367 Kb. Our results indicate that NGS combined with multiple variant detection tools can reduce false positives and significantly speed up the mutation discovery rate.     

Unexpected cryptic species diversity in the widespread coral Seriatopora hystrix masks spatial‐genetic patterns of connectivity

Mounting evidence of cryptic species in a wide range of taxa highlights the need for careful analyses of population genetic data sets to unravel within‐species diversity from potential interspecies relationships. Here, we use microsatellite loci and hierarchical clustering analysis to investigate cryptic diversity in sympatric and allopatric (separated by 450 km) populations of the widespread coral Seriatopora hystrix on the Great Barrier Reef. Structure analyses delimited unique genetic clusters that were confirmed by phylogenetic and extensive population‐level analyses. Each of four sympatric yet distinct genetic clusters detected within S. hystrix demonstrated greater genetic cohesion across regional scales than between genetic clusters within regions (<10 km). Moreover, the magnitude of genetic differentiation between different clusters (>0.620 G”_ST) was similar to the difference between S. hystrix clusters and the congener S. caliendrum (mean G”_ST 0.720). Multiple lines of evidence, including differences in habitat specificity, mitochondrial identity, Symbiodinium associations and morphology, corroborate the nuclear genetic evidence that these distinct clusters constitute different species. Hierarchical clustering analysis combined with more traditional population genetic methods provides a powerful approach for delimiting species and should be regularly applied to ensure that ecological and evolutionary patterns interpreted for single species are not confounded by the presence of cryptic species.     

Expanding the population genetic perspective of cnidarian‐Symbiodinium symbioses

The modern synthesis was a seminal period in the biological sciences, establishing many of the core principles of evolutionary biology that we know today. Significant catalysts were the contributions of R.A. Fisher, J.B.S. Haldane and Sewall Wright (and others) developing the theoretical underpinning of population genetics, thus demonstrating adaptive evolution resulted from the interplay of forces such as natural selection and mutation within groups of individuals occupying the same space and time (i.e. a population). Given its importance, it is surprising that detailed population genetic data remain lacking for numerous organisms vital to many ecosystems. For example, the coral reef ecosystem is well recognized for its high biodiversity and productivity, numerous ecological services and significant economic and societal values (Moberg & Folke 1999; Cinner 2014). Many coral reef invertebrates form symbiotic relationships with single‐celled dinoflagellates within the genus Symbiodinium Freudenthal (Taylor 1974), with hosts providing these (typically) intracellular symbionts with by‐products of metabolism and in turn receiving photosynthetically fixed carbon capable of meeting hosts’ respiratory demands (Falkowski et al. 1984; Muscatine et al. 1984). Unfortunately, the health and integrity of the coral reef ecosystem has been significantly and negatively impacted by onslaughts like anthropogenic eutrophication and disease in addition to global climate change, with increased incidences of ‘bleaching’ events (characterized as the loss of photosynthetic pigments from the algal cell or massive reduction of Symbiodinium density from hosts’ tissue) and host mortality leading to staggering declines in geographic coverage (Bruno & Selig 2007) that have raised questions on the viability of this ecosystem as we know it (Bellwood et al. 2004; Parmesan 2006). One avenue towards anticipating the future of the coral reef ecosystem is by developing a broader and deeper understanding of the current genotypic diversity encompassed within and between populations of their keystone species, the scleractinian corals and dinoflagellate symbionts, as they potentially possess functional variation (either singularly or in combination) that may come under selection due to the ongoing and rapid environmental changes they are experiencing. However, such studies, especially for members of the genus Symbiodinium, are sparse. In this issue, Baums et al. (2014) provide a significant contribution by documenting the range‐wide population genetics of Symbiodinium ‘fitti’ (Fig. 1 ) in the context of complementary data from its host, the endangered Caribbean elkhorn coral Acropora palmata (Fig. 1 ). Notable results of this study include a single S. ‘fitti’ genotype typically dominates an individual A. palmata colony both spatially and temporally, gene flow among coral host populations is a magnitude higher to that of its symbiont populations, and the partners possess disparate patterns of genetic differentiation across the Greater Caribbean. The implications of such findings are discussed herein.     

Influence of northern limit range on genetic diversity and structure in a widespread North American tree,sugar maple (Acer saccharum Marshall)

Due to climate change, the ranges of many North American tree species are expected to shift northward. Sugar maple (Acer saccharum Marshall) reaches its northern continuous distributional limit in northeastern North America at the transition between boreal mixed‐wood and temperate deciduous forests. We hypothesized that marginal fragmented northern populations from the boreal mixed wood would have a distinct pattern of genetic structure and diversity. We analyzed variation at 18 microsatellite loci from 23 populations distributed along three latitudinal transects (west, central, and east) that encompass the continuous–discontinuous species range. Each transect was divided into two zones, continuous (temperate deciduous) and discontinuous (boreal mixed wood), based on sugar maple stand abundance. Respective positive and negative relationships were found between the distance of each population to the northern limit (D_north), and allelic richness (A_R) and population differentiation (F_ST). These relations were tested for each transect separately; the pattern (discontinuous–continuous) remained significant only for the western transect. structure analysis revealed the presence of four clusters. The most northern populations of each transect were assigned to a distinct group. Asymmetrical gene flow occurred from the southern into the four northernmost populations. Southern populations in Québec may have originated from two different postglacial migration routes. No evidence was found to validate the hypothesis that northern populations were remnants of a larger population that had migrated further north of the species range after the retreat of the ice sheet. The northernmost sugar maple populations possibly originated from long‐distance dispersal.     

Morphologic and genetic variability in the Barbus fishes (Teleostei,Cyprinidae) of Central Italy

Italian freshwaters are highly biodiverse, with species present including the native fishes Barbus plebejus and Barbus tyberinus that are threatened by habitat alteration, fish stocking and invasive fishes, especially European barbel Barbus barbus. In central Italy, native fluvio‐lacustrine barbels are mainly allopatric and so provide an excellent natural system to evaluate the permeability of the Apennine Mountains. Here, the morphologic and genetic distinctiveness was determined for 611 Barbus fishes collected along the Padany–Venetian (Adriatic basins; PV) and Tuscany–Latium (Tyrrhenian basins; TL) districts. Analyses of morphological traits and mitochondrial DNA sequence data explored the natural and anthropogenic factors that have shaped their distribution ranges. Over 100 alien B. barbus were recorded in the Tiber basin (TL district) and Metauro basin (PV district). Comparisons of genetic and morphometric data revealed that morphometric data could identify alien B. barbus from native Barbus, but could not differentiate between B. tyberinus and B. plebejus. Genetic analyses revealed ~50 D‐loop mtDNA haplotypes and identified a distinct Barbus lineage present only in the Vomano River at the southern boundary of PV district. Demographic expansion and molecular variance analyses revealed a lack of geographic structuring across the sampling regions. While the contemporary B. plebejus distribution has been driven primarily by anthropogenic fish translocations, the dispersal of B. tyberinus has been via natural dispersion, including their crossing of the Apennine Mountains via temporary river connectivity. The results also revealed that the Barbus fishes of the mid‐Adriatic region of Europe have a complex pattern of local endemism. To conserve these patterns of genetic uniqueness, especially in the mid‐Adriatic basins, Barbus fishes should be managed by treating them as unique evolutionary units and ceasing translocations of all Barbus fishes between river basins.     

The influence of breeding phenology on the genetic structure of four pond‐breeding salamanders

Understanding metapopulation dynamics requires knowledge about local population dynamics and movement in both space and time. Most genetic metapopulation studies use one or two study species across the same landscape to infer population dynamics; however, using multiple co‐occurring species allows for testing of hypotheses related to different life history strategies. We used genetic data to study dispersal, as measured by gene flow, in three ambystomatid salamanders (Ambystoma annulatum , A. maculatum , and A. opacum ) and the Central Newt (Notophthalmus viridescens louisianensis ) on the same landscape in Missouri, USA . While all four salamander species are forest dependent organisms that require fishless ponds to reproduce, they differ in breeding phenology and spatial distribution on the landscape. We use these differences in life history and distribution to address the following questions: (1) Are there species‐level differences in the observed patterns of genetic diversity and genetic structure? and (2) Is dispersal influenced by landscape resistance? We detected two genetic clusters in A. annulatum and A. opacum on our landscape; both species breed in the fall and larvae overwinter in ponds. In contrast, no structure was evident in A. maculatum and N. v. louisianensis , species that breed during the spring. Tests for isolation by distance were significant for the three ambystomatids but not for N. v. louisianensis . Landscape resistance also contributed to genetic differentiation for all four species. Our results suggest species‐level differences in dispersal ability and breeding phenology are driving observed patterns of genetic differentiation. From an evolutionary standpoint, the observed differences in dispersal distances and genetic structure between fall breeding and spring breeding species may be a result of the trade‐off between larval period length and size at metamorphosis which in turn may influence the long‐term viability of the metapopulation. Thus, it is important to consider life history differences among closely related and ecologically similar species when making management decisions.     

Endemic species may have complex histories: within‐refugium phylogeography of an endangered Iberian vole

Glacial refugia protected and promoted biodiversity during the Pleistocene, not only at a broader scale, but also for many endemics that contracted and expanded their ranges within refugial areas. Understanding the evolutionary history of refugial endemics is especially important in the case of endangered species to recognize the origins of their genetic structure and thus produce better informed conservation practices. The Iberian Peninsula is an important European glacial refugium, rich in endemics of conservation concern, including small mammals, such as the Cabrera vole (Microtus cabrerae). This near‐threatened rodent is characterized by an unusual suite of genetic, life history and ecological traits, being restricted to isolated geographic nuclei in fast‐disappearing Mediterranean subhumid herbaceous habitats. To reconstruct the evolutionary history of the Cabrera vole, we studied sequence variation at mitochondrial, autosomal and sex‐linked loci, using invasive and noninvasive samples. Despite low overall mitochondrial and nuclear nucleotide diversities, we observed two main well‐supported mitochondrial lineages, west and east. Phylogeographic modelling in the context of the Cabrera vole's detailed fossil record supports a demographic scenario of isolation of two populations during the Last Glacial Maximum from a single focus in the southern part of the Iberian Peninsula. In addition, our data suggest subsequent divergence within the east, and secondary contact and introgression of the expanding western population, during the late Holocene. This work emphasizes that refugial endemics may have a phylogeographic history as rich as that of more widespread species, and conservation of such endemics includes the preservation of that genetic legacy.     

Genetic delineation between and within the widespread coccolithophore morpho‐species Emiliania huxleyi and Gephyrocapsa oceanica (Haptophyta)

Emiliania huxleyi and Gephyrocapsa oceanica are abundant coccolithophore morpho‐species that play key roles in ocean carbon cycling due to their importance as both primary producers and cal‐cifiers. Global change processes such as ocean acidification impact these key calcifying species. The physiology of E. huxleyi, a developing model species, has been widely studied, but its genetic delineation from G. oceanica remains unclear due to a lack of resolution in classical genetic markers. Using nuclear (18S rDNA and 28S rDNA), mitochondrial (cox1, cox2, cox3, rpl16, and dam), and plastidial (16S rDNA, rbcL, tufA, and petA) DNA markers from 99 E. huxleyi and 44 G. oceanica strains, we conducted a multigene/multistrain survey to compare the suitability of different markers for resolving phylogenetic patterns within and between these two morpho‐species. The nuclear genes tested did not provide sufficient resolution to discriminate between the two morpho‐species that diverged only 291Kya. Typical patterns of incomplete lineage sorting were generated in phylogenetic analyses using plastidial genes. In contrast, full morpho‐species delineation was achieved with mitochondrial markers and common intra‐morpho‐species phylogenetic patterns were observed despite differing rates of DNA substitution. Mitochondrial genes are thus promising barcodes for distinguishing these coccolithophore morpho‐species, in particular in the context of environmental monitoring.     

Population genetic structure,genetic diversity,and natural history of the South American species of Nothofagus subgenus Lophozonia (Nothofagaceae) inferred from nuclear microsatellite data

The effect of glaciation on the levels and patterns of genetic variation has been well studied in the Northern Hemisphere. However, although glaciation has undoubtedly shaped the genetic structure of plants in the Southern Hemisphere, fewer studies have characterized the effect, and almost none of them using microsatellites. Particularly, complex patterns of genetic structure might be expected in areas such as the Andes, where both latitudinal and altitudinal glacial advance and retreat have molded modern plant communities. We therefore studied the population genetics of three closely related, hybridizing species of Nothofagus (N. obliqua, N. alpina, and N. glauca, all of subgenus Lophozonia; Nothofagaceae) from Chile. To estimate population genetic parameters and infer the influence of the last ice age on the spatial and genetic distribution of these species, we examined and analyzed genetic variability at seven polymorphic microsatellite DNA loci in 640 individuals from 40 populations covering most of the ranges of these species in Chile. Populations showed no significant inbreeding and exhibited relatively high levels of genetic diversity (H_E = 0.502–0.662) and slight, but significant, genetic structure (R_ST = 8.7–16.0%). However, in N. obliqua, the small amount of genetic structure was spatially organized into three well‐defined latitudinal groups. Our data may also suggest some introgression of N. alpina genes into N. obliqua in the northern populations. These results allowed us to reconstruct the influence of the last ice age on the genetic structure of these species, suggesting several centers of genetic diversity for N. obliqua and N. alpina, in agreement with the multiple refugia hypothesis.