Molecular evolution and nucleotide diversity of nuclear plastid phosphoglycerate kinase (PGK) gene in Triticeae (Poaceae)

Levels of nucleotide divergence provide key evidence in the evolution of polyploids. The nucleotide diversity of 226 sequences of pgk1 gene in Triticeae species was characterized. Phylogenetic analyses based on the pgk1 gene were carried out to determine the diploid origin of polyploids within the tribe in relation to their A^u, B, D, St, Ns, P, and H haplomes. Sequences from the Ns genome represented the highest nucleotide diversity values for both polyploid and diploid species with π = 0.03343 and θ = 0.03536 for polyploid Ns genome sequences and π = 0.03886 and θ = 0.03886 for diploid Psathyrostachys sequences, while Triticum urartu represented the lowest diversity among diploid species at π = 0.0011 and θ = 0.0011. Nucleotide variation of diploid Aegilops speltoides (π = 0.2441, presumed the B genome donor of Triticum species) is five times higher than that (π = 0.00483) of B genome in polyploid species. Significant negative Tajima's D values for the St, A^u, and D genomes along with high rates of polymorphisms and low sequence diversity were observed. Origins of the A^u, B, and D genomes were linked to T. urartu, A. speltoides, and A. tauschii, respectively. Putative St genome donor was Pseudoroegneria, while Ns and P donors were Psathyrostachys and Agropyron. H genome diploid donor is Hordeum.     

Phylogeny and molecular evolution of the Acc1 gene within the StH genome species in Triticeae (Poaceae)

To estimate the phylogeny and molecular evolution of a single-copy gene encoding plastid acetyl-CoA carboxylase (Acc1) within the StH genome species, two Acc1 homoeologous sequences were isolated from nearly all the sampled StH genome species and were analyzed with those from 35 diploid taxa representing 19 basic genomes in Triticeae. Sequence diversity patterns and genealogical analysis suggested that (1) the StH genome species from the same areas or neighboring geographic regions are closely related to each other; (2) the Acc1 gene sequences of the StH genome species from North America and Eurasia are evolutionarily distinct; (3) Dasypyrum has contributed to the nuclear genome of Elymus repens and Elymus mutabilis; (4) the StH genome polyploids have higher levels of sequence diversity in the H genome homoeolog than the St genome homoeolog; and (5) the Acc1 sequence may evolve faster in the polyploid species than in the diploids. Our result provides some insight on evolutionary dynamics of duplicate Acc1 gene, the polyploidy speciation and phylogeny of the StH genome species.     

Phylogenetic and molecular dating analysis of Taiwan Blue Pheasant (Lophura swinhoii)

The Swinhoe's Pheasant (Lophura swinhoii) is an endemic and most endangered species to Taiwan, China. It belongs to the genus Lophura, family Phasianidae. To further investigate the evolutionary history of L. swinhoii, we determined its complete mitochondrial genome and reconstructed a single, robust phylogenetic tree. Our results showed that L. swinhoii is clustered with Lophura nycthemera and forms a sister group of Lophura ignita. The genus Lophura is strongly supported as the sister taxon of the genus Crossoptilon. The molecular clock analysis showed that the genetic divergence of L. swinhoii occurred in 2.71 (1.31–4.22) Mya. The most common ancestor of L. swinhoii might have migrated from mainland of South East Asia to Taiwan Island by the land bridge at 2.71 Mya ago. Taiwan Island is separated from the mainland by the sea (Taiwan Strait) and formed a separate island at around 2.5 Mya because of the transgression and regression. Therefore, geographical isolation and climate change may accelerate the evolution of L. swinhoii. In this study, we propose a biogeographic hypothesis for speciation of the L. swinhoii based on known events of the geographic and geological history of South East Asia and southeast China, which would benefit the understanding of evolutionary history of L. swinhoii as well as other galliform birds.     

Molecular characterization of Indian Sugarcane streak mosaic virus isolates reveals recombination and negative selection in the P1 gene

