Phylogenomics reveals extreme gene tree discordance in a lineage of dominant trees: hybridization, introgression, and incomplete lineage sorting blur deep evolutionary relationships despite clear species groupings in Eucalyptus subgenus Eudesmia.

Eucalypts are a large and ecologically important group of plants on the Australian continent, and understanding their evolution is important in understanding evolution of the unique Australian flora. Previous phylogenies using plastome DNA, nuclear-ribosomal DNA, or random genome-wide SNPs, have been confounded by limited genetic sampling or by idiosyncratic biological features of the eucalypts, including widespread plastome introgression. Here we present phylogenetic analyses of Eucalyptus subgenus Eudesmia (22 species from western, northern, central and eastern Australia), in the first study to apply a target-capture sequencing approach using custom, eucalypt-specific baits (of 568 genes) to a lineage of Eucalyptus. Multiple accessions of all species were included, and target-capture data were supplemented by separate analyses of plastome genes (average of 63 genes per sample). Analyses revealed a complex evolutionary history likely shaped by incomplete lineage sorting and hybridization. Gene tree discordance generally increased with phylogenetic depth. Species, or groups of species, toward the tips of the tree are mostly supported, and three major clades are identified, but the branching order of these clades cannot be confirmed with confidence. Multiple approaches to filtering the nuclear dataset, by removing genes or samples, could not reduce gene tree conflict or resolve these relationships. Despite inherent complexities in eucalypt evolution, the custom bait kit devised for this research will be a powerful tool for investigating the evolutionary history of eucalypts more broadly.

A box on the river: The phylogenetics and phylogeography of Eucalyptus baueriana (Eucalyptus sect. Adnataria ser. Heterophloiae).

We present a phylogeographic study of the tree species Eucalyptus baueriana Schauer, which occurs in disjunct areas on the near coastal plains and ranges of the south-east Australian mainland. DArTseq data are used to build a phylogeny including E. baueriana and closely related taxa to test its monophyly, test the genetic distinctness of the three subspecies of E. baueriana, and investigate relationships between its disjunct populations. Additionally, we use population structure analysis to investigate the genetic distinctness of populations, and MaxEnt to investigate the environmental factors potentially influencing the species' distribution. We show E. baueriana is monophyletic and most closely related to three other Blue Box eucalypt species: E. conica H.Deane & Maiden, E. dalveenica T.L.Collins, R.L.Andrew & J.J.Bruhl and E. magnificata L.A.S.Johnson & K.D.Hill, with some evidence for genetic introgression between these taxa. Within E. baueriana, the deepest genetic breaks do not correspond with the subspecies classification as the two geographically restricted subspecies, together with samples of the more widespread E. baueriana subsp. baueriana from west of the Gippsland lowlands, form a south-western clade with that is sister to other populations of subsp. baueriana. The oldest genetic break in the species occurs in far eastern Gippsland (Victoria), corresponding to one of the shortest geographic disjunctions in the species' distribution. Genetic breaks in other species have been observed in this region which is broadly referred to as the southern transition zone. Both total annual rainfall and the seasonality of this rainfall are hypothesised to affect the species' distribution; gaps in its distribution are in areas of higher rainfall that support closed forest and in regions with more winter dominated rainfall.

A genome resource for Acacia, Australia's largest plant genus.

Acacia (Leguminosae, Caesalpinioideae, mimosoid clade) is the largest and most widespread genus of plants in the Australian flora, occupying and dominating a diverse range of environments, with an equally diverse range of forms. For a genus of its size and importance, Acacia currently has surprisingly few genomic resources. Acacia pycnantha, the golden wattle, is a woody shrub or tree occurring in south-eastern Australia and is the country's floral emblem. To assemble a genome for A. pycnantha, we generated long-read sequences using Oxford Nanopore Technology, 10x Genomics Chromium linked reads, and short-read Illumina sequences, and produced an assembly spanning 814 Mb, with a scaffold N50 of 2.8 Mb, and 98.3% of complete Embryophyta BUSCOs. Genome annotation predicted 47,624 protein-coding genes, with 62.3% of the genome predicted to comprise transposable elements. Evolutionary analyses indicated a shared genome duplication event in the Caesalpinioideae, and conflict in the relationships between Cercis (subfamily Cercidoideae) and subfamilies Caesalpinioideae and Papilionoideae (pea-flowered legumes). Comparative genomics identified a suite of expanded and contracted gene families in A. pycnantha, and these were annotated with both GO terms and KEGG functional categories. One expanded gene family of particular interest is involved in flowering time and may be associated with the characteristic synchronous flowering of Acacia. This genome assembly and annotation will be a valuable resource for all studies involving Acacia, including the evolution, conservation, breeding, invasiveness, and physiology of the genus, and for comparative studies of legumes.

