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.

Genetic control of the operculum and capsule morphology of Eucalyptus globulus.

BACKGROUND AND AIMS: The petaline operculum that covers the inner whorls until anthesis and the woody capsule that develops after fertilization are reproductive structures of eucalypts that protect the flower and seeds. Although they are distinct organs, they both develop from flower buds and this common ontogeny suggests shared genetic control. In Eucalyptus globulus their morphology is variable and we aimed to identify the quantitative trait loci (QTL) underlying this variation and determine whether there is common genetic control of these ecologically and taxonomically important reproductive structures. METHODS: Samples of opercula and capsules were collected from 206 trees that belong to a large outcrossed F2E. globulus mapping population. The morphological variation in these structures was characterized by measuring six operculum and five capsule traits. QTL analysis was performed using these data and a linkage map consisting of 480 markers. KEY RESULTS: A total of 27 QTL were detected for operculum traits and 28 for capsule traits, with the logarithm of odds ranging from 2.8 to 11.8. There were many co-located QTL associated with operculum or capsule traits, generally reflecting allometric relationships. A key finding was five genomic regions where co-located QTL affected both operculum and capsule morphology, and the overall trend for these QTL was to affect elongation of both organs. Some of these QTL appear to have a significant effect on the phenotype, with the strongest QTL explaining 26.4 % of the variation in operculum shape and 16.4 % in capsule shape. Flower bud measurements suggest the expression of these QTL starts during bud development. Several candidate genes were found associated with the QTL and their putative function is discussed. CONCLUSIONS: Variation in both operculum and capsule traits in E. globulus is under strong genetic control. Our results suggest that these reproductive structures share a common genetic pathway during flower bud development.

Chemical Traits that Predict Susceptibility of Pinus radiata to Marsupial Bark Stripping.

Bark stripping by mammals is a major problem in managed conifer forests worldwide. In Australia, bark stripping in the exotic plantations of Pinus radiata is mainly caused by native marsupials and results in reduced survival, growth, and in extreme cases death of trees. Herbivory is influenced by a balance between primary metabolites that are sources of nutrition and secondary metabolites that act as defences. Identifying the compounds that influence herbivory may be a useful tool in the management of forest systems. This study aimed to detect and identify both constitutive and induced compounds that are associated with genetic differences in susceptibility of two-year-old P. radiata trees to bark stripping by marsupials. An untargeted profiling of 83 primary and secondary compounds of the needles and bark samples from 21 susceptible and 21 resistant families was undertaken. These were among the most and least damaged families, respectively, screened in a trial of 74 families that were exposed to natural field bark stripping by marsupials. Experimental plants were in the same field trial but protected from bark stripping and a subset were subjected to artificial bark stripping to examine induced and constitutive chemistry differences between resistant and susceptible families. Machine learning (random forest), partial least squares plus discriminant analysis (PLS-DA), and principal components analysis with discriminant analysis (PCA-DA), as well as univariate methods were used to identify the most important totals by compound group and individual compounds differentiating the resistant and susceptible families. In the bark, the constitutive amount of two sesquiterpenoids - bicyclogermacrene and an unknown sesquiterpenoid alcohol -were shown to be of higher levels in the resistant families, whereas the constitutive sugars, fructose, and glucose, as well individual phenolics, were higher in the more susceptible families. The chemistry of the needles was not useful in differentiating the resistant and susceptible families to marsupial bark stripping. After artificial bark stripping, the terpenes, sugars, and phenolics responded in both the resistant and susceptible families by increasing or reducing amounts, which leveled the differences in the amounts of the compounds between the different resistant and susceptible classes observed at the constitutive level. Overall, based on the families with extreme values for less and more susceptibility, differences in the amounts of secondary compounds were subtle and susceptibility due to sugars may outweigh defence as the cause of the genetic variation in bark stripping observed in this non-native tree herbivory system.

Additive genetic variation in Pinus radiata bark chemistry and the chemical traits associated with variation in mammalian bark stripping.

