MitoPhAST, a new automated mitogenomic phylogeny tool in the post-genomic era with a case study of 89 decapod mitogenomes including eight new freshwater crayfish mitogenomes

Genome of the early spider-orchid Ophrys sphegodes provides insights into sexual deception and pollinator adaptation

Pollinator-driven evolution of floral traits is thought to be a major driver of angiosperm speciation and diversification. Ophrys orchids mimic female insects to lure male pollinators into pseudocopulation. This strategy, called sexual deception, is species-specific, thereby providing strong premating reproductive isolation. Identifying the genomic architecture underlying pollinator adaptation and speciation may shed light on the mechanisms of angiosperm diversification. Here, we report the 5.2 Gb chromosome-scale genome sequence of Ophrys sphegodes. We find evidence for transposable element expansion that preceded the radiation of the O. sphegodes group, and for gene duplication having contributed to the evolution of chemical mimicry. We report a highly differentiated genomic candidate region for pollinator-mediated evolution on chromosome 2. The Ophrys genome will prove useful for investigations into the repeated evolution of sexual deception, pollinator adaptation and the genomic architectures that facilitate evolutionary radiations.

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

The mitochondrial genome of the mountain wooly tapir, Tapirus pinchaque and a formal test of the effect of altitude on the adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates

BACKGROUND

The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments.

RESULTS

The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4.

CONCLUSION

The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments.

Comparative mitochondrial genomics of American nectar-feeding long-nosed bats Leptonycteris spp. with insights into the phylogeny of the family Phyllostomidae

Among leaf-nosed bats (family Phyllostomidae), the genus Leptonycteris (subfamily Glossophaginae), contains three migratory and obligate nectar-dwelling species of great ecological and economic importance; the greater long-nosed bat L. nivalis, the lesser long-nosed bat L. yerbabuenae, and the southern long-nosed bat L. curasoae. According to the IUCN, the three species are categorized as 'vulnerable', 'endangered', and 'near threatened', respectively. In this study, we assembled and characterized in detail the mitochondrial genome of Leptonycteris spp. and examined the phylogenetic position of this genus in the family Phyllostomidae based on protein coding genes (PCGs). The mitogenomes of L. nivalis, L. curasoae, and L. yerbabuenae are 16,708, 16,758, and 16,729 bp in length and each encode 13 PCGs, 22 transfer RNA genes, 2 rRNA genes, and a putative control region (CR). Mitochondrial gene order is identical to that reported before for the family Phyllostomidae. All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that is missing the DHU arm in the three species. All PCGs are exposed to purifying selection with atp8 experiencing the most relaxed purifying selection as the ω ratio was higher for this gene compared to the other PCGs in each species. The CR of each species contains three functional domains: extended termination associated sequence (ETAS), Central, and conserved sequence block domain (CSB). A phylomitogenomic analysis revealed that Leptonycteris is monophyletic and most closely related to the genus Glossophaga. The analysis also supported the monophyly of the family Glossophaginae in the speciose family Phyllostomidae. The mitochondria characterization of these species provides relevant information to develop molecular markers for conservation purposes.

Comparative analysis of mitochondrial genomes reveals family-specific architectures and molecular features in scorpions (Arthropoda: Arachnida: Scorpiones)

Scorpions are a group of arachnids with great evolutionary success that comprise more than 2,000 described species. Mitochondrial genomes have been little studied in this clade. We describe and compare different scorpion mitochondrial genomes and analyze their architecture and molecular characteristics. We assembled eight new scorpion mitochondrial genomes from transcriptomic datasets, annotated them, predicted the secondary structures of tRNAs, and compared the nucleotide composition, codon usage, and relative synonymous codon usage of 16 complete scorpion mitochondrial genomes. Lastly, we provided a phylogeny based on all mitochondrial protein coding genes. We characterized the mitogenomes in detail and reported particularities such as dissimilar synteny in the family Buthidae compared to other scorpions, unusual tRNA secondary structures, and unconventional start and stop codons in all scorpions. Our comparative analysis revealed that scorpion mitochondrial genomes exhibit different architectures and features depending on taxonomic identity. We highlight the parvorder Buthida, particularly the family Buthidae, as it invariably exhibited different mitogenome features such as synteny, codon usage, and AT-skew compared to the parvorder Iurida that included the rest of the scorpion families we analyzed in this study. Our results provide a better understanding of the evolution of mitogenome features and phylogenetic relationships in scorpions.

Mitochondrial genomes assembled from non-invasive eDNA metagenomic scat samples in the endangered Amur tiger Panthera tigris altaica

