The use of transcriptomic next-generation sequencing data to assemble mitochondrial genomes of Ancistrus spp. (Loricariidae)

In silico prospecting of the mtDNA of Macrobrachium amazonicum from transcriptome data

BACKGROUND

Macrobrachium amazonicum is a freshwater prawn widely distributed in South America that is undergoing speciation, so the denomination "M. amazonicum complex" is used for it. The mitochondrial cytochrome c oxidase subunit I (COI) gene has been used to elucidate this speciation, but heteroplasmies and pseudogenes have been recorded, making separation difficult. Obtaining genes from cDNA (RNA) rather than genomic DNA is an effective tool to mitigate those two types of occurrences. The aim of this study was to assemble in silico the mitochondrial DNA (mtDNA) of the Amazonian coastal population of M. amazonicum inhabiting the state of Pará.

RESULTS

Sequences were obtained from the prawn's transcriptome using the de novo approach. Six libraries of cDNA from the androgen gland, hepatopancreas, and muscle tissue were used. The mtDNA of M. amazonicum was 14,960 bp in length. It contained 13 protein-coding genes, 21 complete transfer RNAs, and the 12S and 16S subunits of ribosomal RNA. All regions were found on the light strand except tRNAGln, which was on the heavy strand. The control region (D-loop) was not recovered, making for a gap of 793 bp. The cladogram showed the formation of the well-defined Macrobrachium clade, with high support value in the established branches (91-100). The three-dimensional spatial conformation of the mtDNA-encoded proteins showed that most of them were mainly composed of major α-helices that typically shows in those proteins inserted in the membrane (mitochondrial).

CONCLUSIONS

It was possible to assemble a large part of the mitochondrial genome of M. amazonicum in silico using data from other genomes deposited in GenBank and to validate it through the similarities between its COI and 16S genes and those from animals of the same region deposited in GenBank. Depositing the M. amazonicum mtDNA sequences in GenBank may help solve the taxonomic problems recorded for the species, in addition to providing complete sequences of candidate coding genes for use as biomarkers in ecological studies.

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.

Complete mitochondrial genome of Phoneutria depilata (Araneae, Ctenidae): New Insights into the Phylogeny and Evolution of Spiders

Spiders (Araneae) are the most abundant terrestrial predators and megadiverse on earth. In recent years, the mitochondrial genome of a great diversity of species has been sequenced, mainly for ecological and commercial purposes. These studies have uncovered the existence of a variety of mitochondrial genome rearrangements. However, there is poor genetic information in several taxonomic families of spiders. We have sequenced the complete genome of Phoneutria depilata (Ctenidae) and, based on this, extract the mitogenomes of other ctenid species from published transcriptomes to perform a comparative study among spider species to determine the relationship between the level of mitochondrial rearrangements and its possible relationship with molecular variability in spiders. Complete mitochondrial genomes of eighteen spiders (including eight Ctenidae species) were obtained by two different methodologies (sequencing and transcriptome extraction). Fifty-eight spider mitochondrial genomes were downloaded from the NCBI database for gene order analysis. After verifying the annotation of each mitochondrial gene, a phylogenetic and a gene order analysis from 76 spider mitochondrial genomes were carried out. Our results show a high rate of annotation error in the published spider mitochondrial genomes, which could lead to errors in phylogenetic inference. Moreover, to provide new mitochondrial genomes in spiders by two different methodologies to obtain them, our analysis identifies six different mitochondrial architectures among all spiders. Translocation or tandem duplication random loss (TDRL) events in tRNA genes were identified to explain the evolution of the spider mitochondrial genome. In addition, our findings provide new insights into spider mitochondrial evolution.

Extension of Mitogenome Enrichment Based on Single Long-Range PCR: mtDNAs and Putative Mitochondrial-Derived Peptides of Five Rodent Hibernators

