Genetic structure among continental and island populations of gyrfalcons

Little is known about the possible influence that past glacial events have had on the phylogeography and population structure of avian predators in the Arctic and sub-Arctic. In this study, we use microsatellite and mitochondrial control region DNA variation to investigate the population genetic structure of gyrfalcons (Falco rusticolus) throughout a large portion of their circumpolar distribution. In most locations sampled, the mtDNA data revealed little geographic structure; however, five out of eight mtDNA haplotypes were unique to a particular geographic area (Greenland, Iceland, or Alaska) and the Iceland population differed from others based on haplotype frequency differences (F(ST)). With the microsatellite results, significant population structure (F(ST), principal components analysis, and cluster analysis) was observed identifying Greenland and Iceland as separate populations, while Norway, Alaska and Canada were identified as a single population consistent with contemporary gene flow across Russia. Within Greenland, differing levels of gene flow between western and eastern sampling locations was indicated with apparent asymmetric dispersal in western Greenland from north to south. This dispersal bias is in agreement with the distribution of plumage colour variants with white gyrfalcons in much higher proportion in northern Greenland. Lastly, because the mtDNA control region sequence differed by only one to four nucleotides from a common haplotype among all gyrfalcons, we infer that the observed microsatellite population genetic structure has developed since the last glacial maximum. This conclusion is further supported by our finding that a closely related species, the saker falcon (Falco cherrug), has greater genetic heterogeneity, including mtDNA haplotypes differing by 1-16 nucleotide substitutions from a common gyrfalcon haplotype. This is consistent with gyrfalcons having expanded rapidly from a single glacial-age refugium to their current circumpolar distribution. Additional sampling of gyrfalcons from Fennoscandia and Russia throughout Siberia is necessary to test putative gene flow between Norway and Alaska and Canada as suggested by this study.     

Phylogeography of Varicorhinus barbatulus (Cyprinidae) in Taiwan based on nucleotide variation of mtDNA and allozymes

The phylogeographical patterns and population genetic structures of Varicorhinus barbatulus in Taiwan were investigated based on genetic diversity of 34 allozyme loci and nucleotide sequences of 3' end of the cytochrome b gene, tRNA genes, D-loop control region, and the 5' end of the 12S rRNA of mtDNA. Allozyme and mtDNA analyses revealing evident geographical structuring suggest limited gene flow between populations (F(ST)=0.511 and 0.791, respectively). Low genetic variability within populations (P=5.56%; He=0.018) based on allozymes and significantly negative Tajima's D statistics based on mtDNA suggest that most populations in Taiwan may have originated from a small number of founders followed by demographic expansion. The gene genealogy of mtDNA identified six lineages corresponding to major drainages that were separated by the geological barriers due to vicariant events. A minimum spanning network based on nucleotide substitutions reflects divergence from populations of the Miao-li Plateau to northern and southern regions of the island. In contrast to a previous hypothesis that suggests an early invasion to eastern part of Taiwan prior to the lifting of central mountain range some one million years ago, the mtDNA genealogy and molecular dating reveal very recent colonization of the eastern population. Nested clade analyses revealing significant associations between genetic structure and geographical division identify past fragmentation and range expansion as major phylogeographical events that shaped the geographical distribution of this species in Taiwan.     

Intraspecific variation in Podostemum ceratophyllum (Podostemaceae): evidence of refugia and colonization since the last glacial maximum

? Premise of the Study: Intraspecific variation among 20 populations of Podostemum ceratophyllum Michx. was investigated to test the hypothesis of range expansion from southern refugia since the last glacial maximum. ? Methods: Six noncoding regions of chloroplast DNA were sequenced in 60 individuals. Populations were divided into two groups, north and south of the glacial boundary, in addition to isolated populations in Arkansas and Honduras. Variation in populations north of the boundary was compared with variation in populations to the south and in the isolated populations. ? Key Results: Nucleotide diversity was an order of magnitude lower in populations north of the glacial boundary than in those to the south. The Arkansas and Honduras populations showed no variation. The predominant haplotype in northern populations was also found in a Virginia population. ? Conclusions: Reduced variation north of the glacial boundary suggests a founder event associated with range expansion since the last glacial maximum. Colonization probably occurred from populations in refugia located several hundred kilometers south of the glacial boundary. The results provide insight into the effects of past and current climate change on patterns of geographic distribution and genetic variation in aquatic plants.     

