Genetic resources of common ash (Fraxinus excelsior L.) in Poland

BACKGROUND

Knowledge of genetic structure and the factors that shape it has an impact on forest management practices. European ash (Fraxinus excelsior L.) has declined dramatically throughout its range as a result of a disease caused by the fungus Hymenoscyphus fraxineus. Despite the need for conservation and restoration of the species, genetic data required to guide these efforts at the country level are scarce. Thereofore, we studied the chloroplast and nuclear genetic diversity of 26 natural common ash populations (1269 trees) in Poland.

RESULTS

Chloroplast polymorphisms grouped the populations into two geographically structured phylogenetic lineages ascribed to different glacial refugia (the Balkans and the Eastern Alps). However, the populations demonstrated high genetic diversity (mean AR = 12.35; mean Ho = 0.769; mean He = 0.542) but low differentiation based on nuclear microsatellites (FST = 0.045). Significant spatial genetic structure, consistent with models of isolation by distance, was detected in 14 out of 23 populations. Estimated effective population size was moderate-to-high, with a harmonic mean of 57.5 individuals per population.

CONCLUSIONS

Genetic diversity was not homogeneously distributed among populations within phylogenetic gene pools, indicating that ash populations are not equal as potential sources of reproductive material. Genetic differences among populations could be related to their histories, including founder effects or gene flow between evolutionary lineages (admixture). Our results suggest that ash stands across Poland could be treated as two main management units (seed zones). Therefore, despite the homogenizing effect of pollen gene flow known for this species, the genetic structure should be taken into account in the management of the genetic resources of the common ash. Although ash dieback poses an additional challenge for the management of genetic resources, efforts should be directed towards protecting populations with high genetic diversity within defined phylogenetic units, as they may be an important source of adaptive variation for future stands.

Conserving a threatened North American walnut: a chromosome-scale reference genome for butternut (Juglans cinerea)

With the advent of affordable and more accurate third-generation sequencing technologies, and the associated bioinformatic tools, it is now possible to sequence, assemble, and annotate more species of conservation concern than ever before. Juglans cinerea, commonly known as butternut or white walnut, is a member of the walnut family, native to the Eastern United States and Southeastern Canada. The species is currently listed as Endangered on the IUCN Red List due to decline from an invasive fungus known as Ophiognomonia clavigignenti-juglandacearum (Oc-j) that causes butternut canker. Oc-j creates visible sores on the trunks of the tree which essentially starves and slowly kills the tree. Natural resistance to this pathogen is rare. Conserving butternut is of utmost priority due to its critical ecosystem role and cultural significance. As part of an integrated undergraduate and graduate student training program in biodiversity and conservation genomics, the first reference genome for Juglans cinerea is described here. This chromosome-scale 539 Mb assembly was generated from over 100 × coverage of Oxford Nanopore long reads and scaffolded with the Juglans mandshurica genome. Scaffolding with a closely related species oriented and ordered the sequences in a manner more representative of the structure of the genome without altering the sequence. Comparisons with sequenced Juglandaceae revealed high levels of synteny and further supported J. cinerea's recent phylogenetic placement. Comparative assessment of gene family evolution revealed a significant number of contracting families, including several associated with biotic stress response.

Domestication and selection footprints in Persian walnuts (Juglans regia)

Walnut (Juglans) species are economically important hardwood trees cultivated worldwide for both edible nuts and high-quality wood. Broad-scale assessments of species diversity, evolutionary history, and domestication are needed to improve walnut breeding. In this study, we sequenced 309 walnut accessions from around the world, including 55 Juglans relatives, 98 wild Persian walnuts (J. regia), 70 J. regia landraces, and 86 J. regia cultivars. The phylogenetic tree indicated that J. regia samples (section Dioscaryon) were monophyletic within Juglans. The core areas of genetic diversity of J. regia germplasm were southwestern China and southern Asia near the Qinghai-Tibet Plateau and the Himalayas, and the uplift of the Himalayas was speculated to be the main factor leading to the current population dynamics of Persian walnut. The pattern of genomic variation in terms of nucleotide diversity, linkage disequilibrium, single nucleotide polymorphisms, and insertions/deletions revealed the domestication and selection footprints in Persian walnut. Selective sweep analysis, GWAS, and expression analysis further identified two transcription factors, JrbHLH and JrMYB6, that influence the thickness of the nut diaphragm as loci under selection during domestication. Our results elucidate the domestication and selection footprints in Persian walnuts and provide a valuable resource for the genomics-assisted breeding of this important crop.

Genetic diversity in North American Cercis Canadensis reveals an ancient population bottleneck that originated after the last glacial maximum

Understanding of the present-day genetic diversity, population structure, and evolutionary history of tree species can inform resource management and conservation activities, including response to pressures presented by a changing climate. Cercis canadensis (Eastern Redbud) is an economically valuable understory tree species native to the United States (U.S.) that is also important for forest ecosystem and wildlife health. Here, we document and explain the population genetics and evolutionary history of this deciduous tree species across its distributed range. In this study, we used twelve microsatellite markers to investigate 691 wild-type trees sampled at 74 collection sites from 23 Eastern U.S. states. High genetic diversity and limited gene flow were revealed in wild, natural stands of C. canadensis with populations that are explained by two major genetic clusters. These findings indicate that an ancient population bottleneck occurred coinciding with the last glacial maximum (LGM) in North America. The structure in current populations likely originated from an ancient population in the eastern U.S. that survived LGM and then later diverged into two contemporary clusters. Data suggests that populations have expanded since the last glaciation event from one into several post-glacial refugia that now occupy this species’ current geographic range. Our enhanced understanding benchmarks the genetic variation preserved within this species and can direct future efforts in conservation, and resource utilization of adaptively resilient populations that present the greatest genetic and structural diversity.

A high quality reference genome for Fraxinus pennsylvanica for ash species restoration and research

