Phylogenomics and the rise of the angiosperms

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5-7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

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.

Continental-scale distribution and diversity of Ceratobasidium orchid mycorrhizal fungi in Australia

BACKGROUND AND AIMS

Mycorrhizal fungi are a critical component of the ecological niche of most plants and can potentially constrain their geographical range. Unlike other types of mycorrhizal fungi, the distributions of orchid mycorrhizal fungi (OMF) at large spatial scales are not well understood. Here, we investigate the distribution and diversity of Ceratobasidium OMF in orchids and soils across the Australian continent.

METHODS

We sampled 217 Ceratobasidium isolates from 111 orchid species across southern Australia and combined these with 311 Ceratobasidium sequences from GenBank. To estimate the taxonomic diversity of Ceratobasidium associating with orchids, phylogenetic analysis of the ITS sequence locus was undertaken. Sequence data from the continent-wide Australian Microbiome Initiative were used to determine the geographical range of operational taxonomic units (OTUs) detected in orchids, with the distribution and climatic correlates of the two most frequently detected OTUs modelled using MaxEnt.

KEY RESULTS

We identified 23 Ceratobasidium OTUs associating with Australian orchids, primarily from the orchid genera Pterostylis, Prasophyllum, Rhizanthella and Sarcochilus. OTUs isolated from orchids were closely related to, but distinct from, known pathogenic fungi. Data from soils and orchids revealed that ten of these OTUs occur on both east and west sides of the continent, while 13 OTUs were recorded at three locations or fewer. MaxEnt models suggested that the distributions of two widespread OTUs are correlated with temperature and soil moisture of the wettest quarter and far exceeded the distributions of their host orchid species.

CONCLUSIONS

Ceratobasidium OMF with cross-continental distributions are common in Australian soils and frequently have geographical ranges that exceed that of their host orchid species, suggesting these fungi are not limiting the distributions of their host orchids at large spatial scales. Most OTUs were distributed within southern Australia, although several OTUs had distributions extending into central and northern parts of the continent, illustrating their tolerance of an extraordinarily wide range of environmental conditions.

Continental-scale metagenomics, BLAST searches, and herbarium specimens: The Australian Microbiome Initiative and the National Herbarium of Victoria

PREMISE

Motivated to make sensible interpretations of the massive volume of data from the Australian Microbiome Initiative (AusMic), we characterize the soil mycota of Australia. We establish operational taxonomic units (OTUs) from the data and compare these to GenBank and a data set from the National Herbarium of Victoria (MEL), Melbourne, Australia. We also provide visualizations of Agaricomycete diversity, drawn from our analyses of the AusMic sequences and taxonomy.

METHODS

The AusMic internal transcribed spacer (ITS) data were filtered to create OTUs, which were searched against the National Center for Biotechnology Information Nucleotide database and the MEL database. We further characterized a portion of our OTUs by graphing the counts of the families and orders of Agaricomycetes. We also graphed AusMic species determinations for Australian Agaricomycetes against latitude.

RESULTS

Our filtering process generated 192,325 OTUs; for Agaricomycetes, there were 27,730 OTUs. Based on the existing AusMic taxonomy at species level, we inferred the diversity of Australian Agaricomycetes against latitude to be lowest between -20 and -25 decimal degrees.

DISCUSSION

BLAST comparisons provided reciprocal insights between the three data sets, including the detection of unusual root-associated species in the AusMic data, insights into mushroom morphology from the MEL data, and points of comparison for the taxonomic determinations between AusMic, GenBank, and MEL. This study provides a tabulation of Australian fungi, different visual snapshots of a subset of those taxa, and a springboard for future studies.

Characterization of the complete plastid genome of Astelia australiana (J. H. Willis) L. B. Moore (Asteliaceae, Asparagales)

