Evolution underground: A molecular phylogenetic investigation of Australian burrowing freshwater crayfish (Decapoda: Parastacidae) with particular focus on Engaeus Erichson

Correlated evolution of conspicuous colouration and burrowing in crayfish

Conspicuous colours have fascinated biologists for centuries, leading to research on the evolution and functional significance of colour traits. In many cases, research suggests that conspicuous colours are adaptive and serve a function in sexual or aposematic signalling. In other cases, a lack of evidence for the adaptive value of conspicuous colours garners interest from biologists, such as when organisms that live underground and are rarely exposed to the surface are nevertheless colourful. Here, we use phylogenetic comparative methods to investigate colour evolution throughout freshwater crayfishes that vary in burrowing ability. Within the taxa we analysed, conspicuous colours have evolved independently over 50 times, and these colours are more common in semi-terrestrial crayfishes that construct extensive burrows. The intuitive but not evolutionarily justified assumption when presented with these results is to assume that these colours are adaptive. But contrary to this intuition, we discuss the hypothesis that colouration in crayfish is neutral. Supporting these ideas, the small population sizes and reduced gene flow within semi-terrestrial burrowing crayfishes may lead to the fixation of colour-phenotype mutations. Overall, our work brings into question the traditional view of animal colouration as a perfectly adapted phenotype.

Living on the edge: Crayfish as drivers to anoxification of their own shelter microenvironment

Burrowing is a common trait among crayfish thought to help species deal with adverse environmental challenges. However, little is known about the microhabitat ecology of crayfish taxa in relation to their burrows. To fill this knowledge gap, we assessed the availability of oxygen inside the crayfish shelter by series of in-vivo and in-silico modelling experiments. Under modeled condition, we found that, except for the entrance region of the 200 mm, a flooded burrow microenvironment became anoxic within 8 h, on average. Multiple 12-hour day-night cycles, with burrows occupied by crayfish for 12 h and empty for 12 h, were not sufficient for refreshing the burrow microenvironment. We then examined the degree to which crayfish species with different propensities for burrowing are tolerant of self-created anoxia. From these experiments, primary and secondary burrowers showed best and most consistent tolerance-exhibiting ≥ 64% survival to anoxia and 25-91% survival of ≥ 9 h at anoxia, respectively. Tertiary burrowers exhibited little to no tolerance of anoxia with 0-50% survival to anoxia and only one species exhibiting survival (2%) of ≥ 9 h at anoxia. Results suggest that moderate to strongly burrowing crayfish can quickly draw down the dissolved oxygen in burrow water but appear to have conserved a legacy of strong tolerance of anoxia from their monophyletic ancestors-the lobsters-whereas tertiary burrowers have lost (or never evolved) this ability.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

Phylogenetic and morphological discord indicates introgressive hybridisation in two genera of Australian millipedes (Diplopoda, Polydesmida, Paradoxosomatidae)

Discord between molecular and morphological datasets was observed in two pairs of species of Australian millipedes in the family Paradoxosomatidae using morphological and molecular phylogenetic analysis (mitochondrial COI rDNA and 16 rRNA, and nuclear 28S rRNA). Close to the presumed distributional boundary between Pogonosternumnigrovirgatum (Carl, 1912) and Pogonosternumjeekeli Decker, 2017, near Dargo in Central Gippsland, Victoria, Pogonosternum specimens were collected which are phylogenetically closer to P.jeekeli in COI and 16S sequences, but are morphologically closer to P.nigrovirgatum. At Mount Osmond, Adelaide, South Australia, eight morphologically typical Somethuscastaneus (Attems, 1944) specimens were collected are phylogenetically closer to S.castaneus in 28S genealogy, but three of the eight are closer to S.lancearius Jeekel, 2002 in COI genealogy. These two cases are discussed in terms of hybridisation, past introgressive hybridisation events and aberrant morphology.

More evolution underground: Accelerated mitochondrial substitution rate in Australian burrowing freshwater crayfishes (Decapoda: Parastacidae)