Divergent lineages in a semi-arid mallee species, Eucalyptus behriana, correspond to a major geographic break in southeastern Australia.

AIM: To infer relationships between populations of the semi-arid, mallee eucalypt, Eucalyptus behriana, to build hypotheses regarding evolution of major disjunctions in the species' distribution and to expand understanding of the biogeographical history of southeastern Australia. LOCATION: Southeastern Australia. TAXON: Eucalyptus behriana (Myrtaceae, Angiospermae). METHODS: We developed a large dataset of anonymous genomic loci for 97 samples from 11 populations of E. behriana using double digest restriction site-associated DNA sequencing (ddRAD-seq), to determine genetic relationships between the populations. These relationships, along with species distribution models, were used to construct hypotheses regarding environmental processes that have driven fragmentation of the species' distribution. RESULTS: Greatest genetic divergence was between populations on either side of the Lower Murray Basin. Populations west of the Basin showed greater genetic divergence between one another than the eastern populations. The most genetically distinct population in the east (Long Forest) was separated from others by the Great Dividing Range. A close relationship was found between the outlying northernmost population (near West Wyalong) and those in the Victorian Goldfields despite a large disjunction between them. CONCLUSIONS: Patterns of genetic variation are consistent with a history of vicariant differentiation of disjunct populations. We infer that an early disjunction to develop in the species distribution was that across the Lower Murray Basin, an important biogeographical barrier separating many dry sclerophyll plant taxa in southeastern Australia. Additionally, our results suggest that the western populations fragmented earlier than the eastern ones. Fragmentation, both west and east of the Murray Basin, is likely tied to climatic changes associated with glacial-interglacial cycles although it remains possible that major geological events including uplift of the Mount Lofty Ranges and basalt flows in the Newer Volcanics Province also played a role.

Long-read assemblies reveal structural diversity in genomes of organelles - an example with Acacia pycnantha.

Organelle genomes are typically represented as single, static, circular molecules. However, there is evidence that the chloroplast genome exists in two structural haplotypes and that the mitochondrial genome can display multiple circular, linear or branching forms. We sequenced and assembled chloroplast and mitochondrial genomes of the Golden Wattle, Acacia pycnantha, using long reads, iterative baiting to extract organelle-only reads, and several assembly algorithms to explore genomic structure. Using a de novo assembly approach agnostic to previous hypotheses about structure, we found that different assemblies revealed contrasting arrangements of genomic segments; a hypothesis supported by mapped reads spanning alternate paths.

Chloroplast genome analysis of Australian eucalypts--Eucalyptus, Corymbia, Angophora, Allosyncarpia and Stockwellia (Myrtaceae).

We present a phylogenetic analysis and comparison of structural features of chloroplast genomes for 39 species of the eucalypt group (genera Eucalyptus, Corymbia, Angophora, and outgroups Allosyncarpia and Stockwellia). We use 41 complete chloroplast genome sequences, adding 39 finished-quality chloroplast genomes to two previously published genomes. Maximum parsimony and Bayesian analyses, based on >7000 variable nucleotide positions, produced one fully resolved phylogenetic tree (35 supported nodes, 27 with 100% bootstrap support). Eucalyptus and its sister lineage Angophora+Corymbia show a deep divergence. Within Eucalyptus, three lineages are resolved: the 'eudesmid', 'symphyomyrt' and 'monocalypt' groups. Corymbia is paraphyletic with respect to Angophora. Gene content and order do not vary among eucalypt chloroplasts; length mutations, especially frame shifts, are uncommon in protein-coding genes. Some non-synonymous mutations are highly incongruent with the overall phylogenetic signal, notably in rbcL, and may be adaptive. Application of custom informatics pipelines (GYDLE Inc.) enabled direct chloroplast genome assembly, resolving each genome to finished-quality with no need for PCR gap-filling or contig order resolution. Analysis of whole chloroplast genomes resolved major eucalypt clades and revealed variable regions of the genome that will be useful in lower-level genetic studies (including phylogeography and geneflow).