Enriching mitochondrial DNA (mtDNA) for sequencing entire mitochondrial genomes (mitogenomes) can be achieved by single long-range PCR. This avoids interference from the omnipresent nuclear mtDNA sequences (NUMTs). The approach is currently restricted to the use of samples collected from humans and ray-finned fishes. Here, we extended the use of single long-range PCR by introducing back-to-back oligonucleotides that target a sequence of extraordinary homology across vertebrates. The assay was applied to five hibernating rodents, namely alpine marmot, Arctic and European ground squirrels, and common and garden dormice, four of which have not been fully sequenced before. Analysis of the novel mitogenomes focussed on the prediction of mitochondrial-derived peptides (MDPs) providing another level of information encoded by mtDNA. The comparison of MOTS-c, SHLP4 and SHLP6 sequences across vertebrate species identified segments of high homology that argue for future experimentation. In addition, we evaluated four candidate polymorphisms replacing an amino acid in mitochondrially encoded subunits of the oxidative phosphorylation (OXPHOS) system that were reported in relation to cold-adaptation. No obvious pattern was found for the diverse sets of mammalian species that either apply daily or multiday torpor or otherwise cope with cold. In summary, our single long-range PCR assay applying a pair of back-to-back primers that target a consensus sequence motif of Vertebrata has potential to amplify (intact) mitochondrial rings present in templates from a taxonomically diverse range of vertebrates. It could be promising for studying novel mitogenomes, mitotypes of a population and mitochondrial heteroplasmy in a sensitive, straightforward and flexible manner.

MITGARD: an automated pipeline for mitochondrial genome assembly in eukaryotic species using RNA-seq data

MOTIVATION

Over the past decade, the field of next-generation sequencing (NGS) has seen dramatic advances in methods and a decrease in costs. Consequently, a large expansion of data has been generated by NGS, most of which have originated from RNA-sequencing (RNA-seq) experiments. Because mitochondrial genes are expressed in most eukaryotic cells, mitochondrial mRNA sequences are usually co-sequenced within the target transcriptome, generating data that are commonly underused or discarded. Here, we present MITGARD, an automated pipeline that reliably recovers the mitochondrial genome from RNA-seq data from various sources. The pipeline identifies mitochondrial sequence reads based on a phylogenetically related reference, assembles them into contigs, and extracts a complete mtDNA for the target species.

RESULTS

We demonstrate that MITGARD can reconstruct the mitochondrial genomes of several species throughout the tree of life. We noticed that MITGARD can recover the mitogenomes in different sequencing schemes and even in a scenario of low-sequencing depth. Moreover, we showed that the use of references from congeneric species diverging up to 30 million years ago (MYA) from the target species is sufficient to recover the entire mitogenome, whereas the use of species diverging between 30 and 60 MYA allows the recovery of most mitochondrial genes. Additionally, we provide a case study with original data in which we estimate a phylogenetic tree of snakes from the genus Bothrops, further demonstrating that MITGARD is suitable for use on biodiversity projects. MITGARD is then a valuable tool to obtain high-quality information for studies focusing on the phylogenetic and evolutionary aspects of eukaryotes and provides data for easily identifying a sample using barcoding, and to check for cross-contamination using third-party tools.

Complete mitochondrial genomes from transcriptomes: assessing pros and cons of data mining for assembling new mitogenomes

Thousands of eukaryotes transcriptomes have been generated, mainly to investigate nuclear genes expression, and the amount of available data is constantly increasing. A neglected but promising use of this large amount of data is to assemble organelle genomes. To assess the reliability of this approach, we attempted to reconstruct complete mitochondrial genomes from RNA-Seq experiments of Reticulitermes termite species, for which transcriptomes and conspecific mitogenomes are available. We successfully assembled complete molecules, although a few gaps corresponding to tRNAs had to be filled manually. We also reconstructed, for the first time, the mitogenome of Reticulitermes banyulensis. The accuracy and completeness of mitogenomes reconstruction appeared independent from transcriptome size, read length and sequencing design (single/paired end), and using reference genomes from congeneric or intra-familial taxa did not significantly affect the assembly. Transcriptome-derived mitogenomes were found highly similar to the conspecific ones obtained from genome sequencing (nucleotide divergence ranging from 0% to 3.5%) and yielded a congruent phylogenetic tree. Reads from contaminants and nuclear transcripts, although slowing down the process, did not result in chimeric sequence reconstruction. We suggest that the described approach has the potential to increase the number of available mitogenomes by exploiting the rapidly increasing number of transcriptomes.

Leishmania Mitochondrial Genomes: Maxicircle Structure and Heterogeneity of Minicircles