Sugarcane streak mosaic virus (SCSMV), a member of the genus Poacevirus is an important viral pathogen affecting sugarcane production in India. The P1 gene of ten Indian isolates was sequenced and compared with previously reported SCSMV isolates. Comparative sequence analysis revealed a high level of diversity in the P1 gene (83–98% nucleotide sequence identity; 87–100% amino acid sequence identity), and the Indian SCSMV isolates were found to be the most variable (up to 9% diversity at the amino acid level). Phylogenetic tree analysis showed clustering of 17 SCSMV isolates into two groups: group I included isolates from India (except SCSMV-TPT) and Pakistan, and group II consisted of isolates from Japan, Indonesia, Thailand and SCSMV-TPT. The results obtained from phylogenetic study were further supported by the different in silico analysis viz. SNPs (single nucleotide polymorphism), INDELs (insertion and deletion) and evolutionary distance analysis. A significant proportion of recombination sites were observed at the N terminal region of P1 gene. Analysis of selection pressure indicated that the P1 gene of the Indian SCSMV isolates is under strong negative or purifying selection. It is likely that recombination identified in Indian SCSMV isolates, along with strong purifying selection, enhances the speed of elimination of deleterious mutations in the P1 gene. The evolutionary processes (recombination and selection pressure) together contributed to the observed genetic diversity and population structure of Indian SCSMV isolates.     

Analysis of genetic variability and relationships among Mentha L. using the limonene synthase gene, LS

The genus Mentha comprises a group of aromatic plants with worldwide distribution. Because of frequent interspecific hybridization, the genetic relationships within the genus are not clearly understood. Limonene synthase, which catalyses the first committed step in the essential oil monoterpene biosynthetic pathway, is considered to be a possible rate limiting enzyme. With the homology-based cloning method, primers were designed according to cDNA sequence to amplify full-length DNA sequences in 13 Mentha samples from five species, using Perilla as an outgroup. Analyses of gene structure, length variation, GC-content, Ts/Tv ratio and evolutionary diversity were carried out. Consensus phylogenetic trees were obtained using maximum likelihood, neighbor-joining, and maximum parsimony, respectively, based on the full-length genomic DNA sequences, complete ORF coding sequences and predicted amino acid sequences. The results presented here based on the sequence of MhLS provide the first credibly supported genetic relationships for Mentha, which enables a basis for further mint taxonomy, cultivation and breeding.     

A mitogenomic phylogeny and genetic history of sable (Martes zibellina)

We assessed phylogeny of sable (Martes zibellina, Linnaeus, 1758) by sequence analysis of nearly complete, new mitochondrial genomes in 36 specimens from different localities in northern Eurasia (Primorye, Khabarovsk and Krasnoyarsk regions, the Kamchatka Peninsula, the Kuril Islands and the Urals). Phylogenetic analysis of mtDNA sequences demonstrates that two clades, A and BC, radiated about 200–300 thousand years ago (kya) according to results of Bayesian molecular clock and RelTime analyses of different mitogenome alignments (nearly complete mtDNA sequences, protein-coding region, and synonymous sites), while the age estimates of clades A, B and C fall within the Late Pleistocene (~ 50–140 kya). Bayesian skyline plots (BSPs) of sable population size change based on analysis of nearly complete mtDNAs show an expansion around 40 kya in the warm Karganian time, without a decline of population size around the Last Glacial Maximum (21 kya). The BSPs based on synonymous clock rate indicate that M. zibellina experienced demographic expansions later, approximately 22 kya. The A2a clade that colonized Kamchatka ~ 23–50 kya (depending on the mutation rate used) survived the last glaciation there as demonstrated by the BSP analysis. In addition, we have found evidence of positive selection acting at ND4 and cytochrome b genes, thereby suggesting adaptive evolution of the A2a clade in Kamchatka.     

Phylogenetic relationships of the Cobitoidea (Teleostei: Cypriniformes) inferred from mitochondrial and nuclear genes with analyses of gene evolution

