Interisland Mutation of a Novel Phospholipase A<Subscript>2</Subscript> from <Emphasis Type="Italic">Trimeresurus flavoviridis</Emphasis> Venom and Evolution of Crotalinae Group II Phospholipases A<Subscript>2</Subscript>

Trimeresurus flavoviridis (Crotalinae) snakes inhabit the southwestern islands of Japan: Amami-Oshima, Tokunoshima, and Okinawa. Affinity and conventional chromatographies of Amami-Oshima T. flavoviridis venom led to isolation of a novel phospholipase A₂ (PLA₂). This protein was highly homologous (91%) in sequence to trimucrotoxin, a neurotoxic PLA₂, which had been isolated from T. mucrosquamatus (Taiwan) venom, and exhibited weak neurotoxicity. This protein was named PLA-N. Its LD₅₀ for mice was 1.34 µg/g, which is comparable to that of trimucrotoxin. The cDNA encoding PLA-N was isolated from both the Amami-Oshima and the Tokunoshima T. flavoviridis venom-gland cDNA libraries. Screening of the Okinawa T. flavoviridis venom-gland cDNA library with PLA-N cDNA led to isolation of the cDNA encoding one amino acid-substituted PLA-N homologue, named PLA-N(O), suggesting that interisland mutation occurred and that Okinawa island was separated from a former island prior to dissociation of Amami-Oshima and Tokunoshima islands. Construction of a phylogenetic tree of Crotalinae venom group II PLA₂s based on the amino acid sequences revealed that neurotoxic PLA₂s including PLA-N and PLA-N(O) form an independent cluster which is distant from other PLA₂ groups such as PLA₂ type, basic [Asp⁴⁹]PLA₂ type, and [Lys⁴⁹]PLA₂ type. Comparison of the nucleotide sequence of PLA-N cDNA with those of the cDNAs encoding other T. flavoviridis venom PLA₂s showed that they have evolved in an accelerated manner. However, when comparison was made within the cDNAs encoding Crotalinae venom neurotoxic PLA₂s, their evolutionary rates appear to be reduced to a level between accelerated evolution and neutral evolution. It is likely that ancestral genes of neurotoxic PLA₂s evolved in an accelerated manner until they had acquired neurotoxic function and since then they have evolved with less frequent mutation, possibly for functional conservation. The nucleotide sequences reported in this paper are available from the GenBank/EMBL/DDBJ databases under accession numbers AB102728 and AB102729.     

Accelerated evolution of Trimeresurus okinavensis venom gland phospholipase A2 isozyme-encoding genes

Three Trimeresurus okinavensis (To; himehabu snake, Crotalinae) venom gland phospholipase A₂ (PLA₂) isozymeencoding genes, gPLA₂-o1, gPLA₂-o2 and gPLA₂-o3, were isolated from its genomic DNA library. The nucleotide (nt) sequence analysis revealed that two of the three genes (gPLA₂-o2) and (gPLA₂-o3) occasionally have been converted to inactivated genes by introduction of one base insertion or substitution. It was confirmed from Southern blot analysis that the To haploid genome contains only three venom gland PLA₂ isozyme genes herein isolated. Comparison of these genes showed that nonsynonymous nt substitutions have occurred more frequently than synonymous nt substitutions in the protein-coding regions, except for the signal-peptide coding domain, implying that To venom gland PLA₂isozyme genes have evolved via accelerated evolution. Such an evolutionary feature of To venom gland PLA₂ isozyme genes proves the general universality of accelerated evolution previously drawn for venom gland PLA₂ isozyme genes of other crotalinae snakes. The variability in the mature protein-coding regions of three To venom gland PLA₂ isozyme genes appears to have been brought about by natural selection for point mutations.     

The molecular cloning of a phospholipase A2 from Bothrops jararacussu snake venom: Evolution of venom group II phospholipase A2's may imply gene duplications

The sequence coding for a snake venom phospholipase A₂ (PLA₂), BJUPLA₂, has been cloned from a Bothrops jararacussu venom gland cDNA library. The cDNA sequence predicts a precursor containing a 16-residue signal peptide followed by a molecule of 122 amino acid residues with a strong sequence similarity to group II snake venom PLA₂'s. A striking feature of the cDNA is the high sequence conservation of the 5 and 3 untranslated regions in cDNAs coding for PLA₂'s from a number of viper species. The greatest sequence variation was observed between the regions coding for the mature proteins, with most substitutions occurring in nonsynonymous sites. The phylogenetic tree constructed by alignment of the amino acid sequence of BJUPLA₂ with group II PLA₂'s in general groups them according to current taxonomical divisions and/or functional activity. It also suggests that gene duplications may have occurred at a number of different points during the evolution of snake venom group II PLA₂'s.The nucleotide sequence reported in this paper has been submitted to the GenBank/EMBL Data Bank with accession number X76289.Correspondence to: A.M. Moura-da-Silva     

