Assembly and phylogenomic analysis of cotton mitochondrial genomes provide insights into the history of cotton evolution

Cotton is a major crop that provides the most important renewable textile fibers in the world.Studies of the taxonomy and evolution of cotton species have received wide attentions,not only due to cotton’s economic value but also due to the fact that Gossypium is an ideal model system to study the origin,evolution,and cultivation of polyploid species.Previous studies suggested the involvement of mitochondrial genome editing sites and copy number as well as mitochondrial functions in cotton fiber elongation.Whereas,with only a few mitogenomes assembled in the cotton genus Gossypium,our knowledge about their roles in cotton evolution and speciation is still scarce.To close this gap,here we assembled 20 mitogenomes from 15 cotton species spanning all the cotton clades(A–G,K,and AD genomes)and 5 cotton relatives using short and long sequencing reads.Systematic analyses uncovered a high level of mitochondrial gene sequence conservation,abundant sequence repeats and many insertions of foreign sequences,as well as extensive structural variations in cotton mitogenomes.The sequence repeats and foreign sequences caused significant mitogenome size inflation in Gossypium and its close relative Kokia in general,while there is no significant difference between the lint and fuzz cotton mitogenomes in terms of gene content,RNA editing,and gene expression level.Interestingly,we further revealed the specific presence and expression of two novel mitochondrial open reading frames(ORFs)in lint-fiber cotton species.Finally,these structural features and novel ORFs help us gain valuable insights into the history of cotton evolution and polyploidization and the origin of species producing long lint fibers from a mitogenomic perspective.

Molecular phylogenetic relationships based on mitochondrial genomes of novel deep-sea corals(Octocorallia:Alcyonacea):Insights into slow evolution and adaptation to extreme deep-sea environments

A total of 10 specimens of Alcyonacea corals were collected at depths ranging from 905 m to 1633 m by the manned submersible Shenhai Yongshi during two cruises in the South China Sea(SCS).Based on mitochondrial genomic characteristics,morphological examination,and sclerite scanning electron microscopy,the samples were categorized into four suborders(Calcaxonia,Holaxonia,Scleraxonia,and Stolonifera),and identified as 9 possible new cold-water coral species.Assessments of GC-skew dissimilarity,phylogenetic distance,and average nucleotide identity(ANI)revealed a slow evolutionary rate for the octocoral mitochondrial sequences.The nonsynonymous(Ka)to synonymous(Ks)substitution ratio(Ka/Ks)suggested that the 14 protein-coding genes(PCGs)were under purifying selection,likely due to specific deep-sea environmental pressures.Correlation analysis of the median Ka/Ks values of five gene families and environmental factors indicated that the genes encoding cytochrome b(cyt b)and DNA mismatch repair protein(mutS)may be influenced by environmental factors in the context of deep-sea species formation.This study highlights the slow evolutionary pace and adaptive mechanisms of deep-sea corals.

Gene characterization and phylogenetic analysis of four mitochondrial genomes in Caenogastropoda

Caenogastropoda is a highly diverse group,containing~60%of all existing gastropods.Species in this subclass predominantly inhabit marine environments and have a high ecological and economic value.Owing to the increase in relevant phylogenetic studies,our understanding of between species relatedness in Caenogastropoda has improved.However,the biodiversity,taxonomic status,and phylogenetic relationships of this group remain unclear.In the present study,we performed next-generation sequencing of four complete mitochondrial genomes from three families(Buccinidae,Columbellidae,and Cypraeidae)and the four mitogenomes were classical circular structures,with a length of 16177 bp in Volutharpa ampullacea,16244 bp in Mitrella albuginosa,16926bp in Mauritia arabica asiatica and 15422 bp in Erronea errones.Base composition analysis indicated that whole sequences were biased toward A and T.Then compared them with 171 complete mitochondrial genomes of Caenogastropoda.The phylogenetic relationship of Caenogastropoda derived from Maximum Likelihood(ML)and Bayesian Inference(BI)trees constructed based on CDS sequences was consistent with the results of traditional morphological analysis,with all three families showing close relationships.This study supported Caenogastropoda at the molecular level as a separate clade of Mollusca.According to our divergence time estimations,Caenogastropoda was formed during the middle Triassic period(~247.2–237 Ma).Our novel mitochondrial genomes provide evidence for the speciation of Caenogastropoda in addition to elucidating the mitochondrial genomic evolution of this subclass.

