Detection and visualization of spatial genetic structure in continuous Eucalyptus globulus forest

Genome-wide analysis highlights genetic admixture in exotic germplasm resources of Eucalyptus and unexpected ancestral genomic composition of interspecific hybrids

Eucalyptus is an economically important genus comprising more than 890 species in different subgenera and sections. Approximately twenty species of subgenus Symphyomyrtus account for 95% of the world's planted eucalypts. Discrimination of closely related eucalypt taxa is challenging, consistent with their recent phylogenetic divergence and occasional hybridization in nature. Admixture, misclassification or mislabeling of Eucalyptus germplasm resources maintained as exotics have been suggested, although no reports are available. Moreover, hybrids with increased productivity and traits complementarity are planted worldwide, but little is known about their actual genomic ancestry. In this study we examined a set of 440 trees of 16 different Eucalyptus species and 44 interspecific hybrids of multi-species origin conserved in germplasm banks in Brazil. We used genome-wide SNP data to evaluate the agreement between the alleged phylogenetic classification of species and provenances as registered in their historical records, and their observed genetic clustering derived from SNP data. Genetic structure analyses correctly assigned each of the 16 species to a different cluster although the PCA positioning of E. longirostrata was inconsistent with its current taxonomy. Admixture was present for closely related species' materials derived from local germplasm banks, indicating unintended hybridization following germplasm introduction. Provenances could be discriminated for some species, indicating that SNP-based discrimination was directly proportional to geographical distance, consistent with an isolation-by-distance model. SNP-based genomic ancestry analysis showed that the majority of the hybrids displayed realized genomic composition deviating from the expected ones based on their pedigree records, consistent with admixture in their parents and pervasive genome-wide directional selection toward the fast-growing E. grandis genome. SNP data in support of tree breeding provide precise germplasm identity verification, and allow breeders to objectively recognize the actual ancestral origin of superior hybrids to more realistically guide the program toward the development of the desired genetic combinations.

The Genetics and Ecology of Post-Fire Eucalyptus globulus Recruitment in an Isolated Stand in Central Portugal

Eucalyptus globulus Labill. is a widespread forest tree species, reported as naturalized across the introduced range, often showing abundant natural regeneration after wildfires. This paper studies a post-fire cohort of seedlings derived from a small, isolated E. globulus stand in central Portugal. The aim is to better understand the genetic dynamics and dispersal mechanisms of naturally established E. globulus populations in the introduced range. The seedling density at 55 m from adult trees was 12,000 ha−1, the farthest seedling being registered at 101 m. Post-fire expansion occurred in a southward direction, in accordance with predominant wind. Seedlings had significantly lower levels of genetic diversity (Ae = 5.8.; He = 0.8) than adult trees (Ae = 6.5; He = 0.8). Crossings were strongly unbalanced, with only eight trees contributing to the sampled seedlings, and one single tree contributing to 52% of these seedlings. Moreover, the co-ancestry amongst seedlings more than doubled (from θ = 0.03 to 0.076), and the population status number was around one third of the value registered for the adult population (from Ns = 16.2 to 6.6). Despite its isolation, external pollen was detected in the stand, and appears to have travelled at least 700 m, contributing to 22% of the offspring. Overall, the seedling cohort is much less diverse than the parent trees, with expected lower resilience and persistence to environmental stresses.

Climate Adaptation, Drought Susceptibility, and Genomic-Informed Predictions of Future Climate Refugia for the Australian Forest Tree Eucalyptus globulus

