ITS sequence data resolve higher level relationships among the eucalypts

Chloroplast variation is incongruent with classification of the Australian bloodwood eucalypts (genus Corymbia, family Myrtaceae)

Previous molecular phylogenetic analyses have resolved the Australian bloodwood eucalypt genus Corymbia (~100 species) as either monophyletic or paraphyletic with respect to Angophora (9-10 species). Here we assess relationships of Corymbia and Angophora using a large dataset of chloroplast DNA sequences (121,016 base pairs; from 90 accessions representing 55 Corymbia and 8 Angophora species, plus 33 accessions of related genera), skimmed from high throughput sequencing of genomic DNA, and compare results with new analyses of nuclear ITS sequences (119 accessions) from previous studies. Maximum likelihood and maximum parsimony analyses of cpDNA resolve well supported trees with most nodes having >95% bootstrap support. These trees strongly reject monophyly of Corymbia, its two subgenera (Corymbia and Blakella), most taxonomic sections (Abbreviatae, Maculatae, Naviculares, Septentrionales), and several species. ITS trees weakly indicate paraphyly of Corymbia (bootstrap support <50% for maximum likelihood, and 71% for parsimony), but are highly incongruent with the cpDNA analyses, in that they support monophyly of both subgenera and some taxonomic sections of Corymbia. The striking incongruence between cpDNA trees and both morphological taxonomy and ITS trees is attributed largely to chloroplast introgression between taxa, because of geographic sharing of chloroplast clades across taxonomic groups. Such introgression has been widely inferred in studies of the related genus Eucalyptus. This is the first report of its likely prevalence in Corymbia and Angophora, but this is consistent with previous morphological inferences of hybridisation between species. Our findings (based on continent-wide sampling) highlight a need for more focussed studies to assess the extent of hybridisation and introgression in the evolutionary history of these genera, and that critical testing of the classification of Corymbia and Angophora requires additional sequence data from nuclear genomes.

Special trends in CBF and DREB2 groups in Eucalyptus gunnii vs Eucalyptus grandis suggest that CBF are master players in the trade-off between growth and stress resistance

Annotation of the Eucalyptus grandis genome showed a large amplification of the dehydration-responsive element binding 1/C-repeat binding factor (DREB1/CBF) group without recent DREB2 gene duplication compared with other plant species. The present annotation of the CBF and DREB2 genes from a draft of the Eucalyptus gunnii genome sequence reveals at least one additional CBF copy in the E. gunnii genome compared with E. grandis, suggesting that this group is still evolving, unlike the DREB2 group. This study aims to investigate the redundancy/neo- or sub-functionalization of the duplicates and the relative involvement of the two groups in abiotic stress responses in both E. grandis and E. gunnii (lower growth but higher cold resistance). A comprehensive transcriptional analysis using high-throughput quantitative real-time polymerase chain reaction (qRT-PCR) was performed on leaves, stems and roots from the two Eucalyptus species after cold, heat or drought treatment. A large CBF cluster accounted for most of the cold response in all the organs, whereas heat and drought responses mainly involved a small CBF cluster and the DREB2 genes. In addition, CBF putative target genes, known to be involved in plant tolerance and development, were found to be cold-regulated. The higher transcript amounts of both the CBF and target genes in the cold tolerant E. gunnii contrasted with the higher CBF induction rates in the fast growing E. grandis. Altogether, the present results, in agreement with previous data about Eucalyptus transgenic lines over-expressing CBF, suggest that these factors, which promote both stress protection and growth limitation, participate in the trade-off between growth and resistance in this woody species.

Early physiological flood tolerance is followed by slow post-flooding root recovery in the dryland riparian tree Eucalyptus camaldulensis subsp. refulgens

We investigated physiological and morphological responses to flooding and recovery in Eucalyptus camaldulensis subsp. refulgens, a riparian tree species from a dryland region prone to intense episodic floods. Seedlings in soil flooded for 88 d produced extensive adventitious roots, displayed stem hypertrophy (stem diameter increased by 93%) and increased root porosity owing to aerenchyma formation. Net photosynthesis (Pn) and stomatal conductance (gs) were maintained for at least 2 weeks of soil flooding, contrasting with previous studies of other subspecies of E. camaldulensis. Gradual declines followed in both gs (30% less than controls) and Pn (19% less). Total leaf soluble sugars did not differ between flooded and control plants. Root mass did not recover 32 d after flooding ceased, but gs was not lower than controls, suggesting the root system was able to functionally compensate. However, the limited root growth during recovery after flooding was surprising given the importance of extensive root systems in dryland environments. We conclude that early flood tolerance could be an adaptation to capitalize on scarce water resources in a water-limited environment. Overall, our findings highlight the need to assess flooding responses in relation to a species' fitness for particular flood regimes or ecological niches.

