The mitochondrial genome sequence of the Australian crayfish, Euastacus yarraensis, is documented and compared with other Australian crayfish genera. Euastacus yarraensis has a mitogenome of 15,548 base pairs consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a non-coding AT-rich region. The base composition of E. yarraensis mitogenome is 32.39% for T, 22.45% for C, 34.43% for A, and 10.73% for G, with an AT bias of 66.82%. The mitogenome gene order conforms to what is considered the primitive arrangement for parastacid crayfish.
The complete mitochondrial genome of the invasive house crow (Corvus splendens) was sequenced (GenBank accession number: KJ766304) using the MiSeq Personal Sequencer (Illumina, San Diego, CA). The mitochondrial genome is 16,962 bp in length, comprising 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal subunit genes and a non-coding control region. The mitogenome structural organization is identical to that of the other Corvus species and related genera. The overall base composition of C. splendens is 30.65% for A, 29.71% for C, 14.84% for G and 24.80% for T, with an AT content of 55.45%. We propose to use full mitochondrial genome to address taxonomic issues and to study the population genetics of crows.
The mitochondrial genome sequence of the stone crab, Myomenippe fornasinii, second of the superfamily Eriphioidea is documented. Myomenippe fornasinii has a mitogenome of 15,658 base pairs consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs and a non-coding AT-rich region. The base composition of the M. fornasinii mitogenome is 36.10% for T, 18.52% for C, 35.48% for A, and 9.90% for G, with an AT bias of 71.58%. The mitogenome gene order conforms to what is the standard arrangement for brachyuran crabs.
The mitochondrial genome sequence of the Morton Bay bug, Thenus orientalis, is documented, which makes it the second mitogenome for species of the family Scyllaridae and the ninth for members of the superfamily Palinuroidae. Thenus orientalis has a mitogenome of 16,826 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 the T. orientalis mitogenome is 31.31% for T, 23.77% for C, 31.05% for A, and 13.87% for G, with an AT bias of 62.36%. In addition to a duplicated trnS1 and several other tRNA gene rearrangements, the mitogenome gene order has novel protein coding gene order with the nad6 and cob genes translocated as a block to a location downstream of the nad3 gene.
The complete mitogenome of the ray Pastinachus atrus was recovered from a partial genome scan using the HiSeq sequencing system. The P. atrus mitogenome has 18,162 base pairs (61% A + T content) made up of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a 2516 bp non-coding AT-rich region. This mitogenome sequence is the first for a ray from Australian waters, the first for the Genus Pastinachus, and the 6th for the family Dasyatidae.
The complete mitochondrial genome of the swimming crab Thalamita crenata was obtained from a partial genome scan using the MiSeq sequencing system. The Thalamita crenata mitogenome has 15,787 base pairs (70% A+T content) made up of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a putative 897 bp non-coding AT-rich region. This Thalamita mitogenome sequence is the first for the genus and the eighth for the family Portunidae.
The complete mitochondrial genome of the moon crab Ashtoret lunaris was obtained from a partial genome scan using the MiSeq sequencing system. The Ashtoret lunaris mitogenome is 15,807 base pairs in length (70% A + T content) and made up of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a putative 956 bp non-coding AT-rich region. This A. lunaris mitogenome sequence is the first for the genus, as well as the family Matutidae and superfamily Calappoidea.
The complete mitochondrial genome of the commercially important snout otter clam Lutraria rhynchaena was obtained from low-coverage shotgun sequencing data on the MiSeq platform. The L. rhynchaena mitogenome has 16,927 base pairs (69% A + T content) and made up of 12 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a 953 bp non-coding AT-rich region. This is the first mitogenome to be sequenced from the genus Lutraria, and the seventh to be reported for the family Mactridae.
The complete mitochondrial genome of a highland freshwater crayfish, Cherax monticola, was recovered by shotgun sequencing. The mitogenome consists of 15,917 base pairs containing 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs and a non-coding AT-rich region. The base composition of C. monticola is 33.46% for T, 21.48% for C, 33.71% for A and 11.35% for G, with an AT bias of 67.17%.
The commercial freshwater crayfish Cherax quadricarinatus complete mitochondrial genome was recovered from partial genome sequencing using the MiSeq Personal Sequencer. The mitogenome has 15,869 base pairs consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a non-coding AT-rich region. The base composition of C. quadricarinatus is 32.16% for T, 23.39% for C, 33.26% for A, and 11.19% for G, with an AT bias of 65.42%.
The complete mitochondrial genome of the conservationally significant Macquarie perch (Macquaria australasica) was obtained from low-coverage shotgun sequencing using the MiSeq sequencer. The M. australasica mitogenome has 16,496 base pairs (55% A + T content) made up of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a 819 bp non-coding AT-rich region. This is the first mitogenome sequence for the genus Macquaria, and the third to be reported for the family Percichthyidae.
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.
The complete mitochondrial genomes of two jungle crows (Corvus macrorhynchos) were sequenced. DNA was extracted from tissue samples obtained from shed feathers collected in the field in Sri Lanka and sequenced using the Illumina MiSeq Personal Sequencer. Jungle crow mitogenomes have a structural organization typical of the genus Corvus and are 16,927 bp and 17,066 bp in length, both comprising 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal subunit genes, and a non-coding control region. In addition, we complement already available house crow (Corvus spelendens) mitogenome resources by sequencing an individual from Singapore. A phylogenetic tree constructed from Corvidae family mitogenome sequences available on GenBank is presented. We confirm the monophyly of the genus Corvus and propose to use complete mitogenome resources for further intra- and interspecies genetic studies.
The complete mitochondrial genomes of four fish species of the commercially important family Latidae were sequenced using the Illumina MiSeq, thereby significantly increasing the mitogenomic resources for the family. Whole mitogenome-based phylogenetic analysis supports the monophyly of the genus Lates and more generally the family Latidae. The mitogenome sequences from this study will be useful for future assessments of the diversity within and between Lates species and studies of phylogenetic relationships within the diverse and taxonomically challenging perciform fishes.
The marine clam Lutraria rhynchaena is gaining popularity as an aquaculture species in Asia. Lutraria populations are present in the wild throughout Vietnam and several stocks have been established and translocated for breeding and aquaculture grow-out purposes. In this study, we demonstrate the feasibility of utilising Illumina next-generation sequencing technology to streamline the identification and genotyping of microsatellite loci from this clam species. Based on an initial partial genome scan, 48 microsatellite markers with similar melting temperatures were identified and characterised. The 12 most suitable polymorphic loci were then genotyped using 51 individuals from a population in Quang Ninh Province, North Vietnam. Genetic variation was low (mean number of alleles per locus = 2.6; mean expected heterozygosity = 0.41). Two loci showed significant deviation from Hardy-Weinberg equilibrium (HWE) and the presence of null alleles, but there was no evidence of linkage disequilibrium among loci. Three additional populations were screened (n = 7-36) to test the geographic utility of the 12 loci, which revealed 100 % successful genotyping in two populations from central Vietnam (Nha Trang). However, a second population from north Vietnam (Co To) could not be successfully genotyped and morphological evidence and mitochondrial variation suggests that this population represents a cryptic species of Lutraria. Comparisons of the Qang Ninh and Nha Trang populations, excluding the 2 loci out of HWE, revealed statistically significant allelic variation at 4 loci. We reported the first microsatellite loci set for the marine clam Lutraria rhynchaena and demonstrated its potential in differentiating clam populations. Additionally, a cryptic species population of Lutraria rhynchaena was identified during initial loci development, underscoring the overlooked diversity of marine clam species in Vietnam and the need to genetically characterise population representatives prior to microsatellite development. The rapid identification and validation of microsatellite loci using next-generation sequencing technology warrant its integration into future microsatellite loci development for key aquaculture species in Vietnam and more generally, aquaculture countries in the South East Asia region.
Adaptive differences across species' ranges can have important implications for population persistence and conservation management decisions. Despite advances in genomic technologies, detecting adaptive variation in natural populations remains challenging. Key challenges in gene-environment association studies involve distinguishing the effects of drift from those of selection and identifying subtle signatures of polygenic adaptation. We used paired-end restriction site-associated DNA sequencing data (6,605 biallelic single nucleotide polymorphisms; SNPs) to examine population structure and test for signatures of adaptation across the geographic range of an iconic Australian endemic freshwater fish species, the Murray cod Maccullochella peelii. Two univariate gene-association methods identified 61 genomic regions associated with climate variation. We also tested for subtle signatures of polygenic adaptation using a multivariate method (redundancy analysis; RDA). The RDA analysis suggested that climate (temperature- and precipitation-related variables) and geography had similar magnitudes of effect in shaping the distribution of SNP genotypes across the sampled range of Murray cod. Although there was poor agreement among the candidate SNPs identified by the univariate methods, the top 5% of SNPs contributing to significant RDA axes included 67% of the SNPs identified by univariate methods. We discuss the potential implications of our findings for the management of Murray cod and other species generally, particularly in relation to informing conservation actions such as translocations to improve evolutionary resilience of natural populations. Our results highlight the value of using a combination of different approaches, including polygenic methods, when testing for signatures of adaptation in landscape genomic studies.
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 studies designed to understand the evolution of mtDNA rearrangements. We anticipate that our bioinformatics pipeline will substantially help mitogenome-based studies increase the speed, accuracy and efficiency of phylogenetic studies utilising mitogenome information. MitoPhAST is available for download at https://github.com/mht85/MitoPhAST.
The infraorder Anomura consists of a morphologically and ecologically heterogeneous group of decapod crustaceans, and has attracted interest from taxonomists for decades attempting to find some order out of the seemingly chaotic diversity within the group. Species-level diversity within the Anomura runs the gamut from the "hairy" spindly-legged yeti crab found in deep-sea hydrothermal vent environments to the largest known terrestrial invertebrate, the robust coconut or robber crab. Owing to a well-developed capacity for parallel evolution, as evidenced by the occurrence of multiple independent carcinization events, Anomura has long tested the patience and skill of both taxonomists attempting to find order, and phylogeneticists trying to establish stable hypotheses of evolutionary inter-relationships. In this study, we performed genome skimming to recover the mitogenome sequences of 12 anomuran species including the world's largest extant invertebrate, the robber crab (Birgus latro), thereby over doubling these resources for this group, together with 8 new brachyuran mitogenomes. Maximum-likelihood (ML) and Bayesian-inferred (BI) phylogenetic reconstructions based on amino acid sequences from mitogenome protein-coding genes provided strong support for the monophyly of the Anomura and Brachyura and their sister relationship, consistent with previous studies. The majority of relationships within families were supported and were largely consistent with current taxonomic classifications, whereas many relationships at higher taxonomic levels were unresolved. Nevertheless, we have strong support for a polyphyletic Paguroidea and recovered a well-supported clade of a subset of paguroids (Diogenidae + Coenobitidae) basal to all other anomurans, though this requires further testing with greater taxonomic sampling. We also introduce a new feature to the MitoPhAST bioinformatics pipeline (https://github.com/mht85/MitoPhAST) that enables the extraction of mitochondrial gene order (MGO) information directly from GenBank files and clusters groups based on common MGOs. Using this tool, we compared MGOs across the Anomura and Brachyura, identifying Anomura as a taxonomic "hot spot" with high variability in MGOs among congeneric species from multiple families while noting the broad association of highly-rearranged MGOs with several anomuran lineages inhabiting extreme niches. We also demonstrate the value of MGOs as a source of novel synapomorphies for independently reinforcing tree-based relationships and for shedding light on relationships among challenging groups such as the Aegloidea and Lomisoidea that were unresolved in phylogenetic reconstructions. Overall, this study contributes a substantial amount of new genetic material for Anomura and attempts to further resolve anomuran evolutionary relationships where possible based on a combination of sequence and MGO information. The new feature in MitoPhAST adds to the relatively limited number of bioinformatics tools available for MGO analyses, which can be utilized widely across animal groups.
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.
Aquaculture production of the Pacific white shrimp is the largest in the world for crustacean species. Crucial to the sustainable global production of this important seafood species is a fundamental understanding of the shrimp gut microbiota and its relationship to the microbial ecology of shrimp pond. This is especially true, given the recently recognized role of beneficial microbes in promoting shrimp nutrient intake and in conferring resistance against pathogens. Unfortunately, aquaculture-related microbiome studies are scarce in Southeast Asia countries despite the severe impact of early mortality syndrome outbreaks on shrimp production in the region. In this study, we employed the 16S rRNA amplicon (V3-V4 region) sequencing and amplicon sequence variants (ASV) method to investigate the microbial diversity of shrimp guts and pond water samples collected from aquaculture farms located in Malaysia and Vietnam. Substantial differences in the pond microbiota were observed between countries with the presence and absence of several taxa extending to the family level. Microbial diversity of the shrimp gut was found to be generally lower than that of the pond environments with a few ubiquitous genera representing a majority of the shrimp gut microbial diversity such as Vibrio and Photobacterium, indicating host-specific selection of microbial species. Given the high sequence conservation of the 16S rRNA gene, we assessed its veracity at distinguishing Vibrio species based on nucleotide alignment against type strain reference sequences and demonstrated the utility of ASV approach in uncovering a wider diversity of Vibrio species compared to the conventional OTU clustering approach.