The genetic pattern of population threat and loss: a case study of butterflies

Recent decreases in biodiversity in Europe are commonly thought to be due to land use and climate change. However, the genetic diversity of populations is also seen as one essential factor for their fitness. Genetic diversity in species across the continent of Europe has been recognized as being in part a consequence of ice age isolation in southern refugia and postglacial colonization northwards, and these phylogeographical patterns may themselves affect the adaptability of populations. Recent work on butterfly species with different refugia, colonization paths and genetic structures allows this idea to be examined. The 'chalk-hill blue' pattern is one of decreasing genetic diversity from south to north, whereas the 'woodland ringlet' pattern shows greater genetic diversity in eastern than in western lineages. Comparison of population demographic trends in species with these biogeographical patterns reveals higher rates of decrease with lower genetic diversity. This indicates reduced adaptability due to genetic impoverishment as a result of glacial and postglacial range changes. Analysis of phylogeographical pattern may be a useful guide to interpreting demographic trends and in conservation planning.     

Mitochondrial DNA sequence analysis reveals multiple Pleistocene glacial refugia for the Yellow‐spotted mountain newt,Neurergus derjugini (Caudata: Salamandridae) in the mid‐Zagros range in Iran and Iraq

Phylogeography is often used to investigate the effects of glacial cycles on current genetic structure of various plant and animal species. This approach can also identify the number and location of glacial refugia as well as the recolonization routes from those refugia to the current locations. To identify the location of glacial refugia of the Yellow‐spotted mountain newt, Neurergus derjugini, we employed phylogeography patterns and genetic variability of this species by analyzing partial ND4 sequences (867 bp) of 67 specimens from 15 sampling localities from the whole species range in Iran and Iraq. Phylogenetic trees concordant with haplotype networks showed a clear genetic structure among populations as three groups corresponding to the populations in the north, center, and south. Evolutionary ages of clades north and south ranging from 0.15 to 0.17 Myr, while the oldest clade is the central clade, corresponding to 0.32 Myr. Bayesian skyline plots of population size change through time show a relatively slight increase until about 25 kyr (around the last glacial maximum) and a decline of population size about 2.5 kyr. The presence of geographically structured clades in north, center, and south sections of the species range signifies the disjunct populations that have emerged in three different refugium. This study illustrates the importance of the effect of previous glacial cycles in shaping the genetic structure of mountain species in the Zagros range. These areas are important in terms of long‐term species persistence and therefore valuable areas for conservation of biodiversity.     

Genetic consequences of glacial survival and postglacial colonization in Norway spruce: combined analysis of mitochondrial DNA and fossil pollen

Norway spruce (Picea abies [L.] Karst.) is a broadly distributed European conifer tree whose history has been intensively studied by means of fossil records to infer the location of full‐glacial refugia and the main routes of postglacial colonization. Here we use recently compiled fossil pollen data as a template to examine how past demographic events have influenced the species’ modern genetic diversity. Variation was assessed in the mitochondrial nad1 gene containing two minisatellite regions. Among the 369 populations (4876 trees) assayed, 28 mitochondrial variants were identified. The patterns of population subdivision superimposed on interpolated fossil pollen distributions indicate that survival in separate refugia and postglacial colonization has led to significant structuring of genetic variation in the southern range of the species. The populations in the northern range, on the other hand, showed a shallow genetic structure consistent with the fossil pollen data, suggesting that the vast northern range was colonized from a single refugium. Although the genetic diversity decreased away from the putative refugia, there were large differences between different colonization routes. In the Alps, the diversity decreased over short distances, probably as a result of population bottlenecks caused by the presence of competing tree species. In northern Europe, the diversity was maintained across large areas, corroborating fossil pollen data in suggesting that colonization took place at high population densities. The genetic diversity increased north of the Carpathians, probably as a result of admixture of expanding populations from two separate refugia.     

Phylogeography of the montane caddisfly Drusus discolor: evidence for multiple refugia and periglacial survival