The superfamily Cobitoidea of the order Cypriniformes is a diverse group of fishes, inhabiting freshwater ecosystems across Eurasia and North Africa. The phylogenetic relationships of this well-corroborated natural group and diverse clade are critical to not only informing scientific communities of the phylogeny of the order Cypriniformes, the world's largest freshwater fish order, but are key to every area of comparative biology examining the evolution of traits, functional structures, and breeding behaviors to their biogeographic histories, speciation, anagenetic divergence, and divergence time estimates. In the present study, two mitochondrial gene sequences (COI, ND4+5) and four single-copy nuclear gene segments (RH1, RAG1, EGR2B, IRBP) were used to infer the phylogenetic relationships of the Cobitoidea as reconstructed from maximum likelihood (ML) and partitioned Bayesian Analysis (BA). Analyses of the combined mitochondrial/nuclear gene datasets revealed five strongly supported monophyletic Cobitoidea families and their sister-group relationships: Botiidae+(Vaillantellidae+(Cobitidae+(Nemacheilidae+Balitoridae))). These recovered relationships are in agreement with previous systematic studies on the order Cypriniformes and/or those focusing on the superfamily Cobitoidea. Using these relationships, our analyses revealed pattern lineage- or ecological-group-specific evolution of these genes for the Cobitoidea. These observations and results corroborate the hypothesis that these group-specific-ancestral ecological characters have contributed in the diversification and/or adaptations within these groups. Positive selections were detected in RH1 of nemacheilids and in RAG1 of nemacheilids and genus Vaillantella, which indicated that evolution of RH1 (related to eye's optic sense) and RAG1 (related to immunity) genes appeared to be important for the diversification of these groups. The balitorid lineage (those species inhabiting fast-flowing riverine habitats) had, as compared with other cobitoid lineages, significantly different dN/dS, dN and dS values for ND4 and IRBP genes. These significant differences are usually indicative of weaker selection pressure, and lineage-specific evolution on genes along the balitorid lineage. Furthermore, within Cobitoidea, excluding balitorids, species living in subtropics had significantly higher dN/dS values in RAG1 and IRBP genes than those living in temperate and tropical zones. Among tropical cobitoids, genes COI, ND5, EGR2B, IRBP and RH1, had a significantly higher mean dS value than those species in subtropical and temperate groups. These findings suggest that the evolution of these genes could also be ecological-group-specific and may have played an important role in the adaptive evolution and diversification of these groups. Thus, we hypothesize that the genes included in the present study were actively involved in lineage- and/or ecological-group-specific evolutionary processes of the highly diverse Cobitoidea. These two evolutionary patterns, both subject to further testing, are hypothesized as integral in the diversification with this major clade of the world's most diverse group of freshwater fishes.     

Effects of garlic on Schistosoma mansoni harbored in albino mice: Molecular characterization of the host and parasite

Garlic has been used for its health benefits for thousands of years. Modern research confirmed many of the healing properties of garlic, including its antiparasitic activity. This study was designed to evaluate the antischistosomal action of garlic through detecting the changes in DNA profile of Schistosoma mansoni worms and the infected mouse. Forty mice were subcutaneously infected with ~ 200 Schistosoma mansoni cercariae/mouse. Infected mice were divided into four equal groups: non-treated, prophylactic, therapeutic, and continuously-treated. Non-infected control and garlic-treated groups were assigned for the sake of comparison. Garlic extract (50 mg/kg bw/mouse) was given orally, day after day, at a fixed daytime. Seven weeks post-infection, adult schistosomes were recovered by perfusion and the livers of the mice were excised out and were processed for DNA extraction and Random Amplification of Polymorphic DNA-Polymerase Chain Reaction (RAPD-PCR). The results showed that garlic exerted no major changes in the genome of schistosomes. Nevertheless, that schistosomal infection induced genetic alterations in the DNA of mice, and garlic was able to ameliorate such alterations to a great extent.     

The complete mitochondrial genome sequence analysis of Tibetan argali (Ovis ammon hodgsoni): Implications of Tibetan argali and Gansu argali as the same subspecies

The genus Ovis (Bovidae, Artiodactyla) includes six species, i.e. Ovis ammon, Ovis aries, Ovis canadensis, Ovis dalli, Ovis nivicola and Ovis vignei. Based on morphology, geographical location, habitat, etc., the species O. ammon is divided into nine subspecies. The near threatened Tibetan argali is distributed across the Tibetan Plateau and its peripheral mountains, and believed to be one of the O. ammon subspecies (O. a. hodgsoni). However, considering its morphological features and distributions, a question has been proposed by some researchers about the subspecies status of Tibetan argali. In this study, we employed complete mitochondrial DNA (mtDNA) to explore the phylogenetic relationship and population genetic structure of Tibetan argali. The results revealed that the nucleotide composition, gene arrangement and codon usage pattern of the mitochondrial genome of Tibetan argali are similar to those of other caprines. Phylogenetic analyses showed that Tibetan argali was clustered with O. ammon. Interestingly, five Tibetan argali individuals and one of the three Gansu argali (O. a. dalailamae) individuals were clustered in the same branch, which is a sister group to other two Gansu argali individuals. Together with morphological characteristics, our results suggested that Tibetan argali and Gansu argali may belong to the same subspecies (O. a. hodgsoni) of O. ammon, rather than two different subspecies.     