Astelia australiana is a robust understorey plant with a highly restricted distribution in southeastern Australia. Here we report its complete plastid genome. The genome was 157,943 bp in length and comprises a pair of inverted repeats (IRs) of 27,028 bp separated by a large single-copy region (LSC) of 85,699 bp, and a small single-copy region (SSC) of 18,188 bp. The GC content was 37.7%. In total, 132 genes were annotated including 81 protein-coding genes (PCGs), 38 tRNA genes, and 8 rRNA genes. Phylogenetic analysis of the PCGs from A. australiana aligned with those from 10 Asparagales representatives confirms that, based on these taxa, A. australiana is sister to A. pumila and sits within the Asteliaceae.

Testing efficacy of distance and tree-based methods for DNA barcoding of grasses (Poaceae tribe Poeae) in Australia

In Australia, Poaceae tribe Poeae are represented by 19 genera and 99 species, including economically and environmentally important native and introduced pasture grasses [e.g. Poa (Tussock-grasses) and Lolium (Ryegrasses)]. We used this tribe, which are well characterised in regards to morphological diversity and evolutionary relationships, to test the efficacy of DNA barcoding methods. A reference library was generated that included 93.9% of species in Australia (408 individuals, [Formula: see text] = 3.7 individuals per species). Molecular data were generated for official plant barcoding markers (rbcL, matK) and the nuclear ribosomal internal transcribed spacer (ITS) region. We investigated accuracy of specimen identifications using distance- (nearest neighbour, best-close match, and threshold identification) and tree-based (maximum likelihood, Bayesian inference) methods and applied species discovery methods (automatic barcode gap discovery, Poisson tree processes) based on molecular data to assess congruence with recognised species. Across all methods, success rate for specimen identification of genera was high (87.5-99.5%) and of species was low (25.6-44.6%). Distance- and tree-based methods were equally ineffective in providing accurate identifications for specimens to species rank (26.1-44.6% and 25.6-31.3%, respectively). The ITS marker achieved the highest success rate for specimen identification at both generic and species ranks across the majority of methods. For distance-based analyses the best-close match method provided the greatest accuracy for identification of individuals with a high percentage of "correct" (97.6%) and a low percentage of "incorrect" (0.3%) generic identifications, based on the ITS marker. For tribe Poeae, and likely for other grass lineages, sequence data in the standard DNA barcode markers are not variable enough for accurate identification of specimens to species rank. For recently diverged grass species similar challenges are encountered in the application of genetic and morphological data to species delimitations, with taxonomic signal limited by extensive infra-specific variation and shared polymorphisms among species in both data types.

Using Transcriptomics to Identify Differential Gene Expression in Response to Salinity among Australian Phragmites australis Clones

Common Reed (Phragmites australis) is a frequent component of inland and coastal wetlands in temperate zones worldwide. Ongoing environmental changes have resulted in the decline of this species in many areas and invasive expansion in others. In the Gippsland Lakes coastal waterway system in south-eastern Australia, increasing salinity is thought to have contributed to the loss of fringing P. australis reed beds leading to increased shoreline erosion. A major goal of restoration in this waterway is to address the effect of salinity by planting a genetically diverse range of salt-tolerant P. australis plants. This has prompted an interest in examining the variation in salinity tolerance among clones and the underlying basis of this variation. Transcriptomics is an approach for identifying variation in genes and their expression levels associated with the exposure of plants to environmental stressors. In this paper we present initial results of the first comparative culm transcriptome analysis of P. australis clones. After sampling plants from sites of varied surface water salinity across the Gippsland Lakes, replicates from three clones from highly saline sites (>18 g L(-1) TDS) and three from low salinity sites (<6 g L(-1)) were grown in containers irrigated with either fresh (<0.1 g L(-1)) or saline water (16 g L(-1)). An RNA-Seq protocol was used to generate sequence data from culm tissues from the 12 samples allowing an analysis of differential gene expression. Among the key findings, we identified several genes uniquely up- or down-regulated in clones from highly saline sites when irrigated with saline water relative to clones from low salinity sites. These included the higher relative expression levels of genes associated with photosynthesis and lignan biosynthesis indicative of a greater ability of these clones to maintain growth under saline conditions. Combined with growth data from a parallel study, our data suggests local adaptation of certain clones to salinity and provides a basis for more detailed studies.