To further understand the evolutionary history and mitogenomic features of Australia's highly distinctive freshwater crayfish fauna, we utilized a recently described rapid mitogenome sequencing pipeline to generate 24 new crayfish mitogenomes including a diversity of burrowing crayfish species and the first for Astacopsis gouldi, the world's largest freshwater invertebrate. Whole mitogenome-based phylogeny estimates using both Bayesian and Maximum Likelihood methods substantially strengthen existing hypotheses for systematic relationships among Australian freshwater crayfish with evidence of pervasive diversifying selection and accelerated mitochondrial substitution rate among the members of the clade representing strongly burrowing crayfish that may reflect selection pressures for increased energy requirement for adaptation to terrestrial environment and a burrowing lifestyle. Further, gene rearrangements are prevalent in the burrowing crayfish mitogenomes involving both tRNA and protein coding genes. In addition, duplicated control regions were observed in two closely related Engaeus species, together with evidence for concerted evolution. This study significantly adds to the understanding of Australian freshwater crayfish evolutionary relationships and suggests a link between mitogenome evolution and adaptation to terrestrial environments and a burrowing lifestyle in freshwater crayfish.

A novel genus and cryptic species harboured within the monotypic freshwater crayfish genus Tenuibranchiurus Riek, 1951 (Decapoda: Parastacidae)

Identifying species groups is an important yet difficult task, with there being no single accepted definition as to what constitutes a species, nor a set of criteria by which they should be delineated. Employing the General Lineage Concept somewhat circumvents these issues, as this concept allows multiple concordant lines of evidence to be used as support for species delimitation, where a species is defined as any independently evolving lineage. Genetically diverse groups have previously been identified within the monotypic parastacid genus Tenuibranchiurus Riek, 1951, but no further investigation of this diversity has previously been undertaken. Analysis of two mitochondrial DNA gene regions has previously identified two highly divergent groups within this taxon, representing populations from Queensland (Qld) and New South Wales (NSW), respectively. Additional testing within this study of both mitochondrial and nuclear DNA through species discovery analyses identified genetically diverse groups within these regions, which were further supported by lineage validation methods. The degree of genetic differentiation between Qld and NSW populations supports the recognition of two genera; with Qld retaining the original genus name Tenuibranchiurus, and NSW designated as Gen. nov. until a formal description is completed. Concordance between the species discovery and lineage validation methods supports the presence of six species within Tenuibranchiurus and two within Gen. nov. The recognition of additional species removes the monotypy of the genus, and the methods used can improve species identification within groups of organisms with taxonomic problems and cryptic diversity.

Phylogenetic analysis of the Australian trans-Bass Strait millipede genus Pogonosternum (Carl, 1912) (Diplopoda, Polydesmida, Paradoxosomatidae) indicates multiple glacial refugia in southeastern Australia

This study documents the first detailed phylogenetic analysis of an Australian paradoxosomatid millipede genus. Two mitochondrial genes (partial COI and 16S) as well as partial nuclear 28S rDNA were amplified and sequenced for 41 individuals of the southeastern Australian genus Pogonosternum Jeekel, 1965. The analysis indicates that five species groups of Pogonosternum occur across New South Wales, Victoria and Tasmania: Pogonosternumnigrovirgatum (Carl, 1912), Pogonosternumadrianae Jeekel, 1982, Pogonosternumlaetificum Jeekel, 1982 and two undescribed species. Pogonosternumconiferum (Jeekel, 1965) specimens cluster within Pogonosternumnigrovirgatum. Most of these five species groups exhibit a pattern of high intraspecific genetic variability and highly localized haplotypes, suggesting that they were confined to multiple Pleistocene refugia on the southeastern Australian mainland. The phylogenetic data also show that northwestern Tasmania was colonized by Pogonosternumnigrovirgatum, probably from central Victoria, and northeastern Tasmania by an as yet undescribed species from eastern Victoria.

A synthetic phylogeny of freshwater crayfish: insights for conservation

Phylogenetic systematics is heading for a renaissance where we shift from considering our phylogenetic estimates as a static image in a published paper and taxonomies as a hardcopy checklist to treating both the phylogenetic estimate and dynamic taxonomies as metadata for further analyses. The Open Tree of Life project (opentreeoflife.org) is developing synthesis tools for harnessing the power of phylogenetic inference and robust taxonomy to develop a synthetic tree of life. We capitalize on this approach to estimate a synthesis tree for the freshwater crayfish. The crayfish make an exceptional group to demonstrate the utility of the synthesis approach, as there recently have been a number of phylogenetic studies on the crayfishes along with a robust underlying taxonomic framework. Importantly, the crayfish have also been extensively assessed by an IUCN Red List team and therefore have accurate and up-to-date area and conservation status data available for analysis within a phylogenetic context. Here, we develop a synthesis phylogeny for the world's freshwater crayfish and examine the phylogenetic distribution of threat. We also estimate a molecular phylogeny based on all available GenBank crayfish sequences and use this tree to estimate divergence times and test for divergence rate variation. Finally, we conduct EDGE and HEDGE analyses and identify a number of species of freshwater crayfish of highest priority in conservation efforts.