Structures of genes encoding phospholipase A2 inhibitors from the serum of Trimeresurus flavoviridis snake

Inhibitors (PLIs) against snake venom gland phospholipases A₂ (PLA₂s) have been found in their sera. A cDNA encoding a PLI from Trimeresurus flavoviridis (Tf, habu snake, Crotalinae) serum, cPLI-A, was isolated from the Tf liver cDNA library and sequenced. Northern blot analysis with cPLI-A showed that PLIs are expressed only in liver. Genes for PLIs, gPLI-A and gPLI-B, were isolated from the Tf genomic DNA library and their nucleotide (nt) sequences were determined. The genes consisted of four exons and three introns, and exon 4 encoded the carbohydrate recognition domain (CRD)-like motif. Comparison of the nt sequences between gPLI-A and gPLI-B showed that these genes are highly homologous, including introns, except that exon 3 is rich in nonsynonymous nt substitutions which are almost four times as frequent as synonymous nt substitutions. This evolutionary feature of PLI genes is different from that of venom gland PLA₂ isozyme genes in which nonsynonymous nt substitutions are spread over the entire mature protein-coding region.     

Plant low-molecular-weight phospholipase A2s (PLA2s) are structurally related to the animal secretory PLA2s and are present as a family of isoforms in rice (Oryza sativa)

Recently, we purified to homogeneity and characterized a low-molecular-weight calcium-dependent phospholipase A₂ (PLA₂) from developing elm seed endosperm. This represented the first purified and characterized PLA₂ from a plant tissue. The full sequences of two distinct but homologous rice (Oryza sativa) cDNAs are given here. These encode mature proteins of 119 amino acids (PLA₂-I, preceded by a 19 amino acid signal peptide) and 128 amino acids (PLA₂-II, preceded by a 25 amino acid signal peptide), and were derived from four expressed sequence tag (EST) clones. Both proteins were homologous to the N-terminal amino acid sequence of the elm PLA₂. They contained twelve conserved cysteine residues and sequences that are likely to represent the Ca²⁺-binding loop and active-site motif, which are characteristic of animal secretory PLA₂s. A soluble PLA₂ activity was purified 145 000-fold from green rice shoots. This had the same biochemical characteristics as the elm and animal secretory PLA₂s. The purified rice PLA₂ consisted of two proteins, with a molecular weight of 12 440 and 12 920, that had identical N-terminal amino acid sequences. This sequence was different from but homologous to the PLA₂-I and PLA₂-II sequences. Taken together, the results suggest that at least three different low-molecular-weight PLA₂s are expressed in green rice shoots. Southern blot analysis suggested that multiple copies of such genes are likely to occur in the rice and in other plant genomes.     

cDNA cloning and sequencing of phospholipase A2 from the pyloric ceca of the starfish Asterina pectinifera

Three cDNA from the pyloric ceca of the starfish Asterina pectinifera, (namely, cDNA 1, 2, and 3), encoding phospholipase A₂ (PLA₂), were isolated and sequenced. These cDNAs were composed of 415 bp with an open reading frame of 414 bp at nucleotide positions 1–414, which encodes 138 amino acids including N-terminal Met derived from the PCR primer. The amino acid sequence deduced from the cDNA 1 was completely consistent with the sequence determined with the starfish PLA₂ protein, while those deduced from cDNA 2 and cDNA 3 differed at one and twelve amino acid residual positions, respectively, from the sequence of the PLA₂ protein, suggesting the presence of multiple forms in the starfish PLA₂. All of the sequences deduced from cDNA 1, 2, and 3 required two amino acid deletions in pancreatic loop region, and sixteen insertions and three deletions in β-wing region when aligned with the sequence of mammalian pancreatic PLA₂. In phylogenetic tree, the starfish PLA₂ should be classified into an independent group, but hardly to the established groups IA and IB. The characteristic structure in the pancreatic loop and β-wing regions may account for the specific properties of the starfish PLA₂, e.g. the higher activity and characteristic substrate specificity compared with commercially available PLA₂ from porcine pancreas.     