The Amur or Siberian tiger Panthera tigris altaica (Temminck, 1844) is currently restricted to a small region of its original geographical range in northwestern Asia and is considered 'endangered' by the IUCN Red List of Threatened Species. This solitary, territorial, and large top predator is in major need of genomic resources to inform conservation management strategies. This study formally tested if complete mitochondrial genomes of P. tigris altaica can be assembled from non-enriched metagenomic libraries generated from scat eDNA samples using the Illumina sequencing platform and open-access bioinformatics pipelines. The mitogenome of P. tigris altaica was assembled and circularized using the pipeline GetOrganelle with a coverage ranging from 322.7x to 17.6x in four different scat eDNA samples. A nearly complete mitochondrial genome (101x) was retrieved from a fifth scat eDNA sample. The complete or nearly complete mitochondrial genomes of P. tigris altaica were AT-rich and composed of 13 protein coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a putative control region. Synteny observed in all assembled mitogenomes was identical to that reported before for P. tigris altaica and other felids. A phylogenomic analysis based on all PCGs demonstrated that the mitochondrial genomes assembled from scat eDNA reliably identify the sequenced samples as belonging to P. tigris and distinguished the same samples from closely and distantly related congeneric species. This study demonstrates that it is viable to retrieve accurate whole and nearly complete mitochondrial genomes of P. tigris altaica (and probably other felids) from scat eDNA samples without library enrichment protocols and using open-access bioinformatics workflows. This new genomic resource represents a new tool to support conservation strategies (bio-prospecting and bio-monitoring) in this iconic cat.

The complete mitochondrial genome of the roosterfish Nematistius pectoralis Gill 1862: purifying selection in protein coding genes, organization of the control region, and insights into family-level phylogenomic relationships in the recently erected order Carangiformes

The roosterfish Nematistius pectoralis is considered as one of the most magnificent sportfishes worldwide. This study developed the first genomic resource for this trophy-fish that is heavily targeted by the fly-fishing industry. The complete mitochondrial genome of N. pectoralis was assembled using short read sequences and analyzed in detail. The mitochondrial genome of N. pectoralis is 16,537 bp in length and comprises 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (12S and 16S), and 22 transfer RNA genes. A long intergenic space 770 bp in length was assumed to be the D-loop or Control Region (CR). Most of the PCGs and tRNA genes are encoded in the L-strand. All PCGs are under purifying selection and atp8 and nad6 experienced the least selective pressure. All tRNAs exhibit a cloverleaf secondary structure except tRNA-Serine 1 that lacked the D-arm loop. The D-loop of N. pectoralis exhibits three domains commonly described in other fishes; extended terminal associated sequences (ETAS), central, and conserved sequence block (CSB). A ML phylogenetic reconstruction of the newly recognized order Carangiformes based on all 13 mitochondrial PCGs did not support the monophyly of this clade but recognized several families as monophyletic, including Bothidae, Carangidae, Istiophoridae, Latidae, Paralichthyidae, Polynemidae, and Rhombosoleidae. Nematistius pectoralis was sister to a clade composed of Toxotes chatareus (fam. Toxotidae) + Lactarius lactarius (fam. Lactariidae). This genomic resource developed for N. pectoralis will aid in improving our understanding of the population genomics of and strengthen conservation and management strategies in this remarkable trophy-fish.