Biogeography and speciation of terrestrial fauna in the south-western Australian biodiversity hotspot

The south-western land division of Western Australia (SWWA), bordering the temperate Southern and Indian Oceans, is the only global biodiversity hotspot recognised in Australia. Renowned for its extraordinary diversity of endemic plants, and for some of the largest and most botanically significant temperate heathlands and woodlands on Earth, SWWA has long fascinated biogeographers. Its flat, highly weathered topography and the apparent absence of major geographic factors usually implicated in biotic diversification have challenged attempts to explain patterns of biogeography and mechanisms of speciation in the region. Botanical studies have always been central to understanding the biodiversity values of SWWA, although surprisingly few quantitative botanical analyses have allowed for an understanding of historical biogeographic processes in both space and time. Faunistic studies, by contrast, have played little or no role in defining hotspot concepts, despite several decades of accumulating quantitative research on the phylogeny and phylogeography of multiple lineages. In this review we critically analyse datasets with explicit supporting phylogenetic data and estimates of the time since divergence for all available elements of the terrestrial fauna, and compare these datasets to those available for plants. In situ speciation has played more of a role in shaping the south-western Australian fauna than has long been supposed, and has occurred in numerous endemic lineages of freshwater fish, frogs, reptiles, snails and less-vagile arthropods. By contrast, relatively low levels of endemism are found in birds, mammals and highly dispersive insects, and in situ speciation has played a negligible role in generating local endemism in birds and mammals. Quantitative studies provide evidence for at least four mechanisms driving patterns of endemism in south-western Australian animals, including: (i) relictualism of ancient Gondwanan or Pangaean taxa in the High Rainfall Province; (ii) vicariant isolation of lineages west of the Nullarbor divide; (iii) in situ speciation; and (iv) recent population subdivision. From dated quantitative studies we derive four testable models of historical biogeography for animal taxa in SWWA, each explicit in providing a spatial, temporal and topological perspective on patterns of speciation or divergence. For each model we also propose candidate lineages that may be worthy of further study, given what we know of their taxonomy, distributions or relationships. These models formalise four of the strongest patterns seen in many animal taxa from SWWA, although other models are clearly required to explain particular, idiosyncratic patterns. Generating numerous new datasets for suites of co-occurring lineages in SWWA will help refine our understanding of the historical biogeography of the region, highlight gaps in our knowledge, and allow us to derive general postulates from quantitative (rather than qualitative) results. For animals, this process has now begun in earnest, as has the process of taxonomically documenting many of the more diverse invertebrate lineages. The latter remains central to any attempt to appreciate holistically biogeographic patterns and processes in SWWA, and molecular phylogenetic studies should - where possible - also lead to tangible taxonomic outcomes.

The complete mitogenome of the Australian crayfish Geocharax gracilis Clark 1936 (Crustacea: Decapoda: Parastacidae)

The mitogenome of the black yabby, Geocharax gracilis, was sequenced using the MiSeq Personal Sequencer. It has 15,924 base pairs consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 23 transfer RNAs, and a non-coding AT-rich region. The base composition of G. gracilis mitogenome is 32.18% for T, 22.32% for C, 34.83% for A, and 10.68% for G, with an AT bias of 67.01%. The mitogenome gene order is typical for that of parastacid crayfish with the exception of some minor rearrangements involving tRNA genes.

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

BACKGROUND AND AIMS

While molecular approaches can often accurately reconstruct species relationships, taxa that are incompletely differentiated pose a challenge even with extensive data. Such taxa are functionally differentiated, but may be genetically differentiated only at small and/or patchy regions of the genome. This issue is considered here in Poa tussock grass species that dominate grassland and herbfields in the Australian alpine zone.

METHODS

Previously reported tetraploidy was confirmed in all species by sequencing seven nuclear regions and five microsatellite markers. A Bayesian approach was used to co-estimate nuclear and chloroplast gene trees with an overall dated species tree. The resulting species tree was used to examine species structure and recent hybridization, and intertaxon fertility was tested by experimental crosses.