We studied the genetic population structure and phylogeography of the montane caddisfly Drusus discolor across its entire range in central and southern Europe. The species is restricted to mountain regions and exhibits an insular distribution across the major mountain ranges. Mitochondrial sequence data (COI) of 254 individuals from the entire species range is analysed to reveal population genetic structure. The data show little molecular variation within populations and regions, but distinct genetic differentiation between mountain ranges. Most populations are significantly differentiated based on F(ST) and exact tests of population differentiation and most haplotypes are unique to a single mountain range. Phylogenetic analyses reveal deep divergence between geographically isolated lineages. Combined, these results suggest that past fragmentation is the prominent process structuring the populations across Europe. We use tests of selective neutrality and mismatch distributions, to study the demographic population history of regions with haplotype overlap. The high level of genetic differentiation between mountain ranges and estimates of demographic history provide evidence for the existence of multiple glacial refugia, including several in central Europe. The study shows that these aquatic organisms reacted differently to Pleistocene cooling than many terrestrial species. They persisted in numerous refugia over multiple glacial cycles, allowing many local endemic clades to form.     

Life on the edge: the long-term persistence and contrasting spatial genetic structure of distinct brown trout life histories at their ecological limits

The Atlantic Iberian brown trout is at the southwestern limit of its distribution. At this ecological edge, which was once a glacial refugia, anadromy becomes less common as increased water temperatures restricted populations closer to the headwaters. We examined 847 individuals from 20 populations from throughout this region and assessed spatial genetic structure using 11 protein and four microsatellite loci. The higher levels of heterozygosity and allelic diversity north of the southernmost limit of anadromy (SLA), as well as an isolation-by-distance model of population structure, likely influenced by the anadromous forms, suggest that more stable demographic conditions existed over time in this region. Populations south of the SLA were highly differentiated given the restricted size of the area (protein F(ST) = 0.16 in the north and 0.63 in the south of the SLA; microsatellite F(ST) = 0.18 in the north and 0.70 in the south of the SLA). The low levels of heterozygosity and the pattern of southward allele depletion in resident populations is indicative of fragmentation, caused by stressful ecological conditions that reduced the anadromy (restricting gene flow) and the effective population sizes (higher genetic drift), which, in combination, dramatically decreased within-population genetic variation and increased among-population genetic variation. The higher gene diversity north of the SLA does not reflect ancestry but rather the signature of a population size expansion, as evidence suggest the persistence of older populations (with several private alleles) south of the SLA. These data support a scenario that demonstrates how contemporary events (critical ecological conditions) can moderate historical influences, suggesting that careful interpretation of the evolutionary history of glacial refugia is necessary, especially where populations persisted for a long time but not always with optimal ecological conditions. These peripheral populations are of high conservation value and should be preserved to help conserve the future potential of the species.     

Diversification in subtropical mountains: phylogeography, Pleistocene demographic expansion, and evolution of polyphenic mandibles in Taiwanese stag beetle, Lucanus formosanus

Pleistocene glacial oscillations have had profound impacts on the historical population dynamics of extant species. However, the genetic consequences of past climatic changes depend largely on the latitude and topography of the regions in question. This study investigates the effect of Pleistocene glacial periods and the Central Mountain Range on the phylogeography, historical demography, and phenotypic differentiation of a montane forest-dwelling stag beetle, Lucanus formosanus (Coleoptera: Lucanidae), which exhibits extensive mandible variations across mountain ranges in subtropical Taiwan. Analyses of mitochondrial (cox1) and nuclear (wg) loci reveal that L. formosanus originated nearly 1.6 million years ago (Mya) in the early Pleistocene period and consisted of geographically overlapping Alishan and Widespread clades. A drastic population expansion starting approximately 0.2 Mya in the Widespread clade likely resulted from altitudinal range shift of the temperate forests, which was closely tied to the arrival of the Riss glacial period in the late Middle Pleistocene. A ring-like pattern of historical gene flow among neighboring populations in the vicinity of the Central Mountain Range indicates that the mountains constitute a strong vicariant barrier to the east-west gene flow of L. formosanus populations. A geographic cline of decreasing mandible size from central to north and south, and onto southeast of Taiwan is inconsistent with the low overall phylogeographic structures. The degree of mandible variation does not correlate with the expected pattern of neutral evolution, indicating that the evolutionary diversification of this morphological weapon is most likely subject to sexual or natural selection. We hypothesize that the adaptive evolution of mandibles in L. formosanus is shaped largely by the habitat heterogeneity.     

Sex-biased gene flow and colonization in the Formosan lesser horseshoe bat: inference from nuclear and mitochondrial markers