Green ash (Fraxinus pennsylvanica) is the most widely distributed ash tree in North America. Once common, it has experienced high mortality from the non‐native invasive emerald ash borer (EAB; Agrilus planipennis). A small percentage of native green ash trees that remain healthy in long‐infested areas, termed “lingering ash,” display partial resistance to the insect, indicating that breeding and propagating populations with higher resistance to EAB may be possible. To assist in ash breeding, ecology and evolution studies, we report the first chromosome‐level assembly from the genus Fraxinus for F. pennsylvanica with over 99% of bases anchored to 23 haploid chromosomes, spanning 757 Mb in total, composed of 49.43% repetitive DNA, and containing 35,470 high‐confidence gene models assigned to 22,976 Asterid orthogroups. We also present results of range‐wide genetic variation studies, the identification of candidate genes for important traits including potential EAB‐resistance genes, and an investigation of comparative genome organization among Asterids based on this reference genome platform. Residual duplicated regions within the genome probably resulting from a recent whole genome duplication event in Oleaceae were visualized in relation to wild olive (Olea europaea var. sylvestris). We used our F. pennsylvanica chromosome assembly to construct reference‐guided assemblies of 27 previously sequenced Fraxinus taxa, including F. excelsior. Thus, we present a significant step forward in genomic resources for research and protection of Fraxinus species.

Evolution of Castanea in North America: restriction-site-associated DNA sequencing and ecological modeling reveal a history of radiation, range shifts, and disease

PREMISE

Although chestnuts and chinquapins are some of the best known and most widely loved of any plants in North America, relatively little genomic sequencing has been done, and much is still unknown about their evolution.

METHODS

We used double-digest restriction-site-associated DNA (ddRAD) sequencing data to infer the species-level phylogeny for Castanea and assess the phylogeography of the North American species using samples collected from populations that span the full extent of the species' ranges. We also constructed species distribution models using digitized herbarium specimens and observational data from field surveys.

RESULTS

We identified strong population structure within Castanea dentata (American chestnut) that reflects a stepwise northern migration since the last glacial maximum. Our species distribution models further confirmed this scenario and matched closely with the Castanea fossil pollen record. We also found significant structure within the Castanea pumila lineage, most notably a genetic cluster that corresponds to the frequently recognized Castanea pumila var. ozarkensis.

CONCLUSIONS

The two North American Castanea species have contrasting patterns of population structure, but each is typical of plant phylogeography in North America. Within the C. pumila complex, we found novel genetic structure that provides new insights about C. pumila taxonomy. Our results also identified a series of distinctive populations that will be valuable in ongoing efforts to conserve and restore chestnuts and chinquapins in North America.

Comparing Genetic Diversity in Three Threatened Oaks

Genetic diversity is a critical resource for species’ survival during times of environmental change. Conserving and sustainably managing genetic diversity requires understanding the distribution and amount of genetic diversity (in situ and ex situ) across multiple species. This paper focuses on three emblematic and IUCN Red List threatened oaks (Quercus, Fagaceae), a highly speciose tree genus that contains numerous rare species and poses challenges for ex situ conservation. We compare the genetic diversity of three rare oak species—Quercus georgiana, Q. oglethorpensis, and Q. boyntonii—to common oaks; investigate the correlation of range size, population size, and the abiotic environment with genetic diversity within and among populations in situ; and test how well genetic diversity preserved in botanic gardens correlates with geographic range size. Our main findings are: (1) these three rare species generally have lower genetic diversity than more abundant oaks; (2) in some cases, small population size and geographic range correlate with genetic diversity and differentiation; and (3) genetic diversity currently protected in botanic gardens is inadequately predicted by geographic range size and number of samples preserved, suggesting non-random sampling of populations for conservation collections. Our results highlight that most populations of these three rare oaks have managed to avoid severe genetic erosion, but their small size will likely necessitate genetic management going forward.

Weak founder effects but significant spatial genetic imprint of recent contraction and expansion of European beech populations

Understanding the ecological and evolutionary processes occurring during species range shifts is important in the current context of global change. Here, we investigate the interplay between recent expansion, gene flow and genetic drift, and their consequences for genetic diversity and structure at landscape and local scales in European beech (Fagus sylvatica L.) On Mont Ventoux, South-Eastern France, we located beech forest refugia at the time of the most recent population minimum, ~150 years ago, and sampled 71 populations (2042 trees) in both refugia and expanding populations over an area of 15,000 ha. We inferred patterns of gene flow and genetic structure using 12 microsatellite markers. We identified six plots as originating from planting, rather than natural establishment, mostly from local genetic material. Comparing genetic diversity and structure in refugia versus recent populations did not support the existence of founder effects: heterozygosity (He = 0.667) and allelic richness (Ar = 4.298) were similar, and FST was low (0.031 overall). Still, significant spatial evidence of colonization was detected, with He increasing along the expansion front, while genetic differentiation from the entire pool (βWT) decreased. Isolation by distance was found in refugia but not in recently expanding populations. Our study indicates that beech capacities for colonization and gene flow were sufficient to preserve genetic diversity despite recent forest contraction and expansion. Because beech has long distance pollen and seed dispersal, these results illustrate a 'best case scenario' for the maintenance of high genetic diversity and adaptive potential under climate-change-related range change.

Whole-exome sequencing reveals a long-term decline in effective population size of red spruce (Picea rubens)