Mitochondrial phylogenomics reveal the origin and adaptive evolution of the deep-sea caridean shrimps(Decapoda:Caridea)

The deep-sea is considered as the most extensive ecosystem on the Earth.It is meaningful for elucidating the life origins by exploring the origin and adaptive genetic mechanisms of the large deepsea organisms.Caridean shrimps have colonized and successfully adapted to deep-sea environments.They provide an ideal model to analyze the origin and adaptive evolution of modern deep-sea fauna.Here,we conducted the phylogenetic analyses of mitocho ndrial genomes(mitogenomes)from carideans,including 11 newly sequences reported in this investigation to explore the habitat origins,divergence times,and adaptive evolution of deep-sea(seamounts and hydrothermal vents)caridean shrimps.The results showed that the species of deep-sea Caridea formed a monophyletic group.Phylogenetic analysis supported that the deepsea caridean shrimps may originated from shallow sea.The hydrothermal vents alvinocaridid shrimps and Lebbeus shinkaiae from Thoridae underwent a second range expansion from seamounts to vent ecosystems.Estimates of divergence time showed that the caridean shrimps invaded into deep-sea at 147.75 Ma.The divergence of most of the modern seamount and hydrothermal vent species are in the late Cretaceous/early Tertiary.This may associate with the geological events of the Western Pacific,the climate change,and the global deep-water anoxic/dysoxic events during this period.Twenty-two potentially important adaptive residues were detected in the deep-sea shrimp lineage,which were located in atp6,atp8,cox1,cox3,cytb,nad2,nad4 l,and nad5.This investigation adds our understanding of the evolutionary history of deep-sea caridean shrimps,and provides insights into the mitochondrial genetic basis of deep-sea adaptation in this group.

Mitochondrial phylogenomics provides insights into the phylogeny and evolution of spiders(Arthropoda:Araneae)

Spiders are among the most varied terrestrial predators,with highly diverse morphology,ecology,and behavior.Morphological and molecular data have greatly contributed to advances in the phylogeny and evolutionary dynamics of spiders.Here,we performed comprehensive mitochondrial phylogenomics analysis on 78 mitochondrial genomes(mitogenomes)representing 29 families;of these,23 species from eight families were newly generated.Mesothelae retained the same gene arrangement as the arthropod ancestor(Limulus polyphemus),while Opisthothelae showed extensive rearrangement,with 12 rearrangement types in transfer RNAs(tRNAs)and control region.Most spider tRNAs were extremely truncated and lacked typical dihydrouridine or TΨC arms,showing high tRNA structural diversity;in particular,trnS1 exhibited anticodon diversity across the phylogeny.The evolutionary rates of mitochondrial genes were potentially associated with gene rearrangement or truncated tRNAs.Both mitogenomic sequences and rearrangements possessed phylogenetic characteristics,providing a robust backbone for spider phylogeny,as previously reported.The monophyly of suborder,infraorder,retrolateral tibial apophysis clade,and families(except for Pisauridae)was separately supported,and high-level relationships were resolved as(Mesothelae,(Mygalomorphae,(Entelegynae,(Synspermiata,Hypochilidae)))).The phylogenetic positions of several families were also resolved(e.g.,Eresidae,Oecobiidae and Titanoecidae).Two reconstructions of ancestral web type obtained almost identical results,indicating that the common ancestor of spiders likely foraged using a silk-lined burrow.This study,the largest mitochondrial phylogenomics analysis of spiders to date,highlights the usefulness of mitogenomic data not only for providing efficient phylogenetic signals for spider phylogeny,but also for characterizing trait diversification in spider evolution.