Secondary metabolites are suggested as a major mechanism explaining genetic variation in herbivory levels in Pinus radiata. The potential to incorporate these chemical traits into breeding/deployment programmes partly depends on the presence of additive genetic variation for the relevant chemical traits. In this study, near-infrared spectroscopy was used to quantify the constitutive and induced levels of 54 compounds in the bark of trees from 74 P. radiata full-sib families. The trees sampled for chemistry were protected from browsing and induced levels were obtained by subjecting half of the trees to artificial bark stripping. The treatment effect on bark chemistry was assessed along with narrow-sense heritability, the significance of non-additive genetic effects and the additive genetic correlations of compounds with bark stripping by mammalian herbivores that was observed in unprotected replicates of the field trial. The results indicated: (i) significant additive genetic variation, with low-moderate narrow-sense heritability estimates for most compounds; (ii) while significant induced effects were detected for some chemicals, no significant genetic variation in inducibility was detected; and (iii) sugars, fatty acids and a diterpenoid positively genetically correlated while a sesquiterpenoid negatively genetically correlated with bark stripping by the mammalian herbivore, the Bennett's wallaby (Macropus rufogriseus). At the onset of browsing, a trade-off with height was detected for selecting higher amounts of this sesquiterpenoid. However, overall, results showed potential to incorporate chemical traits into breeding/deployment programmes. The quantitative genetic analyses of the near infrared predicted chemical traits produced associations with mammalian bark stripping that mostly conform with those obtained using standard wet chemistry.

Consistent community genetic effects in the context of strong environmental and temporal variation in Eucalyptus.

Provenance translocations of tree species are promoted in forestry, conservation, and restoration in response to global climate change. While this option is driven by adaptive considerations, less is known of the effects translocations can have on dependent communities. We investigated the relative importance and consistency of extended genetic effects in Eucalyptus using two species-E. globulus and E. pauciflora. In E. globulus, the dependent arthropod and pathogen canopy communities were quantified based on the abundance of 49 symptoms from 722 progeny from 13 geographic sub-races across 2 common gardens. For E. pauciflora, 6 symptoms were quantified over 2 years from 238 progeny from 16 provenances across 2 common gardens. Genetic effects significantly influenced communities in both species. However, site and year effects outweighed genetic effects with site explaining approximately 3 times the variation in community traits in E. globulus and site and year explaining approximately 6 times the variation in E. pauciflora. While the genetic effect interaction terms were significant in some community traits, broad trends in community traits associated with variation in home-site latitude for E. globulus and home-site altitude for E. pauciflora were evident. These broad-scale trends were consistent with patterns of adaptive differentiation within each species, suggesting there may be extended consequences of local adaptation. While small in comparison to site and year, the consistency of genetic effects highlights the importance of provenance choice in tree species, such as Eucalyptus, as adaptive divergence among provenances may have significant long-term effects on biotic communities.

Independent genetic control of drought resistance, recovery, and growth of Eucalyptus globulus seedlings.

Drought is a major stress impacting forest ecosystems worldwide. We utilized quantitative trait loci (QTL) analysis to study the genetic basis of variation in (a) drought resistance and recovery and (b) candidate traits that may be associated with this variation in the forest tree Eucalyptus globulus. QTL analysis was performed using a large outcrossed F2 mapping population from which 300 trees were phenotyped based on the mean performance of their open-pollinated F3 progeny. Progenies were grown in a glasshouse in a randomized complete block design. A subset of seedlings was subjected to a drought treatment after which they were rewatered and scored for damage and growth postdrought. Nondroughted seedlings were assessed for growth traits as well as lignotuber size and resprouting following severe damage to the main stem. QTL were detected for most traits. Importantly, independent QTL were detected for (a) drought damage and plant size, (b) drought damage and growth recovery, and (c) lignotuber size and resprouting capacity. Such independence argues that trade-offs are unlikely to be a major limitation to the response to selection and at the early life history stage studied; there are opportunities to improve resilience to drought without adverse effects on productivity.

Comparison of host susceptibilities to native and exotic pathogens provides evidence for pathogen-imposed selection in forest trees.