Vaporization of Inorganic Mercury by Cell-free Extracts of Drug Resistant Escherichia coli

The cDNAs encoding venom phospholipase A₂ (PLA₂) inhibitors (PLIs), named Protobothrops elegans (Pe)γPLI-A, PeγPLI-B, PeαPLI-A, and PeαPLI-B, were cloned from the P. elegans liver cDNA library. They were further divided into several constituents due to nucleotide substitutions in their open reading frames. For PeαPLI-A, two constituents, PeαPLI-A_a and PeαPLI-A_b, were identified due to three nonsynonymous substitutions in exon 3. Far-western blot and mass-spectrometry analysis of the P. elegans serum proteins showed the presence of γPLIs, and αPLIs, which can bind venom PLA₂s. In αPLIs from Protobothrops sera, A or B subtype-specific amino acid substitutions are concentrated only in exon 3. A comparison of γPLIs showed that γPLI-As are conserved and γPLI-Bs diversified. Mathematical analysis of the nucleotide sequences of Protobothrops γPLI-B cDNAs revealed that the particular loops in the three-finger motifs diversified by accelerated evolution. Such evolutionary features should have made serum PLIs acquire their respective inhibitory activities to adapt to venom PLA₂ isozymes.     

Isolation, characterization and expression of cDNA clones encoding a mitochondrial malate translocator from Panicum miliaceum L.

Three cDNA clones that hybridize to a partial rice cDNA that show similarity to bovine mitochondrial 2-oxoglutarate/malate translocator were isolated from leaves of Panicum miliaceum L. (proso millet), an NAD-malic enzyme-type C₄ plant. The nucleotide sequences of the clones resemble each other, and some of the isolated cDNAs contained extra sequences that seemed to be introns. The predicted proteins encoded by the cDNAs have 302 amino acids and molecular weights of 32211 and 32150. The hydrophobic profile of the amino acid sequence predicted the existence of six transmembrane -helices that is a common property of members in the mitochondrial transporter family. The predicted amino acid sequence showed the highest similarity with that of the 2-oxoglutarate/malate translocator from mammalian mitochondria. An expression plasmid containing the coding region of the cDNAs was used to over-express recombinant protein with a C-terminal histidine tag Escherichia coli, which was affinity-purified. The antibody against the recombinant protein cross-reacted with proteins of 31–32 kDa in the membrane fraction from P. miliaceum mitochondria, but not with the chloroplast fraction. The recombinant protein reconstituted in liposomes efficiently transported malate, citrate, and 2-oxoglutarate.     

Rapid Evolution by Positive Selection and Gene Gain and Loss: PLA<Subscript>2</Subscript> Venom Genes in Closely Related <Emphasis Type="Italic">Sistrurus</Emphasis> Rattlesnakes with Divergent Diets

Rapid evolution of snake venom genes by positive selection has been reported previously but key features of this process such as the targets of selection, rates of gene turnover, and functional diversity of toxins generated remain unclear. This is especially true for closely related species with divergent diets. We describe the evolution of PLA₂ gene sequences isolated from genomic DNA from four taxa of Sistrurus rattlesnakes which feed on different prey. We identified four to seven distinct PLA₂ sequences in each taxon and phylogenetic analyses suggest that these sequences represent a rapidly evolving gene family consisting of both paralogous and homologous loci with high rates of gene gain and loss. Strong positive selection was implicated as a driving force in the evolution of these protein coding sequences. Exons coding for amino acids that make up mature proteins have levels of variation two to three times greater than those of the surrounding noncoding intronic sequences. Maximum likelihood models of coding sequence evolution reveal that a high proportion (∼30%) of all codons in the mature protein fall into a class of codons with an estimated d _N /d _S (ω) ratio of at least 2.8. An analysis of selection on individual codons identified nine residues as being under strong (p < 0.01) positive selection, with a disproportionately high proportion of these residues found in two functional regions of the PLA₂ protein (surface residues and putative anticoagulant region). This is direct evidence that diversifying selection has led to high levels of functional diversity due to structural differences in proteins among these snakes. Overall, our results demonstrate that both gene gain and loss and protein sequence evolution via positive selection are important evolutionary forces driving adaptive divergence in venom proteins in closely related species of venomous snakes.     

Interisland Evolution of <Emphasis Type="Italic">Trimeresurus flavoviridis</Emphasis> Venom Phospholipase A<Subscript>2</Subscript> Isozymes