Plant degreening: evolution and expression of tomato (Solanum lycopersicum) dephytylation enzymes

Chlorophyll is the most abundant pigment on earth and even though it is known that its high photo-excitability necessitates a tight regulation of its degradation pathway, to date there are still several steps in chlorophyll breakdown that remain obscure. In order to better understand the ‘degreening’ processes that accompany leaf senescence and fruit ripening, we characterized the enzyme-encoding genes involved in dephytylation from tomato (Solanum lycopersicum). A single pheophytinase (PPH) gene and four chlorophyllase (CLH) genes were identified in the tomato genome. A phenetic analysis revealed two groups of CLHs in eudicot species and further evolutionary analysis indicated that these enzymes are under diverse selection pressures. A comprehensive expression profile analysis also suggested functional specificity for these dephytylating enzymes. The integrated analysis allows us to propose three general roles for chlorophyll dephytylation: i) PPH, which is under high selective constraint, is responsible for chlorophyll degradation during developmentally programed physiological processes; ii) Group I CLHs, which are under relaxed selection constraint, respond to environmental and hormonal stimuli and play a role in plant adaptation plasticity; and iii) Group II CLHs, which are also under high selective constraint, are mostly involved in chlorophyll recycling.     

Complete nucleotide sequence and gene organization of the mitochondrial genome of Paa spinosa (Anura: Ranoidae)

The mt genome of Paa spinosa (Anura: Ranoidae) is a circular molecule of 18,012 bp in length, containing 38 genes (including an extra copy of tRNA-Met gene). This mt genome is characterized by three distinctive features: a cluster of rearranged tRNA genes (LTPF tRNA gene cluster), a tandem duplication of tRNA-Met gene (Met1 and Met2), and distinct repeat regions at both 5′ and 3′-sides in the control region. Comparing the locations and the sequences of all tRNA-Met genes among Ranoidae, and constructing NJ tree of the nucleotide of those tRNA-Met genes, we suggested a tandem duplication of tRNA-Met gene can be regarded as a synapomorphy of Dicroglossinae. To further investigate the phylogenetic relationships of anurans, phylogenetic analyses (BI, ML and MP) based on the nucleotide dataset and the corresponding amino acid dataset of 11 protein-coding genes (except ND5 and ATP8) arrived at the similar topology.     

The mitochondrial genome of Ifremeria nautilei and the phylogenetic position of the enigmatic deep-sea Abyssochrysoidea (Mollusca: Gastropoda)

The complete nucleotide sequence of the mitochondrial (mt) genome of the deep-sea vent snail Ifremeria nautilei (Gastropoda: Abyssochrysoidea) was determined. The double stranded circular molecule is 15,664 pb in length and encodes for the typical 37 metazoan mitochondrial genes. The gene arrangement of the Ifremeria mt genome is most similar to genome organization of caenogastropods and differs only on the relative position of the trnW gene. The deduced amino acid sequences of the mt protein coding genes of Ifremeria mt genome were aligned with orthologous sequences from representatives of the main lineages of gastropods and phylogenetic relationships were inferred. The reconstructed phylogeny supports that Ifremeria belongs to Caenogastropoda and that it is closely related to hypsogastropod superfamilies. Results were compared with a reconstructed nuclear-based phylogeny. Moreover, a relaxed molecular-clock timetree calibrated with fossils dated the divergence of Abyssochrysoidea in the Late Jurassic–Early Cretaceous indicating a relatively modern colonization of deep-sea environments by these snails.     

Hemocyanin gene family evolution in spiders (Araneae), with implications for phylogenetic relationships and divergence times in the infraorder Mygalomorphae

Hemocyanins are multimeric copper-containing hemolymph proteins involved in oxygen binding and transport in all major arthropod lineages. Most arachnids have seven primary subunits (encoded by paralogous genes a–g), which combine to form a 24-mer (4 × 6) quaternary structure. Within some spider lineages, however, hemocyanin evolution has been a dynamic process with extensive paralog duplication and loss. We have obtained hemocyanin gene sequences from numerous representatives of the spider infraorders Mygalomorphae and Araneomorphae in order to infer the evolution of the hemocyanin gene family and estimate spider relationships using these conserved loci. Our hemocyanin gene tree is largely consistent with the previous hypotheses of paralog relationships based on immunological studies, but reveals some discrepancies in which paralog types have been lost or duplicated in specific spider lineages. Analyses of concatenated hemocyanin sequences resolved deep nodes in the spider phylogeny and recovered a number of clades that are supported by other molecular studies, particularly for mygalomorph taxa. The concatenated data set is also used to estimate dates of higher-level spider divergences and suggests that the diversification of extant mygalomorphs preceded that of extant araneomorphs. Spiders are diverse in behavior and respiratory morphology, and our results are beneficial for comparative analyses of spider respiration. Lastly, the conserved hemocyanin sequences allow for the inference of spider relationships and ancient divergence dates.     