The mitochondrial DNA (mtDNA), which is present in almost all eukaryotic organisms, is a useful marker for phylogenetic studies due to its relative high conservation and its inheritance manner. In Leishmania and other trypanosomatids, the mtDNA (also referred to as kinetoplast DNA or kDNA) is composed of thousands of minicircles and a few maxicircles, catenated together into a complex network. Maxicircles are functionally similar to other eukaryotic mtDNAs, whereas minicircles are involved in RNA editing of some maxicircle-encoded transcripts. Next-generation sequencing (NGS) is increasingly used for assembling nuclear genomes and, currently, a large number of genomic sequences are available. However, most of the time, the mitochondrial genome is ignored in the genome assembly processes. The aim of this study was to develop a pipeline to assemble Leishmania minicircles and maxicircle DNA molecules, exploiting the raw data generated in the NGS projects. As a result, the maxicircle molecules and the plethora of minicircle classes for Leishmania major, Leishmania infantum and Leishmania braziliensis have been characterized. We have observed that whereas the heterogeneity of minicircle sequences existing in a single cell hampers their use for Leishmania typing and classification, maxicircles emerge as an extremely robust genetic marker for taxonomic studies within the clade of kinetoplastids.

Transcriptome-based molecular systematics: Rhodnius montenegrensis (Triatominae) and its position within the Rhodnius prolixus-Rhodnius robustus cryptic-species complex

BACKGROUND

Rhodnius montenegrensis (Triatominae), a potential vector of Chagas disease, was described after R. robustus-like bugs from southwestern Amazonia. Mitochondrial cytb sequence near-identity with sympatric R. robustus (genotype II) raised doubts about the taxonomic status of R. montenegrensis, but comparative studies have reported fairly clear morphological and genetic differences between R. montenegrensis and laboratory stocks identified as R. robustus. Here, we use a transcriptome-based approach to investigate this apparent paradox.

RESULTS

We retrieved publicly-available transcriptome sequence-reads from R. montenegrensis and from the R. robustus stocks used as the taxonomic benchmark in comparative studies. We (i) aligned transcriptome sequence-reads to mitochondrial (cytb) and nuclear (ITS2, D2-28S and AmpG) query sequences (47 overall) from members of the R. prolixus-R. robustus cryptic-species complex and related taxa; (ii) computed breadth- and depth-coverage for the 259 consensus sequences generated by these alignments; and, for each locus, (iii) appraised query sequences and full-breadth-coverage consensus sequences in terms of nucleotide-sequence polymorphism and phylogenetic relations. We found evidence confirming that R. montenegrensis and R. robustus genotype II are genetically indistinguishable and, hence, implying that they are, in all likelihood, the same species. Furthermore, we found compelling genetic evidence that the benchmark 'R. robustus' stocks used in R. montenegrensis description and in later transcriptome-based comparisons are in fact R. prolixus, although likely mixed to some degree with R. robustus (probably genotype II, a.k.a. R. montenegrensis).

CONCLUSIONS

We illustrate how public-domain genetic/transcriptomic data can help address challenging issues in disease-vector systematics. In our case-study, taxonomic confusion apparently stemmed from the misinterpretation of sequence-data analyses and misidentification of taxonomic-benchmark stocks. More generally, and together with previous reports of mixed and/or misidentified Rhodnius spp. laboratory colonies, our results call into question the conclusions of many studies (on morphology, genetics, physiology, behavior, bionomics or interactions with microorganisms including trypanosomes) based on non-genotyped 'R. prolixus' or 'R. robustus' stocks. Correct species identification is a prerequisite for investigating the factors that underlie the physiological, behavioral or ecological differences between primary domestic vectors of Chagas disease, such as R. prolixus, and their sylvatic, medically less-relevant relatives such as R. robustus (s.l.) including R. montenegrensis.

Reducing the information gap on Loricarioidei (Siluriformes) mitochondrial genomics

Background

The genetic diversity of Neotropical fish fauna is underrepresented in public databases. This distortion is evident for the order Siluriformes, in which the suborders Siluroidei and Loricarioidei share equivalent proportion of species, although far less is known about the genetics of the latter clade, endemic to the Neotropical Region. Recently, this information gap was evident in a study about the structural diversity of fish mitochondrial genomes, and hampered a precise chronological resolution of Siluriformes. It has also prevented molecular ecology investigations about these catfishes, their interactions with the environment, responses to anthropogenic changes and potential uses.

Results

