Molecular insights into the dynamics of species invasion by hybridisation in Tasmanian eucalypts.

In plants where seed dispersal is limited compared with pollen dispersal, hybridisation may enhance gene exchange and species dispersal. We provide genetic evidence of hybridisation contributing to the expansion of the rare Eucalyptus risdonii into the range of the widespread Eucalyptus amygdalina. These closely related tree species are morphologically distinct, and observations suggest that natural hybrids occur along their distribution boundaries and as isolated trees or in small patches within the range of E. amygdalina. Hybrid phenotypes occur outside the range of normal dispersal for E. risdonii seed, yet in some hybrid patches small individuals resembling E. risdonii occur and are hypothesised to be a result of backcrossing. Using 3362 genome-wide SNPs assessed from 97 individuals of E. risdonii and E. amygdalina and 171 hybrid trees, we show that (i) isolated hybrids match the genotypes expected of F1 /F2 hybrids, (ii) there is a continuum in the genetic composition among the isolated hybrid patches from patches dominated by F1 /F2 -like genotypes to those dominated by E. risdonii-backcross genotypes, and (iii) the E. risdonii-like phenotypes in the isolated hybrid patches are most-closely related to proximal larger hybrids. These results suggest that the E. risdonii phenotype has been resurrected in isolated hybrid patches established from pollen dispersal, providing the first steps in its invasion of suitable habitat by long-distance pollen dispersal and complete introgressive displacement of E. amygdalina. Such expansion accords with the population demographics, common garden performance data, and climate modelling which favours E. risdonii and highlights a role of interspecific hybridisation in climate change adaptation and species expansion.

Leaf Economic and Hydraulic Traits Signal Disparate Climate Adaptation Patterns in Two Co-Occurring Woodland Eucalypts.

With climate change impacting trees worldwide, enhancing adaptation capacity has become an important goal of provenance translocation strategies for forestry, ecological renovation, and biodiversity conservation. Given that not every species can be studied in detail, it is important to understand the extent to which climate adaptation patterns can be generalised across species, in terms of the selective agents and traits involved. We here compare patterns of genetic-based population (co)variation in leaf economic and hydraulic traits, climate-trait associations, and genomic differentiation of two widespread tree species (Eucalyptus pauciflora and E. ovata). We studied 2-year-old trees growing in a common-garden trial established with progeny from populations of both species, pair-sampled from 22 localities across their overlapping native distribution in Tasmania, Australia. Despite originating from the same climatic gradients, the species differed in their levels of population variance and trait covariance, patterns of population variation within each species were uncorrelated, and the species had different climate-trait associations. Further, the pattern of genomic differentiation among populations was uncorrelated between species, and population differentiation in leaf traits was mostly uncorrelated with genomic differentiation. We discuss hypotheses to explain this decoupling of patterns and propose that the choice of seed provenances for climate-based plantings needs to account for multiple dimensions of climate change unless species-specific information is available.

Genomic divergence in sympatry indicates strong reproductive barriers and cryptic species within Eucalyptus salubris.

Genetic studies are increasingly detecting cryptic taxa that likely represent a significant component of global biodiversity. However, cryptic taxa are often criticized because they are typically detected serendipitously and may not receive the follow-up study required to verify their geographic or evolutionary limits. Here, we follow-up a study of Eucalyptus salubris that unexpectedly detected two divergent lineages but was not sampled sufficiently to make clear interpretations. We undertook comprehensive sampling for an independent genomic analysis (3,605 SNPs) to investigate whether the two purported lineages remain discrete genetic entities or if they intergrade throughout the species' range. We also assessed morphological and ecological traits, and sequenced chloroplast DNA. SNP results showed strong genome-wide divergence (F ST = 0.252) between two discrete lineages: one dominated the north and one the southern regions of the species' range. Within lineages, gene flow was high, with low differentiation (mean F ST = 0.056) spanning hundreds of kilometers. In the central region, the lineages were interspersed but maintained their genomic distinctiveness: an indirect demonstration of reproductive isolation. Populations of the southern lineage exhibited significantly lower specific leaf area and occurred on soils with lower phosphorus relative to the northern lineage. Finally, two major chloroplast haplotypes were associated with each lineage but were shared between lineages in the central distribution. Together, these results suggest that these lineages have non-contemporary origins and that ecotypic adaptive processes strengthened their divergence more recently. We conclude that these lineages warrant taxonomic recognition as separate species and provide fascinating insight into eucalypt speciation.

Regarding the F-word: The effects of data filtering on inferred genotype-environment associations.