Sex-biased behaviours are expected to play an important role in partitioning genetic variance in animal populations. Comparing genetic structure at markers with different modes of inheritance provides a means of detecting these behaviours and their consequences for population genetic structure. In colonially breeding mammals, the common combination of female philopatry and male vagility can promote contrasting patterns of genetic differentiation between the sexes, both via their effects on recurrent gene flow and on colonization. We examined sex differences in gene flow and structure by comparing maternally inherited mitochondrial DNA (mtDNA) and biparentally inherited autosomal loci in the Formosan lesser horseshoe bat Rhinolophus monoceros . We found that genetic partitioning was higher at mtDNA than autosomal markers in both sexes, indicative of female-biased philopatry and male-biased dispersal. Across Taiwan, isolation-by-distance was detected for all sex/marker combinations but was steeper for mtDNA than for nuclear markers. We suggest that isolation-by-distance shown from mtDNA at large scales is likely to reflect the stepwise founding of new breeding colonies by females during colonization. In contrast, no isolation-by-distance was found at smaller distances of up to 100 km, indicating that gene flow and/or recent shared ancestry homogenises genetic structure among nearby sites. Our results highlight the value of an indirect genetic approach to understanding sex-biased behaviours and their consequences in a little-studied species.     

Foundations for conservation of intraspecific genetic diversity revealed by analysis of phylogeographical structure in the endangered endemic lizard Podarcis lilfordi

Aim  To describe and analyse phylogeographical patterns in the endangered endemic lizard Podarcis lilfordi from across its remaining range and thereby establish baseline information on genetic diversity that will help determine conservation priorities and assist future reintroduction programs.
Location  Balearic Islands, Spain.
Methods  We analysed mitochondrial DNA (2382 bp sequence from eight genes) from 118 individuals and characterized the relationships among haplotypes using parsimony networks, as well as phylogenetic inference. Analyses of historical gene flow and population growth were used to provide further insights into population histories.
Results  Four unconnected parsimony networks were obtained that mirrored the main clades in the phylogenetic tree: (I) all Menorcan populations, (II) Dragonera, Malgrats and Toro islands (Western Mallorca) (III and IV) and the remaining populations from Cabrera and Mallorca. Two major haplotype groups were detected in Menorca (I) and these provided signatures of a demographic expansion and asymmetrical historical gene flow, respectively, concordant with the expected direction of colonization from south to north of the island. Populations from western Mallorca (II) showed evidence of historical allopatric fragmentation events following isolation around the start of the Pleistocene. In networks III and IV, Cabreran populations appear to have become isolated from north and south Mallorca quite recently, with asymmetric gene flow indicating a northwards dispersal direction.
Main conclusions  P. lilfordi is a genetically diverse species that shows substantial mtDNA structuring both between regions and, at a finer scale, between some islet populations within regions. The precarious state of some islet populations shown here to be quite divergent (e.g. Toro island in western Mallorca) means that conservation of this intraspecific biodiversity requires urgent action.     

Hierarchical genetic structure and gene flow in three sympatric species of Amazonian rodents

The population genetic structure of three species of Amazonian rodents ( Oligoryzomys microtis, Oryzomys capito , and Mesomys hispidus ) is examined for mtDNA sequence haplotypes of the cytochrome b gene by hierarchical analysis of variance and gene flow estimates based on fixation indices ( N _ST) and coalescence methods. Species samples are from the same localities along 1000 km of the Rio Juruá in western Amazonian Brazil, but each species differs in important life history traits such as population size and reproductive rate. Average haplotype differentiation, hierarchical haplotype apportionment, and gene flow estimates are contrasted in discussing the current and past population structure. Two species exhibit isolation by distance patterns wherein gene flow is largely limited to geographically adjacent localities. Mesomys exhibits this pattern throughout its range along the river. More than 75% of haplotype variation is apportioned among localities and regions, and estimates of Nm for pair-wise comparisons are nearly always less than 1. Oligoryzomys shows weak isolation by distance, but only over the largest geographical distances. Nm values for this species are nearly always above 1 and most (about 80%) of haplotype variation is contained within local populations. In contrast, Oryzomys exhibits no genetic structure throughout its entire distribution; Nm values average 17 and nearly 90% of the total haplotype variance is contained within local populations. Although gene flow estimates are high, the pattern of Nm as a function of geographical distance suggests that this species experienced a more recent invasion of the region and is still in genetic disequilibrium under its current demographic conditions.     

