Rapid Homogenization of Multiple Sources: Genetic Structure of a Recolonizing Population of Fishers

Abstract: Fishers (Martes pennanti) were extirpated from much of southern Ontario, Canada, prior to the 1950s. We hypothesised that the recent recolonization of this area originated from an expansion of the population in Algonquin Provincial Park, which historically served as a refuge for fishers. To test this hypothesis, we created a sampling lattice to encompass Algonquin and the surrounding area, and we collected contemporaneous DNA samples. We sampled fishers from each of 35 sites and genotyped them at 16 microsatellite loci. Using a Bayesian assignment approach, with no a priori geographic information, we inferred 5 discrete genetic populations and used genetic population assignment as a means to cluster sites together. We concluded that the Algonquin Park fisher population has not been a substantial source for recolonization and expansion, which has instead occurred from a number of remnant populations within Ontario, Quebec, and most recently from the Adirondacks in New York, USA. The genetic structure among sampling sites across the entire area revealed a pattern of isolation-by-distance (IBD). However, an examination of the distribution of genetic structure (F_ST/1^- F_ST) at different distances showed higher rates of gene flow than predicted under a strict IBD model at small distances (40 km) within clusters and at larger distances up to 100 km among clusters. This pattern of genetic structure suggests increased migration and gene flow among expanding reproductive fronts.     

Is genetic structure of the southern pygmy mouse Baiomys musculus (Cricetidae) related to human-induced spatial landscape heterogeneity in a tropical dry forest?

Tropical dry forests are biologically important biomes sustaining a high rate of endemic species. However, these forests are highly threatened by human activities that negatively impact them on distinct levels, including the genetic diversity. Within the framework of landscape genetics (that seeks to evaluate the relation of characteristics of the environmental matrix with population genetics), we used ISSR markers to evaluate the relationship between the alteration of tropical dry forest by human-induced activities (conserved vs. disturbed) with the genetic structure of four breeding sites of the southern pygmy mouse Baiomys musculus. Averaging among the 105 loci used, the unbiased heterozygosis per population (0.247–0.305) was statistically similar among the four sites, as well as between conserved and disturbed conditions (2-way ANOVA F _(3,16)?=?1.47, P?=?0.1984). The genetic differentiation among the four breeding sites was high (F _ST?=?0.1122; 95?% CI 0.082–0.146) considering the geographical scale evaluated (ca. 4?km). Both, a factorial correspondence analysis and a model-based clustering analysis showed the existence of four genetic groups (one per breeding site). However, no association of genetic structure with disturbance conditions was revealed by these analyses. The absence of differences in mean genetic diversity and the lack of association of genetic structure with habitat transformation suggest that B. musculus has enough behavioral plasticity and enough genetic diversity to respond to environmental heterogeneity caused by human activities. We discussed conditions in which habitat transformation could indeed favor B. musculus.     

Small-scale genetic structure in an endangered wetland specialist: possible effects of landscape change and population recovery

The effects of anthropogenic landscape change on genetic population structure are well studied, but the temporal and spatial scales at which genetic structure can develop, especially in taxa with high dispersal capabilities like birds, are less well understood. We investigated population structure in the Hawaiian gallinule (Gallinula galeata sandvicensis), an endangered wetland specialist bird on the island of O`ahu (Hawai`i, USA). Hawaiian gallinules have experienced a gradual population recovery from near extinction in the 1950s, and have recolonized wetlands on O`ahu in the context of a rapidly urbanizing landscape. We genotyped 152 Hawaiian gallinules at 12 microsatellite loci and sequenced a 520 base-pair fragment of the ND2 region of mitochondrial DNA (mtDNA) from individuals captured at 13 wetland locations on O`ahu in 2014–2016. We observed moderate to high genetic structuring (overall microsatellite F_ST?=?0.098, mtDNA F_ST?=?0.248) among populations of Hawaiian gallinules occupying wetlands at very close geographic proximity (e.g., 1.5–55 km). Asymmetry in gene flow estimates suggests that Hawaiian gallinules may have persisted in 2–3 strongholds which served as source populations that recolonized more recently restored habitats currently supporting large numbers of birds. Our results highlight that genetic structure can develop in taxa that are expanding their range after severe population decline, and that biologically significant structuring can occur over small geographic distances, even in avian taxa.     