Genotype-environment association (GEA) methods have become part of the standard landscape genomics toolkit, yet, we know little about how to best filter genotype-by-sequencing data to provide robust inferences for environmental adaptation. In many cases, default filtering thresholds for minor allele frequency and missing data are applied regardless of sample size, having unknown impacts on the results, negatively affecting management strategies. Here, we investigate the effects of filtering on GEA results and the potential implications for assessment of adaptation to environment. We use empirical and simulated data sets derived from two widespread tree species to assess the effects of filtering on GEA outputs. Critically, we find that the level of filtering of missing data and minor allele frequency affect the identification of true positives. Even slight adjustments to these thresholds can change the rate of true positive detection. Using conservative thresholds for missing data and minor allele frequency substantially reduces the size of the data set, lessening the power to detect adaptive variants (i.e., simulated true positives) with strong and weak strengths of selection. Regardless, strength of selection was a good predictor for GEA detection, but even some SNPs under strong selection went undetected. False positive rates varied depending on the species and GEA method, and filtering significantly impacted the predictions of adaptive capacity in downstream analyses. We make several recommendations regarding filtering for GEA methods. Ultimately, there is no filtering panacea, but some choices are better than others, depending on the study system, availability of genomic resources, and desired objectives.

The potential of genomics for restoring ecosystems and biodiversity.

Billions of hectares of natural ecosystems have been degraded through human actions. The global community has agreed on targets to halt and reverse these declines, and the restoration sector faces the important but arduous task of implementing programmes to meet these objectives. Existing and emerging genomics tools offer the potential to improve the odds of achieving these targets. These tools include population genomics that can improve seed sourcing, meta-omics that can improve assessment and monitoring of restoration outcomes, and genome editing that can generate novel genotypes for restoring challenging environments. We identify barriers to adopting these tools in a restoration context and emphasize that regulatory and ethical frameworks are required to guide their use.

Microsatellite analysis of population structure in Eucalyptus globulus.

Eucalyptus globulus subsp. globulus Labill. (Tasmanian Blue Gum), native to southeast Australia, is a benchmark species for the pulp and paper industry. We genotyped 397 trees from 16 populations of E. globulus representing the native diversity in Australia using 24 microsatellite loci. Eight genetically distinct groups were detected, consistent with genetic groupings detected in previous quantitative and molecular studies. A sample of 29 Portuguese individuals was added to help clarify the origin of the Portuguese landrace. The results suggest a southern and eastern Tasmania origin for the Portuguese landrace. This genetic framework will enable researchers to investigate the provenance of individuals of unknown pedigree and assess the levels of representation of E. globulus natural variation in the Portuguese landrace.

Bioclimatic transect networks: Powerful observatories of ecological change.

Transects that traverse substantial climate gradients are important tools for climate change research and allow questions on the extent to which phenotypic variation associates with climate, the link between climate and species distributions, and variation in sensitivity to climate change among biomes to be addressed. However, the potential limitations of individual transect studies have recently been highlighted. Here, we argue that replicating and networking transects, along with the introduction of experimental treatments, addresses these concerns. Transect networks provide cost-effective and robust insights into ecological and evolutionary adaptation and improve forecasting of ecosystem change. We draw on the experience and research facilitated by the Australian Transect Network to demonstrate our case, with examples, to clarify how population- and community-level studies can be integrated with observations from multiple transects, manipulative experiments, genomics, and ecological modeling to gain novel insights into how species and systems respond to climate change. This integration can provide a spatiotemporal understanding of past and future climate-induced changes, which will inform effective management actions for promoting biodiversity resilience.

Genomic Scans across Three Eucalypts Suggest that Adaptation to Aridity is a Genome-Wide Phenomenon.

Widespread species spanning strong environmental (e.g., climatic) gradients frequently display morphological and physiological adaptations to local conditions. Some adaptations are common to different species that occupy similar environments. However, the genomic architecture underlying such convergent traits may not be the same between species. Using genomic data from previous studies of three widespread eucalypt species that grow along rainfall gradients in southern Australia, our probabilistic approach provides evidence that adaptation to aridity is a genome-wide phenomenon, likely to involve multiple and diverse genes, gene families and regulatory regions that affect a multitude of complex genetic and biochemical processes.

High density, genome-wide markers and intra-specific replication yield an unprecedented phylogenetic reconstruction of a globally significant, speciose lineage of Eucalyptus.