Population genetics and phylogeography of freshwater mussels in North America, Elliptio dilatata and Actinonaias ligamentina (Bivalvia: Unionidae)

Extrinsic and intrinsic forces combined shape the population structure of every species differently. Freshwater mussels are obligate parasites to a host fish during a juvenile stage (glochidia). Elliptio dilatata (ED) and Actinonaias ligamentina (AL) are co-occurring freshwater mussel taxa with similar North American distribution and share some potential host fish. Using mitochondrial DNA, we determined the genotypes of 190 + individuals from collection sites in at least two tributaries in the Lake Erie and Ohio River watersheds, along with the Ouachita and Strawberry rivers in the southeast. Both species had followed a stepping-stone model of dispersal, with greater pairwise genetic structure among collection sites of ED. Also, phylogeographical analysis for ED found significant geographical structuring of haplotype diversity. Overall, within-population variation increased significantly from north to south, with low genetic diversity in the Strawberry River. We calculated significant among-population structure for both species (ED: Phi(ST) = 0.62, P < 0.001; AL: Phi(ST) = 0.16, P < 0.001). Genetic analysis identified the Ouachita River as an area of significant reproductive isolation for both species. Results for AL indicated dispersal into northern areas from two genetically distinct glacial refugia, where results for ED indicated dispersal followed by low gene flow in northern areas. The conservation strategies for mussels that co-occur in the same 'bed' could be species specific. Species such as ED have management units on the population scale, where AL has a more homogeneous genetic structure across its range.     

Genetic diversity, structure, and demographic change in tanoak, Lithocarpus densiflorus (Fagaceae), the most susceptible species to sudden oak death in California

Knowledge of population genetic structure of tanoak (Lithocarpus densiflorus) is of interest to pathologists seeking natural variation in resistance to sudden oak death disease, to resource managers who need indications of conservation priorities in this species now threatened by the introduced pathogen (Phytophthora ramorum), and to biologists with interests in demographic processes that have shaped plant populations. We investigated population genetic structure using nuclear and chloroplast DNA (cpDNA) and inferred the effects of past population demographic processes and contemporary gene flow. Our cpDNA results revealed a strong pattern of differentiation of four regional groups (coastal California, southern Oregon, Klamath mountains, and Sierra Nevada). The chloroplast haplotype phylogeny suggests relatively deep divergence of Sierra Nevada and Klamath populations from those of coastal California and southern Oregon. A widespread coastal California haplotype may have resulted from multiple refugial sites during the Last Glacial Maximum or from rapid recolonization from few refugia. Analysis of nuclear microsatellites suggests two major groups: (1) central coastal California and (2) Sierra Nevada/Klamath/southern Oregon and an area of admixture in north coastal California. The low level of nuclear differentiation is likely to be due to pollen gene flow among populations during postglacial range expansion.     

利用近似贝氏计算推论台湾海峡沿岸秋茄种群的拓殖路线

由于地理关系,台湾海峡两岸的红树植物组成具有高度的相似性,都以耐寒性较强的秋茄为优势种。中国台湾(以下简称"台湾")与大陆仅一水之隔,因此台湾的秋茄种群来源最有可能来自东南沿海种群,然而台湾南、北红树植物种群的拓殖路线以及与大陆东南沿海种群的遗传关系的研究至今仍未见报道。通过SSR分子标记,利用近似贝氏计算(Approximate Bayesian Computation)推测海峡两岸4个分布区域秋茄的起源及其拓殖路线。结果表明4个区域的种群出现明显分化,大陆东南北部种群与其他种群间分化程度最高。通过推测台湾北部种群起源可追溯到29000—48400a前,早于末次冰期时间,且台湾北部种群遗传结构与大陆东南南部种群最相近,推测它们可能共同起源于南方祖先。大陆东南沿海南北种群的溯祖时间约为15.1万年至25.2a年前,约为更新世中期末,则意味东南沿海南、北种群的遗传分化可能受到更新世后期气候变化与海侵海退的影响而出现隔离,或东南沿海南、北种群可能来自不同的起源。而台湾南部种群与台湾北部种群的相似性,表明台湾南部种群是由北部种群拓殖而来,近似贝氏计算亦支持这个假说。因而,可以推测海峡两岸秋茄的拓殖路线是从大陆东南南方种群随黑潮迁移至台湾北部,再从北部拓殖到台湾南部。利用近似贝氏计算推论台湾海峡两岸红树林种群起源及拓殖路线,为未来我国东南沿海红树林植物的生物地理研究提供参考。     

