Search

There was a change in the author names in the published article. The new list should read: Kathryn E. Holt1, Johanna J. Kenyon2, Mohammad Hamidian3, Mark B. Schultz4, Derek J. Pickard5, Gordon Dougan5, Ruth M. Hall3 1,Department of Biochemistry & Molecular Biology, The University of Melbourne, Royal Parade, Parkville, Victoria, Australia 2, School of Biomedical Science, Queensland University of Technology, Brisbane, Queensland, Australia 3, School of Molecular Bioscience, The University of Sydney, Sydney, New South Wales, Australia 4, Centre for Systems Genomics, The University of Melbourne, Parkville, Victoria, Australia 5, Wellcome Sanger Trust Institute, Hinxton, Cambridge, UK The author apologizes for any inconvenience caused.

Public health agencies are increasingly relying on genomics during Legionnaires' disease investigations. However, the causative bacterium (Legionella pneumophila) has an unusual population structure, with extreme temporal and spatial genome sequence conservation. Furthermore, Legionnaires' disease outbreaks can be caused by multiple L. pneumophila genotypes in a single source. These factors can confound cluster identification using standard phylogenomic methods. Here, we show that a statistical learning approach based on L. pneumophila core genome single nucleotide polymorphism (SNP) comparisons eliminates ambiguity for defining outbreak clusters and accurately predicts exposure sources for clinical cases. We illustrate the performance of our method by genome comparisons of 234 L. pneumophila isolates obtained from patients and cooling towers in Melbourne, Australia, between 1994 and 2014. This collection included one of the largest reported Legionnaires' disease outbreaks, which involved 125 cases at an aquarium. Using only sequence data from L. pneumophila cooling tower isolates and including all core genome variation, we built a multivariate model using discriminant analysis of principal components (DAPC) to find cooling tower-specific genomic signatures and then used it to predict the origin of clinical isolates. Model assignments were 93% congruent with epidemiological data, including the aquarium Legionnaires' disease outbreak and three other unrelated outbreak investigations. We applied the same approach to a recently described investigation of Legionnaires' disease within a UK hospital and observed a model predictive ability of 86%. We have developed a promising means to breach L. pneumophila genetic diversity extremes and provide objective source attribution data for outbreak investigations.

Phylogenetic relationships and species boundaries of Australian burrowing freshwater crayfish belonging to the genera Engaeus, Engaewa, Geocharax, Gramastacus and Tenuibranchiurus are investigated using combined mitochondrial and nuclear DNA sequence data and Bayesian and Maximum Parsimony methods. Phylogenies are statistically compared to previously published hypotheses. Engaeus, Engaewa, Geocharax, Gramastacus and Tenuibranchiurus form a strongly supported monophyletic clade. This grouping is independently supported by morphology but unites geographically highly disjunct lineages. Our data show two cryptic species in Geocharax, one cryptic species in Gramastacus and two cryptic species within the highly divergent Engaeus lyelli lineage. Using a Bayesian relaxed molecular clock method, the 16S rDNA data show generic-level diversification coinciding with the transition from a wet to arid palaeoclimate near the mid Miocene

Historical sea levels have been influential in shaping the phylogeography of freshwater limited taxa via palaeodrainage and palaeoshoreline connections. In this study, we demonstrate an approach to phylogeographic analysis incorporating historical sea-level information in a nested clade phylogeographic analysis (NCPA) framework, using burrowing freshwater crayfish as the model organism. Our study area focuses on the Bass Strait region of southeastern Australia, which is marine region encompassing a shallow seabed that has emerged as a land bridge during glacial cycles connecting mainland Australia and Tasmania. Bathymetric data were analysed using Geographical Information Systems (GIS) to delineate a palaeodrainage model when the palaeocoastline was 150 m below present-day sea level. Such sea levels occurred at least twice in the past 500 000 years, perhaps more often or of larger magnitude within the last 10 million years, linking Victoria and Tasmania. Inter-locality distance measures confined to the palaeodrainage network were incorporated into an NCPA of crayfish (Engaeus sericatus Clark 1936) mitochondrial 16S rDNA haplotypes. The results were then compared to NCPAs using present-day river drainages and traditional great-circle distance measures. NCPA inferences were cross-examined using frequentist and Bayesian procedures in the context of geomorphological and historical sea-level data. We found distribution of present-day genetic variation in E. sericatus to be partly explained not only by connectivity through palaeodrainages but also via present-day drainages or overland (great circle) routes. We recommend that future studies consider all three of these distance measures, especially for studies of coastally distributed species.