Abstract Trimeresurus flavoviridis snakes inhabit the southwestern islands of Japan. A phospholipase A₂ (PLA₂), named PL-Y, was isolated from Okinawa T. flavoviridis venom and its amino acid sequence was determined from both protein and cDNA. PL-Y was unable to induce edema. In contrast, PLA-B, a PLA₂ from Tokunoshima T. flavoviridis venom, which is different at only three positions from PL-Y, is known to induce edema. A new PLA₂, named PLA-B′, which is similar to PLA-B, was cloned from Amami-Oshima T. flavoviridis venom gland. Three T. flavoviridis venom basic [Asp⁴⁹]PLA₂ isozymes, PL-Y (Okinawa), PLA-B (Tokunoshima), and PLA-B′ (Amami-Oshima), are identical in the N-terminal half but have one to four amino acid substitutions in the β1-sheet and its vicinity. Such interisland sequence diversities among them are due to isolation in the different environments over 1 to 2 million years and appear to have been brought about by natural selection for point mutation in their genes. Otherwise, a major PLA₂, named PLA2, ubiquitously exists in the venoms of T. flavoviridis snakes from the three islands with one to three synonymous substitutions in their cDNAs. It is assumed that the PLA2 gene is a prototype among T. flavoviridis venom PLA₂ isozyme genes and has hardly undergone nonsynonymous mutation as a principal toxic component. Phylogenetic analysis based on the amino acid sequences revealed that T. flavoviridis PLA₂ isozymes are clearly separated into three groups, PLA2 type, basic [Asp⁴⁹]PLA₂ type, and [Lys⁴⁹]PLA₂ type. Basic [Asp⁴⁹]PLA₂-type isozymes may manifest their own particular toxic functions different from those of the isozymes of the PLA2 type and [Lys⁴⁹]PLA₂ type.     

Accelerated evolution in the protein-coding region of galectin cDNAs, congerin I and congerin II, from skin mucus of conger eel (Conger myriaster).

Two cDNAs encoding galectins named congerins I and II from the skin mucus of conger eel (Conger myriaster) were isolated and sequenced. Comparison of the nucleotide sequences of congerins I and II showed that the sequence similarities of the 5' and 3' untranslated regions (86 and 88%, respectively) were much higher than those of the protein-coding region (73%). The numbers of nucleotide substitutions per site (KN) for the untranslated regions are smaller than the numbers of nucleotide substitutions per synonymous site (KS) for the protein coding region. Furthermore, nonsynonymous nucleotide substitutions have accelerated more frequently than synonymous nucleotide substitutions in the protein coding region (KA/KS = 2.57). These results suggest that accelerated substitutions have occurred in the protein-coding regions of galectin genes to generate diverse galectins with different molecular properties. Northern blot analysis showed that both congerins were expressed not only in the skin tissues but also in the stomach of conger eel.     

Characterization of a stress-induced,developmentally regulated gene family from soybean

We describe a family of stress-induced, developmentally regulated soybean genes for which cDNAs have been obtained from two different cultivars (Glycine max cv. Mandarin and Glycine max cv. Williams). The mRNAs corresponding to these cDNAs, called SAM22 and H4, respectively, accumulate predominantly in the roots of soybean seedlings but are present at high levels in the roots and leaves of mature plants. SAM22 accumulation is especially dramatic in senescent leaves. In addition, SAM22 accumulation can be induced in young leaves by wounding or by transpiration-mediated uptake of salicylic acid, methyl viologen, fungal elicitor, hydrogen peroxide or sodium phosphate (pH 6.9). Taken together, these data indicate that the genes corresponding to SAM22 and H4 are induced by various stresses and developmental cues. Southern blot analysis indicates that multiple copies of sequences related to SAM22 exist in the soybean genome. We also show that the nucleotide sequences of the cDNAs corresponding to SAM22 and H4 are 86% identical at the nucleotide level to each other and 70% identical at the amino acid level to the disease resistance response proteins of Pisum sativum.     

Identification of novel phospholipase A2 group IX members in metazoans