Shifts in species interactions due to the evolution of functional differences between endemics and non-endemics: an endemic syndrome hypothesis

Species ranges have been shifting since the Pleistocene, whereby fragmentation, isolation, and the subsequent reduction in gene flow have resulted in local adaptation of novel genotypes and the repeated evolution of endemic species. While there is a wide body of literature focused on understanding endemic species, very few studies empirically test whether or not the evolution of endemics results in unique function or ecological differences relative to their widespread congeners; in particular while controlling for environmental variation. Using a common garden composed of 15 Eucalyptus species within the subgenus Symphyomyrtus (9 endemic to Tasmania, 6 non-endemic), here we hypothesize and show that endemic species are functionally and ecologically different from non-endemics. Compared to non-endemics, endemic Eucalyptus species have a unique suite of functional plant traits that have extended effects on herbivores. We found that while endemics occupy many diverse habitats, they share similar functional traits potentially resulting in an endemic syndrome of traits. This study provides one of the first empirical datasets analyzing the functional differences between endemics and non-endemics in a common garden setting, and establishes a foundation for additional studies of endemic/non-endemic dynamics that will be essential for understanding global biodiversity in the midst of rapid species extinctions and range shifts as a consequence of global change.

Leaf hydraulic vulnerability influences species' bioclimatic limits in a diverse group of woody angiosperms

The ability of plants to maintain water flow through leaves under water stress-induced tension (assessed as the leaf hydraulic vulnerability; P50(leaf)) is intimately linked with survival. We examined the significance of P50(leaf) as an adaptive trait in influencing the dry-end distributional limits of cool temperate woody angiosperm species. We also examined differences in within-site variability in P50(leaf) between two high-rainfall montane rainforest sites in Tasmania and Peru, respectively. A significant relationship between P50(leaf) and the 5th percentile of mean annual rainfall across each species distribution was found in Tasmania, suggesting that P50(leaf) influences species climatic limits. Furthermore, a strong correlation between P50(leaf) and the minimum rainfall availability was found using five phylogenetically independent species pairs in wet and dry evergreen tree species, suggesting that rainfall is an important selective agent in the evolution of leaf hydraulic vulnerability. Greater within-site variability in P50(leaf) was found among dominant montane rainforest species in Tasmania than in Peru and this result is discussed within the context of differences in spatial and temporal environmental heterogeneity and parochial historical ecology.

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.

An AFLP marker approach to lower-level systematics in Eucalyptus (Myrtaceae)

Genus Eucalyptus, with over 700 species, presents a number of systematic difficulties including taxa that hybridize or intergrade across environmental gradients. To date, no DNA marker has been found capable of resolving phylogeny below the sectional level in the major subgenera. Molecular markers are needed to support taxonomic revision, assess the extent of genetic divergence at lower taxonomic levels, and inform conservation efforts. We examined the utility of 930 amplified fragment length polymorphisms (AFLPs) for analyzing relationships among Tasmanian taxa of subgenus Symphyomyrtus section Maidenaria. Phenetic and cladistic analyses resolved species into clusters demonstrating significant genetic partitioning, largely concordant with series defined in the most recent taxonomic revision of Eucalyptus. Some departures from current taxonomy were noted, indicating possible cases of morphological convergence and character reversion. Although the resolution obtained using AFLP was greatly superior to that of single sequence markers, the data demonstrated high homoplasy and incomplete resolution of closely related species. The results of this study and others are consistent with recent speciation and reticulate evolution in Maidenaria. We conclude that a combination of phylogenetic and population genetic approaches using multiple molecular markers offers the best prospects for understanding taxonomic relationships below the sectional level in Eucalyptus.

Nuclear ribosomal pseudogenes resolve a corroborated monophyly of the eucalypt genus Corymbia despite misleading hypotheses at functional ITS paralogs