Nucleotide sequence data from the mitochondrial 16S rDNA region were utilised to investigate phylogenetic relationships and species boundaries among Australian freshwater crayfish belonging to the genera Engaeus Erichson, 1846, Geocharax Clark, 1936 and Gramastacus Riek, 1972. Geocharax and Gramastacus were found to be monophyletic genera but one species currently assigned to Engaeus may belong to another genus. Relationships between the three existing genera were not resolved. Analysis of species boundaries within Geocharax suggests that there are an additional two species in this genus, and our analysis of Gramastacus indicates that undescribed populations from central New South Wales may comprise a second species. The data provide at least one instance of a taxon crossing the Great Dividing Range and provide confirmation of previously proposed hypotheses seeking to explain trans-Bass Strait distributions of species.

The increased rate at which complete mitogenomes are being sequenced and their increasing use for phylogenetic studies have resulted in a bioinformatic bottleneck in preparing and utilising such data for phylogenetic analysis. Hence, we present MitoPhAST, an automated tool that (1) identifies annotated protein-coding gene features and generates a standardised, concatenated and partitioned amino acid alignment directly from complete/partial GenBank/EMBL-format mitogenome flat files, (2) generates a maximum likelihood phylogenetic tree using optimised protein models and (3) reports various mitochondrial genes and sequence information in a table format. To demonstrate the capacity of MitoPhAST in handling a large dataset, we used 81 publicly available decapod mitogenomes, together with eight new complete mitogenomes of Australian freshwater crayfishes, including the first for the genus Gramastacus, to undertake an updated test of the monophyly of the major groups of the order Decapoda and their phylogenetic relationships. The recovered phylogenetic trees using both Bayesian and ML methods support the results of studies using fragments of mtDNA and nuclear markers and other smaller-scale studies using whole mitogenomes. In comparison to the fragment-based phylogenies, nodal support values are generally higher despite reduced taxon sampling suggesting there is value in utilising more fully mitogenomic data. Additionally, the simple table output from MitoPhAST provides an efficient summary and statistical overview of the mitogenomes under study at the gene level, allowing the identification of missing or duplicated genes and gene rearrangements. The finding of new mtDNA gene rearrangements in several genera of Australian freshwater crayfishes indicates that this group has undergone an unusually high rate of evolutionary change for this organelle compared to other major families of decapod crustaceans. As a result, freshwater crayfishes are likely to be a useful model for

Antarctica contains some of the most challenging environmental conditions on the planet due to freezing temperatures, prolonged winters and lack of liquid water. Whereas 99.7% of Antarctica is permanently covered by ice and snow, some coastal areas and mountain ridges have remained ice-free and are able to sustain populations of microinvertebrates. Tardigrades are one of the more dominant groups of microfauna in soil and limno-terrestrial habitats, but little is known of their diversity and distribution across Antarctica. Here, we examine tardigrades sampled from across an extensive region of continental Antarctica, and analyse and compare their partial mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences with those from the Antarctic Peninsula, maritime and sub-Antarctica, Tierra del Fuego and other worldwide locations in order to recognise operational taxonomic units (OTUs). From 439 new tardigrade COI sequences, we identified 98 unique haplotypes (85 from Antarctica) belonging to Acutuncus, Diphascon, Echiniscus, Macrobiotus, Milnesium and unidentified Parachela. Operational taxonomic units were delimited by Poisson tree processes and general mixed Yule coalescent methods, resulting in 58 and 55 putative species, respectively. Most tardigrades appear to be locally endemic (i.e. restricted to a single geographic region), but some (e.g. Acutuncus antarcticus (Richters, 1904)) are widespread across continental Antarctica. Our molecular results reveal: (i) greater diversity than has previously been appreciated with distinct OTUs that potentially represent undescribed species, and (ii) a lack of connectivity between most OTUs from continental Antarctica and those from other Antarctic geographical zones.

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.