The complete mitochondrial genome of the grey bamboo shark (Chiloscyllium griseum) (Orectolobiformes: Hemiscylliidae): genomic characterization and phylogenetic application

The complete mitochondrial genome of the grey bamboo shark, Chiloscyllium griseum （Orectolobiformes： HemiscyUiidae） was obtained, with a total of 16 755 base pairs （bp）. The mitogenome encoded 13 protein- coding genes, two ribosomal RNAs, 22 transfer RNAs and a noncoding control region, a gene arrangement identical to the most common type found in vertebrates. All tRNA sequences were folded into typical clover- leaf secondary structure with the exception of tRNA-Ser2, in which the DHU arm stem was replaced with 12 unpaired nucleotides and formed a simple loop. In the control region, the putative termination-associated sequences （TAS） with hairpin-loop structure and the conserved sequence blocks （CSB） I-III were identified and considered to be associated with the replication and transcription of mtDNA. No repeat motifs were found in the mitogenome of C. griseum. On the basis of the mitogenomic data available in the Class Chon- drichthyes from both Genbank and this study, and the phylogenetic analyses using maximum likelihood （ML） and Bayesian inference （BI） methods, the results strongly support the basal division between batoids and sharks and the monophyly of the Superorder Galeomorphii, and confirm the phylogenetic position of C. griseum in the Subclass Elasmobranchii. This study suggests that the mitogenomic data are more robust for reconstructing phylogenetic relationships than individual genes in the Elasmobranchii.

Assembly and analysis of the Populus deltoides mitochondrial genome:the first report of a multicircular mitochondrial conformation for the genus Populus

Genomics research of Populus deltoides,an important timber species that is widely planted worldwide,is an important part of poplar breeding.Currently,the nuclear and chloroplast genome of P.deltoides have been sequenced,but its mitochondrial genome(mitogenome)has not been reported.To further explore the evolution and phylogeny of P.deltoides,the mitogenome of P.deltoides I-69 was assembled using reads from Nanopore and Illumina sequencing platforms and found to consist of 802,637 bp and three circular chromosomes(336,205,280,841,and 185,591 bp)containing 58 genes(34 protein-coding genes,21 tRNA genes,and 3 rRNA genes).RNA analysis in combination with several species showed signifi cantly fewer RNA editingsites in the mitogenomes of poplar and other angiosperms than in gymnosperms.Sequence transfer analysis showed extensive mitogenome rearrangements in Populus species,and with evolution from lower to higher plants,tRNA transfer from chloroplasts to mitochondria became increasingly frequent.In a phylogenetic analysis,the evolutionary status of P.deltoides was determined,and the section Populus was supported.Our results based on the fi rst report of a multicircular conformation of the Populus mitogenome provide a basis for further study of the evolution and genetics of P.deltoides and other Populus species and for breeding programs.

Evolutionary timescale of chalcidoid wasps inferred from over one hundred mitochondrial genomes