Spatial Scales of Pollen and Seed-Mediated Gene Flow in Tropical Rain Forest Trees

Gene flow via seed and pollen is a primary determinant of genetic and species diversity in plant communities at different spatial scales. This paper reviews studies of gene flow and population genetic structure in tropical rain forest trees and places them in ecological and biogeographic context. Although much pollination is among nearest neighbors, an increasing number of genetic studies report pollination ranging from 0.5–14 km for canopy tree species, resulting in extensive breeding areas in disturbed and undisturbed rain forest. Direct genetic measures of seed dispersal are still rare; however, studies of fine scale spatial genetic structure (SGS) indicate that the bulk of effective seed dispersal occurs at local scales, and we found no difference in SGS (Sp statistic) between temperate (N?=?24 species) and tropical forest trees (N?=?15). Our analysis did find significantly higher genetic differentiation in tropical trees (F _ST?=?0.177; N?=?42) than in temperate forest trees (F _ST?=?0.116; N?=?82). This may be due to the fact that tropical trees experience low but significant rates of self-fertilization and bi-parental inbreeding, whereas half of the temperate tree species in our survey are wind pollinated and are more strictly allogamous. Genetic drift may also be more pronounced in tropical trees due to the low population densities of most species.     

Conservation genetics of the yellow-blotched sawback Graptemys flavimaculata (Testudines: Emydidae)

Turtles within the genus Graptemys are increasingly becoming a group of conservation priority due to (1) the number of species endemic to single Gulf of Mexico rivers and (2) human alterations of riverine habitat. The yellow-blotched sawback (Graptemys flavimaculata) is a federally threatened species that is endemic to the Pascagoula River system of southeastern Mississippi, USA. Currently, there is nothing known about genetic variation and population structure in G. flavimaculata. We used microsatellite data to analyze population genetic structure, assess historical demography, and determine effective population size at six sites throughout the Pascagoula River system. Considerable genetic diversity was found within each site (mean allelic richness: 6.65–8.08) and two analyses found no evidence of genetic bottlenecks. All of the pairwise F _ST values, while low (average = 0.026), were significant, with most sites possessing one or more private alleles. Pairwise F _ST values with the Escatawpa River site were larger (0.030–0.047), which likely reflect its historical isolation. Genetic distance was correlated to geographic distance between sites, with the exception of the Escatawpa River site; a similar pattern was also found with estimates of recent rates of migration among sites. While an analysis of molecular variance indicated that most variation was partitioned within rather than among sites, STRUCTURE analysis strongly supported the recognition of two distinct groups (mainstem Pascagoula vs. Escatawpa), with the possibility of additional substructure within the mainstem Pascagoula.     

Unspecific histological and hematological alterations in anadromous and resident <Emphasis Type="Italic">Salvelinus malma</Emphasis> induced by volcanogenic pollution

Knowledge of larval dispersal and connectivity in coral reef species is crucial for understanding population dynamics, resilience, and evolution of species. Here, we use ten microsatellites and one mitochondrial marker (cytochrome b) to investigate the genetic population structure, genetic diversity, and historical demography of the powder-blue tang Acanthurus leucosternon across more than 1000 km of the scarcely studied Eastern African region. The global AMOVA results based on microsatellites reveal a low but significant F _ST value (F _ST = 0.00252 P < 0.001; D _EST = 0.025 P = 0.0018) for the 336 specimens sampled at ten sample sites, while no significant differentiation could be found in the mitochondrial cytochrome b dataset. On the other hand, pairwise F _ST, PCOA, and hierarchical analysis failed to identify any genetic breaks among the Eastern African populations, supporting the hypothesis of genetic homogeneity. The observed genetic homogeneity among Eastern African sample sites can be explained by the lengthy post-larval stage of A. leucosternon, which can potentiate long-distance dispersal. Tests of neutrality and mismatch distribution signal a population expansion during the mid-Pleistocene period.     