Understanding the factors influencing the current distribution of genetic diversity across a species range is one of the main questions of evolutionary biology, especially given the increasing threat to biodiversity posed by climate change. Historical demographic processes such as population expansion or bottlenecks and decline are known to exert a predominant influence on past and current levels of genetic diversity, and revealing this demo-genetic history can have immediate conservation implications. We used a whole-exome capture sequencing approach to analyze polymorphism across the gene space of red spruce (Picea rubens Sarg.), an endemic and emblematic tree species of eastern North America high-elevation forests that are facing the combined threat of global warming and increasing human activities. We sampled a total of 340 individuals, including populations from the current core of the range in northeastern USA and southeastern Canada and from the southern portions of its range along the Appalachian Mountains, where populations occur as highly fragmented mountaintop "sky islands." Exome capture baits were designed from the closely relative white spruce (P. glauca Voss) transcriptome, and sequencing successfully captured most regions on or near our target genes, resulting in the generation of a new and expansive genomic resource for studying standing genetic variation in red spruce applicable to its conservation. Our results, based on over 2 million exome-derived variants, indicate that red spruce is structured into three distinct ancestry groups that occupy different geographic regions of its highly fragmented range. Moreover, these groups show small Ne , with a temporal history of sustained population decline that has been ongoing for thousands (or even hundreds of thousands) of years. These results demonstrate the broad potential of genomic studies for revealing details of the demographic history that can inform management and conservation efforts of nonmodel species with active restoration programs, such as red spruce.

Molecular bases of responses to abiotic stress in trees

Trees are constantly exposed to climate fluctuations, which vary with both time and geographic location. Environmental changes that are outside of the physiological favorable range usually negatively affect plant performance and trigger responses to abiotic stress. Long-living trees in particular have evolved a wide spectrum of molecular mechanisms to coordinate growth and development under stressful conditions, thus minimizing fitness costs. The ongoing development of techniques directed at quantifying abiotic stress has significantly increased our knowledge of physiological responses in woody plants. However, it is only within recent years that advances in next-generation sequencing and biochemical approaches have enabled us to begin to understand the complexity of the molecular systems that underlie these responses. Here, we review recent progress in our understanding of the molecular bases of drought and temperature stresses in trees, with a focus on functional, transcriptomic, epigenetic, and population genomic studies. In addition, we highlight topics that will contribute to progress in our understanding of the plastic and adaptive responses of woody plants to drought and temperature in a context of global climate change.

Is Niagara Falls a barrier to gene flow in riverine fishes? A test using genome-wide SNP data from seven native species

Since the early Holocene, fish population genetics in the Laurentian Great Lakes have been shaped by the dual influences of habitat structure and post-glacial dispersal. Riverscape genetics theory predicts that longitudinal habitat corridors and unidirectional downstream water-flow drive the downstream accumulation of genetic diversity, whereas post-glacial dispersal theory predicts that fish genetic diversity should decrease with increasing distance from glacial refugia. This study examines populations of seven native fish species codistributed above and below the 58 m high Niagara Falls - a hypothesized barrier to gene flow in aquatic species. A better understanding of Niagara Falls' role as a barrier to gene flow and dispersal is needed to identify drivers of Great Lakes genetic diversity and guide strategies to limit exotic species invasions. We used genome-wide SNPs and coalescent models to test whether populations are: (a) genetically distinct, consistent with the Niagara Falls barrier hypothesis; (b) more genetically diverse upstream, consistent with post-glacial expansion theory, or downstream, consistent with the riverscape habitat theory; and (c) have migrated either upstream or downstream past Niagara Falls. We found that genetic diversity is consistently greater below Niagara Falls and the falls are an effective barrier to migration, but two species have probably dispersed upstream past the falls after glacial retreat yet before opening of the Welland Canal. Models restricting migration to after opening of the Welland Canal were generally rejected. These results help explain how river habitat features affect aquatic species' genetic diversity and highlight the need to better understand post-glacial dispersal pathways.

Ecological Status of Juglans cinerea in New Brunswick

Butternut (Juglans cinerea L.), an early successional riparian hardwood species native to Canada and the United States, is under serious threat from a nonnative fungal pathogen, Ophiognomonia clavigignenti-juglandacearum. Since it was first reported in Canada in 1990, this fungal pathogen has spread rapidly and established in New Brunswick in 1997. Apart from the first report in 1997 and another in 2004, no surveys have been conducted to assess the spread of the pathogen in the province. The purpose of this research was to survey butternut throughout its range in New Brunswick, evaluate disease occurrence as well as tree health, and investigate the impact of different topographic and tree health factors on canker incidence. Results showed that the disease has spread throughout the range of butternut in New Brunswick. The disease likely only recently (2007) spread to the northeastern-most populations, given that lower rates of canker occurrence and higher health ratings are found further away from the point of initial occurrence of the disease. Although canker incidence is high throughout the province, tree dieback is minimal, and trees still producing nuts could support opportunities for ex situ conservation. Because of the rate of pathogen spread in the province, implementing a tree improvement strategy might be the only means for maintaining the butternut genome on the landscape.

A rangewide herbarium-derived dataset indicates high levels of gene flow in black cherry (Prunus serotina)

Isolation by Distance (IBD) is a genetic pattern in which populations geographically closer to one another are more genetically similar to each other than populations which are farther apart. Black cherry (Prunus serotina Ehrh.) (Rosaceae) is a forest tree species widespread in eastern North America, and found sporadically in the southwestern United States, Mexico, and Guatemala. IBD has been studied in relatively few North American plant taxa, and no study has rigorously sampled across the range of such a widespread species. In this study, IBD and overall genetic structure were assessed in eastern black cherry (P. serotina Ehrh. var. serotina), the widespread variety of eastern North America. Eastern North America. Prunus serotina Ehrh. var. serotina (Rosaceae). Dense sampling across the entire range of eastern black cherry was made possible by genotyping 15 microsatellite loci in 439 herbarium samples from all portions of the range. Mantel tests and STRUCTURE analyses were performed to evaluate the hypothesis of IBD and genetic structure. Mantel tests demonstrated significant but weak IBD, while STRUCTURE analyses revealed no clear geographic pattern of genetic groups. The modest geographic/genetic structure across the eastern black cherry range suggests widespread gene flow in this taxon. This is consistent with P. serotina's status as a disturbance-associated species. Further studies should similarly evaluate IBD in species characteristic of low-disturbance forests.

Next-generation conservation genetics and biodiversity monitoring