KEY RESULTS

Species tree estimation revealed Poa gunnii, a Tasmanian endemic species, as sister to the rest of the Australian alpine Poa. The taxa have radiated in the last 0·5-1·2 million years and the non-gunnii taxa are not supported as genetically distinct. Recent hybridization following past species divergence was also not supported. Ongoing gene flow is suggested, with some broad-scale geographic structure within the group.

CONCLUSIONS

The Australian alpine Poa species are not genetically distinct despite being distinguishable phenotypically, suggesting recent adaptive divergence with ongoing intertaxon gene flow. This highlights challenges in using conventional molecular taxonomy to infer species relationships in recent, rapid radiations.

Integrated shotgun sequencing and bioinformatics pipeline allows ultra-fast mitogenome recovery and confirms substantial gene rearrangements in Australian freshwater crayfishes

BACKGROUND

Although it is possible to recover the complete mitogenome directly from shotgun sequencing data, currently reported methods and pipelines are still relatively time consuming and costly. Using a sample of the Australian freshwater crayfish Engaeus lengana, we demonstrate that it is possible to achieve three-day turnaround time (four hours hands-on time) from tissue sample to NCBI-ready submission file through the integration of MiSeq sequencing platform, Nextera sample preparation protocol, MITObim assembly algorithm and MITOS annotation pipeline.

RESULTS

The complete mitochondrial genome of the parastacid freshwater crayfish, Engaeus lengana, was recovered by modest shotgun sequencing (1.2 giga bases) using the Illumina MiSeq benchtop sequencing platform. Genome assembly using the MITObim mitogenome assembler recovered the mitochondrial genome as a single contig with a 97-fold mean coverage (min. = 17; max. = 138). The mitogenome consists of 15,934 base pairs and contains the typical 37 mitochondrial genes and a non-coding AT-rich region. The genome arrangement is similar to the only other published parastacid mitogenome from the Australian genus Cherax.

CONCLUSIONS

We infer that the gene order arrangement found in Cherax destructor is common to Australian crayfish and may be a derived feature of the southern hemisphere family Parastacidae. Further, we report to our knowledge, the simplest and fastest protocol for the recovery and assembly of complete mitochondrial genomes using the MiSeq benchtop sequencer.

A next-generation sequencing method for overcoming the multiple gene copy problem in polyploid phylogenetics, applied to Poa grasses

BACKGROUND

Polyploidy is important from a phylogenetic perspective because of its immense past impact on evolution and its potential future impact on diversification, survival and adaptation, especially in plants. Molecular population genetics studies of polyploid organisms have been difficult because of problems in sequencing multiple-copy nuclear genes using Sanger sequencing. This paper describes a method for sequencing a barcoded mixture of targeted gene regions using next-generation sequencing methods to overcome these problems.

RESULTS

Using 64 3-bp barcodes, we successfully sequenced three chloroplast and two nuclear gene regions (each of which contained two gene copies with up to two alleles per individual) in a total of 60 individuals across 11 species of Australian Poa grasses. This method had high replicability, a low sequencing error rate (after appropriate quality control) and a low rate of missing data. Eighty-eight percent of the 320 gene/individual combinations produced sequence reads, and >80% of individuals produced sufficient reads to detect all four possible nuclear alleles of the homeologous nuclear loci with 95% probability.We applied this method to a group of sympatric Australian alpine Poa species, which we discovered to share an allopolyploid ancestor with a group of American Poa species. All markers revealed extensive allele sharing among the Australian species and so we recommend that the current taxonomy be re-examined. We also detected hypermutation in the trnH-psbA marker, suggesting it should not be used as a land plant barcode region. Some markers indicated differentiation between Tasmanian and mainland samples. Significant positive spatial genetic structure was detected at <100 km with chloroplast but not nuclear markers, which may be a result of restricted seed flow and long-distance pollen flow in this wind-pollinated group.

CONCLUSIONS

Our results demonstrate that 454 sequencing of barcoded amplicon mixtures can be used to reliably sample all alleles of homeologous loci in polyploid species and successfully investigate phylogenetic relationships among species, as well as to investigate phylogeographic hypotheses. This next-generation sequencing method is more affordable than and at least as reliable as bacterial cloning. It could be applied to any experiment involving sequencing of amplicon mixtures.