Inference of Potential Genetic Risks Associated with Large-Scale Releases of Red Sea Bream in Kanagawa Prefecture, Japan Based on Nuclear and Mitochondrial DNA Analysis

Since 1978, millions of hatchery-reared red sea bream (Pagrus major) juveniles have been released in Sagami Bay and Tokyo Bay in Kanagawa Prefecture, Japan. The stock enhancement program has contributed to total catch; however, no information regarding the genetic interactions with wild counterparts is available. Here, we combined 15 microsatellite loci and mitochondrial D-loop sequencing to characterize the genetic resources of red sea bream in Sagami Bay and Tokyo Bay and to elucidate the potential harmful genetic effects associated with fish releases. Both types of markers evidenced higher levels of genetic diversity in wild samples (SB and TB) compared with offspring before stocking (H07 and H08) as well as a hatchery-released sample recaptured in Sagami Bay (HR). Microsatellite F _ST estimates and Bayesian clustering analysis found significant genetic differences among samples (F _ST?=?0.013–0.054), except for the two wild samples (F _ST?=?0.002) and HR vs. H07 (F _ST?=?0.007). On the other hand, mitochondrial-based Ф _ST suggested haplotypic similarity between SB, H07, and HR. The low effective number of females contributing to the offspring over multiple generations may be responsible for the lack of haplotypic differentiation. Moreover, the putative hatchery origin to three fish (8 %) without deformity in the inter-nostril epidermis was inferred for the first time. Our results showed the usefulness of combining nuclear and mitochondrial markers to elucidate genetic interactions between hatchery-released and wild red sea bream and warned about potential harmful genetic effects should interbreeding takes place.     

High genetic diversity in an endangered medicinal plant, <Emphasis Type="Italic">Saussurea involucrata</Emphasis> (<Emphasis Type="Italic">Saussurea</Emphasis>, Asteraceae), in western Tianshan Mountains,China

Saussurea involucrata (Asteraceae) is a medicinal and second-degree national priority endangered plant that is mainly distributed in the high latitude region of the western Tianshan Mountains. The population is fragmented and isolated, and extensive human impact merits a suitable and specific conservation strategy, which can be compiled based on the genetic diversity, population structure, and demographic history. Phylogeographic studies were conducted on a total of five natural populations and 150 individuals were sampled. Data from three cpDNA intergenic spacer regions (trnL-F, matK, and ndhF-rpl32) and nrDNA ITS sequences showed that twelve haplotypes in cpDNA and five haplotypes in nrDNA indicated high genetic diversity among populations sampled (H _T?=?0.820 and 0.756) and within populations sampled (H _S?=?0.792 and 0.721). Additionally, the high genetic diversity did not mirror genetic structure in either cpDNA (F _ST?=?0.03153, G _ST?>?N _ST, p?<?0.05) or nrDNA (F _ST?=?0.03666, meaningless G _ST and N _ST). Two groups (north and south) were determined for a SAMOVA analysis. Based on this analysis, the demographic history was conducted with a Bayesian Skyline Plot and Isolation with Migration analysis, which showed sustainable and stable extension without a marked bottleneck. Divergence time was indicated at c. 6.25 Mya (90%HPD: 15.30–0.22 Mya) in the Miocene, which is consistent with the formation of the Kelasu section of Tianshan. The southern populations in the Bayanbulak and Gonglu regions require additional attention and transplanting would be an effective way to restore rare cpDNA haplotypes, increase effective population size, and migration rate. Our results suggested that in situ conservation of S. involucrata in western Tianshan should be the main strategy for protection and recovery of the species.     

Genetic differentiation in Elaeocarpus photiniifolia (Elaeocarpaceae) associated with geographic distribution and habitat variation in the Bonin (Ogasawara) Islands