We used genome-wide markers and an unprecedented scale of sampling to construct a phylogeny for a globally significant Eucalyptus lineage that has been impacted by hybridisation, recent radiation and morphological convergence. Our approach, using 3109 DArT markers distributed throughout the genome and 540 samples covering 185 terminal taxa in sections Maidenaria, Exsertaria, Latoangulatae and related smaller sections, with multiple geographically widespread samples per terminal taxon, produced a phylogeny that largely matched the morphological treatment of sections, though sections Exsertaria and Latoangulatae were polyphyletic. At lower levels there were numerous inconsistencies between the morphological treatment and the molecular phylogeny, and taxa within the three main sections were generally not monophyletic at the series (at least 62% polyphyly) or species (at least 52% polyphyly) level. Some of the discrepancies appear to be the result of morphological convergence or misclassifications, and we propose some taxonomic reassessments to address this. However, many inconsistencies appear to be the products of incomplete speciation and/or hybridisation. Our analysis represents a significant advance on previous phylogenies of these important eucalypt sections (which have mainly used single samples to represent each species), thus providing a robust phylogenetic framework for evolutionary and ecological studies.

Patterns of Reproductive Isolation in Eucalyptus-A Phylogenetic Perspective.

We assess phylogenetic patterns of hybridization in the speciose, ecologically and economically important genus Eucalyptus, in order to better understand the evolution of reproductive isolation. Eucalyptus globulus pollen was applied to 99 eucalypt species, mainly from the large commercially important subgenus, Symphyomyrtus. In the 64 species that produce seeds, hybrid compatibility was assessed at two stages, hybrid-production (at approximately 1 month) and hybrid-survival (at 9 months), and compared with phylogenies based on 8,350 genome-wide DArT (diversity arrays technology) markers. Model fitting was used to assess the relationship between compatibility and genetic distance, and whether or not the strength of incompatibility "snowballs" with divergence. There was a decline in compatibility with increasing genetic distance between species. Hybridization was common within two closely related clades (one including E. globulus), but rare between E. globulus and species in two phylogenetically distant clades. Of three alternative models tested (linear, slowdown, and snowball), we found consistent support for a snowball model, indicating that the strength of incompatibility accelerates relative to genetic distance. Although we can only speculate about the genetic basis of this pattern, it is consistent with a Dobzhansky-Muller-model prediction that incompatibilities should snowball with divergence due to negative epistasis. Different rates of compatibility decline in the hybrid-production and hybrid-survival measures suggest that early-acting postmating barriers developed first and are stronger than later-acting barriers. We estimated that complete reproductive isolation can take up to 21-31 My in Eucalyptus. Practical implications for hybrid eucalypt breeding and genetic risk assessment in Australia are discussed.

Genome-wide scans detect adaptation to aridity in a widespread forest tree species.

Patterns of adaptive variation within plant species are best studied through common garden experiments, but these are costly and time-consuming, especially for trees that have long generation times. We explored whether genome-wide scanning technology combined with outlier marker detection could be used to detect adaptation to climate and provide an alternative to common garden experiments. As a case study, we sampled nine provenances of the widespread forest tree species, Eucalyptus tricarpa, across an aridity gradient in southeastern Australia. Using a Bayesian analysis, we identified a suite of 94 putatively adaptive (outlying) sequence-tagged markers across the genome. Population-level allele frequencies of these outlier markers were strongly correlated with temperature and moisture availability at the site of origin, and with population differences in functional traits measured in two common gardens. Using the output from a canonical analysis of principal coordinates, we devised a metric that provides a holistic measure of genomic adaptation to aridity that could be used to guide assisted migration or genetic augmentation.

Plasticity of functional traits varies clinally along a rainfall gradient in Eucalyptus tricarpa.

Widespread species often occur across a range of climatic conditions, through a combination of local genetic adaptations and phenotypic plasticity. Species with greater phenotypic plasticity are likely to be better positioned to cope with rapid anthropogenic climate changes, while those displaying strong local adaptations might benefit from translocations to assist the movement of adaptive genes as the climate changes. Eucalyptus tricarpa occurs across a climatic gradient in south-eastern Australia, a region of increasing aridity, and we hypothesized that this species would display local adaptation to climate. We measured morphological and physiological traits reflecting climate responses in nine provenances from sites of 460 to 1040 mm annual rainfall, in their natural habitat and in common gardens near each end of the gradient. Local adaptation was evident in functional traits and differential growth rates in the common gardens. Some traits displayed complex combinations of plasticity and genetic divergence among provenances, including clinal variation in plasticity itself. Provenances from drier locations were more plastic in leaf thickness, whereas leaf size was more plastic in provenances from higher rainfall locations. Leaf density and stomatal physiology (as indicated by δ(13)C and δ(18)O) were highly and uniformly plastic. In addition to variation in mean trait values, genetic variation in trait plasticity may play a role in climate adaptation.

Novel distances for dollo data.