This special issue of Evolutionary Applications consists of 10 publications investigating the use of next-generation tools and techniques in population genetic analyses and biodiversity assessment. The special issue stems from a 2016 Next Generation Genetic Monitoring Workshop, hosted by the National Institute for Mathematical and Biological Synthesis (NIMBioS) in Tennessee, USA. The improved accessibility of next-generation sequencing platforms has allowed molecular ecologists to rapidly produce large amounts of data. However, with the increased availability of new genomic markers and mathematical techniques, care is needed in selecting appropriate study designs, interpreting results in light of conservation concerns, and determining appropriate management actions. This special issue identifies key attributes of successful genetic data analyses in biodiversity evaluation and suggests ways to improve analyses and their application in current population and conservation genetics research.

Factors Influencing the Regional Dynamics of Butternut Canker

Butternut (Juglans cinerea) is an important component of native biodiversity in eastern North America. Of urgent concern is the survival of butternut, whose populations are declining rapidly, in large part due to an exotic pathogen, Ophiognomonia clavigignenti-juglandacearum, that causes butternut canker. The disease presently occurs throughout the range of butternut in North America, causing branch and stem cankers, dieback, and tree mortality. Despite the existential threat posed by O. clavigignenti-juglandacearum to butternut, a detailed understanding of the factors that drive cross-scale disease patterns is lacking. Therefore, we investigated the association of a range of factors, including tree attributes, topography, and weather, with butternut canker spatial dynamics at different scales using data collected in the province of Quebec, Canada. Trunk canker damage and dieback showed distinct geographic patterns. Bark phenotype was not significantly associated with trunk canker damage. Results suggest that open or dominant trees may show less dieback than intermediate or suppressed trees. Probability of the presence of trunk canker and percent dieback were proportional to the tree diameter at breast height. Temperature was positively associated with disease severity at a 1-km² scale. Our results provide strong evidence that multiple factors, notably weather, influence butternut canker epidemiology.

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 (AR) and population differentiation (FST). 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.

Demographically idiosyncratic responses to climate change and rapid Pleistocene diversification of the walnut genus Juglans (Juglandaceae) revealed by whole-genome sequences

Whether species demography and diversification are driven primarily by extrinsic environmental changes such as climatic oscillations in the Quaternary or by intrinsic biological interactions like coevolution between antagonists is a matter of active debate. In fact, their relative importance can be assessed by tracking past population fluctuations over considerable time periods. We applied the pairwise sequentially Markovian coalescent approach on the genomes of 11 temperate Juglans species to estimate trajectories of changes in effective population size (N_e ) and used a Bayesian-coalescent based approach that simultaneously considers multiple genomes (G-PhoCS) to estimate divergence times between lineages. N_e curves of all study species converged 1.0 million yr ago, probably reflecting the time when the walnut genus last shared a common ancestor. This estimate was confirmed by the G-PhoCS estimates of divergence times. But all species did not react similarly to the dramatic climatic oscillations following early Pleistocene cooling, so the timing and amplitude of changes in N_e differed among species and even among conspecific lineages. The population histories of temperate walnut species were not driven by extrinsic environmental changes alone, and a key role was probably played by species-specific factors such as coevolutionary interactions with specialized pathogens.

Ice ages leave genetic diversity 'hotspots' in Europe but not in Eastern North America

After the last glacial cycle, temperate European trees migrated northward, experiencing genetic bottlenecks and founder effects, which left high haplotype endemism in southern populations and clines in genetic diversity northward. These patterns are thought to be ubiquitous across temperate forests, and are therefore used to anticipate the potential genetic consequences of future warming. We compared existing and new phylogeographic data sets (chloroplast DNA) from 14 woody taxa in Eastern North America (ENA) to data sets from 21 ecologically similar European species to test for common impacts of Quaternary climate swings on genetic diversity across diverse taxa and between continents. Unlike their European counterparts, ENA taxa do not share common southern centres of haplotype endemism and they generally maintain high genetic diversity even at their northern range limits. Differences between the genetic impacts of Quaternary climate cycles across continents suggest refined lessons for managing genetic diversity in today's warming world.

Back to America: tracking the origin of European introduced populations of Quercus rubra L

Quercus rubra has been introduced in Europe since the end of the 17th century. It is widely distributed today across this continent and considered invasive in some countries. Here, we investigated the distribution of genetic diversity of both native and introduced populations with the aim of tracing the origin of introduced populations. A large sampling of 883 individuals from 73 native and 38 European locations were genotyped at 69 SNPs. In the natural range, we found a continuous geographic gradient of variation with a predominant latitudinal component. We explored the existence of ancestral populations by performing Bayesian clustering analysis and found support for two or three ancestral genetic clusters. Approximate Bayesian Computations analyses based on these two or three clusters support recent extensive secondary contacts between them, suggesting that present-day continuous genetic variation resulted from recent admixture. In the introduced range, one main genetic cluster was not recovered in Europe, suggesting that source populations were preferentially located in the northern part of the natural distribution. However, our results cannot refute the introduction of populations from the southern states that did not survive in Europe.

Population Differentiation in Common Walnut (Juglans regia L.) across Major Parts of Its Native Range-Insights from Molecular and Morphometric Data

Juglans regia is an economically highly important species for fruit and wood production in the warm temperate and subtropical zones of the Northern Hemisphere. Besides the natural influence of climatic and geomorphological barriers, its genetic structure has been strongly modified by humans and the population history is still unclear. For this reason, we investigated mainly natural walnut populations across the Eurasian continent on a molecular (44 populations, 581 trees) and morphometric level (23 populations, 1391 ripe nuts). Population genetic diversity and differentiation were examined by using 7 microsatellite loci. Morphometric characteristics of the nuts (mainly roundness index and nut density) were used to estimate trait variation and population differentiation. Highest allelic richness Rs12 = 7.05 was observed in a Pakistani and the lowest value Rs12 = 3.04 in a Kyrgyz population. The genetic differentiation among populations was high (FST = 0.217; RST = 0.530) indicating a strong phylogeographic pattern. While variation of the roundness index within single populations was high, this trait neither differentiated geographical regions nor was it associated to genetic clusters. Approximated QST based on this trait equalled FST, while approximated QST based on nut density considerably exceeded FST, indicating selection. Nut density was moderately correlated with altitude, latitude, and longitude, and differentiated populations according to their origin. Pakistani and Indian populations showed highest nut densities. These South Asian populations contain putatively ancestral nut forms, which probably have been lost in other populations as a consequence of human selection.