The extent to which spatial structuring of host resistance in wild plant populations reflects direct pathogen-imposed selection is a subject of debate. To examine this issue, genetic susceptibilities to an exotic and a coevolved native fungal pathogen were compared using two Australian host tree species. Damage to common host germplasm of Corymbia citriodora ssp. variegata (CCV) and Eucalyptus globulus, caused by recently introduced (Austropuccinia psidii) and native (Quambalaria pitereka and Teratosphaeria sp.) pathogens was evaluated in common-garden experiments. There was significant additive genetic variation within host species for susceptibility to both the exotic and native pathogens. However, susceptibility to A. psidii was not genetically correlated with susceptibility to either native pathogen, providing support for pathogen-specific rather than general mechanisms of resistance. Population differentiation (QST ) for susceptibility to the native pathogens was greater than neutral expectations (molecular FST ), arguing for divergent selection. Coupled with lower native, but not exotic, pathogen susceptibility in host populations from areas climatically more prone to fungal proliferation, these findings suggest that pathogen-imposed selection has contributed directly to a geographic mosaic of host resistance to native pathogens.

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.

Annotation of the Corymbia terpene synthase gene family shows broad conservation but dynamic evolution of physical clusters relative to Eucalyptus.

Terpenes are economically and ecologically important phytochemicals. Their synthesis is controlled by the terpene synthase (TPS) gene family, which is highly diversified throughout the plant kingdom. The plant family Myrtaceae are characterised by especially high terpene concentrations, and considerable variation in terpene profiles. Many Myrtaceae are grown commercially for terpene products including the eucalypts Corymbia and Eucalyptus. Eucalyptus grandis has the largest TPS gene family of plants currently sequenced, which is largely conserved in the closely related E. globulus. However, the TPS gene family has been well studied only in these two eucalypt species. The recent assembly of two Corymbia citriodora subsp. variegata genomes presents an opportunity to examine the conservation of this important gene family across more divergent eucalypt lineages. Manual annotation of the TPS gene family in C. citriodora subsp. variegata revealed a similar overall number, and relative subfamily representation, to that previously reported in E. grandis and E. globulus. Many of the TPS genes were in physical clusters that varied considerably between Eucalyptus and Corymbia, with several instances of translocation, expansion/contraction and loss. Notably, there was greater conservation in the subfamilies involved in primary metabolism than those involved in secondary metabolism, likely reflecting different selective constraints. The variation in cluster size within subfamilies and the broad conservation between the eucalypts in the face of this variation are discussed, highlighting the potential contribution of selection, concerted evolution and stochastic processes. These findings provide the foundation to better understand terpene evolution within the ecologically and economically important Myrtaceae.

Evidence that divergent selection shapes a developmental cline in a forest tree species complex.

Background and Aims: Evolutionary change in developmental trajectories (heterochrony) is a major mechanism of adaptation in plants and animals. However, there are few detailed studies of the variation in the timing of developmental events among wild populations. We here aimed to identify the climatic drivers and measure selection shaping a genetic-based developmental cline among populations of an endemic tree species complex on the island of Tasmania. Methods: Seed lots from 38 native provenances encompassing the clinal transition from the heteroblastic Eucalyptus tenuiramis to the homoblastic Eucalyptus risdonii were grown in a common-garden field trial in southern Tasmania for 20 years. We used 27 climatic variables to model the provenance variation in vegetative juvenility as assessed at age 5 years. A phenotypic selection analysis was used to measure the fitness consequences of variation in vegetative juvenility based on its impact on the survival and reproductive capacity of survivors at age 20 years. Key Results: Significant provenance divergence in vegetative juvenility was shown to be associated with home-site aridity, with the retention of juvenile foliage increasing with increasing aridity. Our results indicated that climate change may lead to different directions of selection across the geographic range of the complex, and in our mesic field site demonstrated that total directional selection within phenotypically variable provenances was in favour of reduced vegetative juvenility. Conclusions: We provide evidence that heteroblasty is adaptive and argue that, in assessing the impacts of rapid global change, developmental plasticity and heterochrony are underappreciated processes which can contribute to populations of long-lived organisms, such as trees, persisting and ultimately adapting to environmental change.

Quantitative Trait Loci (QTLs) for Intumescence Severity in Eucalyptus globulus and Validation of QTL Detection Based on Phenotyping Using Open-Pollinated Families of a Mapping Population.