We investigate distances on binary (presence/absence) data in the context of a Dollo process, where a trait can only arise once on a phylogenetic tree but may be lost many times. We introduce a novel distance, the Additive Dollo Distance (ADD), that applies to data generated under a Dollo model and show that it has some useful theoretical properties including an intriguing link to the LogDet/paralinear distance. Simulations of Dollo data are used to compare a number of binary distances including ADD, LogDet, a restriction-site-based distance, and some simple, but to our knowledge previously unstudied, variations on common binary distances. The simulations suggest that ADD outperforms other distances on Dollo data. Interestingly, we found that the LogDet distance performs poorly in the context of a Dollo process; this may have implications for its use in connection with conditioned genome reconstruction. We apply the ADD to two Diversity Arrays Technology data sets, one that broadly covers Eucalyptus species and one that focuses on the Eucalyptus series Adnataria. We also reanalyze gene family presence/absence data from bacterial genomes obtained from the COG database and compare the results with previous phylogenies estimated using the conditioned genome reconstruction approach. The results for these case studies are largely congruent with previous studies, in some cases giving more phylogenetic resolution.

What does population structure analysis reveal about the Pterostylis longifolia complex (Orchidaceae)?

Morphologically similar groups of species are common and pose significant challenges for taxonomists. Differences in approaches to classifying unique species can result in some species being overlooked, whereas others are wrongly conserved. The genetic diversity and population structure of the Pterostylis longifolia complex (Orchidaceae) in Tasmania was investigated to determine if four species, and potential hybrids, could be distinguished through genomic AFLP and chloroplast restriction-fragment-length polymorphism (RFLP) markers. Analysis of molecular variance (AMOVA) results indicated that little genetic variation was present among taxa, whereas PCoA analyses revealed genetic variation at a regional scale irrespective of taxa. Population genetic structure analyses identified three clusters that correspond to regional genetic and single taxon-specific phenotypic variation. The results from this study suggest that "longifolia" species have persisted throughout the last glacial maximum in Tasmania and that the complex may be best treated as a single taxon with several morphotypes. These results could have serious evolutionary and conservation implications as taxonomic changes could result in the instatement of a single, widespread taxon in which rarer morphotypes are not protected.

Genomic characterization of DArT markers based on high-density linkage analysis and physical mapping to the Eucalyptus genome.

Diversity Arrays Technology (DArT) provides a robust, high throughput, cost-effective method to query thousands of sequence polymorphisms in a single assay. Despite the extensive use of this genotyping platform for numerous plant species, little is known regarding the sequence attributes and genome-wide distribution of DArT markers. We investigated the genomic properties of the 7,680 DArT marker probes of a Eucalyptus array, by sequencing them, constructing a high density linkage map and carrying out detailed physical mapping analyses to the Eucalyptus grandis reference genome. A consensus linkage map with 2,274 DArT markers anchored to 210 microsatellites and a framework map, with improved support for ordering, displayed extensive collinearity with the genome sequence. Only 1.4 Mbp of the 75 Mbp of still unplaced scaffold sequence was captured by 45 linkage mapped but physically unaligned markers to the 11 main Eucalyptus pseudochromosomes, providing compelling evidence for the quality and completeness of the current Eucalyptus genome assembly. A highly significant correspondence was found between the locations of DArT markers and predicted gene models, while most of the 89 DArT probes unaligned to the genome correspond to sequences likely absent in E. grandis, consistent with the pan-genomic feature of this multi-Eucalyptus species DArT array. These comprehensive linkage-to-physical mapping analyses provide novel data regarding the genomic attributes of DArT markers in plant genomes in general and for Eucalyptus in particular. DArT markers preferentially target the gene space and display a largely homogeneous distribution across the genome, thereby providing superb coverage for mapping and genome-wide applications in breeding and diversity studies. Data reported on these ubiquitous properties of DArT markers will be particularly valuable to researchers working on less-studied crop species who already count on DArT genotyping arrays but for which no reference genome is yet available to allow such detailed characterization.

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.

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.

A high-density Diversity Arrays Technology (DArT) microarray for genome-wide genotyping in Eucalyptus.