High-quality genetic mapping with ddRADseq in the non-model tree Quercus rubra

Background

Restriction site associated DNA sequencing (RADseq) has the potential to be a broadly applicable, low-cost approach for high-quality genetic linkage mapping in forest trees lacking a reference genome. The statistical inference of linear order must be as accurate as possible for the correct ordering of sequence scaffolds and contigs to chromosomal locations. Accurate maps also facilitate the discovery of chromosome segments containing allelic variants conferring resistance to the biotic and abiotic stresses that threaten forest trees worldwide. We used ddRADseq for genetic mapping in the tree Quercus rubra, with an approach optimized to produce a high-quality map. Our study design also enabled us to model the results we would have obtained with less depth of coverage.

Results

Our sequencing design produced a high sequencing depth in the parents (248×) and a moderate sequencing depth (15×) in the progeny. The digital normalization method of generating a de novo reference and the SAMtools SNP variant caller yielded the most SNP calls (78,725). The major drivers of map inflation were multiple SNPs located within the same sequence (77% of SNPs called). The highest quality map was generated with a low level of missing data (5%) and a genome-wide threshold of 0.025 for deviation from Mendelian expectation. The final map included 849 SNP markers (1.8% of the 78,725 SNPs called). Downsampling the individual FASTQ files to model lower depth of coverage revealed that sequencing the progeny using 96 samples per lane would have yielded too few SNP markers to generate a map, even if we had sequenced the parents at depth 248×.

Conclusions

The ddRADseq technology produced enough high-quality SNP markers to make a moderately dense, high-quality map. The success of this project was due to high depth of coverage of the parents, moderate depth of coverage of the progeny, a good framework map, an optimized bioinformatics pipeline, and rigorous premapping filters. The ddRADseq approach is useful for the construction of high-quality genetic maps in organisms lacking a reference genome if the parents and progeny are sequenced at sufficient depth. Technical improvements in reduced representation sequencing (RRS) approaches are needed to reduce the amount of missing data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3765-8) contains supplementary material, which is available to authorized users.

Genetic Structure in the Northern Range Margins of Common Ash, Fraxinus excelsior L

During post glacial colonization, loss of genetic diversity due to leading edge effects may be attenuated in forest trees because of their prolonged juvenile phase, allowing many migrants to reach the colonizing front before populations become reproductive. The northern range margins of temperate tree taxa in Europe are particularly suitable to study the genetic processes that follow colonization because they have been little affected by northern refugia. Here we examined how post glacial range dynamics have shaped the genetic structure of common ash (Fraxinus excelsior L.) in its northern range compared to its central range in Europe. We used four chloroplast and six nuclear microsatellites to screen 42 populations (1099 trees), half of which corresponded to newly sampled populations in the northern range and half of which represented reference populations from the central range obtained from previously studies. We found that northern range populations of common ash have the same chloroplast haplotypes as south-eastern European populations, suggesting that colonization of the northern range took place along a single migration route, a result confirmed by the structure at the nuclear microsatellites. Along this route, diversity strongly decreased only in the northern range, concomitantly with increasing population differentiation and complex population substructures, a pattern consistent with a leading edge colonization model. Our study highlights that while diversity is maintained in the central range of common ash due to broad colonizing fronts and high levels of gene flow, it profoundly decreases in the northern range, where colonization was unidirectional and probably involved repeated founder events and population fluctuations. Currently, common ash is threatened by ash dieback, and our results on northern populations will be valuable for developing gene conservation strategies.

Phylogeography of postglacial range expansion in Juglans mandshurica (Juglandaceae) reveals no evidence of bottleneck, loss of genetic diversity, or isolation by distance in the leading-edge populations

The past studies of postglacial recolonization patterns in high latitude regions have revealed a significant role of dispersal capacity in shaping the genetic diversity and population structure of temperate trees. However, most of these studies have focused on species with long-distance dispersal followed by exponential population growth and were therefore unable to reveal the patterns in the case of a gradual expansion. Here we studied the impacts of postglacial range expansions on the distribution of genetic diversity in the Manchurian walnut (Juglans mandshurica), a common tree of East Asian cool-temperate deciduous forests that apparently lacks long-distance seed dispersal ability. The genetic diversity and structure of 19 natural walnut populations in Northeast China and the Korean Peninsula were examined using 17 nuclear simple sequence repeat (SSR) loci. Potential habitats under current and past climatic conditions were predicted using the ecological niche modelling (ENM) method. Bayesian clustering analysis revealed three groups, which were inferred to have diverged through multiple glacial-interglacial cycles in multiple refugia during the Quaternary Period. ENM estimated a southward range shift at the LGM, but high suitability scores still occurred in the western parts of the Changbai Mountains (Northeast China), the Korean peninsula and the exposed seafloor of the Yellow Sea. In contrast to most other cool-temperate trees co-occurring in the same region, the Manchurian walnut did not show any evidence of a population bottleneck, loss of genetic diversity or isolation by distance during the postglacial expansion. Our study clearly indicates that current northern populations originated from one glacial lineage and recolonization via a gradually advancing front due to the lack of a long-distance seed dispersal mechanism led to no latitudinal decrease in genetic diversity.

High genetic diversity declines towards the geographic range periphery of Adonis vernalis, a Eurasian dry grassland plant