Intumescence is a nonpathogenic physiological disorder characterized by leaf blistering. This disorder can affect growth and development in glasshouses and growth chambers and may be confused with pathogenic diseases. We used quantitative trait loci (QTL) analysis to examine the genetic basis of variation in intumescence severity in Eucalyptus globulus, and test for colocation with previously detected QTLs for pathogen susceptibility. QTL analysis used the phenotype means of open-pollinated (OP) families of an outcrossed F2 mapping family (OP F3; n = 300) of E. globulus and the linkage map constructed in the F2. We validate this phenotyping approach for QTL analysis by assessing a trait previously used for QTL discovery in the F2 and showing the same major QTL was detected with the OP F3. For intumescence severity, five putative QTLs were detected across four linkage groups. Four of these did not colocate with previously reported QTLs for fungal pathogen susceptibility in Eucalyptus, suggesting the mechanisms underlying susceptibility to intumescence and to the two fungal pathogens are largely independent. This study demonstrates there is a genetic basis for variation in intumescence severity, reports the first QTL for intumescence severity in plants, and provides a robust framework for investigating the potential mechanisms involved.

Phylogeny is a powerful tool for predicting plant biomass responses to nitrogen enrichment.

Increasing rates of anthropogenic nitrogen (N) enrichment to soils often lead to the dominance of nitrophilic plant species and reduce plant diversity in natural ecosystems. Yet, we lack a framework to predict which species will be winners or losers in soil N enrichment scenarios, a framework that current literature suggests should integrate plant phylogeny, functional tradeoffs, and nutrient co-limitation. Using a controlled fertilization experiment, we quantified biomass responses to N enrichment for 23 forest tree species within the genus Eucalyptus that are native to Tasmania, Australia. Based on previous work with these species' responses to global change factors and theory on the evolution of plant resource-use strategies, we hypothesized that (1) growth responses to N enrichment are phylogenetically structured, (2) species with more resource-acquisitive functional traits have greater growth responses to N enrichment, and (3) phosphorus (P) limits growth responses to N enrichment differentially across species, wherein P enrichment increases growth responses to N enrichment more in some species than others. We built a hierarchical Bayesian model estimating effects of functional traits (specific leaf area, specific stem density, and specific root length) and P fertilization on species' biomass responses to N, which we then compared between lineages to determine whether phylogeny explains variation in responses to N. In concordance with literature on N limitation, a majority of species responded strongly and positively to N enrichment. Mean responses ranged three-fold, from 6.21 (E. pulchella) to 16.87 (E. delegatensis) percent increases in biomass per g N·m-2 ·yr-1 added. We identified a strong difference in responses to N between two phylogenetic lineages in the Eucalyptus subgenus Symphyomyrtus, suggesting that shared ancestry explains variation in N limitation. However, our model indicated that after controlling for phylogenetic non-independence, eucalypt responses to N were not associated with functional traits (although post-hoc analyses show a phylogenetic pattern in specific root length similar to that of responses to N), nor were responses differentially limited by P. Overall, our model results suggest that phylogeny is a powerful predictor of winners and losers in anthropogenic N enrichment scenarios in Tasmanian eucalypts, which may have implications for other species.

A water availability gradient reveals the deficit level required to affect traits in potted juvenile Eucalyptus globulus.

Background and aims: Drought leading to soil water deficit can have severe impacts on plants. Water deficit may lead to plant water stress and affect growth and chemical traits. Plant secondary metabolite (PSM) responses to water deficit vary between compounds and studies, with inconsistent reports of changes to PSM concentrations even within a single species. This disparity may result from experimental water deficit variation among studies, and so multiple water deficit treatments are used to fully assess PSM responses in a single species. Methods: Juvenile Eucalyptus globulus were grown for 8 weeks at one of ten water deficit levels based on evapotranspiration from control plants (100 %). Treatments ranged from 90 % of control evapotranspiration (mild water deficit) to 0 % of control evapotranspiration (severe water deficit) in 10 % steps. Plant biomass, foliar abscisic acid (ABA) levels, Ψ leaf , leaf C/N, selected terpenes and phenolics were quantified to assess responses to each level of water deficit relative to a control. Key Results: Withholding ≥30 % water resulted in higher foliar ABA levels and withholding ≥40 % water reduced leaf water content. Ψ leaf became more negative when ≥60 % water was withheld. Plant biomass was lower when ≥80 % water was withheld, and no water for 8 weeks (0 % water) resulted in plant death. The total oil concentration was lower and C/N was higher in dead and desiccated juvenile E. globulus leaves (0 % water). Concentrations of individual phenolic and terpene compounds, along with condensed tannin and total phenolic concentrations, remained stable regardless of water deficit or plant stress level. Conclusions: These juvenile E. globulus became stressed with a moderate reduction in available water, and yet the persistent concentrations of most PSMs in highly stressed or dead plants suggests no PSM re-metabolization and continued ecological roles of foliar PSMs during drought.