Divergent paralogs can create both obstacles and opportunities for phylogenetic reconstruction. Phylogenetic relationships among eucalypt genera have been incongruent among datasets in previous studies, where morphological characters supported monophyly of the genus Corymbia, while intergenic spacers of the nuclear ribosomes (ITS) and chloroplast loci (trnL, trnH, psbA) showed Corymbia as either equivocal or paraphyletic. Ribosomal DNA occurs in multiple copies in a genome. We cloned and sequenced the nrITS to investigate if gene duplication was the cause of incongruence among trees in the eucalypts. Three ITS riboforms, two of them widespread, were recovered within some genomes. One of the ITS riboforms recovered a robust phylogeny showing Corymbia as a monophyletic genus, corroborating the evidence from morphology, fossil data, a recent ITS/ETS dataset and microsatellites (SSRs). Compelling evidence suggested that this divergent riboform is a pseudogene, i.e., non-functional paralog: comparatively lower GC content suggesting lower structural stability, deamination-like mutations at potential methylation sites, lack of conserved helices and hairpins and conspicuously lower thermodynamic stability in secondary structures. Phylogenies from the apparently functional riboform retained Corymbia as paraphyletic. We show here that pseudogenes can recover a well-corroborated phylogeny whereas their functional paralogs show misleading hypotheses. We explain that phylogenetic signals may be obscured when functional constraints in ITS necessitate compensatory mutations in the secondary structure helices involved in RNA transcription, whereas pseudogenes mutate under neutrality.

The impact of intragenic recombination on phylogenetic reconstruction at the sectional level in Eucalyptus when using a single copy nuclear gene (cinnamoyl CoA reductase)

Low copy number nuclear genes have been found to be useful for phylogenetic reconstruction at different taxonomic levels. This study investigated the utility of a single copy gene, cinnamoyl CoA reductase (CCR), for resolving phylogenetic relationships at the sectional level within Eucalyptus. The monophyly of sections Exsertaria and Latoangulatae was explored, using section Maidenaria as an outgroup, and the impact of intragenic recombination on phylogenetic reconstruction examined. Phylogenetic analysis did not resolve monophyletic groups. Latoangulatae and Maidenaria were polyphyletic or paraphyletic. Exsertaria species formed a clade but included a single Latoangulatae species (E. major). Recombination analysis identified two intragenic recombination events that involved species from different sections, which have probably been facilitated by inter-sectional hybridisation. One of the events most likely occurred prior to speciation, with several Latoangulatae species having the recombinant allele. The other event may have occurred after speciation, since only one of two E. globulus samples possessed the recombinant allele. This is the first detailed report of intragenic recombination in both CCR and Eucalyptus, and between species of different sections of a plant genus. The occurrence of intragenic recombination may explain the anomalous positions of some species within the phylogenetic tree, and indicates that phylogenetic analysis of Eucalyptus using nuclear genes will be problematic unless recombination is taken into account.

Genomic research in Eucalyptus

Eucalyptus L'Hérit. is a genus comprised of more than 700 species that is of vital importance ecologically to Australia and to the forestry industry world-wide, being grown in plantations for the production of solid wood products as well as pulp for paper. With the sequencing of the genomes of Arabidopsis thaliana and Oryza sativa and the recent completion of the first tree genome sequence, Populus trichocarpa, attention has turned to the current status of genomic research in Eucalyptus. For several eucalypt species, large segregating families have been established, high-resolution genetic maps constructed and large EST databases generated. Collaborative efforts have been initiated for the integration of diverse genomic projects and will provide the framework for future research including exploiting the sequence of the entire eucalypt genome which is currently being sequenced. This review summarises the current position of genomic research in Eucalyptus and discusses the direction of future research.

The rare silver gum, Eucalyptus cordata, is leaving its trace in the organellar gene pool of Eucalyptus globulus

The process of genetic assimilation of rare species by hybridizing congeners has been documented in a number of plant genera. This raises the possibility that some of the genetic diversity found in phylogeographical studies of widespread species has been acquired through hybridization with species that are now rare or extinct. In this fine-scale phylogeographical analysis, we show that a rare eucalypt species is leaving its trace in the chloroplast genome of a more abundant congener. The heart-leafed silver gum, Eucalyptus cordata, is a rare endemic of south-eastern Tasmania. Its populations are scattered amidst populations of more abundant related species, including the Tasmanian blue gum, Eucalyptus globulus. Using 339 samples from across the full range of both species, we compared chloroplast (cp) DNA haplotype phylogeography in E. globulus and E. cordata. The genealogy and distribution of chloroplast haplotypes suggest that E. globulus has acquired cpDNA from E. cordata in at least four different mixed populations. Shared haplotypes are highest in E. globulus sampled within 2 km of known E. cordata populations and drop to zero at a distance of 25 km from the nearest known E. cordata population. Localized haplotype sharing occurs in the absence of obvious hybrid zones or locally shared nuclear ribosomal DNA sequences. Given that the future loss of E. cordata from some mixed populations is likely, these findings indicate that phylogeographical analyses of organellar DNA should consider the possibility of introgression, even from species that have been eliminated from the sites of interest.