Genetic diversity is important for species' fitness and evolutionary processes but our knowledge on how it varies across a species' distribution range is limited. The abundant centre hypothesis (ACH) predicts that populations become smaller and more isolated towards the geographic range periphery - a pattern that in turn should be associated with decreasing genetic diversity and increasing genetic differentiation. We tested this hypothesis in Adonis vernalis, a dry grassland plant with an extensive Eurasian distribution. Its life-history traits and distribution characteristics suggest a low genetic diversity that decreases and a high genetic differentiation that increases towards the range edge. We analysed AFLP fingerprints in 28 populations along a 4698-km transect from the geographic range core in Russia to the western range periphery in Central and Western Europe. Contrary to our expectation, our analysis revealed high genetic diversity (range of proportion of polymorphic bands = 56-81%, He = 0.168-0.238) and low genetic differentiation across populations (Φ(ST) = 0.18). However, in congruence with the genetic predictions of the ACH, genetic diversity decreased and genetic differentiation increased towards the range periphery. Spanish populations were genetically distinct, suggesting a divergent post-glacial history in this region. The high genetic diversity and low genetic differentiation in the remaining A. vernalis populations is surprising given the species' life-history traits and points to the possibility that the species has been widely distributed in the studied region or that it has migrated from a diverse source in an East-West direction, in the past.

Contrasting patterns of genetic diversity across the ranges of Pinus monticola and P. strobus: a comparison between eastern and western North American postglacial colonization histories

• Premises of the study: Understanding the influence of recent glacial and postglacial periods on species' distributions is key for predicting the effects of future environmental changes. We investigated the influence of two physiographic landscapes on population structure and postglacial colonization of two white pine species of contrasting habitats: P. monticola, which occurs in the highly mountainous region of western North America, and P. strobus, which occurs in a much less mountainous area in eastern North America.•

METHODS

To characterize the patterns of genetic diversity and population structure across the ranges of both species, 158 and 153 single nucleotide polymorphism (SNP) markers derived from expressed genes were genotyped on range-wide samples of 61 P. monticola and 133 P. strobus populations, respectively.•

KEY RESULTS

In P. monticola, a steep latitudinal decrease in genetic diversity likely resulted from postglacial colonization involving rare long-distance dispersal (LDD) events. In contrast, no geographic patterns of diversity were detected in P. strobus, suggesting recolonization via a gradually advancing front or frequent LDD events. For each species, structure analyses identified two distinct southern and northern genetic groups that likely originated from two different glacial lineages. At a finer scale, and for the two species, smaller subgroups were detected that could be remnants of cryptic refugia.•

CONCLUSION

During postglacial colonization, the western and eastern North American landscapes had different impacts on genetic signatures in P. monticola compared with P. strobus. We discuss the importance of our findings for conservation programs and predictions of species' response to climate change.

Holocene re-colonisation, central-marginal distribution and habitat specialisation shape population genetic patterns within an Atlantic European grass species

Corynephorus canescens (L.) P.Beauv. is an outbreeding, short-lived and wind-dispersed grass species, highly specialised on scattered and disturbance-dependent habitats of open sandy sites. Its distribution ranges from the Iberian Peninsula over Atlantic regions of Western and Central Europe, but excludes the two other classical European glacial refuge regions on the Apennine and Balkan Peninsulas. To investigate genetic patterns of this uncommon combination of ecological and biogeographic species characteristics, we analysed AFLP variation among 49 populations throughout the European distribution range, expecting (i) patterns of SW European glacial refugia and post-glacial expansion to the NE; (ii) decreasing genetic diversity from central to marginal populations; and (iii) interacting effects of high gene flow and disturbance-driven genetic drift. Decreasing genetic diversity from SW to NE and distinct gene pool clustering imply refugia on the Iberian Peninsula and in western France, from where range expansion originated towards the NE. High genetic diversity within and moderate genetic differentiation among populations, and a significant pattern of isolation-by-distance indicate a gene flow drift equilibrium within C. canescens, probably due to its restriction to scattered and dynamic habitats and limited dispersal distances. These features, as well as the re-colonisation history, were found to affect genetic diversity gradients from central to marginal populations. Our study emphasises the need for including the specific ecology into analyses of species (re-)colonisation histories and range centre-margin analyses. To account for discontinuous distributions, new indices of marginality were tested for their suitability in studies of centre-periphery gradients.

Genetic structuring of remnant forest patches in an endangered medicinal tree in North-western Ethiopia

Background

Habitat loss and fragmentation may have detrimental impacts on genetic diversity, population structure and overall viability of tropical trees. The response of tropical trees to fragmentation processes may, however, be species, cohort or region-specific. Here we test the hypothesis that forest fragmentation is associated with lower genetic variability and higher genetic differentiation in adult and seedling populations of Prunus africana in North-western Ethiopia. This is a floristically impoverished region where all but a few remnant forest patches have been destroyed, mostly by anthropogenic means.

Results

Genetic diversity (based on allelic richness) was significantly greater in large and less-isolated forest patches as well as in adults than seedlings. Nearly all pairwise F_ST comparisons showed evidence for significant population genetic differentiation. Mean F_ST values were significantly greater in seedlings than adults, even after correction for within population diversity, but varied little with patch size or isolation.

Conclusions

Analysis of long-lived adult trees suggests the formerly contiguous forest in North-western Ethiopia probably exhibited strong spatial patterns of genetic structure. This means that protecting a range of patches including small and isolated ones is needed to conserve the extant genetic resources of the valuable forests in this region. However, given the high livelihood dependence of the local community and the high impact of foreign investors on forest resources of this region, in situ conservation efforts alone may not be helpful. Therefore, these efforts should be supported with ex situ gene conservation actions.

The number of markers and samples needed for detecting bottlenecks under realistic scenarios, with and without recovery: a simulation-based study