Gene flow between populations in different environmental conditions can be limited due to divergent natural selection, thus promoting genetic differentiation. Elaeocarpus photiniifolia, an endemic tree species in the Bonin Islands, is distributed in two types of habitats, dry scrubs and mesic forests. We aim to elucidate the genetic differentiation in E. photiniifolia within and between islands and between the habitat types. We investigated genotypes of 639 individuals from 19 populations of E. photiniifolia and its closely-related E. sylvestris at 24 microsatellite loci derived from expressed sequence tags. The data revealed genetic differentiation (1) between E. photiniifolia and E. sylvestris (0.307 ≤ F _ST ≤ 0.470), (2) between the E. photiniifolia populations of the Chichijima and Hahajima Island Groups in the Bonin Islands (0.033 ≤ F _ST ≤ 0.121) and (3) between E. photiniifolia populations associated with dry scrubs and mesic forests in the Chichijima Island Group (0.005 ≤ F _ST ≤ 0.071). Principal coordinate analysis and Bayesian clustering analysis also showed that genetically distinct groups were associated with the habitat types, and isolation by distance was not responsible for the genetic differentiation. These findings suggest that E. photiniifolia is divided into genetically differentiated groups associated with different environmental conditions in the Bonin Islands.     

Linking seed dispersal and genetic structure of trees: a biogeographical approach

Aim Natural and human‐induced differences in frugivore assemblages can influence the seed dispersal distances of trees. An important issue in seed dispersal systems is to understand whether differences in seed dispersal distances also affect the genetic structure of mature trees. One possible approach to test for a relationship between seed dispersal and the genetic structure of mature trees is to compare the genetic structure of two closely related tree species between two biogeographical regions that differ in frugivore assemblages and seed dispersal distances. Previous studies on two Commiphora species revealed that Commiphora guillauminii in Madagascar has a much lower seed dispersal distance than Commiphora harveyi in South Africa. We tested whether the lower seed dispersal distance might have caused decreased gene flow, resulting in a stronger genetic structure in Madagascar than in South Africa. Location Madagascar and South Africa. Methods Using amplified fragment length polymorphism markers we investigated the genetic structure of 134 trees in Madagascar and 158 trees in South Africa at a local and a regional spatial scale. Results In concordance with our hypothesis, kinship analysis suggests that gene flow was restricted mostly to 3 km in Madagascar and to 30 km in South Africa. At the local spatial scale, the genetic differentiation among groups of trees within sample sites was marginally significantly higher in Madagascar (F_ST = 0.069) than in South Africa (F_ST = 0.021). However, at a regional spatial scale genetic differentiation was lower in Madagascar (F_ST = 0.053) than in South Africa (F_ST = 0.163). Main conclusions Our results show that lower seed dispersal distances of trees were linked to higher genetic differentiation of trees only at a local spatial scale. This suggests that seed dispersal affects the genetic population structure of trees at a local, but not at a regional, spatial scale.     

Testing the depth-differentiation hypothesis in a deepwater octocoral

The depth-differentiation hypothesis proposes that the bathyal region is a source of genetic diversity and an area where there is a high rate of species formation. Genetic differentiation should thus occur over relatively small vertical distances, particularly along the upper continental slope (200–1000 m) where oceanography varies greatly over small differences in depth. To test whether genetic differentiation within deepwater octocorals is greater over vertical rather than geographical distances, Callogorgia delta was targeted. This species commonly occurs throughout the northern Gulf of Mexico at depths ranging from 400 to 900 m. We found significant genetic differentiation (F_ST = 0.042) across seven sites spanning 400 km of distance and 400 m of depth. A pattern of isolation by depth emerged, but geographical distance between sites may further limit gene flow. Water mass boundaries may serve to isolate populations across depth; however, adaptive divergence with depth is also a possible scenario. Microsatellite markers also revealed significant genetic differentiation (F_ST = 0.434) between C. delta and a closely related species, Callogorgia americana, demonstrating the utility of microsatellites in species delimitation of octocorals. Results provided support for the depth-differentiation hypothesis, strengthening the notion that factors covarying with depth serve as isolation mechanisms in deep-sea populations.     