Molecular characterization of two y-type high molecular weight glutenin subunit alleles 1Ay12 and 1Ay8 from cultivated einkorn wheat (Triticum monococcum ssp. monococcum)

Two y-type high molecular weight glutenin subunits (HMW-GSs) 1Ay12^? and 1Ay8^? from the two accessions PI560720 and PI345186 of cultivated einkorn wheat (Triticum monococcum ssp. monococcum, AA, 2n = 2x = 14), were identified by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). The mobility of 1Ay12^? and 1Ay8^? was similar to that of 1Dy12 and 1By8 from common wheat Chinese Spring, respectively. Their ORFs respectively consisted of 1812 bp and 1935 bp, encoding 602 and 643 amino acid residues with the four typical structural domains of HMW-GS including signal peptide, conserved N-, and C-terminal and central repetitive domains. Compared with the most similar active 1Ay alleles previous published, there were a total of 15 SNPs and 2 InDels in them. Their encoding functions were confirmed by successful heterogeneous expression. The two novel 1Ay alleles were named as 1Ay12^? and 1Ay8^? with the accession No. JQ318694 and JQ318695 in GenBank, respectively. The two alleles were classed into the two distinct groups, Phe-type and Cys-type, which might be relevant to the differentiation of Glu-A1-2 alleles. Of which, 1Ay8^? belonged to Cys-type group, and its protein possessed an additional conserved cysteine residue in central repetitive region besides the six common ones in N- and C-terminal regions of Phe-type group, and was the second longest in all the known active 1Ay alleles. These results suggested that the subunit 1Ay8^? of cultivated einkorn wheat accession PI345186 might have a potential ability to strengthen the gluten polymer interactions and be a valuable genetic resource for wheat quality improvement.     

The complete mitochondrial genome of Sasakia funebris (Leech) (Lepidoptera: Nymphalidae) and comparison with other Apaturinae insects

Sasakia funebris, a member of the lepidopteran family, Nymphalidae (superfamily Papilionoidea) is a rare species and is found only in some areas of South China. In this study, the 15,233 bp long complete mitochondrial genome of S. funebris was determined, and harbors the gene arrangement identical to all other sequenced lepidopteran insects. The nucleotide composition of the genome is highly A + T biased, accounting for 81.2%. All protein-coding genes (PCGs) start with typical ATN codons, except for COI which begins with the CGA codon. All tRNAs have a typical clover-leaf secondary structure, except for tRNA^Ser(AGN), the dihydrouridine (DHU) arm of which forms a simple loop. The S. funebris A + T-rich region of 370 bp contains several features common to the Lepidoptera insects, including the motif ATAGA followed by a 19 bp poly-T stretch, and two tandem repeats consisting of 18 bp repeat units and 14 bp repeat units. The phylogenetic analyses of Apaturinae based on mitogenome sequences showed: (S. funebris + Sasakia charonda) + (Apatura metis + Apatura ilia). This result is consistent with the morphological classification.     

Phylogenetic analysis of questionable tetraploid species in Roegneria and Pseudoroegneria (Poaceae: Triticeae) inferred from a gene encoding plastid acety1-CoA carboxylase

To evaluate the phylogenetic relationships of questionable tetraploid species Roegneria alashanica Keng, Roegneria magnicaespes (D.F. Cui) L.B. Cai, Roegneria elytrigioides C. Yen et J.L. Yang, Roegneria grandis Keng and Pseudoroegneria geniculata (Trin.) Á. Löve, the single copy sequences of the plastid acetyl-CoA carboxylase gene (Acc1) were analyzed among the five species and the related diploid and tetraploid species. The results indicated that: (a) R. alashanica contained one set of modified St genome which was closely related to the E^e genome, and the other set of genome was closely related to the P genome; (b) R. magnicaespes contained one set of St genome, the other set of genome might be closely related to the P genome. There are close affinities between R. magnicaespes and R. alashanica; (c) R. elytrigioides contained two sets of St genomes, and it is reasonable to be treated as Pseudoreogneria elytrigioides (C. Yen et J.L. Yang) B.R. Lu; (d) the genome of R. grandis should be designed as St^gY. The St^g genome was a differentiated form of the St genome in Pseudoroegneria and was homoeologous with the Y genome in Roegneria; (e) the genomic constitution of P. geniculata was similar to that of R. magnicaespes and R. alashanica and distinctly related to P. geniculata ssp. scythica (E^eSt). They should be treated as different species in different genera; and (f) the Y genome was possibly originated from the St genome, and was sister to the St, E^e, E^b and W genomes.     