Diverse processes shape deep phylogeographical divergence in Cobitis sinensis (Teleostei: Cobitidae) in East Asia

Little has been known about the impacts of past vicariance events on the phylogeography and population structure of freshwater fishes in East Asia. The aims of this study are to assess the genetic variability with extensive sampling throughout the range of Chinese spiny loach, Cobitis sinensis, and to infer the genetic structure and evolutionary history of populations. Cobitis sinensis in China may have initiated from two ancestral populations, namely Yangtze and Pearl Rivers, which diverged about 7.24 MYA likely due to drainage systems alteration. In the phylogroup I, a southward dispersal event occurred from East China (Yangtze River) to south ZheMin and Hainan subregions, followed by eastward dispersal from ZheMin to south Taiwan. In the phylogroup II, eastward colonization took place from Pearl River to north Taiwan in the late Pliocene, coupled with loss of genetic diversity in the island populations. This study showed that Cenozoic tectonic movements and climatic and sea‐level fluctuations may have shaped the genetic structure of C. sinensis in concert. Highly diverged mtDNA sequences suggest existence of cryptic species in morphospecies C. sinensis.     

Multilocus phylogeography and population structure of common eiders breeding in North America and Scandinavia

Aim Glacial refugia during the Pleistocene had major impacts on the levels and spatial apportionment of genetic diversity of species in northern latitude ecosystems. We characterized patterns of population subdivision, and tested hypotheses associated with locations of potential Pleistocene refugia and the relative contribution of these refugia to the post‐glacial colonization of North America and Scandinavia by common eiders (Somateria mollissima). Specifically, we evaluated localities hypothesized as ice‐free areas or glacial refugia for other Arctic vertebrates, including Beringia, the High Arctic Canadian Archipelago, Newfoundland Bank, Spitsbergen Bank and north‐west Norway. Location Alaska, Canada, Norway and Sweden. Methods Molecular data from 12 microsatellite loci, the mitochondrial DNA (mtDNA) control region, and two nuclear introns were collected and analysed for 15 populations of common eiders (n = 716) breeding throughout North America and Scandinavia. Population genetic structure, historical population fluctuations and gene flow were inferred using F‐statistics, analyses of molecular variance, and multilocus coalescent analyses. Results Significant inter‐population variation in allelic and haplotypic frequencies were observed (nuclear DNA F_ST = 0.004–0.290; mtDNA Φ_ST = 0.051–0.927). Whereas spatial differentiation in nuclear genes was concordant with subspecific designations, geographic proximity was more predictive of inter‐population variance in mitochondrial DNA haplotype frequency. Inferences of historical population demography were consistent with restriction of common eiders to four geographic areas during the Last Glacial Maximum: Belcher Islands, Newfoundland Bank, northern Alaska and Svalbard. Three of these areas coincide with previously identified glacial refugia: Newfoundland Bank, Beringia and Spitsbergen Bank. Gene‐flow and clustering analyses indicated that the Beringian refugium contributed little to common eider post‐glacial colonization of North America, whereas Canadian, Scandinavian and southern Alaskan post‐glacial colonization is likely to have occurred in a stepwise fashion from the same glacial refugium. Main conclusions Concordance of proposed glacial refugia used by common eiders and other Arctic species indicates that Arctic and subarctic refugia were important reservoirs of genetic diversity during the Pleistocene. Furthermore, suture zones identified at MacKenzie River, western Alaska/Aleutians and Scandinavia coincide with those identified for other Arctic vertebrates, suggesting that these regions were strong geographic barriers limiting dispersal from Pleistocene refugia.     

Exome resequencing reveals signatures of demographic and adaptive processes across the genome and range of black cottonwood (Populus trichocarpa)