Major clades of Australasian Rutoideae (Rutaceae) based on rbcL and atpB sequences

Background

Rutaceae subfamily Rutoideae (46 genera, c. 660 species) is diverse in both rainforests and sclerophyll vegetation of Australasia. Australia and New Caledonia are centres of endemism with a number of genera and species distributed disjunctly between the two regions. Our aim was to generate a high-level molecular phylogeny for the Australasian Rutoideae and identify major clades as a framework for assessing morphological and biogeographic patterns and taxonomy.

Methodology/Principal Findings

Phylogenetic analyses were based on chloroplast genes, rbcL and atpB, for 108 samples (78 new here), including 38 of 46 Australasian genera. Results were integrated with those from other molecular studies to produce a supertree for Rutaceae worldwide, including 115 of 154 genera. Australasian clades are poorly matched with existing tribal classifications, and genera Philotheca and Boronia are not monophyletic. Major sclerophyll lineages in Australia belong to two separate clades, each with an early divergence between rainforest and sclerophyll taxa. Dehiscent fruits with seeds ejected at maturity (often associated with myrmecochory) are inferred as ancestral; derived states include woody capsules with winged seeds, samaras, fleshy drupes, and retention and display of seeds in dehisced fruits (the last two states adaptations to bird dispersal, with multiple origins among rainforest genera). Patterns of relationship and levels of sequence divergence in some taxa, mostly species, with bird-dispersed (Acronychia, Sarcomelicope, Halfordia and Melicope) or winged (Flindersia) seeds are consistent with recent long-distance dispersal between Australia and New Caledonia. Other deeper Australian/New Caledonian divergences, some involving ant-dispersed taxa (e.g., Neoschmidia), suggest older vicariance.

Conclusions/Significance

This comprehensive molecular phylogeny of the Australasian Rutoideae gives a broad overview of the group’s evolutionary and biogeographic history. Deficiencies of infrafamilial classifications of Rutoideae have long been recognised, and our results provide a basis for taxonomic revision and a necessary framework for more focused studies of genera and species.

Limitations to reproductive output and genetic rescue in populations of the rare shrub Grevillea repens (Proteaceae)

BACKGROUND AND AIMS

When conserving rare plant species, managers are often faced with small and/or isolated populations displaying low levels of sexual reproduction and genetic variation. One option for reinvigorating these populations is the introduction of genetic material from other sites, but in some cases fitness may be reduced as a result of outbreeding depression. Here the pollination biology of the rare shrub Grevillea repens is studied across its natural range and reproductive responses following cross-pollination among populations are examined to determine factors that may be limiting sexual reproduction and the potential for genetic rescue.

METHODS

Pollen manipulation treatments (self-, autogamous self-, cross- and open pollination) were applied to flowers to examine the breeding system and fruit and seed production in five populations of G. repens. Pollen production, presentation and viability were investigated and interpopulation crosses of increasing genetic distance performed among the populations.

KEY RESULTS

The study species is self-incompatible and displayed very low natural seed set over two seasons, due partly to low pollen viability in one of the populations. Within-population crossing increased fruit and seed production at some sites, indicating pollinator limitation. Interpopulation crosses further increased reproductive output in one population, suggesting mate limitation, and for this site there was a positive relationship between genetic distance among populations and the size of genetic rescue benefits. However, in other populations there was a decrease in fruit and seed set with increasing genetic distance.

CONCLUSIONS

The results highlight that management strategies involving interpopulation crosses can improve reproductive output in small, isolated populations of rare plants, but guidelines need to be developed on a population by population basis.