First genomic resource for an endangered neotropical mega-herbivore: the complete mitochondrial genome of the forest-dweller (Baird's) tapir (Tapirus bairdii)

Baird's tapir, or the Central American Tapir Tapirus bairdii (family Tapiridae), is one of the largest mammals native to the forests and wetlands of southern North America and Central America, and is categorized as 'endangered' on the 2014 IUCN Red List of Threatened Species. This study reports, for the first time, the complete mitochondrial genome of T. bairdii and examines the phylogenetic position of T. bairdii amongst closely related species in the same family and order to which it belongs using mitochondrial protein-coding genes (PCG's). The circular, double-stranded, A-T rich mitochondrial genome of T. bairdii is 16,697 bp in length consisting of 13 protein-coding genes (PCG's), two ribosomal RNA genes (rrnS (12s ribosomal RNA and rrnL (16s ribosomal RNA)), and 22 transfer RNA (tRNA) genes. A 33 bp long region was identified to be the origin of replication for the light strand (O_L), and a 1,247 bp long control region (CR) contains the origin of replication for the heavy strand (O_H). A majority of the PCG's and tRNA genes are encoded on the positive, or heavy, strand. The gene order in T. baiirdi is identical to that of T. indicus and T. terrestris, the only two other species of extant tapirs with assembled mitochondrial genomes. An analysis of Ka/Ks ratios for all the PCG's show values <1, suggesting that all these PCGs experience strong purifying selection. A maximum-likelihood phylogenetic analysis supports the monophyly of the genus Tapirus and the order Perissodactyla. The complete annotation and analysis of the mitochondrial genome of T. bairdii will contribute to a better understanding of the population genomic diversity and structure of this species, and it will assist in the conservation and protection of its dwindling populations.

The mitochondrial genome of Grapsus albolineatus (Decapoda: Brachyura: Grapsidae) and phylogenetic associations in Brachyura

Complete mitochondrial genomes (mitogenomes) can provide useful information for phylogenetic relationships, gene rearrangement, and evolutionary traits. In this study, we determined the complete mitochondrial DNA sequence of the herbivorous crab Grapsus albolineatus. It is a typical metazoan mitochondrial genome. The total size is 15,583 bp, contains the entire set of 37 genes, and has an AT-rich region. Then, 23 of the 37 genes were encoded by the heavy (+) strand while 14 are encoded by the light (−) strand. Compared with the pan-crustacean ground pattern, two tRNA genes (tRNA-His and tRNA-Gln) were rearranged and the tandem duplication/random loss model was used to explain the observed gene rearrangements. The phylogenetic results showed that all Grapsidae crabs clustered together as a group. Furthermore, the monophyly of each family was well supported, with the exception of Menippidae. In general, the results obtained in this study will contribute to the better understanding of gene rearrangements in Grapsidae crab mitogenomes and provide new insights into the phylogeny of Brachyura.

Characterization of The Complete Mitochondrial Genome of the Bromeliad Crab Metopaulias depressus (Rathbun, 1896) (Crustacea: Decapoda: Brachyura: Sesarmidae)

Metopaulias depressus is a non-marine crab endemic to Jamaica that dwells in rainforest bromeliads and exhibits elaborate active parental care behavior. Current genomic resources on M. depressus are rare, limiting the understanding of its adaptation to terrestrial life in species that evolved from marine ancestors. This study reports the complete mitochondrial genome of M. depressus assembled using Sanger sequencing. The AT-rich mitochondrial genome of M. depressus is 15,765 bp in length and comprises 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 691 bp-long intergenic space is assumed to be the control region (CR) or D-loop. A set of selective pressure analyses indicate that the entirety of the PCGs experience purifying selection. Cox1, cox2, nad5, cox3, and atp6 experience strong purifying selection, and atp8 experiences weak purifying selection compared to the rest of the PCGs. The secondary structures of most tRNA genes exhibit a standard ‘cloverleaf’ structure, with the exception of trnS1, which lacks the dihydroxyuridine (DHU) arm but not the loop, the trnH gene, which lacks the thymine pseudouracil cytosine (T) loop but not the arm, and trnM, which exhibits an overly developed T loop. A maximum likelihood phylogenetic analysis based on all PCGs indicated that M. depressus is more closely related to the genera Clistocoeloma, Nanosesarma, and Parasesarma than to Chiromantes, Geosesarma, and Orisarma. This study contributes to deciphering the phylogenetic relationships within the family Sesarmidae and represents a new genomic resource for this iconic crab species.

The complete mitochondrial genome of the 'Zacatuche' Volcano rabbit (Romerolagus diazi), an endemic and endangered species from the Volcanic Belt of Central Mexico

BACKGROUND

The 'Zacatuche', 'Teporingo', or Volcano rabbit (Romerolagus diazi) belongs to the family Leporidae, is an endemic species restricted to the Central part of the Trans-Mexican Volcanic Belt, and is considered 'endangered' by the IUCN Red List of Threatened Species.

METHODS AND RESULTS

This study reports, for the first time, the complete mitochondrial genome of R. diazi and examined the phylogenetic position of R. diazi among other closely related co-familiar species using mitochondrial protein-coding genes (PCGs). The mitogenome of R. diazi was assembled from short Illumina 150 bp pair-end reads with a coverage of 189x. The AT-rich mitochondrial genome of R. diazi is 17,400 bp in length and is comprised of 13 PCGs, two ribosomal RNA genes, and 22 transfer RNA genes. The gene order observed in the mitochondrial genome of R. diazi is identical to that reported for other leporids. Phylogenetic analyses based on PCGs support the basal position of Romerolagus within the Leporidae, at least when compared to the genera Oryctolagus and Lepus. Nonetheless, additional mitochondrial genomes from species belonging to the genera Bunolagus, Sylvilagus, and Pronolagus, among others, are needed before a more robust conclusion about the derived vs basal placement of Romerolagus within the family Leporidae can be reached based on mitochondrial PCGs.

CONCLUSIONS

This is the first genomic resource developed for R. diazi and it represents a tool to improve our understanding about the ecology and evolutionary biology of this iconic and endangered species.

Evolutionary History of Inversions in Directional Mutational Pressures in Crustacean Mitochondrial Genomes: Implications for Evolutionary Studies

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.

Comparative mitochondrial genomics of sponge-dwelling snapping shrimps in the genus Synalpheus: Exploring differences between eusocial and non-eusocial species and insights into phylogenetic relationships in caridean shrimps

The genus Synalpheus is a cosmopolitan clade of marine shrimps found in most tropical regions. Species in this genus exhibit a range of social organizations, including pair-forming, communal breeding, and eusociality, the latter only known to have evolved within this genus in the marine realm. This study examines the complete mitochondrial genomes of seven species of Synalpheus and explores differences between eusocial and non-eusocial species considering that eusociality has been shown before to affect the strength of purifying selection in mitochondrial protein coding genes. The AT-rich mitochondrial genomes of Synalpheus range from 15,421 bp to 15,782 bp in length and comprise, invariably, 13 protein-coding genes (PCGs), two ribosomal RNA genes, and 22 transfer RNA genes. A 648 bp to 994 bp long intergenic space is assumed to be the D-loop. Mitochondrial gene synteny is identical among the studied shrimps. No major differences occur between eusocial and non-eusocial species in nucleotide composition and codon usage profiles of PCGs and in the secondary structure of tRNA genes. Maximum likelihood phylogenetic analysis of the complete concatenated PCG complement of 90 species supports the monophyly of the genus Synalpheus and its family Alpheidae. Moreover, the monophyletic status of the caridean families Alvinocaridae, Atyidae, Thoridae, Lysmatidae, Palaemonidae, and Pandalidae within caridean shrimps are fully or highly supported by the analysis. We therefore conclude that mitochondrial genomes contain sufficient phylogenetic information to resolve relationships at high taxonomic levels within the Caridea. Our analysis of mitochondrial genomes in the genus Synalpheus contributes to the understanding of the coevolution between genomic architecture and sociality in caridean shrimps and other marine organisms.

Insights into the evolution of Brachyura (Crustacea: Decapoda) from mitochondrial sequences and gene order rearrangements

Brachyura is one of the most species rich and highly derived groups among extant crustaceans, with over 7250 known species. However, brachyuran phylogeny remains controversial and requires further study. Here, we combined 103 brachyuran mitogenomes from GenBank with 10 new mitogenomes to describe gene rearrangement patterns and explore the internal phylogenetic relationships of Brachyura. Most of the 10 novel mitogenomes had the typical 37 genes, except that of Longpotamon depressum, which lacked trnQ. We discovered 15 gene rearrangement patterns among Brachyura and preliminarily determined their rearrangement mechanisms with the help of CREx. We identified seven putative ancestral family gene orders among the 15 rearrangement patterns and expounded systematically upon the mechanisms of their rearrangement. In our phylogenetic analysis, Raninoida shared a sister relationship with an eubrachyuran clade ((Heterotremata [Potamoidea] + Thoracotremata) + Heterotremata) at maximum nodal support rather than Dromiacea, which did not support monophyly of Podotremata. In addition, Potamoidea (Parathelphusidae + Potamidae) retained a close relationship with Thoracotremata rather than their marine relatives in Heterotremata. Our study provides important information for the evolution of Brachyura by using the large taxon sampling currently available for systematic rearrangement and phylogenetic analyses.

The complete mitochondrial genome of the red-jointed brackish-water fiddler crab Minuca minax (LeConte 1855) (Brachyura: Ocypodidae): New family gene order, and purifying selection and phylogenetic informativeness of protein coding genes

Minuca minax is a semi-terrestrial crustacean that commonly lives in low salinity, riverine habitats along the shores of the eastern United States. This study reports, for the first time, the complete mitochondrial genome of M. minax. The AT-rich mitochondrial genome of M. minax is 15,937 bp in length and comprised of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 737 bp long intergenic space is assumed to be the D-loop. Most of the PCGs and tRNA genes are encoded in the L-strand. The gene order observed in the mitochondrial genome of M. minax is new although almost identical to that reported in confamiliar species. In all other confamiliar species to which M. minax is compared, the positions of the trnQ gene and the trnI gene are switched. K_A/K_S ratios calculated for all mitochondrial PCGs show values of <1, indicating that these PCGs are evolving under purifying selection. A maximum likelihood phylogenetic analysis (concatenated PCGs [n = 13], 15 species) supports the monophyly of the subfamilies Ocypodinae and Gelaminidae. Mitochondrial PCGs have enough phylogenetic information to reveal relationships supporting higher taxonomic levels within this family. The knowledge of a complete mitochondrial genome from the red-jointed brackish-water fiddler crab M. minax contributes to the better understanding of meta-population connectivity and the mechanisms involved in the adaptation of marine organisms to near-limnic conditions.

Mitogenome phylogenetics in the genus Palaemon (Crustacea: Decapoda) sheds light on species crypticism in the rockpool shrimp P. elegans

The genus Palaemon comprises worldwide marine and freshwater shrimps and prawns, and some of them are ecologically or commercially important species. Palaemon is not currently a monophyletic group, so phylogenetics and systematics are constantly changing. Species crypticism has been pointed out in several Palaemon species, being the clearest evidence in the European rockpool shrimp P. elegans. Here we sequenced and described seven European Palaemon mitochondrial genomes. The mitochondrial protein-coding genes were used, along with those of three other Palaemon species, to perform mitogenome phylogenetic analyses to clarify the evolutionary relationships within the genus, and particularly to shed light on the cryptic species found within P. elegans. The Messinian Salinity Crisis (5.3-5.9 Ma, late Miocene) was proposed to be the origin of this cryptic species and it was used as aged constraint for calibration analysis. We provide the largest and the first time-calibrated mitogenome phylogeny of the genus Palaemon and mitogenome substitution rate was estimated (1.59% per million years) in Decapoda for the first time. Our results highlighted the need for future systematics changes in Palaemon and crypticism in P. elegans was confirmed. Mitochondrial genome and cox1 (1.41%) substitution rate estimates matched those published elsewhere, arguing that the Messinian Salinity Crisis was a plausible event driving the split between P. elegans and its cryptic species. Molecular dating suggested that Pleistocene glaciations were likely involved in the differentiation between the Atlantic and Mediterranean populations of P. elegans. On the contrary, the divergence between the Atlantic and Mediterranean populations of the common littoral shrimp P. serratus was greater and dated to be much older (4.5-12.3 Ma, Plio-Miocene), so we considered that they could represent two separated species. Therefore, species crypticism in the genus Palaemon seems to be a common phenomenon.

The complete mitochondrial genome of the eusocial sponge-dwelling snapping shrimp Synalpheus microneptunus

In the marine realm, eusociality is only known to have evolved within a clade of sponge-dwelling snapping shrimps in the genus Synalpheus. Deciphering the genomic underpinnings of eusociality in these marine shrimps has been limited by the sparse genomic resources in this genus. Here, we report, for a eusocial shrimp Synalpheus microneptunus, a complete mitochondrial genome (22X coverage) assembled from short Illumina 150 bp paired-end reads. The 15,603 bp long mitochondrial genome of S. microneptunus is AT-rich and includes 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes and an 834 bp intergenic region assumed to be the D-loop. The gene order is identical to that reported for most caridean shrimps and corresponds to the presumed Pancrustacean ground pattern. All PCGs showed signs of purifying selection, with K_A/K_S <<1 across the whole PCGs and most sliding windows within PCGs. Maximum-likelihood and Bayesian inference phylogenetic analyses of 13 PCGs and 68 terminals supports the monophyly of the Caridea and the family Alpheidae. The complete mitochondrial genome of the eusocial shrimp Synalpheus microneptunus will contribute to a better understanding of the selective pressures and rates of molecular evolution in marine eusocial animals.

"Mitochondrial Toolbox" - A Review of Online Resources to Explore Mitochondrial Genomics

Mitochondria play a significant role in many biological systems. There is emerging evidence that differences in the mitochondrial genome may contribute to multiple common diseases, leading to an increasing number of studies exploring mitochondrial genomics. There is often a large amount of complex data generated (for example via next generation sequencing), which requires optimised bioinformatics tools to efficiently and effectively generate robust outcomes from these large datasets. Twenty-four online resources dedicated to mitochondrial genomics were reviewed. This 'mitochondrial toolbox' summary resource will enable researchers to rapidly identify the resource(s) most suitable for their needs. These resources fulfil a variety of functions, with some being highly specialised. No single tool will provide all users with the resources they require; therefore, the most suitable tool will vary between users depending on the nature of the work they aim to carry out. Genetics resources are well established for phylogeny and DNA sequence changes, but further epigenetic and gene expression resources need to be developed for mitochondrial genomics.

Insight into the resistome and quorum sensing system of a divergent Acinetobacter pittii isolate from an untouched site of the Lechuguilla Cave

Acinetobacter are Gram-negative bacteria belonging to the sub-phyla Gammaproteobacteria, commonly associated with soils, animal feeds and water. Some members of the Acinetobacter have been implicated in hospital-acquired infections, with broad-spectrum antibiotic resistance. Here we report the whole-genome sequence of LC510, an Acinetobacter species isolated from deep within a pristine location of the Lechuguilla Cave. Pairwise nucleotide comparison to three type strains within the genus Acinetobacter assigned LC510 as an Acinetobacter pittii isolate. Scanning of the LC510 genome identified two genes coding for b-lactamase resistance, despite the fact that LC510 was isolated from a portion of the cave not previously visited by humans and protected from anthropogenic input. The ability to produce acyl-homoserine lactone (AHL) signal in culture medium, an observation that is consistent with the identification of the luxI and luxR homologues in its genome, suggests that cell-to-cell communication remains important in an isolated cave ecosystem.

The complete mitochondrial genome and description of a new cryptic species of Benedenia Diesing, 1858 (Monogenea: Capsalidae), a major pathogen infecting the yellowtail kingfish Seriola lalandi Valenciennes in the South-East Pacific

Background

The monogenean Benedenia seriolae parasitizes fishes belonging to the genus Seriola, represents a species complex, and causes substantial impact on fish welfare in aquaculture systems worldwide. This study reports, for the first time, the complete mitochondrial genome of B. humboldti n. sp., a new cryptic species from the South-East Pacific (SEP).

Methods

The mitogenome of B. humboldti n. sp. was assembled from short Illumina 150 bp pair-end reads. The phylogenetic position of B. humboldti n. sp. among other closely related congeneric and confamiliar capsalids was examined using mitochondrial protein-coding genes (PCGs). Morphology of B. humboldti n. sp. was examined based on fixed and stained specimens.

Results

The AT-rich mitochondrial genome of B. humboldti is 13,455 bp in length and comprises 12 PCGs (atp8 was absent as in other monogenean genomes), 2 ribosomal RNA genes, and 22 transfer RNA genes. All protein-coding, ribosomal RNA, and transfer RNA genes are encoded on the H-strand. The gene order observed in the mitochondrial genome of B. humboldti n. sp. was identical to that of B. seriolae from Japan but different from that of B. seriolae from Australia. The genetic distance between B. humboldti n. sp. and B. seriolae from Japan was high. Minor but reliable differences in the shape of the penis were observed between Benedenia humboldti n. sp. and congeneric species.

Conclusions

Phylogenetic analyses based on PCGs in association with differences in the shape of the penis permitted us to conclude that the material from the South-East Pacific represents a new species of Benedenia infecting S. lalandi off the coast of Chile. The discovery of this parasite represents the first step to improving our understanding of infestation dynamics and to develop control strategies for this pathogen infecting the farmed yellowtail kingfish, Seriola lalandi, in the South-East Pacific.

Single Primer Enrichment Technology (SPET) for High-Throughput Genotyping in Tomato and Eggplant Germplasm

Single primer enrichment technology (SPET) is a new, robust, and customizable solution for targeted genotyping. Unlike genotyping by sequencing (GBS), and like DNA chips, SPET is a targeted genotyping technology, relying on the sequencing of a region flanking a primer. Its reliance on single primers, rather than on primer pairs, greatly simplifies panel design, and allows higher levels of multiplexing than PCR-based genotyping. Thanks to the sequencing of the regions surrounding the target SNP, SPET allows the discovery of thousands of closely linked, novel SNPs. In order to assess the potential of SPET for high-throughput genotyping in plants, a panel comprising 5k target SNPs, designed both on coding regions and introns/UTRs, was developed for tomato and eggplant. Genotyping of two panels composed of 400 tomato and 422 eggplant accessions, comprising both domesticated material and wild relatives, generated a total of 12,002 and 30,731 high confidence SNPs, respectively, which comprised both target and novel SNPs in an approximate ratio of 1:1.6, and 1:5.5 in tomato and eggplant, respectively. The vast majority of the markers was transferrable to related species that diverged up to 3.4 million years ago (Solanum pennellii for tomato and S. macrocarpon for eggplant). Maximum Likelihood phylogenetic trees and PCA outputs obtained from the whole dataset highlighted genetic relationships among accessions and species which were congruent with what was previously reported in literature. Better discrimination among domesticated accessions was achieved by using the target SNPs, while better discrimination among wild species was achieved using the whole SNP dataset. Our results reveal that SPET genotyping is a robust, high-throughput technology for genetic fingerprinting, with a high degree of cross-transferability between crops and their cultivated and wild relatives, and allows identification of duplicates and mislabeled accessions in genebanks.

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.

Absence of evidence is not evidence of absence: Nanopore sequencing and complete assembly of the European lobster (Homarus gammarus) mitogenome uncovers the missing nad2 and a new major gene cluster duplication

BACKGROUND

The recently published complete mitogenome of the European lobster (Homarus gammarus) that was generated using long-range PCR exhibits unusual gene composition (missing nad2) and gene rearrangements among decapod crustaceans with strong implications in crustacean phylogenetics. Such atypical mitochondrial features will benefit greatly from validation with emerging long read sequencing technologies such as Oxford Nanopore that can more accurately identify structural variation.

RESULTS

We re-sequenced the H. gammarus mitogenome on an Oxford Nanopore Minion flowcell and performed a long-read only assembly, generating a complete mitogenome assembly for H. gammarus. In contrast to previous reporting, we found an intact mitochondrial nad2 gene in the H. gammarus mitogenome and showed that its gene organization is broadly similar to that of the American lobster (H. americanus) except for the presence of a large tandemly duplicated region with evidence of pseudogenization in one of each duplicated protein-coding genes.

CONCLUSIONS

Using the European lobster as an example, we demonstrate the value of Oxford Nanopore long read technology in resolving problematic mitogenome assemblies. The increasing accessibility of Oxford Nanopore technology will make it an attractive and useful tool for evolutionary biologists to verify new and existing unusual mitochondrial gene rearrangements recovered using first and second generation sequencing technologies, particularly those used to make phylogenetic inferences of evolutionary scenarios.

Improved genome of Agrobacterium radiobacter type strain provides new taxonomic insight into Agrobacterium genomospecies 4

The reported Agrobacterium radiobacter DSM 30174^T genome is highly fragmented, hindering robust comparative genomics and genome-based taxonomic analysis. We re-sequenced the Agrobacterium radiobacter type strain, generating a dramatically improved genome with high contiguity. In addition, we sequenced the genome of Agrobacterium tumefaciens B6^T, enabling for the first time, a proper comparative genomics of these contentious Agrobacterium species. We provide concrete evidence that the previously reported Agrobacterium radiobacter type strain genome (Accession Number: ASXY01) is contaminated which explains its abnormally large genome size and fragmented assembly. We propose that Agrobacterium tumefaciens be reclassified as Agrobacterium radiobacter subsp. tumefaciens and that Agrobacterium radiobacter retains it species status with the proposed name of Agrobacterium radiobacter subsp. radiobacter. This proposal is based, first on the high pairwise genome-scale average nucleotide identity supporting the amalgamation of both Agrobacterium radiobacter and Agrobacterium tumefaciens into a single species. Second, maximum likelihood tree construction based on the concatenated alignment of shared genes (core genes) among related strains indicates that Agrobacterium radiobacter NCPPB3001 is sufficiently divergent from Agrobacterium tumefaciens to propose two independent sub-clades. Third, Agrobacterium tumefaciens demonstrates the genomic potential to synthesize the L configuration of fucose in its lipid polysaccharide, fostering its ability to colonize plant cells more effectively than Agrobacterium radiobacter.

Phylomitogenomics reconfirm the phylogenetic position of the genus Metaplax inferred from the two grapsid crabs (Decapoda: Brachyura: Grapsoidea)

Two new complete mitogenomes of the grapsids, Metaplax longipes Stimpson, 1858 and Nanosesarma minutum (De Man, 1887) were sequenced using next-generation sequencing (NGS). By analyzing a combination of 75 Brachyura taxa, our phylomitogenomic inferences suggested that Metaplax crab seperated earlier from the sesarmid crabs and closely related to the varunids with respect to Nanosesarma crab. It reconfirmed that the Metaplax should be removed from the Sesarmidae and assinged to the Varunidae. Additional mitogenomic comparisons including gene rearrangement and genomic organization were conducted among the 33 taxa of Grapsoidea and Ocypodoidea, and a shared rearrangement pattern between Metaplax longipes and the varunids were recovered, which also strongly supported the inference for the phylogenetic position of the Metaplax.

The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus

Panulirus argus is a keystone species and target of the most lucrative fishery in the Caribbean region. This study reports, for the first time, the complete mitochondrial genome of Panulirus argus (average coverage depth nucleotide-1 = 70×) assembled from short Illumina 150 bp PE reads. The AT-rich mitochondrial genome of Panulirus argus was 15 739 bp in length and comprised 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 801 bp long intergenic space was assumed to be the D-loop. Most of the PCGs were encoded on the H-strand. The gene order observed in the mitochondrial genome of Panulirus argus corresponds to the presumed Pancrustacean ground pattern. KA/KS ratios calculated for all mitochondrial PCGs showed values < 1, indicating that all these PCGs are evolving under purifying selection. A maximum likelihood phylogenetic analysis (concatenated PCGs [n = 13], 154 arthropods) supported the monophyly of the Achelata and other infraorders within the Decapoda. Mitochondrial PCGs have enough phylogenetic informativeness to explore high-level genealogical relationships in the Pancrustacea. The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus will contribute to the better understanding of meta-population connectivity in this keystone overexploited species.

Characterization of four new mitogenomes from Ocypodoidea & Grapsoidea, and phylomitogenomic insights into thoracotreme evolution

Four new complete mitochondrial genomes (mitogenomes) from the two superfamilies Ocypodoidea and Grapsoidea were sequenced, which represented Uca (Gelasimus) borealis (Ocypodidae: Ucinae), Dotilla wichmani (Dotillidae), Metopograpsus quadridentatus (Grapsidae: Grapsinae), and Gaetice depressus (Varunidae: Gaeticinae). All of the mitogenomes shared the complete set of 37 mitochondrial genes. Mitogenome lengths were 15,659, 15,600, 15,517, and 16,288 bp, respectively, with A + T contents of 69.41%, 68.46%, 70.30%, and 72.96%, respectively. Comparative genomic analyses suggested that they exhibited different genomic rearrangements. In particular, G. depressus shared a major rearrangement pattern present in Eriocheir crabs, while the remainder shared the brachyuran ground genomic rearrangement patterns. Phylomitogenomic inferences provided new evidence for the strongly supported nesting of Thoracotremata within Heterotremata clades. A close phylogenetic relationship was observed between Varunidae and Macrophthalmidae crabs, and between Dotillidae and Grapsidae crabs, which was consistent with mitochondrial genomic rearrangement similarities. Altogether, these results suggest the presence of reciprocal paraphyly for Ocypodoidea and Grapsoidea.

More limbs on the tree: mitogenome characterisation and systematic position of ‘living fossil’ species Neoglyphea inopinata and Laurentaeglyphea neocaledonica (Decapoda : Glypheidea : Glypheidae)

Glypheids first appeared in the Lower Triassic period and were believed to be extinct until specimens of Neoglyphea inopinata Forest &amp; Saint Laurent and Laurentaeglyphea neocaledonica Richer de Forges were described in 1975 and 2006, respectively. The finding of extant species has meant that molecular data can now be used to complement morphological and fossil-based studies to investigate the relationships of Glypheidea within the Decapoda. However, despite several molecular studies, the placement of this infraorder within the decapod phylogenetic tree is not resolved. One limitation is that molecular resources available for glypheids have been limited to a few nuclear and mitochondrial gene fragments. Many of the more recent large-scale studies of decapod phylogeny have used information from complete mitogenomes, but have excluded the infraorder Glypheidea due to the unavailability of complete mitogenome sequences. Using next-generation sequencing, we successfully sequenced and assembled complete mitogenome sequences from museum specimens of N. inopinata and L. neocaledonica, the only two extant species of glypheids. With these sequences, we constructed the first decapod phylogenetic tree based on whole mitogenome sequences that includes Glypheidea as one of 10 decapod infraorders positioned within the suborder Pleocyemata. From this, the Glypheidea appears to be a relatively derived lineage related to the Polychelida and Astacidea. Also in our study, we conducted a survey on currently available decapod mitogenome resources available on National Center for Biotechnology Information (NCBI) and identified infraorders that would benefit from more strategic and expanded taxonomic sampling.

Complete mitochondrial genomes from two species of Chinese freshwater crabs of the genus Sinopotamon recovered using next-generation sequencing reveal a novel gene order (Brachyura, Potamidae)

Recent morphological and molecular evidence has challenged classical interpretations of eubrachyuran phylogeny and evolution. Complete mitochondrial genomes of two species of potamid freshwater crabs, Sinopotamon yaanense and Sinopotamon yangtsekiense were obtained using next-generation sequencing. The results revealed a novel gene order with translocations of a five-gene block and a tRNA gene in comparison to available brachyuran mitochondrial genomes. DNA sequence comparisons position the Potamidae, a primary freshwater crab family, outside of the clade for the traditional heterotreme families, and closer to the clade that includes the thoracotreme families of grapsoid and ocypodoid crabs. Mitogenomic comparisons using rapid next-generation sequencing and a much wider taxonomic sample are required for a high-resolution examination of the phylogenetic relationships within the Eubrachyura.

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.

The complete mitochondrial genome of the grooved carpet shell, Ruditapes decussatus (Bivalvia, Veneridae)

Despite the large number of animal complete mitochondrial genomes currently available in public databases, knowledge about mitochondrial genomics in invertebrates is uneven. This paper reports, for the first time, the complete mitochondrial genome of the grooved carpet shell, Ruditapes decussatus, also known as the European clam. Ruditapes decussatus is morphologically and ecologically similar to the Manila clam Ruditapes philippinarum, which has been recently introduced for aquaculture in the very same habitats of Ruditapes decussatus, and that is replacing the native species. Currently the production of the European clam is almost insignificant, nonetheless it is considered a high value product, and therefore it is an economically important species, especially in Portugal, Spain and Italy. In this work we: (i) assembled Ruditapes decussatus mitochondrial genome from RNA-Seq data, and validated it by Sanger sequencing; (ii) analyzed and characterized the Ruditapes decussatus mitochondrial genome, comparing its features with those of other venerid bivalves; (iii) assessed mitochondrial sequence polymorphism (SP) and copy number variation (CNV) of tandem repeats across 26 samples. Despite using high-throughput approaches we did not find evidence for the presence of two sex-linked mitochondrial genomes, typical of the doubly uniparental inheritance of mitochondria, a phenomenon known in ∼100 bivalve species. According to our analyses, Ruditapes decussatus is more genetically similar to species of the Genus Paphia than to the congeneric Ruditapes philippinarum, a finding that bolsters the already-proposed need of a taxonomic revision. We also found a quite low genetic variability across the examined samples, with few SPs and little variability of the sequences flanking the control region (Largest Unassigned Regions (LURs). Strikingly, although we found low nucleotide variability along the entire mitochondrial genome, we observed high levels of length polymorphism in the LUR due to CNV of tandem repeats, and even a LUR length heteroplasmy in two samples. It is not clear if the lack of genetic variability in the mitochondrial genome of Ruditapes decussatus is a cause or an effect of the ongoing replacement of Ruditapes decussatus with the invasive Ruditapes philippinarum, and more analyses, especially on nuclear sequences, are required to assess this point.

Digging deeper: new gene order rearrangements and distinct patterns of codons usage in mitochondrial genomes among shrimps from the Axiidea, Gebiidea and Caridea (Crustacea: Decapoda)

BACKGROUND

Whole mitochondrial DNA is being increasingly utilized for comparative genomic and phylogenetic studies at deep and shallow evolutionary levels for a range of taxonomic groups. Although mitogenome sequences are deposited at an increasing rate into public databases, their taxonomic representation is unequal across major taxonomic groups. In the case of decapod crustaceans, several infraorders, including Axiidea (ghost shrimps, sponge shrimps, and mud lobsters) and Caridea (true shrimps) are still under-represented, limiting comprehensive phylogenetic studies that utilize mitogenomic information.

METHODS

Sequence reads from partial genome scans were generated using the Illumina MiSeq platform and mitogenome sequences were assembled from these low coverage reads. In addition to examining phylogenetic relationships within the three infraorders, Axiidea, Gebiidea, and Caridea, we also investigated the diversity and frequency of codon usage bias and mitogenome gene order rearrangements.

RESULTS

We present new mitogenome sequences for five shrimp species from Australia that includes two ghost shrimps, Callianassa ceramica and Trypaea australiensis, along with three caridean shrimps, Macrobrachium bullatum, Alpheus lobidens, and Caridina cf. nilotica. Strong differences in codon usage were discovered among the three infraorders and significant gene order rearrangements were observed. While the gene order rearrangements are congruent with the inferred phylogenetic relationships and consistent with taxonomic classification, they are unevenly distributed within and among the three infraorders.

DISCUSSION

Our findings suggest potential for mitogenome rearrangements to be useful phylogenetic markers for decapod crustaceans and at the same time raise important questions concerning the drivers of mitogenome evolution in different decapod crustacean lineages.