Radiation of the Australian Salicornioideae (Chenopodiaceae)--based on evidence from nuclear and chloroplast DNA sequences

In phylogenetic analyses of nuclear ITS and chloroplast trnL DNA sequences, the mostly endemic Australian genera; Halosarcia, Pachycornia, Sclerostegia, Tecticornia, and Tegicornia of the subfamily Salicornioideae (Chenopodiaceae) together form a monophyletic group, congruent with the hypothesis that they evolved from a common ancestor. However, limited genetic differentiation evident in both nrDNA and cpDNA sequences among these taxa suggests a possible rapid radiation. Based on fossil pollen records and climatic models of other authors, it is hypothesized that the expansion of the Australian endemic Salicornioideae most likely occurred during the Late Miocene to Pliocene, when increasing aridity caused the formation of extensive salt lakes along endorheic paleodrainage channels. Moreover, Australian Sarcocornia representatives were supported as monophyletic, nested within a paraphyletic Sarcocornia clade that also comprised European Salicornia in the ITS analysis. This suggests that Sarcocornia arrived in Australia subsequent to the ancestor of the Australian endemic genera most likely via long-distance dispersal.

EgLFY, the Eucalyptus grandis homolog of the Arabidopsis gene LEAFY is expressed in reproductive and vegetative tissues

The EgLFY gene cloned from Eucalyptus grandis has sequence homology to the floral meristem identity gene LEAFY (LFY) from Arabidopsis and FLORICAULA (FLO) from Antirrhinum. EgLFY is preferentially expressed in the developing eucalypt floral organs in a pattern similar to that described previously for the Arabidopsis LFY. In situ hybridization experiments have shown that EgLFY is strongly expressed in the early floral meristem and then successively in the primordia of sepals, petals, stamens and carpels. It is also expressed in the leaf primordia of adult trees. The expression of the EgLFY coding region under control of the Arabidopsis LFY promoter could complement strong lfy mutations in transgenic Arabidopsis plants. These data suggest that EgLFY plays a similar role to LFY in flower development and that the basic mechanisms involved in flower initiation and development in Eucalyptus may be similar to those occurring in Arabidopsis.

Glacial refugia and reticulate evolution: the case of the Tasmanian eucalypts

Tasmania is a natural laboratory for investigating the evolutionary processes of the Quaternary. It is a large island lying 40-44 degrees S, which was repeatedly glaciated and linked to southeastern continental Australia during the Quaternary. Climate change promoted both the isolation of species in glacial refugia, and an exchange between Tasmanian and mainland floras. Eucalyptus is a complex and diverse genus, which has increased in abundance in Australia over the past 100 kyr, probably in response to higher fire frequency. Morphological evidence suggests that gene flow may have occurred between many eucalypt species after changes in their distribution during the Quaternary. This paper summarizes recent genetic evidence for migration and introgressive hybridization in Tasmanian Eucalyptus. Maternally inherited chloroplast DNA reveals a long-term persistence of eucalypts in southeastern Tasmanian refugia, coupled with introgressive hybridization involving many species. Detailed analysis of the widespread species Eucalyptus globulus suggests that migration from mainland Australia was followed by introgression involving a rare Tasmanian endemic. The data support the hypothesis that changes in distribution of interfertile species during the Quaternary have promoted reticulate evolution in Eucalyptus.

DNA sequence variation in BpMADS2 gene in two populations of Betula pendula

The PISTILLATA (PI) homologue, BpMADS2, was isolated from silver birch (Betula pendula Roth) and used to study nucleotide polymorphism. Two regions (together about 2450 bp) comprising mainly untranslated sequences were sequenced from 10 individuals from each of two populations in Finland. The nucleotide polymorphism was low in the BpMADS2 locus, especially in the coding region. The synonymous site overall nucleotide diversity (pis) was 0.0043 and the nonsynonymous nucleotide diversity (pia) was only 0.000052. For the whole region, the pi values for the two populations were 0.0039 and 0.0045, and for the coding regions, the pi values were only 0 and 0.00066 (for the corresponding coding regions of Arabidopsis thaliana PI world-wide pi was 0.0021). Estimates of pi or theta did not differ significantly between the two populations, and the two populations were not diverged from each other. Two classes of BpMADS2 alleles were present in both populations, suggesting that this gene exhibits allelic dimorphism. In addition to the nucleotide site variation, two microsatellites were also associated within the haplotypes. This allelic dimorphism might be the result of postglacial re-colonization partly from northwestern, partly from southeastern/eastern refugia. The sequence comparison detected five recombination events in the regions studied. The large number of microsatellites in all of the three introns studied suggests that BpMADS2 is a hotspot for microsatellite formation.