Chalcidoidea is one of the most biologically diverse groups among Hymenoptera.Members are characterized by extraordinary parasitic lifestyles and extensive host ranges,among which several species attack plants or serve as pollinators.However,higher-level chalcidoid relationships remain controversial.Here,we performed mitochondrial phylogenomic analyses for major clades(18out of 25 families)of Chalcidoidea based on 139 mitochondrial genomes.The compositional heterogeneity and conflicting backbone relationships in Chalcidoidea were assessed using various datasets and tree inferences.Our phylogenetic results supported the monophyly of 16families and polyphyly of Aphelinidae and Pteromalidae.Our preferred topology recovered the relationship(Mymaridae+(Signiphoridae+Leucospidae)+(Chalcididae+((Perilampidae+Eucharitidae)+remaining Chalcidoidea))).The monophyly of Agaonidae and Sycophaginae was rejected,while the gall-associated((Megastigmidae+Ormyridae)+(Ormocerinae+Eurytomidae))relationship was supported in most results.A six-gene inversion may be a synapomorphy for most families,whereas other derived gene orders may introduce confusion in phylogenetic signals at deeper nodes.Dating estimates suggested that Chalcidoidea arose near the Jurassic/Cretaceous boundary and that two dynamic shifts in diversification occurred during the evolution of Chalcidoidea.We hypothesized that the potential codiversification between chalcidoids and their hosts may be crucial for accelerating the diversification of Chalcidoidea.Ancestral state reconstruction analyses supported the hypothesis that gallinducers were mainly derived from parasitoids of gallinducers,while other gall-inducers were derived from phytophagous groups.Taken together,these findings advance our understanding of mitochondrial genome evolution in the major interfamilial phylogeny of Chalcidoidea.

Characterization of Rheum palmatum mitochondrial genome and comparative analysis among Caryophyllales species

Objective:This work aimed to report the first complete mitochondrial genome(mitogenome)of Rheum palmatum,summarize the features of Caryophyllales mitogenomes,and to reveal the potential of utilizing the mitogenomes of R.palmatum and other Caryophyllales species for inferring phylogenetic relationships and species identification.Methods:Both Illumina short reads and PacBio HiFi reads were utilized to obtain a complete mitogenome of R.palmatum.A variety of bioinformatics tools were employed to characterize the R.palmatum mitogenome,compare the reported mitogenomes in Caryophyllales and conduct phylogenetic analysis.Results:The mitogenome of R.palmatum was assembled into a single master circle of 302,993 bp,encoding 35 known protein-coding genes,18 transfer RNA genes,and three ribosome RNA genes.A total of 249 long repeats and 49 simple sequence repeats were identified in this mitogenome.The sizes of mitogenomes in Caryophyllales varied from 253 kb to 11.3 Mb.Among them,23 mitogenomes were circular molecules,one was linear,and one consisted of relaxed circles,linear molecules,and supercoiled DNA.Out of the total mitogenomes,11 were single-chromosome structure,whereas the remaining 14 were multi-chromosomal organizations.The phylogenetic analysis is consistent with both the Engler system(1964)and the Angiosperm Phylogeny Group III system.Conclusions:We obtained the first mitogenome of R.palmatum,which consists of a master circle.Mitogenomes in Caryophyllales have variable genome sizes and structures even within the same species.Circular molecules are still the dominant pattern in Caryophyllales.Single-chromosome mitogenomes account for nearly a half of all the mitogenomes in Caryophyllales,in contrast to previous studies.It is feasible to utilize mitochondrial genomes for inferring phylogenetic relationships and conducting species identification.

Mitogenomic phylogeny of the Asian colobine genus Trachypithecus with special focus on Trachypithecus phayrei(Blyth, 1847) and description of a new species

Trachypithecus,which currently contains 20 species divided into four groups,is the most speciose and geographically dispersed genus among Asian colobines.Despite several morphological and molecular studies,however,its evolutionary history and phylogeography remain poorly understood.Phayre’s langur(Trachypithecus phayrei) is one of the most widespread members of the genus,but details on its actual distribution and intraspecific taxonomy are limited and controversial.Thus,to elucidate the evolutionary history of Trachypithecus and to clarify the intraspecific taxonomy and distribution of T.phayrei,we sequenced 41 mitochondrial genomes from georeferenced fecal samples and museum specimens,including two holotypes.Phylogenetic analyses revealed a robustly supported phylogeny of Trachypithecus,suggesting that the T.pileatus group branched first,followed by the T.francoisi group,and the T.cristatus and T.obscurus groups most recently.The four species groups diverged from each other 4.5-3.1 million years ago(Ma),while speciation events within these groups occurred much more recently(1.6-0.3 Ma).Within T.phayrei,we found three clades that diverged 1.0-0.9 Ma,indicating the existence of three rather than two taxa.Following the phylogenetic species concept and based on genetic,morphological,and ecological differences,we elevate the T.phayrei subspecies to species level,describe a new species from central Myanmar,and refine the distribution of the three taxa.Overall,our study highlights the importance of museum specimens and provides new insights not only into the evolutionary history of T.phayrei but the entire Trachypithecus genus as well.