Analysis of Genetic Variation in Natural Populations of Medicago truncatula of Southern Tunisian Ecological Areas, Using Morphological Traits and SSR Markers

We used 19 quantitative traits and 14 microsatellite markers (SSRs) to analyze the genetic variation in four natural populations of the model legume Medicago truncatula sampled in southern Tunisia. The greatest genetic variation of quantitative traits and molecular markers occurred within populations (>71%). In contrast to quantitative population differentiation (Q _ST ?=?0.09), a high level of molecular differentiation (F _ST ?=?0.23) was found among populations. The majority of quantitative traits exhibited Q _ST values significantly less than F _ST values, suggesting that selection may be acting to suppress differentiation for these traits. There was no significant correlation between genetic variation of quantitative traits and molecular markers within populations. On the other hand, significant correlations were found between measured quantitative characters and the site-of-origin environmental factors. The eco-geographical factors with the greatest influence on the variation of measured traits among populations were altitude, followed by soil texture, assimilated phosphorus (P₂O₅) and organic matter. Nevertheless, there were no consistent patterns of associations between gene diversity (He) and eco-geographical factors.     

Two decades of genetic consistency in a reproductive population in the face of exploitation: patterns of adult and larval walleye (<Emphasis Type="Italic">Sander vitreus</Emphasis>) from Lake Erie’s Maumee River

Analyses of genetic variability and allelic composition in a species exhibiting reproductive fidelity to natal sites may provide important ecological indication of temporal population dynamics, facilitating understanding responses to past disturbances and future climate change. The walleye is an ecologically and economically valuable species, whose largest fishery centers in Lake Erie of the Laurentian Great Lakes; it exhibits reproductive site fidelity, despite otherwise wide-ranging dispersal. We tested whether genetic composition and diversity have remained temporally stable in Lake Erie’s Maumee River, which is the largest and most highly fished spawning run. This population has experienced over a century of exploitation, habitat alterations, and pollution, which may have affected genetic structure and might influence future sustainability. Fourteen nuclear DNA microsatellite loci were analyzed from 744 spawning run walleye to test genetic patterns across: (1) years (N = 12, spanning 1995–2013), (2) birth year cohorts, (3) the sexes, (4) those reproducing earlier (ages 2–6) versus later (7 or older) in life, and (5) the adults versus larvae. Results indicated stability in genetic diversity levels (mean H _O = 0.76 ± 0.03) and allelic composition across years (F _ST = 0.000–0.006, NS), cohorts (F _ST = 0.000–0.013, NS), sexes (F _ST = 0.000, NS), earlier versus later reproduction (F _ST = 0.000, NS), and between the larvae and adults (F _ST = 0.000–0.004, NS). Number of breeders and effective population size were substantial and consistent. This reproductive population thus has maintained genetic stability and high diversity, despite intensive anthropogenic pressures.     

Isolation by landscape in populations of a prized edible mushroom Tricholoma matsutake

Tricholoma matsutake, a wild edible ectomycorrhizal mushroom, is revered for its distinguished flavor and iconic significance. Here, we test for landscape effects on T. matsutake gene flow and population structure in the Eastern Himalayas. Using single-nucleotide polymorphic (SNP) DNA markers, isolation by distance patterns were tested on eight populations within and between watersheds. We find that high, treeless ridgelines are effective barriers to gene flow, even at distances less than 65 km, whereas populations located within watersheds are structured at greater distances. Mantel tests demonstrated a significant positive correlation between F _st and a “landscape distance” measured as the shortest distance between population pairs below treeline r = 0.574, P = 0.002, whereas strict euclidian distances do not correlate. AMOVA analysis revealed significant partitioning with 91% of the genetic variance found within populations and 7% found between watersheds, indicative of sexually recombining populations with limited gene flow between watersheds. We show that landscape is an important determinant of air-dispersed ectomycorrhizal species population structure in heterogeneous landscapes.     