Genome-wide variation in recombination rate in Eucalyptus.

BACKGROUND: Meiotic recombination is a fundamental evolutionary process. It not only generates diversity, but influences the efficacy of natural selection and genome evolution. There can be significant heterogeneity in recombination rates within and between species, however this variation is not well understood outside of a few model taxa, particularly in forest trees. Eucalypts are forest trees of global economic importance, and dominate many Australian ecosystems. We studied recombination rate in Eucalyptus globulus using genetic linkage maps constructed in 10 unrelated individuals, and markers anchored to the Eucalyptus reference genome. This experimental design provided the replication to study whether recombination rate varied between individuals and chromosomes, and allowed us to study the genomic attributes and population genetic parameters correlated with this variation. RESULTS: Recombination rate varied significantly between individuals (range = 2.71 to 3.51 centimorgans/megabase [cM/Mb]), but was not significantly influenced by sex or cross type (F1 vs. F2). Significant differences in recombination rate between chromosomes were also evident (range = 1.98 to 3.81 cM/Mb), beyond those which were due to variation in chromosome size. Variation in chromosomal recombination rate was significantly correlated with gene density (r = 0.94), GC content (r = 0.90), and the number of tandem duplicated genes (r = -0.72) per chromosome. Notably, chromosome level recombination rate was also negatively correlated with the average genetic diversity across six species from an independent set of samples (r = -0.75). CONCLUSIONS: The correlations with genomic attributes are consistent with findings in other taxa, however, the direction of the correlation between diversity and recombination rate is opposite to that commonly observed. We argue this is likely to reflect the interaction of selection and specific genome architecture of Eucalyptus. Interestingly, the differences amongst chromosomes in recombination rates appear stable across Eucalyptus species. Together with the strong correlations between recombination rate and features of the Eucalyptus reference genome, we maintain these findings provide further evidence for a broad conservation of genome architecture across the globally significant lineages of Eucalyptus.

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.

Genetic control of cuticular wax compounds in Eucalyptus globulus.

Plant cuticular wax compounds perform functions that are essential for the survival of terrestrial plants. Despite their importance, the genetic control of these compounds is poorly understood outside of model taxa. Here we investigate the genetic basis of variation in cuticular compounds in Eucalyptus globulus using quantitative genetic and quantitative trait loci (QTL) analyses. Quantitative genetic analysis was conducted using 246 open-pollinated progeny from 13 native sub-races throughout the geographic range. QTL analysis was conducted using 112 clonally replicated progeny from an outcross F2 population. Nine compounds exhibited significant genetic variation among sub-races with three exhibiting signals of diversifying selection. Fifty-two QTL were found with co-location of QTL for related compounds commonly observed. Notable among these was the QTL for five wax esters, which co-located with a gene from the KCS family, previously implicated in the biosynthesis of cuticular waxes in Arabidopsis. In combination, the QTL and quantitative genetic analyses suggest the variation and differentiation in cuticular wax compounds within E. globulus has a complex genetic origin. Sub-races exhibited independent latitudinal and longitudinal differentiation in cuticular wax compounds, likely reflecting processes such as historic gene flow and diversifying selection acting upon genes that have diverse functions in distinct biochemical pathways.

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.

On the persistence of reproductive barriers in Eucalyptus: the bridging of mechanical barriers to zygote formation by F1 hybrids is counteracted by intrinsic post-zygotic incompatibilities.