Using high-throughput sequencing, we provide the nearly complete mitochondrial genomes for 26 Loricariidae and one Callichthyidae species. Structural features were highly conserved. A notable exception was identified in the monophyletic clade comprising species of the Hemiancistrus, Hypostomini and Peckoltia-clades, a ~60 nucleotide-long deletion encompassing the seven nucleotides at the 3′ end of the Conserved Sequence Block (CSB) D of the control region. The expression of mitochondrial genes followed the usual punctuation pattern. Heteroplasmic sites were identified in most species. The retrieved phylogeny strongly corroborates the currently accepted tree, although bringing to debate the relationship between Schizolecis guntheri and Pareiorhaphis garbei, and highlighting the low genetic variability within the Peckoltia-clade, an eco-morphologically diverse and taxonomically problematic group.

Conclusions

Herein we have launched the use of high-throughput mitochondrial genomics in the studies of the Loricarioidei species. The new genomic resources reduce the information gap on the molecular diversity of Neotropical fish fauna, impacting the capacity to investigate a variety of aspects of the molecular ecology and evolution of these fishes. Additionally, the species showing the partial CSB-D are candidate models to study the replication and transcription of vertebrate mitochondrial genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3709-3) contains supplementary material, which is available to authorized users.

The liver transcriptome of suckermouth armoured catfish (Pterygoplichthys anisitsi, Loricariidae): Identification of expansions in defensome gene families

Pterygoplichthys is a genus of related suckermouth armoured catfishes native to South America, which have invaded tropical and subtropical regions worldwide. Physiological features, including an augmented resistance to organic xenobiotics, may have aided their settlement in foreign habitats. The liver transcriptome of Pterygoplichthys anisitsi was sequenced and used to characterize the diversity of mRNAs potentially involved in the responses to natural and anthropogenic chemicals. In total, 66,642 transcripts were assembled. Among the identified defensome genes, cytochromes P450 (CYP) were the most abundant, followed by sulfotransferases (SULT), nuclear receptors (NR) and ATP binding cassette transporters (ABC). A novel expansion in the CYP2Y subfamily was identified, as well as an independent expansion of the CYP2AAs. Two expansions were also observed among SULT1. Thirty-two transcripts were classified into twelve subfamilies of NR, while 21 encoded ABC transporters. The diversity of defensome transcripts sequenced herein could contribute to this species' resistance to organic xenobiotics.

Next-generation sequencing of mixed genomic DNA allows efficient assembly of rearranged mitochondrial genomes in Amolops chunganensis and Quasipaa boulengeri

Recent improvements in next-generation sequencing (NGS) technologies can facilitate the obtainment of mitochondrial genomes. However, it is not clear whether NGS could be effectively used to reconstruct the mitogenome with high gene rearrangement. These high rearrangements would cause amplification failure, and/or assembly and alignment errors. Here, we choose two frogs with rearranged gene order, Amolops chunganensis and Quasipaa boulengeri, to test whether gene rearrangements affect the mitogenome assembly and alignment by using NGS. The mitogenomes with gene rearrangements are sequenced through Illumina MiSeq genomic sequencing and assembled effectively by Trinity v2.1.0 and SOAPdenovo2. Gene order and contents in the mitogenome of A. chunganensis and Q. boulengeri are typical neobatrachian pattern except for rearrangements at the position of "WANCY" tRNA genes cluster. Further, the mitogenome of Q. boulengeri is characterized with a tandem duplication of trnM. Moreover, we utilize 13 protein-coding genes of A. chunganensis, Q. boulengeri and other neobatrachians to reconstruct the phylogenetic tree for evaluating mitochondrial sequence authenticity of A. chunganensis and Q. boulengeri. In this work, we provide nearly complete mitochondrial genomes of A. chunganensis and Q. boulengeri.

The complete mitochondrial genome of Corydoras nattereri (Callichthyidae: Corydoradinae)

ABSTRACT The complete mitogenome of Corydoras nattereri , a species of mailed catfishes from southeastern Brazil, was reconstructed using next-generation sequencing techniques. The mitogenome was assembled using mitochondrial transcripts from the liver transcriptomes of three individuals, and produced a circular DNA sequence of 16,557 nucleotides encoding 22 tRNA genes, two rRNA genes, 13 protein-coding genes and two noncoding control regions (D-loop, OrigL). Phylogeographic analysis of closely related sequences of Cytochrome Oxydase C subunit I (COI) demonstrates high diversity among morphologically similar populations of C. nattereri . Corydoras nattereri is nested within a complex of populations currently assigned to C. paleatus and C. ehrhardti . Analysis of mitogenome structure demonstrated that an insertion of 21 nucleotides between the ATPase subunit-6 and COIII genes may represent a phylogenetically informative character associated with the evolution of the Corydoradinae.