Detecting bottlenecks is a common task in molecular ecology. While several bottleneck detection method sexist, evaluations of their power have focused only on severe bottlenecks (e.g. to Ne ~10). As a component of a recent review, Peery et al. (2012) analysed the power of two approaches, the M-ratio and heterozygote excess tests, to detect moderate bottlenecks (e.g. to Ne ~100),which is realistic for many conservation situations. In this Comment, we address three important points relevant to but not considered in Peery et al. Under moderate bottleneck scenarios, we test the (i) relative advantage of sampling more markers vs. more individuals, (ii) potential power to detect the bottleneck when utilizing dozens of microsatellites (a realistic possibility for contemporary studies) and (iii) reduction in power when post bottle neck recovery has occurred. For the realistic situations examined,we show that (i) doubling the number of loci shows equal or better power than tripling the number of individuals,(ii) increasing the number of markers (up to 100) results in continued additive gains in power, and (iii)recovery after a moderate amount of time or gradual change in size reduces power, by up to one-half. Our results provide a practical supplement to Peery et al. and encourage the continued use of bottleneck detection methods in the genomic age, but also emphasize that the power under different sampling schemes should be estimated,using simulation modelling, as a routine component of molecular ecology studies.

Phylogeography sheds light on the central-marginal hypothesis in a Mediterranean narrow endemic plant

BACKGROUND AND AIMS

Understanding the factors that shape variation in genetic diversity across the geographic ranges of species is an important challenge in the effort to conserve evolutionary processes sustaining biodiversity. The historical influences leading to a central-marginal organization of genetic diversity have been explored for species whose range is known to have expanded from refugia after glacial events. However, this question has rarely been addressed for Mediterranean endemic plants of azonal habitats such as rocky slopes or screes. In this context, this comprehensive study examined molecular and field data from Arenaria provincialis (Caryophyllaceae), a narrow endemic plant of south-eastern France.

METHODS

Across the whole geographic range, an investigation was made of whether high levels of abundance and genetic diversity (estimated from amplified fragment length polymorphism markers) are centrally distributed, to evaluate the relevance of the central-marginal hypothesis. Phylogeographic patterns inferred from chloroplast DNA (cpDNA) were used, applying Bayesian methods to test the influence of past biogeographic events. Multivariate analysis combining phylogeographic and ecological data was used to reveal the historical and ecological distinctiveness of populations.

KEY RESULTS

Despite the narrow distribution of A. provincialis, a high level of nucleotide variation is found within cpDNA loci, supporting its persistence throughout the Pleistocene period. The area characterized by the highest genetic diversity is centrally located. Structured phylogeography and Bayesian factor analysis supported the hypothesis that the central area of the distribution was the source of both westward and eastward migrations, probably during arid periods of the Pleistocene, and more recently was a crossroads of backward migrations. By contrast, the two areas located today at the range limits are younger, have reduced genetic diversity and are marginal in the ecological gradients.

CONCLUSIONS

This study highlights a case of strong population distinctiveness within a narrow range. Phylogeography sheds light on the historical role of the areas centrally situated in the distribution. The current range size and abundance patterns are not sufficient to predict the organization of genetic diversity.

High variance in reproductive success generates a false signature of a genetic bottleneck in populations of constant size: a simulation study

Background

Demographic bottlenecks can severely reduce the genetic variation of a population or a species. Establishing whether low genetic variation is caused by a bottleneck or a constantly low effective number of individuals is important to understand a species’ ecology and evolution, and it has implications for conservation management. Recent studies have evaluated the power of several statistical methods developed to identify bottlenecks. However, the false positive rate, i.e. the rate with which a bottleneck signal is misidentified in demographically stable populations, has received little attention. We analyse this type of error (type I) in forward computer simulations of stable populations having greater than Poisson variance in reproductive success (i.e., variance in family sizes). The assumption of Poisson variance underlies bottleneck tests, yet it is commonly violated in species with high fecundity.

Results

With large variance in reproductive success (V_k ≥ 40, corresponding to a ratio between effective and census size smaller than 0.1), tests based on allele frequencies, allelic sizes, and DNA sequence polymorphisms (heterozygosity excess, M-ratio, and Tajima’s D test) tend to show erroneous signals of a bottleneck. Similarly, strong evidence of population decline is erroneously detected when ancestral and current population sizes are estimated with the model based method MSVAR.

Conclusions

Our results suggest caution when interpreting the results of bottleneck tests in species showing high variance in reproductive success. Particularly in species with high fecundity, computer simulations are recommended to confirm the occurrence of a population bottleneck.

Do pathogens reduce genetic diversity of their hosts? Variable effects of sylvatic plague in black-tailed prairie dogs

Introduced diseases can cause dramatic declines in-and even the loss of-natural populations. Extirpations may be followed by low recolonization rates, leading to inbreeding and a loss of genetic variation, with consequences on population viability. Conversely, extirpations may create vacant habitat patches that individuals from multiple source populations can colonize, potentially leading to an influx of variation. We tested these alternative hypotheses by sampling 15 colonies in a prairie dog metapopulation during 7 years that encompassed an outbreak of sylvatic plague, providing the opportunity to monitor genetic diversity before, during and after the outbreak. Analysis of nine microsatellite loci revealed that within the metapopulation, there was no change in diversity. However, within extirpated colonies, patterns varied: In half of the colonies, allelic richness after recovery was less than the preplague conditions, and in the other half, richness was greater than the preplague conditions. Finally, analysis of variation within individuals revealed that prairie dogs present in recolonized colonies had higher heterozygosity than those present before plague. We confirmed plague survivorship in six founders; these individuals had significantly higher heterozygosity than expected by chance. Collectively, our results suggest that high immigration rates can maintain genetic variation at a regional scale despite simultaneous extirpations in spatially proximate populations. Thus, virulent diseases may increase genetic diversity of host populations by creating vacant habitats that allow an influx of genetic diversity. Furthermore, even highly virulent diseases may not eliminate individuals randomly; rather, they may selectively remove the most inbred individuals.

Human-impacted landscapes facilitate hybridization between a native and an introduced tree