Sequence homologues of the bacterium Streptomyces violaceoruber and sea anemone Nematostella vectensis PLA₂ pfam09056 members were identified in several bacteria, fungi and metazoans illustrating the evolution of this PLA₂ sub-family. Comparison of their molecular structures revealed that bacteria and fungi members are part of the GXIV of PLA₂s while metazoan representatives are similar with GIX PLA₂ of the marine snail Conus magus. Members of GXIV and GIX PLA₂s show modest overall sequence similarity (21–35%) but considerable motif conservation within the putative Ca²⁺-binding, catalytic sites and cysteine residue positions which are essential for enzyme function. GXIV PLA₂s of bacteria and fungi typically contain four cysteine residues composing two intramolecular disulphide bonds. GIX PLA₂ homologues were identified in cnidarians and molluscs and in a single tunicate but appear to be absent from other metazoan genomes. The mature GIX PLA₂ deduced peptides contain up to ten cysteine residues capable of forming five putative disulphide bonds. Three disulphide bonds were identified in GIX PLA₂s, two of which correspond to those localized in GXIV PLA₂s. Phylogenetic analysis demonstrates that metazoan GIX PLA₂s cluster separate from the bacterial and fungal GXIV PLA₂s and both pfam09056 members form a group separate from the prokaryote and eukaryote GXIIA PLA₂ pfam06951. Duplicate PLA₂ pfam09056 genes were identified in the genomes of sea anemone N. vectensis and oyster Crassostrea gigas suggest that members of this family evolved via species-specific duplication events. These observations indicate that the newly identified metazoan pfam09056 members may be classified as GIX PLA₂s and support the idea of the common evolutionary origin of GXIV and GIX PLA₂ pfam09056 members, which emerged early in bacteria and were maintained in the genomes of fungi and selected extant metazoan taxa.     

Purification and characterization of a potent hemolytic toxin with phospholipase A2 activity from the venom of Indian Russell's viper

A potent toxin with phospholipase A₂ (PLA₂) and hemolytic activity in vitro was purified from the Russell's viper venom of eastern India (RVV-EI). The purified protein (RVV-PFIIc) of 15.3 kDa molecular weight, and a lethal toxicity dose (LD₅₀ i.p.) of 0.1 mg/kg body weight, was the most toxic PLA₂ so far reported from the Indian subcontinent. The material also possessed anticoagulant activity as it enhanced the prothrombin induced plasma clotting time in vitro. The PLA₂ toxin (RVV-PFIIc) was shown to be different from other PLA₂s of RVV in respect to one or more of these parameters e.g. molecular weight, isoelectric pH, in vivo toxicity, specific activity of the enzyme and certain other biological activities. The first 19 amino terminal sequence (NLFQFAEMIVKMTGKEAVH) of RVV-PFIIc showed variable degree of homology (42.1–94.7%) with those of other RVV-PLA₂s described in the literature. Antisera raised against RVV-EI or RVV-PFIIc, though completely neutralized the in vivo lethal toxicity of RVV-EI or RVV-PFIIc, failed to inhibit their PLA₂ activity in vitro thereby suggesting that in vivo toxicity and in vitro activity of the enzyme may not be directly related. Apart from RVV-PFIIc, at least two other PLA₂ isozymes were found to be present in RVV-EI that were distinct from RVV-PFIIc in respect to their molecular, biological as well as serological properties. The significance of these and related data in antivenom therapy is discussed.     

Three extra copies of a C4-related gene in H-2 w7 mice are C4/Slp hybrid genes generated by multiple recombinational events

Mice bearing the H-2 _w7 haplotype have five C4-related genes and constitutively express the Slp antigen. To understand the structure and evolution of the five C4-related genes of the C3H.W7 mouse, we have determined nucleotide sequences of the 5 end region of these genes. A C4/Slp hybrid nature was confirmed for three of five C4-related genes as predicted previously by restriction enzyme analysis. The nucleotide sequences of the 5 flanking regions of these three hybrid genes showed close similarity to that of the C4 gene, while the 3 side of the ninth exon of the three hybrid genes showed close similarity to that of the Slp gene. In contrast, the regions between the first exon and the middle of the ninth exon of the three hybrid genes showed a mosaic structure of C4-like and Slp-like sequences. Moreover, the boundaries of the C4-like and Slp-like sequences were quite different among the three hybrid genes. The pattern of nucleotide sequence diversity in this region among the five C4-related sequences could be mainly explained not by point mutations but by gene conversions or unequal crossovers. These results suggest that multiple genetic recombinational events between two homologous sequences played an important role in the generation and diversification of the extra copies of the C4/Slp gene in the H-2 ^w7 mouse.The nucleotide sequence data reported in this paper have been submitted to the DDBJ, EMBL, and GenBank nucleotide sequence databases and have been assigned the accession numbers D90167-71.     

