Limited genetic divergence among Australian alpine Poa tussock grasses coupled with regional structuring points to ongoing gene flow and taxonomic challenges

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.

Breeding system diversification and evolution in American Poa supersect. Homalopoa (Poaceae: Poeae: Poinae)

BACKGROUND AND AIMS

Poa subgenus Poa supersect. Homalopoa has diversified extensively in the Americas. Over half of the species in the supersection are diclinous; most of these are from the New World, while a few are from South-East Asia. Diclinism in Homalopoa can be divided into three main types: gynomonoecism, gynodioecism and dioecism. Here the sampling of species of New World Homalopoa is expanded to date its origin and diversification in North and South America and examine the evolution and origin of the breeding system diversity.

METHODS

A total of 124 specimens were included in the matrix, of which 89 are species of Poa supersect. Homalopoa sections Acutifoliae, Anthochloa, Brizoides, Dasypoa, Dioicopoa, Dissanthelium, Homalopoa sensu lato (s.l.), Madropoa and Tovarochloa, and the informal Punapoa group. Bayesian and parsimony analyses were conducted on the data sets based on four markers: the nuclear ribosomal internal tanscribed spacer (ITS) and external transcribed spacer (ETS), and plastid trnT-L and trnL-F. Dating analyses were performed on a reduced Poa matrix and enlarged Poaceae outgroup to utilize fossils as calibration points. A relaxed Bayesian molecular clock method was used.

KEY RESULTS

Hermaphroditism appears to be pleisiomorphic in the monophyletic Poa supersect. Homalopoa, which is suggested to have originated in Eurasia 8·4-4·2 million years ago (Mya). The ancestor of Poa supersect. Homalopoa radiated throughout the New World in the Late Miocene-Early Pliocene, with major lineages originating during the Pliocene to Pleistocene (5-2 Mya). Breeding systems are linked to geographic areas, showing an evolutionary pattern associated with different habitats. At least three major pathways from hermaphroditism to diclinism are inferred in New World Homalopoa: two leading to dioecism, one via gynodioecism in South America and another directly from hermaphroditism in North America, a result that needs to be checked with a broader sampling of diclinous species in North America. A third pathway leads from hermaphroditism to gynomonoecism in Andean species of South America, with strictly pistillate species evolving in the highest altitudes.

CONCLUSIONS

Divergence dating provides a temporal context to the evolution of breeding systems in New World Poa supersect. Homalopoa The results are consistent with the infrageneric classification in part; monophyletic sections are confirmed, it is proposed to reclassify species of sect. Acutifoliae, Dasypoa and Homalopoa s.l. and it is acknowledged that revision of the infrageneric taxonomy of the gynomonoecious species is needed.

Biological invasions, climate change and genomics

The rate of biological invasions is expected to increase as the effects of climate change on biological communities become widespread. Climate change enhances habitat disturbance which facilitates the establishment of invasive species, which in turn provides opportunities for hybridization and introgression. These effects influence local biodiversity that can be tracked through genetic and genomic approaches. Metabarcoding and metagenomic approaches provide a way of monitoring some types of communities under climate change for the appearance of invasives. Introgression and hybridization can be followed by the analysis of entire genomes so that rapidly changing areas of the genome are identified and instances of genetic pollution monitored. Genomic markers enable accurate tracking of invasive species' geographic origin well beyond what was previously possible. New genomic tools are promoting fresh insights into classic questions about invading organisms under climate change, such as the role of genetic variation, local adaptation and climate pre-adaptation in successful invasions. These tools are providing managers with often more effective means to identify potential threats, improve surveillance and assess impacts on communities. We provide a framework for the application of genomic techniques within a management context and also indicate some important limitations in what can be achieved.

Strong genetic structure corresponds to small-scale geographic breaks in the Australian alpine grasshopper Kosciuscola tristis

Background

Mountain landscapes are topographically complex, creating discontinuous ‘islands’ of alpine and sub-alpine habitat with a dynamic history. Changing climatic conditions drive their expansion and contraction, leaving signatures on the genetic structure of their flora and fauna. Australia’s high country covers a small, highly fragmented area. Although the area is thought to have experienced periods of relative continuity during Pleistocene glacial periods, small-scale studies suggest deep lineage divergence across low-elevation gaps. Using both DNA sequence data and microsatellite markers, we tested the hypothesis that genetic partitioning reflects observable geographic structuring across Australia’s mainland high country, in the widespread alpine grasshopper Kosciuscola tristis (Sjösted).

Results

We found broadly congruent patterns of regional structure between the DNA sequence and microsatellite datasets, corresponding to strong divergence among isolated mountain regions. Small and isolated mountains in the south of the range were particularly distinct, with well-supported divergence corresponding to climate cycles during the late Pliocene and Pleistocene. We found mixed support, however, for divergence among other mountain regions. Interestingly, within areas of largely contiguous alpine and sub-alpine habitat around Mt Kosciuszko, microsatellite data suggested significant population structure, accompanied by a strong signature of isolation-by-distance.

Conclusions

Consistent patterns of strong lineage divergence among different molecular datasets indicate genetic breaks between populations inhabiting geographically distinct mountain regions. Three primary phylogeographic groups were evident in the highly fragmented Victorian high country, while within-region structure detected with microsatellites may reflect more recent population isolation. Despite the small area of Australia’s alpine and sub-alpine habitats, their low topographic relief and lack of extensive glaciation, divergence among populations was on the same scale as that detected in much more extensive Northern hemisphere mountain systems. The processes driving divergence in the Australian mountains might therefore differ from their Northern hemisphere counterparts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0204-1) contains supplementary material, which is available to authorized users.