Understanding the capacity of forest tree species to adapt to climate change is of increasing importance for managing forest genetic resources. Through a genomics approach, we modelled spatial variation in climate adaptation within the Australian temperate forest tree Eucalyptus globulus, identified putative climate drivers of this genomic variation, and predicted locations of future climate refugia and populations at-risk of future maladaptation. Using 812,158 SNPs across 130 individuals from 30 populations (i.e., localities) spanning the species’ natural range, a gradientForest algorithm found 1177 SNPs associated with locality variation in home-site climate (climate-SNPs), putatively linking them to climate adaptation. Very few climate-SNPs were associated with population-level variation in drought susceptibility, signalling the multi-faceted nature and complexity of climate adaptation. Redundancy analysis (RDA) showed 24% of the climate-SNP variation could be explained by annual precipitation, isothermality, and maximum temperature of the warmest month. Spatial predictions of the RDA climate vectors associated with climate-SNPs allowed mapping of genomically informed climate selective surfaces across the species’ range under contemporary and projected future climates. These surfaces suggest over 50% of the current distribution of E. globulus will be outside the modelled adaptive range by 2070 and at risk of climate maladaptation. Such surfaces present a new integrated approach for natural resource managers to capture adaptive genetic variation and plan translocations in the face of climate change.

Landscape drivers of genomic diversity and divergence in woodland Eucalyptus

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the "eucalypts") are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole-genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome-wide, interspecific differentiation (FST  = 0.15) and intraspecific differentiation between localities (FST  ≈ 0.01-0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.

Life cycle expression of inbreeding depression in Eucalyptus regnans and inter-generational stability of its mixed mating system

BACKGROUND AND AIMS

Many plants exhibit a mixed mating system. Published models suggest that this might be an evolutionarily stable rather than a transitional state despite the presence of inbreeding depression, but there is little empirical evidence. Through field experimentation, we studied the role of inbreeding depression in eliminating inbred progeny from the reproductive cohort of the forest tree Eucalyptus regnans, and demonstrate a stable mixed primary mating system over two successive generations.

METHODS

Two field experiments were conducted using seed from natural populations. We sowed open-pollinated seeds to simulate a natural regeneration event and determined isozyme genotypes of dominant and suppressed individuals over 10 years. We also planted a mixture of open-pollinated, outcross and selfed families with common maternal parentage; monitored survival of cross types over 29 years; and determined the percentage of outcrosses in open-pollinated seed from a sample of reproductively mature trees using microsatellite analysis.

KEY RESULTS

Both experiments demonstrated progressive competitive elimination of inbred plants. By 29 years, the reproductive cohort in the planted experiment consisted only of outcrosses which produced seed which averaged 66 % outcrosses, similar to the estimate for the parental natural population (74 %).

CONCLUSIONS

Selective elimination of inbred genotypes during the intense intra-specific competition characteristic of the pre-reproductive phase of the life cycle of E. regnans results in a fully outcrossed reproductive population, in which self-fertility is comparable with that of its parental generation. The mixed mating system may be viewed as an unavoidable consequence of the species' reproductive ecology, which includes the demonstrated effects of inbreeding depression, rather than a strategy which is actively favoured by natural selection.

A multidisciplinary approach to inform assisted migration of the restricted rainforest tree, Fontainea rostrata

Assisted migration can aid in the conservation of narrowly endemic species affected by habitat loss, fragmentation and climate change. Here, we employ a multidisciplinary approach by examining the population genetic structure of a threatened, dioecious rainforest tree of the subtropical notophyll vine forests of eastern Australia, Fontainea rostrata, and its potential requirements for population enhancement and translocation to withstand the effects of anthropogenic fragmentation and climate change. We used microsatellite markers to gain an understanding of the way genetic diversity is partitioned within and among the nine extant populations of F. rostrata identified in this study. We combined the results with species distribution modelling to identify populations vulnerable to possible future range shifts based on climate change projections. We found regional differences between the species’ main distribution in the south and a disjunct northern population cluster (F_RT = 0.074, F_SR = 0.088, F_ST = 0.155), in mean allelic richness (A_R = 2.77 vs 2.33, p < 0.05), expected heterozygosity (H_E = 0.376 vs 0.328), and inbreeding (F = 0.116 vs 0.219). Species distribution models predicted that while southern populations of F. rostrata are likely to persist for the next 50 years under the RCP6.0 climate change scenario, with potential for a small-scale expansion to the south-east, the more highly inbred and less genetically diverse northern populations will come under increasing pressure to expand southwards as habitat suitability declines. Given the species’ genetic structure and with the aim to enhance genetic diversity and maximise the likelihood of reproductive success, we recommend that plant reintroductions to supplement existing populations should be prioritised over translocation of the species to new sites. However, future conservation efforts should be directed at translocation to establish new sites to increase population connectivity, focussing particularly on habitat areas identified as persisting under conditions of climate change.

Microsatellite resources of Eucalyptus: current status and future perspectives

Eucalyptus is the premier paper pulp, short rotation plantation species grown all over the world. Genetic improvement programs integrating molecular marker tools are in progress in many parts of the globe to increase the productivity. Whole genome sequence and expressed sequence tags (ESTs) of the eucalypts paved way for introduction of molecular genetics and breeding in this genus. Different molecular characterization approaches have been used simultaneously in eucalypts, however, microsatellites or simple sequence repeats (SSRs) with their prolific characteristics could occupy a special niche in Eucalyptus genetic improvement. Further, highly informative SSRs were used for the clonal identity, genetic fidelity and in certification of breeder's rights. Eucalyptus genetic linkage maps generated with microsatellite loci were used successfully to identify quantitative trait loci (QTLs) for various economically important traits. Progressively more numbers of microsatellites are being linked to genes associated with adaptive and functional variations, therefore making their utility broader in genetic applications. Availability of common SSR markers across the species provides an opportunity to validate the expression of QTLs across variable genetic backgrounds and accurately compare the position of QTLs in other species. Recent evidences suggest that the presence of SSRs in micro RNAs of plant species play a role in the quantitative trait expression. Similar studies in eucalypts may provide new insights into the genetic architecture of transcript-level variations and post transcriptional gene regulation. This review on eucalypts microsatellites, highlights the availability and characteristics of genomic and eSSRs and their potential in genetic analysis of natural and breeding populations and also discusses the future prospects in population genetics and marker assisted selection.

Heterosis may result in selection favouring the products of long-distance pollen dispersal in Eucalyptus

Using native trees from near the northern and southern extremities of the relatively continuous eastern distribution of Eucalyptus globulus in Tasmania, we compared the progenies derived from natural open-pollination (OP) with those generated from within-region and long-distance outcrossing. Controlled outcrossing amongst eight parents - with four parents from each of the northern and southern regions - was undertaken using a diallel mating scheme. The progeny were planted in two field trials located within the species native range in southern Tasmania, and their survival and diameter growth were monitored over a 13-year-period. The survival and growth performances of all controlled cross types exceeded those of the OP progenies, consistent with inbreeding depression due to a combination of selfing and bi-parental inbreeding. The poorer survival of the northern regional (♀_N♂_N) outcrosses compared with the local southern regional outcrosses (♀_S♂_S) indicated differential selection against the former. Despite this mal-adaptation of the non-local ♀_N♂_N crosses at both southern sites, the survival of the inter-regional hybrids (♀_N♂_S and ♀_S♂_N) was never significantly different from that of the local ♀_S♂_S crosses. Significant site-dependent heterosis was detected for the growth of the surviving long-distance hybrids. This was expressed as mid-parent heterosis, particularly at the more northern planting site. Heterosis increased with age, while the difference between the regional ♀_N♂_N and ♀_S♂_S crosses remained insignificant at any age at either site. Nevertheless, the results for growth suggest that the fitness of individuals derived from long-distance crossing may be better at the more northern of the planting sites. Our results demonstrate the potential for early-age assessments of pollen dispersal to underestimate realised gene flow, with local inbreeding under natural open-pollination resulting in selection favouring the products of longer-distance pollinations. Indeed, heterosis derived from long-distance pollinations may be sufficient to counter local mal-adaptation, at least in the first generation.

Genetic control of interactions among individuals: contrasting outcomes of indirect genetic effects arising from neighbour disease infection and competition in a forest tree

Indirect genetic effects (IGEs) are heritable effects of individuals on trait values of their conspecifics. IGEs may substantially affect response to selection, but empirical studies on IGEs are sparse and their magnitude and correlation with direct genetic effects are largely unknown in plants. Here we used linear mixed models to estimate genetic (co)variances attributable to direct and indirect effects for growth and foliar disease damage in a large pedigreed population of Eucalyptus globulus. We found significant IGEs for growth and disease damage, which increased with age for growth. The correlation between direct and indirect genetic effects was highly negative for growth, but highly positive for disease damage, consistent with neighbour competition and infection, respectively. IGEs increased heritable variation by 71% for disease damage, but reduced heritable variation by 85% for growth, leaving nonsignificant heritable variation for later age growth. Thus, IGEs are likely to prevent response to selection in growth, despite a considerable ordinary heritability. IGEs change our perspective on the genetic architecture and potential response to selection. Depending on the correlation between direct and indirect genetic effects, IGEs may enhance or diminish the response to natural or artificial selection compared with that predicted from ordinary heritability.

A latitudinal cline in disease resistance of a host tree

The possible drivers and implications of an observed latitudinal cline in disease resistance of a host tree were examined. Mycosphaerella leaf disease (MLD) damage, caused by Teratosphaeria species, was assessed in five Eucalyptus globulus (Tasmanian blue gum) common garden trials containing open-pollinated progeny from 13 native-forest populations. Significant population and family within population variation in MLD resistance was detected, which was relatively stable across different combinations of trial sites, ages, seasons and epidemics. A distinct genetic-based latitudinal cline in MLD damage among host populations was evident. Two lines of evidence argue that the observed genetic-based latitudinal trend was the result of direct pathogen-imposed selection for MLD resistance. First, MLD damage was positively associated with temperature and negatively associated with a prediction of disease risk in the native environment of these populations; and, second, the quantitative inbreeding coefficient (QST) significantly exceeded neutral marker FST at the trial that exhibited the greatest MLD damage, suggesting that diversifying selection contributed to differentiation in MLD resistance among populations. This study highlights the potential for spatial variation in pathogen risk to drive adaptive differentiation across the geographic range of a foundation host tree species.

Multiple evolutionary processes drive the patterns of genetic differentiation in a forest tree species complex

Forest trees frequently form species complexes, complicating taxonomic classification and gene pool management. This is certainly the case in Eucalyptus, and well exemplified by the Eucalyptus globulus complex. This ecologically and economically significant complex comprises four taxa (sspp. bicostata, globulus, maidenii, pseudoglobulus) that are geographically and morphologically distinct, but linked by extensive "intergrade" populations. To resolve their genetic affinities, nine microsatellites were used to genotype 1200 trees from throughout the natural range of the complex in Australia, representing 33 morphological core and intergrade populations. There was significant spatial genetic structure (F(ST) = 0.10), but variation was continuous. High genetic diversity in southern ssp. maidenii indicates that this region is the center of origin. Genetic diversity decreases and population differentiation increases with distance from this area, suggesting that drift is a major evolutionary process. Many of the intergrade populations, along with other populations morphologically classified as ssp. pseudoglobulus or ssp. globulus, belong to a "cryptic genetic entity" that is genetically and geographically intermediate between core ssp. bicostata, ssp. maidenii, and ssp. globulus. Geography, rather than morphology, therefore, is the best predictor of overall genetic affinities within the complex and should be used to classify germplasm into management units for conservation and breeding purposes.

Genetic evidence for a Janzen-Connell recruitment pattern in reproductive offspring of Pinus halepensis trees

Effective seed dispersal, combining both dispersal and postdispersal (establishment) processes, determines population dynamics and colonization ability in plants. According to the Janzen-Connell (JC) model, high mortality near the mother plant shifts the offspring establishment distribution farther away from the mother plant relative to the seed dispersal distribution. Yet, extending this prediction to the distribution of mature (reproductive) offspring remains a challenge for long-living plants. To address this challenge, we selected an isolated natural Aleppo pine (Pinus halepensis) population in Mt. Pithulim (Israel), which expanded from five ancestor trees in the beginning of the 20th century into ∼2000 trees today. Using nine microsatellite markers, we assigned parents to trees established during the early stages of population expansion. To elucidate the effect of the distance from the mother plant on postdispersal survival, we compared the effective seed dispersal kernel, based on the distribution of mother-offspring distances, with the seed dispersal kernel, based on simulations of a mechanistic wind dispersal model. We found that the mode of the effective dispersal kernel is shifted farther away than the mode of the seed dispersal kernel, reflecting increased survival with increasing distance from the mother plant. The parentage analysis demonstrated a highly skewed reproductive success and a strong directionality in effective dispersal corresponding to the wind regime. We thus provide compelling evidence that JC effects act also on offspring that become reproductive and persist as adults for many decades, a key requirement in assessing the role of postdispersal processes in shaping population and community dynamics.

Population genetic analysis and phylogeny reconstruction in Eucalyptus (Myrtaceae) using high-throughput, genome-wide genotyping

A set of over 8000 Diversity Arrays Technology (DArT) markers was tested for its utility in high-resolution population and phylogenetic studies across a range of Eucalyptus taxa. Small-scale population studies of Eucalyptus camaldulensis, Eucalyptus cladocalyx, Eucalyptus globulus, Eucalyptus grandis, Eucalyptus nitens, Eucalyptus pilularis and Eucalyptus urophylla demonstrated the potential of genome-wide genotyping with DArT markers to differentiate species, to identify interspecific hybrids and to resolve biogeographic disjunctions within species. The population genetic studies resolved geographically partitioned clusters in E. camaldulensis, E. cladocalyx, E. globulus and E. urophylla that were congruent with previous molecular studies. A phylogenetic study of 94 eucalypt species provided results that were largely congruent with traditional taxonomy and ITS-based phylogenies, but provided more resolution within major clades than had been obtained previously. Ascertainment bias (the bias introduced in a phylogeny from using markers developed in a small sample of the taxa that are being studied) was not detected. DArT offers an unprecedented level of resolution for population genetic, phylogenetic and evolutionary studies across the full range of Eucalyptus species.

Fine-scale spatial genetic structure and gene dispersal in Silene latifolia

Plants are sessile organisms, often characterized by limited dispersal. Seeds and pollen are the critical stages for gene flow. Here we investigate spatial genetic structure, gene dispersal and the relative contribution of pollen vs seed in the movement of genes in a stable metapopulation of the white campion Silene latifolia within its native range. This short-lived perennial plant is dioecious, has gravity-dispersed seeds and moth-mediated pollination. Direct measures of pollen dispersal suggested that large populations receive more pollen than small isolated populations and that most gene flow occurs within tens of meters. However, these studies were performed in the newly colonized range (North America) where the specialist pollinator is absent. In the native range (Europe), gene dispersal could fall on a different spatial scale. We genotyped 258 individuals from large and small (15) subpopulations along a 60 km, elongated metapopulation in Europe using six highly variable microsatellite markers, two X-linked and four autosomal. We found substantial genetic differentiation among subpopulations (global F(ST)=0.11) and a general pattern of isolation by distance over the whole sampled area. Spatial autocorrelation revealed high relatedness among neighboring individuals over hundreds of meters. Estimates of gene dispersal revealed gene flow at the scale of tens of meters (5-30 m), similar to the newly colonized range. Contrary to expectations, estimates of dispersal based on X and autosomal markers showed very similar ranges, suggesting similar levels of pollen and seed dispersal. This may be explained by stochastic events of extensive seed dispersal in this area and limited pollen dispersal.

Comparison of contemporary mating patterns in continuous and fragmented Eucalyptus globulus native forests

While habitat fragmentation is a central issue in forest conservation studies in the face of broad-scale anthropogenic changes to the environment, its effects on contemporary mating patterns remain controversial. This is partly because of the inherent variation in mating patterns which may exist within species and the fact that few studies have replication at the landscape level. To study the effect of forest fragmentation on contemporary mating patterns, including effective pollen dispersal, we compared four native populations of the Australian forest tree, Eucalyptus globulus. We used six microsatellite markers to genotype 1289 open-pollinated offspring from paired fragmented and continuous populations on the island of Tasmania and in Victoria on mainland Australia. The mating patterns in the two continuous populations were similar, despite large differences in population density. In contrast, the two fragmented populations were variable and idiosyncratic in their mating patterns, particularly in their pollen dispersal kernels. The continuous populations showed relatively high outcrossing rates (86-89%) and low correlated paternity (0.03-0.06) compared with the fragmented populations (65-79% and 0.12-0.20 respectively). A greater proportion of trees contributed to reproduction in the fragmented (de/d>or= 0.5) compared with the continuous populations (de/d = 0.03-0.04). Despite significant inbreeding in the offspring of the fragmented populations, there was little evidence of loss of genetic diversity. It is argued that enhanced medium- and long-distance dispersal in fragmented landscapes may act to partly buffer the remnant populations from the negative effects of inbreeding and drift.

Population fragmentation leads to isolation by distance but not genetic impoverishment in the philopatric Lesser Kestrel: a comparison with the widespread and sympatric Eurasian Kestrel

Population fragmentation is a widespread phenomenon usually associated with human activity. As a result of habitat transformation, the philopatric and steppe-specialist Lesser Kestrel Falco naumanni underwent a severe population decline during the last century that increased population fragmentation throughout its breeding range. In contrast, the ubiquitous Eurasian Kestrel Falco tinnunculus did not suffer such adverse effects, its breeding range still remaining rather continuous. Using microsatellites, we tested the effects of population fragmentation on large-scale spatial patterns of genetic differentiation and diversity by comparing these two sympatric and phylogenetically related species. Our results suggest that habitat fragmentation has increased genetic differentiation between Lesser Kestrel populations, following an isolation-by-distance pattern, while the population of Eurasian Kestrels is panmictic. Contrary to expectations, we did not detect significant evidence of reduced genetic variation or increased inbreeding in Lesser Kestrels. Although this study reports genetic differentiation in a species that has potential for long-distance dispersal but philopatry-limited gene flow, large enough effective population sizes and migration may have been sufficient to mitigate genetic depauperation. A serious reduction of genetic diversity in Lesser Kestrels would, therefore, only be expected after severe population bottlenecks following extreme geographic isolation.

Parallel evolution of dwarf ecotypes in the forest tree Eucalyptus globulus

Three small populations of a dwarf ecotype of the forest tree Eucalyptus globulus are found on exposed granite headlands in south-eastern Australia. These populations are separated by at least 100 km. Here, we used 12 nuclear microsatellites and a chloroplast DNA marker to investigate the genetic affinities of the dwarf populations to one another and to their nearest populations of tall E. globulus. Cape Tourville was studied in greater detail to assess the processes enabling the maintenance of distinct ecotypes in close geographical proximity. The three dwarf populations were not related to one another and were more closely related to adjacent tall trees than to one another. At Cape Tourville the dwarf and tall ecotypes were significantly differentiated in microsatellites and in chloroplast DNA. The dwarf and tall populations differed in flowering time and no evidence of pollen dispersal from the more extensive tall to the dwarf population was found. The three dwarf populations have evolved in parallel from the local tall ecotypes. This study shows that small marginal populations of eucalypts are capable of developing reproductive isolation from nearby larger populations through differences in flowering time and/or minor spatial separation, making parapatric speciation possible.