Molecular cloning and characterization of an Hsp90/70 organizing protein gene from Frankliniella occidentalis (Insecta: Thysanoptera,Thripidae)

The heat shock 90/70 organizing protein (Hop), also known as Sti-1 (stress-induced protein-1), is a co-chaperone that usually mediates the interaction of Hsp90 and Hsp70 and has been extensively characterized in mammals and plants. However, its role in insects remains unknown. In the present study, we isolated and characterized a Hop homologue gene from Frankliniella occidentalis (Fohop). The Fohop contains a 1659 bp ORF encoding a protein of 552 amino acids with a caculated molecular mass of approximately 62.25 kDa, which displays a reasonable degree of identity with the known Hops and shares several canonical motifs, including three tetratricopeptide repeated motif domains (TPR1, TPR2A and TPR2B) and two aspartic acid–proline (DP) repeat motifs (DP1 and DP2). As in other hops, Fohop contains introns, but the number and the position are quite variable. The mRNA expression patterns indicated that Fohop was constitutively expressed throughout the developmental stages, but was obviously upregulated by heat stress both in larvae and adults. Our studies imply that Hop, as in other Hsps, may play an important role in heat shock response of F. occidentalis.     

Complete mitochondrial genome sequences of three bats species and whole genome mitochondrial analyses reveal patterns of codon bias and lend support to a basal split in Chiroptera

Order Chiroptera is a unique group of mammals whose members have attained self-powered flight as their main mode of locomotion. Much speculation persists regarding bat evolution; however, lack of sufficient molecular data hampers evolutionary and conservation studies. Of ~ 1200 species, complete mitochondrial genome sequences are available for only eleven. Additional sequences should be generated if we are to resolve many questions concerning these fascinating mammals. Herein, we describe the complete mitochondrial genomes of three bats: Corynorhinus rafinesquii, Lasiurus borealis and Artibeus lituratus. We also compare the currently available mitochondrial genomes and analyze codon usage in Chiroptera. C. rafinesquii, L. borealis and A. lituratus mitochondrial genomes are 16438 bp, 17048 bp and 16709 bp, respectively. Genome organization and gene arrangements are similar to other bats. Phylogenetic analyses using complete mitochondrial genome sequences support previously established phylogenetic relationships and suggest utility in future studies focusing on the evolutionary aspects of these species. Comprehensive analyses of available bat mitochondrial genomes reveal distinct nucleotide patterns and synonymous codon preferences corresponding to different chiropteran families. These patterns suggest that mutational and selection forces are acting to different extents within Chiroptera and shape their mitochondrial genomes.     

The complete mitochondrial genome sequence of the world's largest fish,the whale shark (Rhincodon typus), and its comparison with those of related shark species

The whale shark (Rhincodon typus) is the largest extant species of fish, belonging to the order Orectolobiformes. It is listed as a “vulnerable” species on the International Union for Conservation of Nature (IUCN)'s Red List of Threatened Species, which makes it an important species for conservation efforts. We report here the first complete sequence of the mitochondrial genome (mitogenome) of the whale shark obtained by next-generation sequencing methods. The assembled mitogenome is a 16,875 bp circle, comprising of 13 protein-coding genes, two rRNA genes, 22 tRNA genes and a control region. We also performed comparative analysis of the whale shark mitogenome to the available mitogenome sequences of 17 other shark species, four from the order Orectolobiformes, five from Lamniformes and eight from Carcharhiniformes. The nucleotide composition, number and arrangement of the genes in whale shark mitogenome are the same as found in the mitogenomes of the other members of the order Orectolobiformes and its closest orders Lamniformes and Carcharhiniformes, although the whale shark mitogenome had a slightly longer control region. The availability of mitogenome sequence of whale shark will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this species.