The complete mitochondrial genome of the Keeled box turtle Pyxidea mouhotii and phylogenetic analysis of major turtle groups

The complete mitochondrial genome （16,837 bp） from the Keeled box turtle （Pyxidea mouhotii） was determined. The genome content, gene order, and base composition conformed to the consensus vertebrate type mtDNA. However, a remarkable feature was found in this molecule： a large number of （ATTATATC） n direct tandem repeats followed by （TA） n microsatellite at the 3＇ end of the control region （D-loop）, which might be useful as molecular markers for studying population genetics and helpful for species identification and conservation. Besides, to review phylogenetic relationships among major turtle lineages, maximum-likelihood （ML） and Bayesian （BI） analyses were conducted based on concatenated sequences of 13 protein-coding genes from 16 taxa. The resultant ML and BI analyses showed homological topologies, which only differed on the exact placement of Platysternon. Nevertheless, the results strongly supported that 1） Pyxidea mouhotii and Cuora aurocapitata formed a monophyletic clade, whereas Cyclemys atripons was not closer to the Pyxidea-Cuora than to Chinemys reevesii, suggesting that Cyclemys and the Cuora group （containing Pyxidea） may have originated from two ancestors; 2） the Geoemydidae with Testudinidae was a sister group rather than with the Emydidae.

The first mitochondrial genome for the butterfly family Riodinidae(Abisara fylloides) and its systematic implications

The Riodinidae is one of the lepidopteran butterfly families. This study describes the complete mitochondrial genome of the butterfly species Abisara fylloides, the first mitochondrial genome of the Riodinidae family. The results show that the entire mitochondrial genome of A. fylloides is 15301 bp in length, and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and a 423 bp A＋T-rich region. The gene content, orientation and order are identical to the majority of other lepidopteran insects. Phylogenetic reconstruction was conducted using the concatenated 13 protein-coding gene （PCG） sequences of 19 available butterfly species covering all the five butterfly families （Papilionidae, Nymphalidae, Peridae, Lycaenidae and Riodinidae）. Both maximum likelihood and Bayesian inference analyses highly supported the monophyly of Lycaenidae＋Riodinidae, which was standing as the sister of Nymphalidae. In addition, we propose that the riodinids be categorized into the family Lycaenidae as a subfamilial taxon.

Characterization of the Complete Mitochondrial Genome of the Hybrid Epinephelus moara♀ × Epinephelus lanceolatus♂,and Phylogenetic Analysis in Subfamily Epinephelinae

This study presents the complete mitochondrial genome of the hybrid Epinephelus moara♀× Epinephelus lanceolatus♂. The genome is 16886 bp in length, and contains 13 protein-coding genes, 2 r RNA genes, 22 t RNA genes, a light-strand replication origin and a control region. Additionally, phylogenetic analysis based on the nucleotide sequences of 13 conserved protein-coding genes using the maximum likelihood method indicated that the mitochondrial genome is maternally inherited. This study presents genomic data for studying phylogenetic relationships and breeding of hybrid Epinephelinae.

Comparative analysis of mitochondrial genome of a deepsea crab Chaceon granulates reveals positive selection and novel genetic features