Analysis of Geographic and Pairwise Genetic Distances Among Sheep Populations

This study assessed the usefulness of geographic and pairwise genetic distances in the characterization of five sheep populations using 15 microsatellite markers. The average F statistics across loci were F _IT = 0.523 ± 0.140, F _ST = 0.363 ± 0.131, and F _IS = 0.263 ± 0.092. The average heterozygosity was 0.716 ± 0.069, polymorphism information content was 0.691 ± 0.070, and effective number of alleles was 3.736 ± 0.998. Sheep populations clustered into group 1 (Hu and Tong breeds) and group 2 (small-tailed Han, Wadi, and Tan breeds). Reynolds’ distance varied from 0.0062 to 0.0499, and the range of gene flow (N _m) was 4.8834–40.0726 among the sheep populations. The results showed that the genetic structure of the five populations was not consistent with their genetic distances, and the population genetic divergence was not linearly related to geographic distance as indicated by a Mantel test (P = 0.7936).     

Population stock structure of leatherback turtles (Dermochelys coriacea) in the Atlantic revealed using mtDNA and microsatellite markers

This study presents a comprehensive genetic analysis of stock structure for leatherback turtles (Dermochelys coriacea), combining 17 microsatellite loci and 763 bp of the mtDNA control region. Recently discovered eastern Atlantic nesting populations of this critically endangered species were absent in a previous survey that found little ocean-wide mtDNA variation. We added rookeries in West Africa and Brazil and generated longer sequences for previously analyzed samples. A total of 1,417 individuals were sampled from nine nesting sites in the Atlantic and SW Indian Ocean. We detected additional mtDNA variation with the longer sequences, identifying ten polymorphic sites that resolved a total of ten haplotypes, including three new variants of haplotypes previously described by shorter sequences. Population differentiation was substantial between all but two adjacent rookery pairs, and F _ST values ranged from 0.034 to 0.676 and 0.004 to 0.205 for mtDNA and microsatellite data respectively, suggesting that male-mediated gene flow is not as widespread as previously assumed. We detected weak (F _ST = 0.008 and 0.006) but significant differentiation with microsatellites between the two population pairs that were indistinguishable with mtDNA data. POWSIM analysis showed that our mtDNA marker had very low statistical power to detect weak structure (F _ST < 0.005), while our microsatellite marker array had high power. We conclude that the weak differentiation detected with microsatellites reflects a fine scale level of demographic independence that warrants recognition, and that all nine of the nesting colonies should be considered as demographically independent populations for conservation. Our findings illustrate the importance of evaluating the power of specific genetic markers to detect structure in order to correctly identify the appropriate population units to conserve.     

High gene flow between populations of Macrotermes michaelseni (Isoptera, Termitidae)

Termite alates are thought to be poor active flyers, and this should lead to considerable genetic differentiation on small spatial scales. However, using four microsatellite loci for the termite Macrotermes michaelseni we found low values of genetic differentiation (F_ST) across a spatial scale of even more than 50 km. Genetic differentiation between populations increased with spatial distance up to 50 km. Furthermore, up to this distance, the scatter around the linear regression of genetic differentiation versus spatial distance increased with spatial distance. This suggests that across such spatial distances gene flow and genetic drift are of about equal importance, and near equilibrium. Using a regional F_ST as well as the distance between populations with non-significant F_ST-values (up to 25 km), gene flow is sufficiently high so that populations may be regarded as panmictic on spatial scales of 25 to 50 km. The apparent contradiction between dispersal distances observed in the field and estimates of gene flow from genetic markers may be due to the masses of swarming alates. Assuming a leptokurtic distribution of dispersal distances, atleast some alates are expected to travel considerable distances, most likely by passive drift. Received 25 January 2005; revised 11 April 2005; accepted 26 April 2005.     

Extreme isolation by distance in a montane frog <Emphasis Type="Italic">Rana cascadae</Emphasis>