Evolution and biogeography of Alectryon (Sapindaceae)

Phylogenetic analysis of nucleotide sequences from four plastid loci (matK, partial trnK-matK introns, rps16 intron) and one nuclear locus (the internal transcribed spacer of rDNA; ITS-1) was conducted for 14 species of Alectryon and five related genera in Sapindaceae. Both matK and rps16 intron provide few informative characters within Alectryon, whereas ITS-1 provides the largest number of parsimony-informative characters and has the greatest sequence divergence between taxa. Support for branches in cladograms produced in PAUP increased markedly upon inclusion of ITS-1 data to matK and rps16 intron data. Analyses of each region alone or combined produced congruent results, suggesting that the regions are complementary. Phylogenetic analysis indicates that there are two main lineages within Alectryon, with A. subcinereus sister to the remaining sampled Alectryon taxa. Two morphological characters, presence/absence of petals and aril patterning, are congruent with the molecular phylogeny. One robustly supported clade is characterized by smooth arils and petals, in contrast to the taxa in the other major clade which have patterned arils and an absence of petals. These analyses also support a number of revised subgeneric groupings for Alectryon. The decision to submerge Heterodendrum in Alectryon is supported, although taxa belonging to Heterodendrum do not form a clade. The majority of the Australian Alectryon appear to belong to the tropical monsoonal/arid flora with species from both lineages being found in representative vine thickets across northern Australia. It appears that the seasonally dry rainforest communities comprise a number of elements that do not share common evolutionary histories within this genus.

Why nuclear ribosomal DNA spacers (ITS) tell different stories in Quercus

The molecular systematics of Quercus (Fagaceae) was recently assessed by two teams using independently generated ITS sequences. Although the results disagreed in several remarkable features, the phylogenetic trees for either hypothesis were highly supported by bootstrap resampling. We have reanalyzed the ITS sequences used by both teams (eight taxa) to reveal the underlying patterns of this divergence. Within species, conspicuous length and G + C% divergence were evident in most sequence comparisons. In addition, a high rate of substitutions and deletions involving highly conserved motifs in both ITS spacers were present in a set of sequences. This was coupled with a less thermodynamic stability in the RNA structure, lacking conserved hairpins that are putatively involved in the processing of the RNA transcripts. Compelling evidence suggests that the divergent ITS alleles reported by one team are pseudogenes, i.e., nonfunctional paralogous loci. The hypothesis that the contrasting phylogenetic histories drawn from Quercus using ITS data are not strictly related to technical differences between laboratories, but that they have rather been generated from the analysis of paralogous sequences, best reconciles the available data. The risk of incorporating ITS paralogues in plant evolutionary studies which can distortion the phylogenetic signal should caution molecular systematists. Without a detailed inspection of some basic features of the sequence, including the integrity of the conserved motifs and the thermodynamic stability of the secondary structures of the RNA transcripts, errors in evolutionary inferences could be easily overlooked.

Chloroplast sharing in the Tasmanian eucalypts

The biogeographic pattern of chloroplast DNA (cpDNA) haplotypes in Eucalyptus on the island of Tasmania is consistent with reticulate evolution, involving at least 12 Tasmanian species from the subgenus Symphyomyrtus. Intraspecific cpDNA polymorphism in 14 of 17 species is coupled with extensive sharing of identical haplotypes across populations of different species in the same geographic area. Haplotype diversity is lowest in central regions of Tasmania formerly occupied by alpine vegetation during glacial intervals and in northern regions that were periodically linked to continental Australia by land bridges. The observed distribution of several cpDNA haplotypes unique to Tasmania coincides with modeled locations of glacial refugia in coastal areas of Tasmania and shows the power of cpDNA in unraveling the complex history of past distributions of Eucalyptus. The results suggest that the model of evolution of the eucalypts should be reassessed to allow for the anastomosing effects of interspecific hybridization and introgression.