Deep-sea organisms survive in an extremely harsh environment.There must be some genetic adaptation mechanisms for them.We systematically characterized and compared the complete mitochondrial genome(mitogenome)of a deep-sea crab(Chaceon granulates)with those of shallow crabs.The mitogenome of the crab was 16126 bp in length,and encoded 37 genes as most of a metazoan mitogenome,including 13 protein-coding genes(PCGs),22 transfer RNA(tRNA)genes,and 2 ribosomal RNA(rRNA)genes.The gene arrangement and orientation was conserved in the crabs.However,a unique mitogenome element regulator,the origin of light-strand replication(O L),was firstly predicted in the present crab mitogenome.In addition,further positive selection analysis showed that two residues(33 S in ND3 and 502 I in ND5)in C.granulates mitogenome were positively selected,indicated the selective evolution of the deep-sea crab.Therefore,the mitogenome of deep-sea C.granulates showed a unique OL element and positive selection.These special features would influence the mitochondrial energy metabolism,and be involved in the adaptation of deepsea environment,such as oxygen deficits and low temperatures.

The mitochondrial genome of Chaeturichthys stigmatias provides novel insight into the interspecific difference with Amblychaeturichthys hexanema

Chaeturichthys stigmatias and Amblychaeturichthys hexanema belong to the family Gobiidae,which are offshore warm fish species and widely distribute in the western Pacific Ocean.In this study,the mitochondrial cytochrome c oxidase subunit I(COI)sequences and 12 S ribosomal RNA(12 S rRNA)sequences were used to analyze the interspecific differences between the two species.The phylogenetic analysis showed that the interspecific distance was significantly higher than the intraspecific genetic distance.The Neighbor-Joining tree showed two separate clusters,without sharing haplotype.The mitochondrial genome sequence of C.stigmatias was also reported.This genome was 17134 bp in size,with a high A+T content of 55.9%.The phylogenetic analysis based on the tandem 13 coding protein genes nucleotide sequences indicated that C.stigmatias showed a close relationship with A.hexanema.This study can provide the basic genetic data for two species and will help for constructing the phylogeny of the Gobiiade.

Assembly and marker analysis of mitochondrial genomes provide insights into origin,evolution and spread of Brassica juncea(L.)Czern.et Coss.

The release of mitochondrial genome sequences provides the basis for characterizing interspecific and intraspecific variation in Brassica mitochondrial genomes.However,few B.juncea(mustard)mitochondrial genomes have been published.We assembled the mitochondrial genomes of three B.juncea subspecies and compared them with previously published genomes.The genomes were phylogenetically classified into A,B,C,and Bna clades.Two variant sites,a transversion(C→A)at nt 79,573 and a 31-bp copy-number variation between nts 65,564 and 65,596,were identified.Based on these variant sites,mitotype-specific sequence markers were developed to characterize the variation among worldwide 558 B.juncea accessions.Three mitochondrial genome types(mitotypes MT1–MT3)were identified.In terms of geographical distribution,MT1 and MT2 accessions were distributed mainly to the north and MT3 to the south of 34°N.Root mustards carried only MT1,leaf and stem mustards carried mainly MT3,and seed mustards carried all three mitotypes,implying that the mitotypes underwent selection during B.juncea domestication.A new form of oil mustard evolved by hybridization between two gene pools in southwest China.

The cotton mitochondrial chimeric gene orf610a causes male sterility by disturbing the dynamic balance of ATP synthesis and ROS burst

Plant cytoplasmic male sterility(CMS)is maternally inherited and often manifested as aborted pollen development,but the molecular basis of abortion remains to be identified.To facilitate an investigation of CMS in cotton,the complete sequence of cotton mitochondrial(mt)genome for CMS-D2 line ZBA was determined.The mt genome was assembled as a single circular molecule with 634,036 bp in length.A total of 194 ORFs,36 protein-coding genes,six r RNAs,and 24 t RNAs were identified.Several chimeric genes encoding hypothetical proteins with transmembrane domains were identified.Among them,a previously unknown chimeric gene,orf610a,which is composed of atp1 and a 485-bp downstream sequence of unknown nature,was identified.RT-PCR and q RT-PCR validation indicated that orf610a was expressed specifically in a sterile line.Ectopic expression of orf610a in yeast resulted in excessive accumulation of reactive oxygen species and reduction in ATP content,in addition to inhibition of cellular growth.Transgenic A.thaliana overexpressing orf610a fused with a mitochondrial targeting peptide displayed partial male sterility.Interaction between ORF610a and the nuclear-encoded protein RD22 indicated an association between ORF610a and pollen abortion.Positive feedback during transcriptional regulation between nuclear regulatory factors and the mt CMS gene may account for the male sterility of ZBA.

