Reticulate evolution in the natural range of the invasive wetland tree species Melaleuca quinquenervia

Characterization of terpene biosynthesis in Melaleuca quinquenervia and ecological consequences of terpene accumulation during myrtle rust infection

Plants use a wide array of secondary metabolites including terpenes as defense against herbivore and pathogen attack, which can be constitutively expressed or induced. Here, we investigated aspects of the chemical and molecular basis of resistance against the exotic rust fungus Austropuccinia psidii in Melaleuca quinquenervia, with a focus on terpenes. Foliar terpenes of resistant and susceptible plants were quantified, and we assessed whether chemotypic variation contributed to resistance to infection by A. psidii. We found that chemotypes did not contribute to the resistance and susceptibility of M. quinquenervia. However, in one of the chemotypes (Chemotype 2), susceptible plants showed higher concentrations of several terpenes including α-pinene, limonene, 1,8-cineole, and viridiflorol compared with resistant plants. Transcriptome profiling of these plants showed that several TPS genes were strongly induced in response to infection by A. psidii. Functional characterization of these TPS showed them to be mono- and sesquiterpene synthases producing compounds including 1,8-cineole, β-caryophyllene, viridiflorol and nerolidol. The expression of these TPS genes correlated with metabolite data in a susceptible plant. These results suggest the complexity of resistance mechanism regulated by M. quinquenervia and that modulation of terpenes may be one of the components that contribute to resistance against A. psidii.

Identifying genetic markers for a range of phylogenetic utility-From species to family level

Resolving the phylogenetic relationships of closely related species using a small set of loci is challenging as sufficient information may not be captured from a limited sample of the genome. Relying on few loci can also be problematic when conflict between gene-trees arises from incomplete lineage sorting and/or ongoing hybridization, problems especially likely in recently diverged lineages. Here, we developed a method using limited genomic resources that allows identification of many low copy candidate loci from across the nuclear and chloroplast genomes, design probes for target capture and sequence the captured loci. To validate our method we present data from Eucalyptus and Melaleuca, two large and phylogenetically problematic genera within the Myrtaceae family. With one annotated genome, one transcriptome and two whole-genome shotgun sequences of one Eucalyptus and four Melaleuca species, respectively, we identified 212 loci representing 263 kbp for targeted sequence capture and sequencing. Of these, 209 were successfully tested from 47 samples across five related genera of Myrtaceae. The average percentage of reads mapped back to the reference was 57.6% with coverage of more than 20 reads per position across 83.5% of the data. The methods developed here should be applicable across a large range of taxa across all kingdoms. The core methods are very flexible, providing a platform for various genomic resource availabilities and are useful from shallow to deep phylogenies.

Interpreting the modern distribution of Myrtaceae using a dated molecular phylogeny

The angiosperm family Myrtaceae has extant and fossil taxa from all southern continents and is assumed to be of Gondwanan origin. Many modern groups contain sister taxa that have disjunct transoceanic distributions, which can be interpreted as a result of either vicariance or long-distance dispersal and establishment (LDDE). Further, some Myrtaceae groups occur on Pacific islands with enigmatic geological histories. We tested hypotheses of vicariance and LDDE by estimating divergence times using a relaxed molecular clock calibrated with 12 fossils. In total, 88 genera and 202 species were sampled, representing both subfamilies and all tribes of Myrtaceae. We reconstructed the family as Gondwanan in origin. Of the 22 geographically disjunct sister groups in our study, up to six are potentially explained as the product of vicariance, three resulting from overland dispersal via new land connections, and 13 due to LDDE events. Nine of the 13 hypothesized LDDE events occurred in fleshy-fruited taxa. Our results indicate that most of the transoceanic distribution patterns in Myrtaceae have occurred since the Miocene due to LDDE, whereas inferred vicariance events all occurred before the Late Eocene. There are many instances of sister relationships between species-poor and species-rich groups in Myrtaceae, and at least three occurrences of geographically isolated taxa on long branches of the phylogeny (Arillastrum, Myrtus, and Tepualia), whose modern-day distributions are difficult to explain without additional fossil or geological evidence.

Clock model makes a large difference to age estimates of long-stemmed clades with no internal calibration: a test using Australian grasstrees

BACKGROUND

Estimating divergence times in phylogenies using a molecular clock depends on accurate modeling of nucleotide substitution rates in DNA sequences. Rate heterogeneity among lineages is likely to affect estimates, especially in lineages with long stems and short crowns ("broom" clades) and no internal calibration. We evaluate the performance of the random local clocks model (RLC) and the more routinely employed uncorrelated lognormal relaxed clock model (UCLN) in situations in which a significant rate shift occurs on the stem branch of a broom clade. We compare the results of simulations to empirical results from analyses of a real rate-heterogeneous taxon - Australian grass trees (Xanthorrhoea) - whose substitution rate is slower than in its sister groups, as determined by relative rate tests.

RESULTS

In the simulated datasets, the RLC model performed much better than UCLN: RLC correctly estimated the age of the crown node of slow-rate broom clades, whereas UCLN estimates were consistently too young. Similarly, in the Xanthorrhoea dataset, UCLN returned significantly younger crown ages than RLC (mean estimates respectively 3-6 Ma versus 25-35 Ma). In both real and simulated datasets, Bayes Factor tests strongly favored the RLC model over the UCLN model.

CONCLUSIONS

The choice of an unsuitable molecular clock model can strongly bias divergence time estimates. In particular, for data predicted to have more rate variation among than within clades, dating with RLC is much more likely to be accurate than with UCLN. The choice of clocks should be informed by the biology of the study group (e.g., life-form) or assessed with relative rate tests and post-hoc model comparisons.

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.