Genomic structure and chromosomal location of the mouse pre-T-cell receptor alpha gene

The mouse pre-T-cell receptor alpha (pT) chain is a 33 000 M _r glycoprotein expressed on the surface of immature thymocytes as a disulfide-linked heterodimer with the T-cell receptor beta (TCR) chain, and in association with proteins of the CD3 complex. The cDNA for pT, isolated previously, encodes a type I transmembrane protein that is a member of the immunoglobulin (Ig) superfamily. Here we report the complete nucleotide sequence, the exon/intron structure, and the chromosomal location of the pTa gene. The gene spans about 8.4 kilobases (kb) and consists of four exons. Exon 1 encodes the 5 untranslated region, the leader peptide, and the first three amino acids of the mature protein. This exon is followed by a relatively long intron of 4.9 kb that contains many short interspersed repeats (SINEs) of the B1 and B2 family. The second exon encodes the extracellular Ig-like domain and exon 3 with just 45 base pairs the connecting peptide (CP), including the cysteine required for heterodimer formation. A similar exon/intron structure encoding corresponding parts of the mature polypeptide is found both in the Tcra and Tcrd constant region genes. The last exon encodes the transmembrane portion, the cytoplasmic tail, and about 540 nucleotides of 3 untranslated sequence, including a B2 repetitive element. In situ hybridization maps the pTa gene to the D/E1 region of mouse chromosome 17.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number U27268     

Lysosomotropic agents activate the capacity for calcium dependent pinocytosis in starvedAmoeba proteus: evidence for a mechanism involving phospholipase A2

Summary Pinocytosis induced by Na⁺ was assayed by phase contrast microscopy in 8–12 days starvedAmoeba proteus. These cultures were inactive with respect to calcium-dependent Na⁺-induced pinocytosis, but treatment with amino acid methyl and ethyl esters increased their capacity for pinocytosis. Besides promoting pinocytosis these compounds also stimulated calcium-sensitive secretion of lysosomal enzymes from normal, 2–3 days starved, cells. Only uncharged 1-forms of the amino acid esters were effective. Also other lysosomotropic compounds including monodansylcadaverine, glycine-phenylalanine-2-naphthylamide, NH₄Cl, and the ionophores monensin and A23187 activated starved cells. The effect of these agents (except A23187) was inhibited by the drug dantrolene suggesting that activation is a consequence of release of Ca²⁺ from intracellular stores. Several of the lysosomotropic agents also lost their activating effect in the presence of phospholipase A₂ (PLA₂) inhibitors. To investigate whether or not PLA₂ activity in the cell culture could imitate the effect of the lysosomotropic agents, we incubated starved cells with snake venom PLA₂s. These enzymes caused rapid, dantrolene-sensitive activation of the cells. Measurement of endogenous PLA₂in normal cells revealed significant cellular activity but no significant secretion of the enzyme into the culture medium was observed. Together the studies with enzyme inhibitors and dantrolene suggest that the process by which lysosomotropic agents affect pinocytosis involves activation of PLA₂ and release of Ca²⁺ from intracellular stores.Abbreviations AnBOMe amino-n-butyric acid methylester - Et ethylester - GPN glycine-1-phenylalanine-2-naphthylamide - MDC monodansylcadaverine - MDTC monodansylthiacadaverine - Me methylester - pBPB p-bromo phenacylbromide - PLA₂ phospholipase A₂     

The complete mitochondrial genome of a threatened loach (<Emphasis Type="Italic">Sinibotia reevesae</Emphasis>) and its phylogeny

In present study, the complete mitochondrial genome of Sinibotia reevesae was first sequenced using the next-generation sequencing technology and annotated using bioinformatic tools. The circular mitochondrial genome was 16,572 bp in length, and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and 1 displacement loop locus. It presents a typical gene organization and order for completely sequenced cypriniformes mitogenomes. The control region could be divided into three parts included the extended termination associated sequence domain, the central conserved domain and the conserved sequence block. Interestingly, two stem-loop domains were found in control region and O_L region, respectively. Furthermore, phylogenetic analyses using concatenated amino acid and nucleotide sequences of the 13 protein-coding genes with two different methods (Maximum likelihood and Bayesian analysis) both highly supported the close relationship of S. reevesae and Sinibotia superciliaris, which was in line with the previous classifications based on morphological and molecular studies. These data provide useful information for a better understanding of the mitogenomic diversities and evolution in fish as well as novel genetic markers for studying population genetics and species identification.