Complete mitochondrial genome of the Thai Red Junglefowl (Gallus gallus) and phylogenetic analysis

DEAR EDITOR,In this study, we sequenced the complete mitochondrial genome （mitogenome） of the Thai Red Junglefowl （RJF; Gallus gallus） using the next-generation sequencing （NGS） platform of the Ion Torrent PGM. Samples were taken from Mae Wang District, Chiang Mai Province, northern Thailand Our data showed the complete mitogenome to be 16 785 bp in length, composed by 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region. The genome nucleotide composition was 30.3% A, 23.7% T, 32.5% C, and 13.5% G, resulting in a high percentage of A＋T （50.4%）. Phylogenetic analysis revealed that the mitogenome belonged to haplogroup X, whereas those of all domestic chickens belong to haplogroups A to G. This newly released mitogenome sequence will advance further evolutionary and population genetics study of the RJF and domestic chicken The availability of the G. gallus mitogenome will also contribute to further conservation genetics research of a unique species, listed as ＇data deficient＇ in Thailand.

Three New Ranidae Mitogenomes and the Evolution of Mitochondrial Gene Rearrangements among Ranidae Species

Various types of gene rearrangements have been discovered in the mitogenoes of the frog family Ranidae. In this study, we determined the complete mitogenome sequence of three Rana frogs. By combining the available mitogenomic data sets from GenBank, we evaluated the phylogenetic relationships of Ranidae at the mitogenome level and analyzed mitogenome rearrangement cases within Ranidae. The three frogs shared an identical mitogenome organization that was extremely similar to the typical Neobatrachian-type arrangement. Except for the genus Babina, the monophyly of each genus was well supported. The genus Amnirana occupied the most basal position among the Ranidae. The [Lithobates ＋ Rana] was the closest sister group of Odorrana. The diversity of mitochondrial gene arrangements in ranid species was unexpectedly high, with 47 mitogenomes from 40 ranids being classified into 10 different gene rearrangement types. Some taxa owned their unique gene rearrangement characteristics, which had significant implication for their phylogeny analysis. All rearrangement events discovered in the Ranidae mitogenomes can be explained by the duplication and random loss model.

A Monophyletic Status of Axis Genus in Subfamily Cervinae Supported by the Complete Mitochondrial Genome of Chinese Hog Deer(Axis porcinus)

Hog deer(Axis porcinus)is a small mammal and listed in the International Union for Conservation of Nature.However,phylogenetic position of hog deer within Axis genus has remained controversial.In the present study,we first assembled complete mitochondrial genome of Chinese hog deer reared in Chengdu Zoo,Sichuan,by the second-generation sequencing technology.This newly assembled mitochondrial genome of hog deer is 16376 bp in length and consists of 13 protein-encoding genes,23 transfer RNA genes and 2 ribosomal RNA genes.Phylogenetic analyses based on complete mitochondrial genome and cytochrome b gene sequences revealed that hog deer is closely clustered together and placed with sister taxon of spotted deer(A.Axis),which therefore supported monophyletic statue of Axis genus.Furthermore,considerable genetic differentiation,up to 139 mutations of complete mitochondrial genome was revealed between geographical populations of hog deer in France and Southeast Asia.However,only six variable sites(nucleotide diversity of 0.00007)and four haplotypes(haplotype diversity of 0.533)were totally detected among ten newly sequenced Chinese hog deer.The results provide a better understanding on the phylogeny of hog deer.