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.
Bactrocera latifrons is a serious pest of solanaceous fruits and Bactrocera umbrosa is a pest of Artocarpus fruits, while Bactrocera melastomatos infests the fruit of Melastomataceae. They are members of the subgenus Bactrocera. We report here the complete mitochondrial genome of these fruit flies determined by next-generation sequencing and their phylogeny with other taxa of the subgenus Bactrocera. The whole mitogenomes of these three species possessed 37 genes namely, 13 protein-coding genes (PCGs), 2 rRNA and 22 tRNA genes. The mitogenome of B. latifrons (15,977 bp) was longer than those of B. melastomatos (15,954 bp) and B. umbrosa (15,898 bp). This difference can be attributed to the size of the intergenic spacers (283 bp in B. latifrons, 261 bp in B. melastomatos, and 211 bp in B. umbrosa). Most of the PCGs in the three species have an identical start codon, except for atp8 (adenosine triphosphate synthase protein 8), which had an ATG instead of GTG in B. umbrosa, whilst the nad3 (NADH dehydrogenase subunit 3) and nad6 (NADH dehydrogenase subunit 6) genes were characterized by an ATC instead of ATT in B. melastomatos. The three species had identical stop codon for the respective PCGs. In B. latifrons and B. melastomatos, the TΨC (thymidine-pseudouridine-cytidine)-loop was absent in trnF (phenylalanine) and DHU (dihydrouracil)-loop was absent in trnS1 (serine S1). In B. umbrosa, trnN (asparagine), trnC (cysteine) and trnF lacked the TψC-loop, while trnS1 lacked the DHU-stem. Molecular phylogeny based on 13 PCGs was in general concordant with 15 mitochondrial genes (13 PCGs and 2 rRNA genes), with B. latifrons and B. umbrosa forming a sister group basal to the other species of the subgenus Bactrocera which was monophyletic. The whole mitogenomes will serve as a useful dataset for studying the genetics, systematics and phylogenetic relationships of the many species of Bactrocera genus in particular, and tephritid fruit flies in general.
In this study, we analyzed the complete mitochondrial genome of the cavity-nesting honeybee, A. koschevnikovi. The mitochondrial genome of A. koschevnikovi was observed to be a circular molecule of 15,278 bp and was similar to that of the other cavity-nesting honeybee species. The average AT content in the A. koschevnikovi mitochondrial genome was 84%. It was predicted to contain 13 protein-coding, 24 tRNA and two rRNA genes, along with one A + T-rich control region, besides three tRNA-Met repeats.
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.
We sequenced and assembled three whole mitogenome sequences of the commercially important snakeskin gourami Trichopodus pectoralis isolated from Malaysia (introduced), Viet Nam (native) and Thailand (native). The mitogenome length range from 16,397 to 16,420 bp. The final partitioned nucleotide alignment consists of 14,002 bp and supports the monophyly of the genus Trichopodus (95% ultrafast bootstrap support) with T. trichopterus forming a sister group with the members of T. pectoralis.
Genetic variation in mitochondrial genes could underlie metabolic adaptations because mitochondrially encoded proteins are directly involved in a pathway supplying energy to metabolism. Macquarie perch from river basins exposed to different climates differ in size and growth rate, suggesting potential presence of adaptive metabolic differences. We used complete mitochondrial genome sequences to build a phylogeny, estimate lineage divergence times and identify signatures of purifying and positive selection acting on mitochondrial genes for 25 Macquarie perch from three basins: Murray-Darling Basin (MDB), Hawkesbury-Nepean Basin (HNB) and Shoalhaven Basin (SB). Phylogenetic analysis resolved basin-level clades, supporting incipient speciation previously inferred from differentiation in allozymes, microsatellites and mitochondrial control region. The estimated time of lineage divergence suggested an early- to mid-Pleistocene split between SB and the common ancestor of HNB+MDB, followed by mid-to-late Pleistocene splitting between HNB and MDB. These divergence estimates are more recent than previous ones. Our analyses suggested that evolutionary drivers differed between inland MDB and coastal HNB. In the cooler and more climatically variable MDB, mitogenomes evolved under strong purifying selection, whereas in the warmer and more climatically stable HNB, purifying selection was relaxed. Evidence for relaxed selection in the HNB includes elevated transfer RNA and 16S ribosomal RNA polymorphism, presence of potentially mildly deleterious mutations and a codon (ATP6113) displaying signatures of positive selection (ratio of nonsynonymous to synonymous substitution rates (dN/dS) >1, radical change of an amino-acid property and phylogenetic conservation across the Percichthyidae). In addition, the difference could be because of stronger genetic drift in the smaller and historically more subdivided HNB with low per-population effective population sizes.
The complete mitochondrial genome of the cavity-nesting honeybee Apis cerana from Sabah on Borneo Island was analyzed using next-generation sequencing. The mitochondrial genome of A. cerana was a circular molecule of 15,884 bp and was similar to that of the other cavity-nesting honeybee species. The average AT content in the A. cerana mitochondrial genome was 84.4%. It was predicted to contain 13 protein-coding, 22 tRNA, and two rRNA genes, along with one A + T-rich control region.
Southeast (SE) Asia is 1 of the most biodiverse regions in the world, and it holds approximately 20% of all mammal species. Despite this, the majority of SE Asia's genetic diversity is still poorly characterized. The growing interest in using environmental DNA to assess and monitor SE Asian species, in particular threatened mammals-has created the urgent need to expand the available reference database of mitochondrial barcode and complete mitogenome sequences. We have partially addressed this need by generating 72 new mitogenome sequences reconstructed from DNA isolated from a range of historical and modern tissue samples. Approximately 55 gigabases of raw sequence were generated. From this data, we assembled 72 complete mitogenome sequences, with an average depth of coverage of ×102.9 and ×55.2 for modern samples and historical samples, respectively. This dataset represents 52 species, of which 30 species had no previous mitogenome data available. The mitogenomes were geotagged to their sampling location, where known, to display a detailed geographical distribution of the species. Our new database of 52 taxa will strongly enhance the utility of environmental DNA approaches for monitoring mammals in SE Asia as it greatly increases the likelihoods that identification of metabarcoding sequencing reads can be assigned to reference sequences. This magnifies the confidence in species detections and thus allows more robust surveys and monitoring programmes of SE Asia's threatened mammal biodiversity. The extensive collections of historical samples from SE Asia in western and SE Asian museums should serve as additional valuable material to further enrich this reference database.
Currently, the effects of chronic, continuous low dose environmental irradiation on the mitochondrial genome of resident small mammals are unknown. Using the bank vole (Myodes glareolus) as a model system, we tested the hypothesis that approximately 50 generations of exposure to the Chernobyl environment has significantly altered genetic diversity of the mitochondrial genome. Using deep sequencing, we compared mitochondrial genomes from 131 individuals from reference sites with radioactive contamination comparable to that present in northern Ukraine before the 26 April 1986 meltdown, to populations where substantial fallout was deposited following the nuclear accident. Population genetic variables revealed significant differences among populations from contaminated and uncontaminated localities. Therefore, we rejected the null hypothesis of no significant genetic effect from 50 generations of exposure to the environment created by the Chernobyl meltdown. Samples from contaminated localities exhibited significantly higher numbers of haplotypes and polymorphic loci, elevated genetic diversity, and a significantly higher average number of substitutions per site across mitochondrial gene regions. Observed genetic variation was dominated by synonymous mutations, which may indicate a history of purify selection against nonsynonymous or insertion/deletion mutations. These significant differences were not attributable to sample size artifacts. The observed increase in mitochondrial genomic diversity in voles from radioactive sites is consistent with the possibility that chronic, continuous irradiation resulting from the Chernobyl disaster has produced an accelerated mutation rate in this species over the last 25 years. Our results, being the first to demonstrate this phenomenon in a wild mammalian species, are important for understanding genetic consequences of exposure to low-dose radiation sources.
The cavity-nesting honeybee Apis nuluensis inhabits only the highlands of Mount Kinabalu of Sabah, Borneo Island. The mitochondrial genome is a circular molecule of approximately 1.6 kb that includes 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and one AT-rich control region. The average AT content was 84.5%. The start codons ATC, ATG, and ATT were found in one, three, and nine genes, respectively, whereas the stop codon TAA was observed in all genes. The phylogenetic relationship, inferred using 13 PCGs, was consistent with that reported in previous studies that predicted a sister taxon relationship between A. nuluensis and A. cerana.
Although the role of nuclear-encoded gene alterations has been well documented in brain tumor development, the involvement of the mitochondrial genome in brain tumorigenesis has not yet been fully elucidated and remains controversial. The present study aimed to identify mutations in the mitochondrial DNA (mtDNA) control region D-loop in patients with brain tumors in Malaysia. A mutation analysis was performed in which DNA was extracted from paired tumor tissue and blood samples obtained from 49 patients with brain tumors. The D-loop region DNA was amplified using the PCR technique, and genetic data from DNA sequencing analyses were compared with the published revised Cambridge sequence to identify somatic mutations. Among the 49 brain tumor tissue samples evaluated, 25 cases (51%) had somatic mutations of the mtDNA D-loop, with a total of 48 mutations. Novel mutations that had not previously been identified in the D-loop region (176 A-deletion, 476 C>A, 566 C>A and 16405 A-deletion) were also classified. No significant associations between the D-loop mutation status and the clinicopathological parameters were observed. To the best of our knowledge, the current study presents the first evidence of alterations in the mtDNA D-loop regions in the brain tumors of Malaysian patients. These results may provide an overview and data regarding the incidence of mitochondrial genome alterations in Malaysian patients with brain tumors. In addition to nuclear genome aberrations, these specific mitochondrial genome alterations may also be considered as potential cancer biomarkers for the diagnosis and staging of brain cancers.
The role of mitochondria in tumorigenesis has regained much attention as it could dysregulate cellular energetics, oxidative stress and apoptosis. However, the role of mitochondria in different grade gliomasis still unknown. This study aimed to identify mitochondrial DNA (mtDNA) sequence variations that could possibly affect the mitochondrial functions and also the oxidative stress status. Three different grades of human glioma cell lines and a normal human astrocyte cell line were cultured in-vitro and tested for oxidative stress biomarkers. Relative oxidative stress level, mitochondria activity, and mitochondrial mass were determined by live cell imaging with confocal laser scanning microscope using CM-H2DCFDA, MitoTracker Green, and MitoTracker Orange stains. The entire mitochondrial genome was sequenced using the AffymetrixGeneChip Human Mitochondrial Resequencing Array 2.0. The mitochondrial sequence variations were subjected to phylogenetic haplogroup assessment and pathogenicity of the mutations were predicted using pMUT and PolyPhen2. The Grade II astrocytoma cells showed increased oxidative stress wherea high level of 8-OHdG and oxidative stress indicator were observed. Simultaneously, Grade II and III glioma cells showed relatively poor mitochondria functions and increased number of mutations in the coding region of the mtDNA which could be due to high levels of oxidative stress in these cells. These non-synonymous mtDNA sequence variations were predicted to be pathogenic and could possibly lead to protein dysfunction, leading to oxidative phosphorylation (OXPHOS) impairment, mitochondria dysfunction and could create a vicious cycle of oxidative stress. The Grade IV cells had no missense mutation but preserved intact mitochondria and excellent antioxidant defense mechanisms thus ensuring better survival. In conclusion, Grade II and III glioma cells demonstrated coding region mtDNA mutations, leading to mitochondrial dysfunction and higher oxidative stress.
In industry, the yeast Rhodotorula mucilaginosa is commonly used for the production of carotenoids. The production of carotenoids is important because they are used as natural colorants in food and some carotenoids are precursors of retinol (vitamin A). However, the identification and molecular characterization of the carotenoid pathway/s in species belonging to the genus Rhodotorula is scarce due to the lack of genomic information thus potentially impeding effective metabolic engineering of these yeast strains for improved carotenoid production. In this study, we report the isolation, identification, characterization and the whole nuclear genome and mitogenome sequence of the endophyte R. mucilaginosa RIT389 isolated from Distemonanthus benthamianus, a plant known for its anti-fungal and antibacterial properties and commonly used as chewing sticks. The assembled genome of R. mucilaginosa RIT389 is 19 Mbp in length with an estimated genomic heterozygosity of 9.29%. Whole genome phylogeny supports the species designation of strain RIT389 within the genus in addition to supporting the monophyly of the currently sequenced Rhodotorula species. Further, we report for the first time, the recovery of the complete mitochondrial genome of R. mucilaginosa using the genome skimming approach. The assembled mitogenome is at least 7,000 bases larger than that of Rhodotorula taiwanensis which is largely attributed to the presence of large intronic regions containing open reading frames coding for homing endonuclease from the LAGLIDADG and GIY-YIG families. Furthermore, genomic regions containing the key genes for carotenoid production were identified in R. mucilaginosa RIT389, revealing differences in gene synteny that may play a role in the regulation of the biotechnologically important carotenoid synthesis pathways in yeasts.
The tephritid fruit fly Zeugodacus tau (Walker) is a polyphagous fruit pest of economic importance in Asia. Studies based on genetic markers indicate that it forms a species complex. We report here (1) the complete mitogenome of Z. tau from Malaysia and comparison with that of China as well as the mitogenome of other congeners, and (2) the relationship of Z. tau taxa from different geographical regions based on sequences of cytochrome c oxidase subunit I gene. The complete mitogenome of Z. tau had a total length of 15631 bp for the Malaysian specimen (ZT3) and 15835 bp for the China specimen (ZT1), with similar gene order comprising 37 genes (13 protein-coding genes-PCGs, 2 rRNA genes, and 22 tRNA genes) and a non-coding A + T-rich control region (D-loop). Based on 13 PCGs and 15 mt-genes, Z. tau NC_027290 (China) and Z. tau ZT1 (China) formed a sister group in the lineage containing also Z. tau ZT3 (Malaysia). Phylogenetic analysis based on partial sequences of cox1 gene indicates that the taxa from China, Japan, Laos, Malaysia, Bangladesh, India, Sri Lanka, and Z. tau sp. A from Thailand belong to Z. tau sensu stricto. A complete cox1 gene (or 13 PCGs or 15 mt-genes) instead of partial sequence is more appropriate for determining phylogenetic relationship.
Diversifying selection between populations that inhabit different environments can promote lineage divergence within species and ultimately drive speciation. The mitochondrial genome (mitogenome) encodes essential proteins of the oxidative phosphorylation (OXPHOS) system and can be a strong target for climate-driven selection (i.e., associated with inhabiting different climates). We investigated whether Pleistocene climate changes drove mitochondrial selection and evolution within Australian birds. First, using phylogeographic analyses of the mitochondrial ND2 gene for 17 songbird species, we identified mitochondrial clades (mitolineages). Second, using distance-based redundancy analyses, we tested whether climate predicts variation in intraspecific genetic divergence beyond that explained by geographic distances and geographic position. Third, we analysed 41 complete mitogenome sequences representing each mitolineage of 17 species using codon models in a phylogenetic framework and a biochemical approach to identify signals of selection on OXPHOS protein-coding genes and test for parallel selection in mitolineages of different species existing in similar climates. Of 17 species examined, 13 had multiple mitolineages (range: 2-6). Climate was a significant predictor of mitochondrial variation in eight species. At least two amino acid replacements in OXPHOS complex I could have evolved under positive selection in specific mitolineages of two species. Protein homology modelling showed one of these to be in the loop region of the ND6 protein channel and the other in the functionally critical helix HL region of ND5. These findings call for direct tests of the functional and evolutionary significance of mitochondrial protein candidates for climate-associated selection.
The complete mitochondrial genome plays an important role in the research on phylogenetic relationship. Here, we reported the first complete mitochondrial genome sequence of Varuna yui Hwang & Takeda, 1986 (Varunidae). The complete mtDNA (15,915 bp in length) consisted of 13 protein-coding genes, 22 tRNAs, two rRNA genes, and a control region. The gene arrangement was identical to those observed in the Varunidae species. The phylogenetic analysis suggested that V. yui had close relationship with other Varunidae species (Helicetient sinensis, Eriocher sinesis, etc.). The newly described genome may facilitate further comparative mitogenomic analysis within Varunidae species.
In this study, we sequenced and analyzed the whole mitochondrial genome of Metopograpsus frontalis Miers, 1880 (Decapoda, Grapsidae). The circular genome is 15,587 bp in length, consisting of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, as well as a control region. Both atp8/atp6 and nad4L/nad4 share 7 nucleotides in their adjacent overlapping region, which is identical to those observed in other Grapsidae crabs. The genome composition and gene order follow a classic crab-type arrangement regulation. The phylogenetic analysis suggested that Grapsidae crabs formed a solid monophyletic group. The newly described mitochondrial genome may provide genetic marker for studies on phylogeny of the grapsid crabs.
The complete mitochondrial genome sequence of Atergatis integerrimus from China has been amplified and sequenced in this study. The mitogenome assembly was found to be 15,924 bp in length with base composition of A (32.88%), G (10.58%), C (20.87%), T (35.66%), A + T (68.54%), and G + C (31.46%). It contained 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and a control region. The phylogenetic position was constructed and the A. integerrimus was closely clustered with Pseudocarcinus gigas and Leptodius sanguineus. The complete mitochondrial genome sequence would be useful for further understanding the evolution of A. integerrimus.
The House Crow (Corvus splendens) is a useful study system for investigating the genetic basis of adaptations underpinning successful range expansion. The species originates from the Indian subcontinent, but has successfully spread through a variety of thermal environments across Asia, Africa and Europe. Here, population mitogenomics was used to investigate the colonisation history and to test for signals of molecular selection on the mitochondrial genome. We sequenced the mitogenomes of 89 House Crows spanning four native and five invasive populations. A Bayesian dated phylogeny, based on the 13 mitochondrial protein-coding genes, supports a mid-Pleistocene (~630,000 years ago) divergence between the most distant genetic lineages. Phylogeographic patterns suggest that northern South Asia is the likely centre of origin for the species. Codon-based analyses of selection and assessments of changes in amino acid properties provide evidence of positive selection on the ND2 and ND5 genes against a background of purifying selection across the mitogenome. Protein homology modelling suggests that four amino acid substitutions inferred to be under positive selection may modulate coupling efficiency and proton translocation mediated by OXPHOS complex I. The identified substitutions are found within native House Crow lineages and ecological niche modelling predicts suitable climatic areas for the establishment of crow populations within the invasive range. Mitogenomic patterns in the invasive range of the species are more strongly associated with introduction history than climate. We speculate that invasions of the House Crow have been facilitated by standing genetic variation that accumulated due to diversifying selection within the native range.
Angiostrongylus mackerrasae is a parasitic nematode of rats found in Australia. When first reported, it was referred to as A. cantonensis. Recent molecular studies, including the mitochondrial genome, indicate that it is highly similar to A. cantonensis. These studies did not include A. malaysiensis, another member of the A. cantonensis species complex, for comparison. The present study examined the genetic distance and phylogenetic relationship between the component taxa (A. cantonensis, A. mackerrasae and A. malaysiensis) of the A. cantonensis species complex, based on the 12 protein-coding genes (PCGs) of their mitochondrial genome. Both the nucleotide and amino acid sequences were analysed. Angiostrongylus mackerrasae and A. cantonensis are members of the same genetic lineage and both are genetically distinct from A. malaysiensis. The genetic distance based on concatenated nucleotide sequences of 12 mt-PCGs between A. mackerrasae and A. cantonensis from Thailand is p = 1.73%, while that between the Thai and Chinese taxa of A. cantonensis is p = 3.52%; the genetic distance between A. mackerrasae and A. cantonensis from China is p = 3.70%. The results indicate that A. mackerrasae and A. cantonensis belong to the same genetic lineage, and that A. mackerrasae may be conspecific with A. cantonensis. It remains to be resolved whether A. mackerrasae is conspecific with A. cantonensis or undergoing incipient speciation.