Given the recent interest in declining amphibian populations, it is surprising that there are so few data on genetic drift and gene flow in anuran species. We used seven microsatellite loci to investigate genetic structure and diversity at both large and small geographic scales, and to estimate gene flow in the Cascades frog, Rana cascadae. We sampled 18 sites in a hierarchical design (inter-population distances ranging from 1–670 km) to test for isolation by distance and to determine the geographic scale over which substantial gene flow occurs. Eleven of these sites were sampled as three fine-scale clusters of three, three, and five sites separated by pairwise distances of 1–23 km to estimate number of migrants exchanged per generation via F _ST and by a coalescent approach. We found R. cascadae exhibits a strong pattern of isolation by distance over the entire species range, and that there is a sharp drop in migrants exchanged between sites separated by greater than 10 km. These data, in conjunction with results of other recent studies, suggest that montane habitats promote unusually strong genetic isolation among frog populations. We discuss our results in light of future management and conservation of R. cascadae.     

Pronounced differences in genetic structure despite overall ecological similarity for two Ambystoma salamanders in the same landscape

Studies linking genetic structure in amphibian species with ecological characteristics have focused on large differences in dispersal capabilities. Here, we test whether two species with similar dispersal potential but subtle differences in other ecological characteristics also exhibit strong differences in genetic structure in the same landscape. We examined eight microsatellites in marbled salamanders (Ambystoma opacum) from 29 seasonal ponds and spotted salamanders (Ambystoma maculatum) from 19 seasonal ponds in a single geographic region in west-central Massachusetts. Despite overall similarity in ecological characteristics of spotted and marbled salamanders, we observed clear differences in the genetic structure of these two species. For marbled salamanders, we observed strong overall genetic differentiation (F _ST = 0.091, F′ _ST = 0.375), three population-level clusters of populations (K = 3), a strong pattern of isolation by distance (r = 0.58), and marked variation in family-level structure (from 1 to 23 full-sibling families per site). For spotted salamanders, overall genetic differentiation was weaker (F _ST = 0.025, F′ _ST = 0.102), there was no evidence of population-level clustering (K = 1), the pattern of isolation by distance (r = 0.17) was much weaker compared to marbled salamanders, and there was less variation in family-level structure (from 10 to 36 full-sibling families per site). We suspect that a combination of breeding site fidelity, effective population size, and generation interval is responsible for these marked differences. Our results suggest that marbled salamanders, compared to spotted salamanders, are more sensitive to fragmentation from various land-use activities and would be less likely to recolonize extirpated sites on an ecologically and conservation-relevant time frame.     

Comparative phylogenetic relationships and genetic structure of the caterpillar fungus Ophiocordyceps sinensis and its host insects inferred from multiple gene sequences

Ophiocordyceps sinensis (Ascomycota: Ophiocordycipitaceae) is a native fungal parasite of Hepialidae caterpillars and one of the most economically important medicinal caterpillar fungi in China. However, little is known about the phylogenetic and evolutionary relationships between O. sinensis and its host insects. In this study, nuclear ITS and β-tubulin sequences from O. sinensis and mitochondrial COI, COII, and Cytb sequences from its hosts were analyzed across 33 populations sampled from five regions in China. Phylogenetically, both O. sinensis and its hosts were divided into three geographically correlated clades, and their phylogenies were congruent. Analysis of molecular variance and calculated coefficients of genetic differentiation revealed significant genetic divergence among the clades within both O. sinensis (F_ST= 0.878, N_ST=0.842) and its hosts (F_ST=0.861, N_ST=0.816). Estimated gene flow was very low for O. sinensis (Nm=0.04) and the host insects (Nm=0.04) among these three clades. Mantel tests demonstrated a significant correlation (P<0.01) between the genetic distances for O. sinensis and its hosts, as well as a significant association (P<0.05) between geographic and genetic distances in both. The similar phylogenetic relationships, geographic distributions, and genetic structure and differentiation between O. sinensis and its hosts imply that they have coevolved.