Extant variation in temperate and boreal plant species has been influenced by both demographic histories associated with Pleistocene glacial cycles and adaptation to local climate. We used sequence capture to investigate the role of these neutral and adaptive processes in shaping diversity in black cottonwood (Populus trichocarpa). Nucleotide diversity and Tajima's D were lowest at replacement sites and highest at intergenic sites, while LD showed the opposite pattern. With samples grouped into three populations arrayed latitudinally, effective population size was highest in the north, followed by south and centre, and LD was highest in the south followed by the north and centre, suggesting a possible northern glacial refuge. F_ST outlier analysis revealed that promoter, 5′‐UTR and intronic sites were enriched for outliers compared with coding regions, while no outliers were found among intergenic sites. Codon usage bias was evident, and genes with synonymous outliers had 30% higher average expression compared with genes containing replacement outliers. These results suggest divergent selection related to regulation of gene expression is important to local adaptation in P. trichocarpa. Finally, within‐population selective sweeps were much more pronounced in the central population than in putative northern and southern refugia, which may reflect the different demographic histories of the populations and concomitant effects on signatures of genetic hitchhiking from standing variation.     

乔氏新银鱼基于细胞色素b序列的种群遗传结构和种群历史

为促进乔氏新银鱼(Neosalanx jordani)种质资源的保护, 采集了长江流域和淮河流域5个乔氏新银鱼地理种群计129个样本, 利用线粒体细胞色素b基因(Cyt b)全序列作为分子标记, 初步分析了乔氏新银鱼种群的遗传多样性、遗传结构及种群历史动态。研究结果共检测到18个Cyt b单倍型, 发现和其他鱼类相比, 乔氏新银鱼具有较高的单倍型多样性(h, 0.590±0.047), 但核苷酸多样性较低(π, 0.00088±0.00011)。分子变异分析(AMOVA)表明, 乔氏新银鱼5个地理种群内个体间和流域内种群间均存在显著的遗传差异, 而长江流域与淮河流域之间遗传分化不明显, 显示乔氏新银鱼遗传分化与当前水系的分布格局不吻合。结果表明乔氏新银鱼目前的遗传格局主要是由于长距离独立的建群事件、基因流限制以及种群的持续扩张的共同作用而形成的。种群历史动态分析结果显示乔氏新银鱼种群为近期扩张种群, 其大部分变异发生在1.897万年之内, 与最后一次冰期全面消退、海平面上升、长江流域及淮河流域中下游大量湖泊(适宜生境)形成的时间相吻合。建议对现存的乔氏新银鱼种群, 特别是遗传多样性较高的那些种群(如鄱阳湖、太湖及洪泽湖种群)分别保护。     

Identification of glacial refugia in south-eastern North America by phylogeographical analyses of a forest understorey plant, Trillium cuneatum

Aim We examine several hypotheses emerging from biogeographical and fossil records regarding glacial refugia of a southern thermophilic plant species. Specifically, we investigated the glacial history and post‐glacial colonization of a forest understorey species, Trillium cuneatum. We focused on the following questions: (1) Did T. cuneatum survive the Last Glacial Maximum (LGM) in multiple refugia, and (if so) where were they located, and is the modern genetic structure congruent with the fossil record‐based reconstruction of refugia for mesic deciduous forests? (2) What are the post‐glacial colonization patterns in the present geographical range? Location South‐eastern North America. Methods We sampled 45 populations of T. cuneatum throughout its current range. We conducted phylogeographical analyses based on maternally inherited chloroplast DNA (cpDNA haplotypes) and used TCS software to reconstruct intraspecific phylogeny. Results We detected six cpDNA haplotypes, geographically highly structured into non‐overlapping areas. With one exception, none of the populations had mixed haplotype composition. TCS analysis resulted in two intraspecific cpDNA lineages, with one clade subdivided further by shallower diversification. Main conclusions Our investigation revealed that T. cuneatum survived the LGM in multiple refugia, belonging to two (western, eastern) genealogical lineages geographically structured across south‐eastern North America. The western clade is confined to the south‐western corner of T. cuneatum’s modern range along the Lower Mississippi Valley, where fossil records document a major refugium of mesic deciduous forest. For the eastern clade, modern patterns of cpDNA haplotype distribution suggest cryptic vicariance, in the form of forest contractions and subsequent expansions associated with Pleistocene glacial cycles, rather than simple southern survival and subsequent northward colonization. The north–south partitioning of cpDNA haplotypes was unexpected, suggesting that populations of this rather southern thermophilic species may have survived in more northern locations than initially expected based on LGM climate reconstruction, and that the Appalachian Mountains functioned as a barrier to the dispersal of propagules originating in more southern refugia. Furthermore, our results reveal south‐west to north‐east directionality in historical migration through the Valley and Ridge region of north‐west Georgia.