Testing the impact of calibration on molecular divergence times using a fossil-rich group: the case of Nothofagus (Fagales).

Although temporal calibration is widely recognized as critical for obtaining accurate divergence-time estimates using molecular dating methods, few studies have evaluated the variation resulting from different calibration strategies. Depending on the information available, researchers have often used primary calibrations from the fossil record or secondary calibrations from previous molecular dating studies. In analyses of flowering plants, primary calibration data can be obtained from macro- and mesofossils (e.g., leaves, flowers, and fruits) or microfossils (e.g., pollen). Fossil data can vary substantially in accuracy and precision, presenting a difficult choice when selecting appropriate calibrations. Here, we test the impact of eight plausible calibration scenarios for Nothofagus (Nothofagaceae, Fagales), a plant genus with a particularly rich and well-studied fossil record. To do so, we reviewed the phylogenetic placement and geochronology of 38 fossil taxa of Nothofagus and other Fagales, and we identified minimum age constraints for up to 18 nodes of the phylogeny of Fagales. Molecular dating analyses were conducted for each scenario using maximum likelihood (RAxML + r8s) and Bayesian (BEAST) approaches on sequence data from six regions of the chloroplast and nuclear genomes. Using either ingroup or outgroup constraints, or both, led to similar age estimates, except near strongly influential calibration nodes. Using "early but risky" fossil constraints in addition to "safe but late" constraints, or using assumptions of vicariance instead of fossil constraints, led to older age estimates. In contrast, using secondary calibration points yielded drastically younger age estimates. This empirical study highlights the critical influence of calibration on molecular dating analyses. Even in a best-case situation, with many thoroughly vetted fossils available, substantial uncertainties can remain in the estimates of divergence times. For example, our estimates for the crown group age of Nothofagus varied from 13 to 113 Ma across our full range of calibration scenarios. We suggest that increased background research should be made at all stages of the calibration process to reduce errors wherever possible, from verifying the geochronological data on the fossils to critical reassessment of their phylogenetic position.

Historical biogeographical patterns in continental Australia: congruence among areas of endemism of two major clades of eucalypts.

Published phylogenies of two eucalypt clades, red bloodwoods Corymbia subgenus Corymbia and eudesmids Eucalyptus subgenus Eudesmia (Myrtaceae), were combined for an analysis of historical biogeographical area relationships within continental Australia. The method of paralogy-free subtree analysis was used to eliminate geographical paralogy; the paralogy-free subtrees were coded as characters for parsimony analysis to find the minimal and area cladogram, which proved to be informative of a continent-wide pattern. The eucalypt fossil record and molecular dating studies allow an interpretation of the biogeographical history in terms of major vicariance events that date from the early Paleogene. The summary area cladogram shows the wet jarrah forest region of South-West Western Australia, a region of high endemism, as the earliest to differentiate from all other areas, isolated by marine inundation across southern Australia and climatic cooling in the Late Eocene-Early Oligocene. From about this time, regionalization continued, with warmer conditions and monsoonal climate developing in central and northern Australia, and cooling in the south-east. Northern and eastern humid and semi-humid areas were related as a track, but with increased aridity in the interior of the continent, the monsoonal climate contracted northwards. The Australian Monsoon Tropics (AMT: Kimberley, Top End, Arnhem, Cape York and inland north-east Queensland) differentiated from eastern areas (Queensland wet tropics to McPherson-Macleay). Our results also show all arid and semi-arid regions as related, suggestive of a historically cohesive interior biota rather than repeated colonizations of the interior from the periphery of the continent. Climate largely differentiates hot arid areas in the north (Pilbara, Northern and Central deserts) from arid areas in the south (south-west interzone, Wheatbelt, Goldfields and Great Victoria Desert). © The Willi Hennig Society 2010.

Description and affinities of a new sequestrate fungus, Barcheria willisiana gen. et sp. nov. (Agaricales) from Australia.

A new sequestrate fungus, Barcheria willisiana gen. et sp. nov., is described and its affinities evaluated using nLSU rDNA sequence data. This unusual fungus has several characters that are reminiscent of species of Agaricus and Lepiota, but with a very reduced basidiome form. The nLSU rDNA of four Australian taxa, Barcheria willisiana, Agaricus xanthodermus, Leucoagaricus naucinus, and Lepiota discolorata, was sequenced for this study. Parsimony analysis of the sequences placed Barcheria within an Agaricus clade.