Spatial and temporal dynamics of hybridization, in particular the influence of local environmental conditions, are well studied for sympatric species but less is known for native-introduced systems, especially for long-lived species. We used microsatellite and chloroplast DNA markers to characterize the influence of anthropogenic landscapes on the extent, direction, and spatial distribution of hybridization between a native North American tree Juglans cinerea (butternut) and an introduced tree Juglans ailantifolia (Japanese walnut) for 1363 trees at 48 locations across the native range of butternut. Remarkably, admixture in anthropogenic sites reached nearly 70%, while fragmented and continuous forests showed minimal admixture (<8%). Furthermore, more hybrids in anthropogenic sites had J. ailantifolia seed parents (95%) than hybrids in fragmented and continuous forests (69% and 59%, respectively). Our results show a strong influence of landscape type on rate and direction of realized gene flow. While hybrids are common in anthropogenic landscapes, our results suggest that even small forested landscapes serve as substantial barriers to hybrid establishment, a key consideration for butternut conservation planning, a species already exhibiting severe decline, and for other North American forest trees that hybridize with introduced congeners.

Population structure of the butternut canker fungus, Ophiognomonia clavigignenti-juglandacearum, in North American forests

The occurrence of multiple introduction events, or sudden emergence from a host jump, of forest pathogens may be an important factor in successful establishment in a novel environment or on a new host; however, few studies have focused on the introduction and emergence of fungal pathogens in forest ecosystems. While Ophiognomonia clavigignenti-juglandacearum (Oc-j), the butternut canker fungus, has caused range-wide mortality of butternut trees in North America since its first observation in 1967, the history of its emergence and spread across the United States and Canada remains unresolved. Using 17 single nucleotide polymorphic loci, we investigated the genetic population structure of 101 isolates of Oc-j from across North America. Clustering analysis revealed that the Oc-j population in North America is made up of three differentiated genetic clusters of isolates, and these genetic clusters were found to have a strong clonal structure. These results, in combination with the geographic distribution of the populations, suggest that Oc-j was introduced or has emerged in North America on more than one occasion, and these clonal lineages have since proliferated across much of the range of butternut. No evidence of genetic recombination was observed in the linkage analysis, and conservation of the distinct genetic clusters in regions where isolates from two or more genetic clusters are present, would indicate a very minimal or non-existent role of sexual recombination in populations of Oc-j in North America.

Genetic consequences of fragmentation in "arbor vitae," eastern white cedar (Thuja occidentalis L.), toward the northern limit of its distribution range

We tested the hypothesis that marginal fragmented populations of eastern white cedar (EWC) are genetically isolated due to reduced pollen and gene flow. In accordance with the central-marginal model, we predicted a decrease in population genetic diversity and an increase in differentiation along the latitudinal gradient from the boreal mixed-wood to northern coniferous forest. A total of 24 eastern white cedar populations were sampled along the north-south latitudinal gradient for microsatellite genotyping analysis. Positive F_is values and heterozygote deficiency were observed in populations from the marginal (F_is = 0.244; P_HW = 0.0042) and discontinuous zones (F_is = 0.166; P_HW = 0.0042). However, populations from the continuous zone were in HW equilibrium (F_is = −0.007; P_HW = 0.3625). There were no significant latitudinal effects on gene diversity (H_s), allelic richness (AR), or population differentiation (F_st). Bayesian and NJT (neighbor-joining tree) analyses demonstrated the presence of a population structure that was partly consistent with the geographic origins of the populations. The impact of population fragmentation on the genetic structure of EWC is to create a positive inbreeding coefficient, which was two to three times higher on average than that of a population from the continuous zone. This result indicated a higher occurrence of selfing within fragmented EWC populations coupled with a higher degree of gene exchange among near-neighbor relatives, thereby leading to significant inbreeding. Increased population isolation was apparently not correlated with a detectable effect on genetic diversity. Overall, the fragmented populations of EWC appear well-buffered against effects of inbreeding on genetic erosion.

Current and historical drivers of landscape genetic structure differ in core and peripheral salamander populations

With predicted decreases in genetic diversity and greater genetic differentiation at range peripheries relative to their cores, it can be difficult to distinguish between the roles of current disturbance versus historic processes in shaping contemporary genetic patterns. To address this problem, we test for differences in historic demography and landscape genetic structure of coastal giant salamanders (Dicamptodon tenebrosus) in two core regions (Washington State, United States) versus the species' northern peripheral region (British Columbia, Canada) where the species is listed as threatened. Coalescent-based demographic simulations were consistent with a pattern of post-glacial range expansion, with both ancestral and current estimates of effective population size being much larger within the core region relative to the periphery. However, contrary to predictions of recent human-induced population decline in the less genetically diverse peripheral region, there was no genetic signature of population size change. Effects of current demographic processes on genetic structure were evident using a resistance-based landscape genetics approach. Among core populations, genetic structure was best explained by length of the growing season and isolation by resistance (i.e. a ‘flat’ landscape), but at the periphery, topography (slope and elevation) had the greatest influence on genetic structure. Although reduced genetic variation at the range periphery of D. tenebrosus appears to be largely the result of biogeographical history rather than recent impacts, our analyses suggest that inherent landscape features act to alter dispersal pathways uniquely in different parts of the species' geographic range, with implications for habitat management.

Consequences of range contractions and range shifts on molecular diversity

Due to past and current climatic changes, range contractions and range shifts are essential stages in the history of a species. However, unlike range expansions, the molecular consequences of these processes have been little investigated. In order to fill this gap, we simulated patterns of molecular diversity within and between populations for various types of range contractions and range shifts. We show that range contractions tend to decrease genetic diversity as compared with population with stable ranges but quite counterintuitively fast range contractions preserve higher levels of diversity and induce lower levels of genetic differentiation among refuge areas than slow contractions. Contrastingly, fast range shifts lead to lower levels of diversity than slow range shifts. At odds with our expectations, we find that species actively migrating toward refuge areas can only preserve higher levels of diversity in refugia if the contraction is rapid. Under slow range contraction or slow range shift, active migration toward refugia lead to a larger loss of diversity as compared with scenarios with isotropic migration and may thus not be a good evolutionary strategy. These results suggest that the levels of diversity preserved after a climate change both within and between refuge areas will not only depend on the dispersal abilities of a species but also on the speed of the change. It also implies that a given episode of climatic change will impact differently species with different generation times.