BACKGROUND: A number of molecular marker technologies have allowed important advances in the understanding of the genetics and evolution of Eucalyptus, a genus that includes over 700 species, some of which are used worldwide in plantation forestry. Nevertheless, the average marker density achieved with current technologies remains at the level of a few hundred markers per population. Furthermore, the transferability of markers produced with most existing technology across species and pedigrees is usually very limited. High throughput, combined with wide genome coverage and high transferability are necessary to increase the resolution, speed and utility of molecular marker technology in eucalypts. We report the development of a high-density DArT genome profiling resource and demonstrate its potential for genome-wide diversity analysis and linkage mapping in several species of Eucalyptus. FINDINGS: After testing several genome complexity reduction methods we identified the PstI/TaqI method as the most effective for Eucalyptus and developed 18 genomic libraries from PstI/TaqI representations of 64 different Eucalyptus species. A total of 23,808 cloned DNA fragments were screened and 13,300 (56%) were found to be polymorphic among 284 individuals. After a redundancy analysis, 6,528 markers were selected for the operational array and these were supplemented with 1,152 additional clones taken from a library made from the E. grandis tree whose genome has been sequenced. Performance validation for diversity studies revealed 4,752 polymorphic markers among 174 individuals. Additionally, 5,013 markers showed segregation when screened using six inter-specific mapping pedigrees, with an average of 2,211 polymorphic markers per pedigree and a minimum of 859 polymorphic markers that were shared between any two pedigrees. CONCLUSIONS: This operational DArT array will deliver 1,000-2,000 polymorphic markers for linkage mapping in most eucalypt pedigrees and thus provide high genome coverage. This array will also provide a high-throughput platform for population genetics and phylogenetics in Eucalyptus. The transferability of DArT across species and pedigrees is particularly valuable for a large genus such as Eucalyptus and will facilitate the transfer of information between different studies. Furthermore, the DArT marker array will provide a high-resolution link between phenotypes in populations and the Eucalyptus reference genome, which will soon be completed.

A geographic mosaic of genetic variation within a foundation tree species and its community-level consequences.

Knowledge of the manner in which genetic variation within a tree species affects associated communities and ecosystem processes across its entire range is important for understanding how geographic mosaics of genetic interactions might develop and support different communities. While numerous studies have investigated the community and ecosystem consequences of genetic variation at the hybrid cross type or genotype level within a species, none has investigated the community-level effects of intraspecific genetic variation across the geographic range of a widespread species. This is the scale at which geographic mosaics of coevolution are hypothesized to exist. Studies at this level are particularly important for foundation tree species, which typically support numerous microbial, fungal, plant, and animal communities. We studied genetic variation across eight geographical races of the forest tree Eucalyptus globulus representing its natural distribution across southeastern Australia. The study was conducted in a 15-year-old common garden trial based on families derived from single-tree open-pollinated seed collections from the wild. Neutral molecular genetic variation within E. globulus was also assessed and compared with genetic divergence in the phenotypic and community traits. Three major findings emerged. First, we found significant genetically based, hierarchical variation in associated communities corresponding to geographical races of E. globulus and families within races. Second, divergence in foliar communities at the racial level was associated with genetically based divergence in specific leaf morphological and chemical traits that have known defensive functions. Third, significant positive correlations between canopy community dissimilarity and both neutral molecular genetic and leaf quantitative genetic dissimilarity at the race level supported a genetic similarity rule. Our results argue that genetic variation within foundation tree species has the potential to be a significant driver of the geographical mosaics of variation typical of forest communities, which could have important ecological and evolutionary implications.

Biodiversity consequences of genetic variation in bark characteristics within a foundation tree species.

The developing field of community genetics has the potential to broaden the contribution of genetics to conservation biology by demonstrating that genetic variation within foundation plant species can act to structure associated communities of microorganisms, invertebrates, and vertebrates. We assessed the biodiversity consequences of natural patterns of intraspecific genetic variation within the widely distributed Australian forest tree, Eucalyptus globulus. We assessed genetic variation among geographic races of E. globulus (i.e., provenances, seed zones) in the characteristics of tree-trunk bark in a 17-year-old common garden and the associated response of a dependent macroarthropod community. In total, 180 macroarthropod taxa were identified following a collection from 100 trees of five races. We found substantial genetically based variation within E. globulus in the quantity and type of decorticating bark. In the community of organisms associated with this bark, significant variation existed among trees of different races in composition, and there was a two-fold difference in species richness (7-14 species) and abundance (22-55 individuals) among races. This community variation was tightly linked with genetically based variation in bark, with 60% of variation in community composition driven by bark characteristics. No detectable correlation was found, however, with neutral molecular markers. These community-level effects of tree genetics are expected to extend to higher trophic levels because of the extensive use of tree trunks as foraging zones by birds and marsupials. Our results demonstrate the potential biodiversity benefits that may be gained through conservation of intraspecific genetic variation within broadly distributed foundation species. The opportunities for enhancing biodiversity values of forestry and restoration plantings are also highlighted because such planted forests are increasingly becoming the dominant forest type in many areas of the world.