BACKGROUND AND AIMS: Many previous studies conclude that pre-zygotic barriers such as mechanical isolation account for most reproductive isolation between pairs of taxa. However, the inheritance and persistence of barriers such as these after the first generation of hybridization is rarely quantified, even though it is a vital consideration in understanding gene flow potential. There is an asymmetrical pre-zygotic mechanical barrier to hybridization between Eucalyptus nitens and Eucalyptus globulus, which completely prevents small-flowered E. nitens pollen from mating with large E. globulus flowers, while the reverse cross is possible. We aimed to determine the relative importance of pre- and post-zygotic barriers in preventing gene flow following secondary contact between E. nitens and E. globulus, including the inheritance of barriers in advanced-generation hybrids. METHODS: Experimental crossing was used to produce outcrossed E. nitens, E. globulus and their F1, F2, BCg and BCn hybrids. The strength and inheritance of a suite of pre- and post-zygotic barriers were assessed, including 20-year survival, growth and reproductive capacity. KEY RESULTS: The mechanical barrier to hybridization was lost or greatly reduced in the F1 hybrid. In contrast, intrinsic post-zygotic barriers were strong and persistent. Line-cross analysis indicated that the outbreeding depression in the hybrids was best explained by epistatic loss. CONCLUSIONS: The removal of strong mechanical barriers between E. nitens and E. globulus allows F1 hybrids to act as a bridge for bi-directional gene flow between these species. However, strong and persistent post-zygotic barriers exist, meaning that wherever F1 hybridization does occur, intrinsic post-zygotic barriers will be responsible for most reproductive isolation in this system. This potential transient nature of mechanical barriers to zygote formation due to additive inheritance in hybrids appears under-appreciated, and highlights the often important role that intrinsic post-mating barriers play in maintaining species boundaries at zones of secondary contact.

Phylogeny Explains Variation in The Root Chemistry of Eucalyptus Species.

Plants are dependent on their root systems for survival, and thus are defended from belowground enemies by a range of strategies, including plant secondary metabolites (PSMs). These compounds vary among species, and an understanding of this variation may provide generality in predicting the susceptibility of forest trees to belowground enemies and the quality of their organic matter input to soil. Here, we investigated phylogenetic patterns in the root chemistry of species within the genus Eucalyptus. Given the known diversity of PSMs in eucalypt foliage, we hypothesized that (i) the range and concentrations of PSMs and carbohydrates in roots vary among Eucalyptus species, and (ii) that phylogenetic relationships explain a significant component of this variation. To test for interspecific variation in root chemistry and the influence of tree phylogeny, we grew 24 Eucalyptus species representing two subgenera (Eucalyptus and Symphyomyrtus) in a common garden for two years. Fine root samples were collected from each species and analyzed for total phenolics, condensed tannins, carbohydrates, terpenes, and formylated phloroglucinol compounds. Compounds displaying significant interspecific variation were mapped onto a molecular phylogeny and tested for phylogenetic signal. Although all targeted groups of compounds were present, we found that phenolics dominated root defenses and that all phenolic traits displayed significant interspecific variation. Further, these compounds displayed a significant phylogenetic signal. Overall, our results suggest that within these representatives of genus Eucalyptus, more closely related species have more similar root chemistry, which may influence their susceptibility to belowground enemies and soil organic matter accrual.

Genomic patterns of species diversity and divergence in Eucalyptus.

We examined genome-wide patterns of DNA sequence diversity and divergence among six species of the important tree genus Eucalyptus and investigated their relationship with genomic architecture. Using c. 90 range-wide individuals of each Eucalyptus species (E. grandis, E. urophylla, E. globulus, E. nitens, E. dunnii and E. camaldulensis), genetic diversity and divergence were estimated from 2840 polymorphic diversity arrays technology markers covering the 11 chromosomes. Species differentiating markers (SDMs) identified in each of 15 pairwise species comparisons, along with species diversity (HHW ) and divergence (FST ), were projected onto the E. grandis reference genome. Across all species comparisons, SDMs totalled 1.1-5.3% of markers and were widely distributed throughout the genome. Marker divergence (FST and SDMs) and diversity differed among and within chromosomes. Patterns of diversity and divergence were broadly conserved across species and significantly associated with genomic features, including the proximity of markers to genes, the relative number of clusters of tandem duplications, and gene density within or among chromosomes. These results suggest that genomic architecture influences patterns of species diversity and divergence in the genus. This influence is evident across the six species, encompassing diverse phylogenetic lineages, geography and ecology.