Protein-coding Gene Sequences Research Articles

According to mitochondrial DNA evidence, Parascaris equorum and Parascaris univalens may represent the same species.

Parascarosis is caused mainly by parasitic infections with Parascaris equorum and Parascaris univalens, the most common ascarid nematodes, in the small intestine of equines. Parascarosis often causes severe illness and even death in foals and yearlings. In this study, we obtained the complete sequence of the P. equorum mitochondrial (mt) genome and compared its organization and structure with that of P. equorum Japan isolate (nearly complete), and the complete mtDNA sequences of P. univalens Switzerland and USA isolates. The complete mtDNA genome of P. equorum China isolate is 13,899 base pairs (bp), making it the smallest of the four genomes. All four Parascaris mt genomes are circular, and all genes are transcribed in the same direction. The P. equorum mtDNA genome consists of 12 protein-coding genes, two ribosomal RNA genes, 22 transfer (t) RNA genes and one non-coding region, which is consistent with P. equorum Japan isolate and P. univalens Switzerland isolate but distinct from P. univalens USA isolate, which has 20 tRNA genes. Differences in nucleotide sequences of the four entire mt genomes range from 0.1-0.9%, and differences in total amino acid sequences of protein-coding genes are 0.2-2.1%. Phylogenetic analyses showed that the four Parascaris species clustered in a clade, indicating that P. equorum and P. univalens are very closely related. These mt genome datasets provide genetic evidence that P. equorum and P. univalens may represent the same species, which will be of use in further studies of the taxonomy, systematics and population genetics of ascarids and other nematodes.

Assessing the functional association of intronic miRNAs with their host genes.

In human, nearly half of the known microRNAs (miRNAs) are encoded within the introns of protein-coding genes. The embedment of these miRNA genes within the sequences of protein-coding genes alludes to a possible functional relationship between intronic miRNAs and their hosting genes. Several studies, using predicted targets, suggested that intronic miRNAs influence their hosts’ function either antagonistically or synergistically. New experimental data of miRNA expression patterns and targets enable exploring this putative association by relying on actual data rather than on predictions. Here, our analysis based on currently available experimental data implies that the potential functional association between intronic miRNAs and their hosting genes is limited. For host-miRNA examples where functional associations were detected, it was manifested by either autoregulation, common targets of the miRNA and hosting gene, or through the targeting of transcripts participating in pathways in which the host gene is involved. This low prevalence of functional association is consistent with our observation that many intronic miRNAs have independent transcription start sites and are not coexpressed with the hosting gene. Yet, the intronic miRNAs that do show functional association with their hosts were found to be more evolutionarily conserved compared to other intronic miRNAs. This might suggest a selective pressure to maintain this architecture when it has a functional consequence.

Extensive gene rearrangements in the mitochondrial genomes of two egg parasitoids, Trichogramma japonicum and Trichogramma ostriniae (Hymenoptera: Chalcidoidea: Trichogrammatidae)

Animal mitochondrial genomes usually exhibit conserved gene arrangement across major lineages, while those in the Hymenoptera are known to possess frequent rearrangements, as are those of several other orders of insects. Here, we sequenced two complete mitochondrial genomes of Trichogramma japonicum and Trichogramma ostriniae (Hymenoptera: Chalcidoidea: Trichogrammatidae). In total, 37 mitochondrial genes were identified in both species. The same gene arrangement pattern was found in the two species, with extensive gene rearrangement compared with the ancestral insect mitochondrial genome. Most tRNA genes and all protein-coding genes were encoded on the minority strand. In total, 15 tRNA genes and seven protein-coding genes were rearranged. The rearrangements of cox1 and nad2 as well as most tRNA genes were novel. Phylogenetic analysis based on nucleotide sequences of protein-coding genes and on gene arrangement patterns produced identical topologies that support the relationship of (Agaonidae + Pteromalidae) + Trichogrammatidae in Chalcidoidea. CREx analysis revealed eight rearrangement operations occurred from presumed ancestral gene order of Chalcidoidea to form the derived gene order of Trichogramma. Our study shows that gene rearrangement information in Chalcidoidea can potentially contribute to the phylogeny of Chalcidoidea when more mitochondrial genome sequences are available.

Complete mitochondrial genome of the Northern Long-eared Owl (Asio otus Linnaeus, 1758) determined using next-generation sequencing

In this study, the mitogenome of Asio otus, the Northern Long-eared Owl, was analysed using Illumina next-generation sequencing. The mitogenome was found to be a circular molecule, 17,735 bp long with a slight AT bias (53.0%). The gene arrangement pattern was the same as that of a typical vertebrate, containing 37 genes (13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a non-coding control region). In the putative control region (1984 bp), there were two types of tandem repeats at the end of the region. A phylogenetic tree was constructed using the 13 PCG sequences discovered in this study and those of that have been previously published of other Strigidae species and revealed a close relationship between A. otus and A. flammeus. The newly generated mitogenome from this study enriches the genomic resources available for future evolutionary studies and promotes conservation genetics of this species.

Organellar phylogenomics inform systematics in the green algal family Hydrodictyaceae (Chlorophyceae) and provide clues to the complex evolutionary history of plastid genomes in the green algal tree of life.

Phylogenomic analyses across the green algae are resolving relationships at the class, order, and family levels and highlighting dynamic patterns of evolution in organellar genomes. Here we present a within-family phylogenomic study to resolve genera and species relationships in the family Hydrodictyaceae (Chlorophyceae), for which poor resolution in previous phylogenetic studies, along with divergent morphological traits, have precluded taxonomic revisions. Complete plastome sequences and mitochondrial protein-coding gene sequences were acquired from representatives of the Hydrodictyaceae using next-generation sequencing methods. Plastomes were characterized, and gene order and content were compared with plastomes spanning the Sphaeropleales. Single-gene and concatenated-gene phylogenetic analyses of plastid and mitochondrial genes were performed. The Hydrodictyaceae contain the largest sphaeroplealean plastomes thus far fully sequenced. Conservation of plastome gene order within Hydrodictyaceae is striking compared with more dynamic patterns revealed across Sphaeropleales. Phylogenetic analyses resolve Hydrodictyon sister to a monophyletic Pediastrum, though the morphologically distinct P. angulosum and P. duplex continue to be polyphyletic. Analyses of plastid data supported the neochloridacean genus Chlorotetraëdron as sister to Hydrodictyaceae, while conflicting signal was found in the mitochondrial data. A phylogenomic approach resolved within-family relationships not obtainable with previous phylogenetic analyses. Denser taxon sampling across Sphaeropleales is necessary to capture patterns in plastome evolution, and further taxa and studies are needed to fully resolve the sister lineage to Hydrodictyaceae and polyphyly of Pediastrum angulosum and P. duplex.

Complete chloroplast genome of the medicinal plant Amomum compactum: gene organization, comparative analysis, and phylogenetic relationships within Zingiberales

BackgroundAmomum compactum is one of the basic species of the traditional herbal medicine amomi fructus rotundus, with great pharmacology effect. The system position of A. compactum is not clear yet, and the introduction of this plant has been hindered by many plant diseases. However, the correlational molecular studies are relatively scarce.MethodsThe total chloroplast (cp) DNA was extracted according to previous studies, and then sequenced by 454 GS FLX Titanium platform. Sequence assembly was complished by Newbler. Genome annotation was preformed by CPGAVAS and tRNA-SCAN. Then, general characteristics of the A. compactum cp genome and genome comparsion with three Zingiberaceae species was analyzed by corresponding softwares. Additionally, phylogenetical trees were reconstructed, based on the shared protein-coding gene sequences among 15 plant taxa by maximum parsimony (MP) and maximum likelihood (ML) methods.ResultsThe A. compactum cp genome with a classic quadripartite structure, consisting of a pair of reverse complement repeat regions (IRa/IRb) of 29,824 bp, a large single copy (LSC, 88,535 bp) region as well as a small single copy (SSC, 15,370 bp) region, is 163,553 bp in total size. The total GC content of this cp genome is 36.0%. The A. compactum cp genome owns 135 functional genes, that 113 genes are unique, containing eighty protein-coding genes, twenty-nine tRNA (transfer RNA) genes and four rRNA (ribosomal RNA) genes. Codon usage of the A. compactum cp genome is biased toward codons ending with A/T. Total 58 SSR loci and 24 large repeats are detected in the A. compactum cp genome. Relative to three other Zingiberaceae cp genomes, the A. compactum cp genome exhibits an obvious expansion in the IR regions. In A. compactum cp genome, the ycf1 pseudogene is 2969 bp away from the IRa/SSC border, whereas in other Zingiberaceae species, it is only 4–5 bp away from the IRa/SSC border. Comparative cp genome sequences analysis of A. compactum with other Zingiberaceae reveals that the gene order and gene content differ slightly among Zingiberaceae species. The phylogenetic analysis based on 67 protein-coding gene sequences supports the phylogenetic position of A. compactum.ConclusionsThe study has identified unique features of the A. compactum cp genome which would be helpful for us to understand the cp genome evolution and offer useful information for phylogenetics and further studies of this traditional medicinal plant.

The complete mitochondrial genome of Thryssa hamiltonii (Engraulinae, Engranlidae, Clupeoidei) and phylogenetic studies of Clupeoidei

The complete mitochondrial genome of Thryssa hamiltonii has been determined. The whole sequence was 16,894 bp in length and included 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and one control region (D-loop). The overall base composition is A 30.77%, C 27.97%, G 16.25%, T 25.01%, with a slightly A + T bias of 55.78%. With the exception of ND6 and eight tRNA genes, all other mitochondrial genes are encoded on the heavy strand. Three tandem repeat sequences were observed in the control region. Phylogenetic tree was constructed based on 13 protein-coding genes sequences of 21 clupeoidei species, and the result showed that T. hamiltonii is most closely related to T. dussumieri. These mitogenome sequence data will be useful for population genetics and phylogenetic analysis of the Clupeoidei.

The complete mitochondrial genome of Gadus chalcogramma and phylogenetic analysis

Here, we report the complete mitochondrial genome sequence of Alaska Pollock, Gadus chalcogramma. The genome sequence was obtained via long PCR reactions using universal primer sets for the fish mitochondrial genome. The total length of mitochondrial genome was 16,571 bp, consisting of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and one control region (D-loop). Except for ND6 and eight tRNA genes, all of the other mitochondrial genes were encoded on the heavy strand. The NJ tree of the combined 13 protein-coding gene sequences of Gadus chalcogramma possesses a relatively closer relationship with the Okhotsk Sea Pollock (Gadus chalcogrammus). Our complete mitochondrial genome will be valuable resources for phylogeny, genetics, and conservation of the Gadus chalcogramma in Korea.

The complete mitochondrial genome of Thryssa setirostris (Engraulinae, Engranlidae, Clupeoidei) and phylogenetic studies of Clupeoidei

The complete mitochondrial genome of Thryssa setirostris has been determined. The whole sequence was 16,926 bp in length and included 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 control region (D-loop). The overall base composition is A 29.88%, C 29.23%, G 16.77%, T 24.02%, with a slightly A + T bias of 54%. With the exception of ND6 and eight tRNA genes, all other mitochondrial genes are encoded on the heavy strand. Phylogenetic tree was constructed based on 13 protein-coding genes sequences of 21 clupeoidei species, and the result showed that T. setirostris is most closely related to T. dussumieri. These mitogenome sequence data will be useful for population genetics and phylogenetic analysis of the Clupeoidei.

The complete mitochondrial genome of Epinephelus Chlorostigma (Serranidae; Epinephelus) with phylogenetic consideration of Epinephelus

The complete mitochondrial genome of Epinephelus chlorostigma has been determined. It contains a typically conserved structure including 13 protein-coding genes, 22 tRNA genes, two rRNA genes and one control region, and the whole sequence was 16,894 bp in length. The overall base composition is A 28.70%, C 27.99%, G 16.28%, T 27.03%. Except ND6 and eight tRNA genes, all other mitochondrial genes are encoded on the heavy strand. Phylogenetic tree was constructed based on 12 protein-coding genes sequences of 18 Epinephelus species, four species belong to Cephalopholis as outgroup, the result showed that E. chlorostigma is most closely related to Epinephelus areolatus. We sequenced the complete mitochondrial genome of E. chlorostigma to enrich the resource of molecular markers for examination of phylogenetic relationships in Epinephelus.

Light in the darkness: New perspective on lanternfish relationships and classification using genomic and morphological data

Massive parallel sequencing allows scientists to gather DNA sequences composed of millions of base pairs that can be combined into large datasets and analyzed to infer organismal relationships at a genome-wide scale in non-model organisms. Although the use of these large datasets is becoming more widespread, little to no work has been done in estimating phylogenetic relationships using UCEs in deep-sea fishes. Among deep-sea animals, the 257 species of lanternfishes (Myctophiformes) are among the most important open-ocean lineages, representing half of all mesopelagic vertebrate biomass. With this relative abundance, they are key members of the midwater food web where they feed on smaller invertebrates and fishes in addition to being a primary prey item for other open-ocean animals. Understanding the evolution and relationships of midwater organisms generally, and this dominant group of fishes in particular, is necessary for understanding and preserving the underexplored deep-sea ecosystem. Despite substantial congruence in the evolutionary relationships among deep-sea lanternfishes at higher classification levels in previous studies, the relationships among tribes, genera, and species within Myctophidae often conflict across phylogenetic studies or lack resolution and support. Herein we provide the first genome-scale phylogenetic analysis of lanternfishes, and we integrate these data from across the nuclear genome with additional protein-coding gene sequences and morphological data to further test evolutionary relationships among lanternfishes. Our phylogenetic hypotheses of relationships among lanternfishes are entirely congruent across a diversity of analyses that vary in methods, taxonomic sampling, and data analyzed. Within the Myctophiformes, the Neoscopelidae is inferred to be monophyletic and sister to a monophyletic Myctophidae. The current classification of lanternfishes is incongruent with our phylogenetic tree, so we recommend revisions that retain much of the traditional tribal structure and recognize five subfamilies instead of the traditional two subfamilies. The revised monophyletic taxonomy of myctophids includes the elevation of three former lampanyctine tribes to subfamilies. A restricted Lampanyctinae was recovered sister to Notolychninae. These two clades together were recovered as the sister group to the Gymnoscopelinae. Combined, these three subfamilies were recovered as the sister group to a clade composed of a monophyletic Diaphinae sister to the traditional Myctophinae. Our results corroborate recent multilocus molecular studies that infer a polyphyletic Myctophum in Myctophinae, and a para- or polyphyletic Lampanyctus and Nannobrachium within Lampanyctinae. We resurrect Dasyscopelus and Ctenoscopelus for the independent clades traditionally classified as species of Myctophum, and we place Nannobrachium into the synonymy of Lampanyctus.

The complete mitochondrial genome sequence and phylogenetic analysis of white-lipped deer (Cervus albirostris)

The white-lipped deer (Cervus albirostris) is a rare and endangered species endemic to the Qinghai–Tibet Plateau in China. In this study, we determined the complete mitochondrial genome sequence of the white-lipped deer. The total length of the mitogenome is 16473 bp. It has the typical structure of a vertebrate mitochondrial genome including 13 protein-coding genes(PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a control region (D-loop). tRNASer lacks the stem of the dihydrouridine (DHU) arm and fails to form a typical cloverleaf secondary structure. A phylogenetic tree based on the 13 PCG sequences of nine species closely related to white-lipped deer shows the classification status of the species within the family Cervidae. These results provide new molecular biology information to further understand the diversity of the white-lipped deer, which will help to protect this population.

The complete mitochondrial genome of the alvinocaridid shrimp Shinkaicaris leurokolos (Decapoda, Caridea): Insight into the mitochondrial genetic basis of deep-sea hydrothermal vent adaptation in the shrimp

Deep-sea hydrothermal vent is one of the most extreme environments on Earth with low oxygen and high levels of toxins. Decapod species from the family Alvinocarididae have colonized and successfully adapted to this extremely harsh environment. Mitochondria plays a vital role in oxygen usage and energy metabolism, thus it may be under selection in the adaptive evolution of the hydrothermal vent shrimps. In this study, the mitochondrial genome (mitogenome) of alvinocaridid shrimp Shinkaicaris leurokolos (Kikuchi & Hashimoto, 2000) was determined through Illumina sequencing. The mitogenome of S. leurokolos was 15,903bp in length, containing 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. The gene order and orientation were identical to those of sequenced alvinocaridids. It has the longest concatenated sequences of protein-coding genes, tRNAs and shortest pooled rRNAs among the alvinocaridids. The control regions (CRs) of alvinocaridid were significantly longer (p<0.01) than those of the other caridaen. The alignment of the alvinocaridid CRs revealed two conserved sequence blocks (CSBs), and each of the CSBs included a noncanonical open reading frame (ORF), which may be involved in adjusting mitochondrial energy metabolism to adapt to the hydrothermal environment. Phylogenetic analysis supported that the deep-sea hydrothermal vent shrimps may have originated from those living in shallow area. Positive selection analysis reveals the evidence of adaptive change in the mitogenome of Alvinocarididae. Thirty potentially important adaptive residues were identified, which were located in atp6, cox1, cox3, cytb and nad1–5. This study explores the mitochondrial genetic basis of hydrothermal vent adaptation in alvinocaridid for the first time, and provides valuable clues regarding the adaptation.

Molecular phylogeny of Bradyrhizobium bacteria isolated from root nodules of tribe Genisteae plants growing in southeast Poland

The phylogeny of 16 isolates from root nodules of Genista germanica, Genista tinctoria, Cytisus ratisbonensis, and Cytisus scoparius growing in southeast Poland was estimated by comparative sequence analysis of core (16S rDNA, atpD, glnII, recA) and symbiosis-related (nodC, nodZ, nifH) genes. All the sequences analyzed placed the studied rhizobia in the genus Bradyrhizobium. Phylogenetic analysis of individual and concatenated housekeeping genes showed that the Genisteae microsymbionts form a homogeneous group with Bradyrhizobium japonicum strains. The phylogeny of nodulation and nitrogen fixation genes indicated a close relationship of the examined rhizobia with B. japonicum, Bradyrhizobium canariense, Bradyrhizobium cytisi, Bradyrhizobium rifense and Bradyrhizobium lupini strains infecting other plants of the tribe Genisteae. For the first time, the taxonomic position of G. germanica and C. ratisbonensis rhizobia, inferred from multigenic analysis, is described. The results of the phylogenetic analysis based on the protein-coding gene sequences presented in this study also indicate potential pitfalls concerning the choice of marker and reference strains, which may lead to conflicting conclusions in species delineation.

Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops.

The use of transcriptome data in the study of the population genetics of a species can capture faint signals of both genetic variation and expression variation and can provide a broad picture of a species' genomic response to environmental conditions. In this study, we characterized the genetic and expression diversity of Miscanthus lutarioriparius by comparing more than 16,225 transcripts obtained from 78 individuals, belonging to 10 populations distributed across the species' entire geographic range. We only observed a low level of nucleotide diversity (π=0.000434) among the transcriptome data of these populations, which is consistent with highly conserved sequences of functional elements and protein-coding genes captured with this method. Tests of population divergence using the transcriptome data were consistent with previous microsatellite data but proved to be more sensitive, particularly if gene expression variation was considered as well. For example, the analysis of expression data showed that genes involved in photosynthetic processes and responses to temperature or reactive oxygen species stimuli were significantly enriched in certain populations. This differential gene expression was primarily observed among populations and not within populations. Interestingly, nucleotide diversity was significantly negatively correlated with expression diversity within populations, while this correlation was positive among populations. This suggests that genetic and expression variation play separate roles in adaptation and population persistence. Combining analyses of genetic and gene expression variation represents a promising approach for studying the population genetics of wild species and may uncover both adaptive and nonadaptive processes.

Optimizing genomic medicine in epilepsy through a gene-customized approach to missense variant interpretation.

Gene panel and exome sequencing have revealed a high rate of molecular diagnoses among diseases where the genetic architecture has proven suitable for sequencing approaches, with a large number of distinct and highly penetrant causal variants identified among a growing list of disease genes. The challenge is, given the DNA sequence of a new patient, to distinguish disease-causing from benign variants. Large samples of human standing variation data highlight regional variation in the tolerance to missense variation within the protein-coding sequence of genes. This information is not well captured by existing bioinformatic tools, but is effective in improving variant interpretation. To address this limitation in existing tools, we introduce the missense tolerance ratio (MTR), which summarizes available human standing variation data within genes to encapsulate population level genetic variation. We find that patient-ascertained pathogenic variants preferentially cluster in low MTR regions (P < 0.005) of well-informed genes. By evaluating 20 publicly available predictive tools across genes linked to epilepsy, we also highlight the importance of understanding the empirical null distribution of existing prediction tools, as these vary across genes. Subsequently integrating the MTR with the empirically selected bioinformatic tools in a gene-specific approach demonstrates a clear improvement in the ability to predict pathogenic missense variants from background missense variation in disease genes. Among an independent test sample of case and control missense variants, case variants (0.83 median score) consistently achieve higher pathogenicity prediction probabilities than control variants (0.02 median score; Mann-Whitney U test, P < 1 × 10−16). We focus on the application to epilepsy genes; however, the framework is applicable to disease genes beyond epilepsy.

The complete mitochondrial genome of the causative agent of the human cercarial dermatitis, the visceral bird schistosome species Trichobilharzia szidati (platyhelminthes: Trematoda: Schistosomatidae)

We report the first complete mitochondrial genome of visceral bird schistosome Trichobilharzia szidati (Platyhelminthes, Trematoda: Schistosomatidae). The circular genome is 14293 bp in length and contains 12 protein-coding genes, 12S and 16S rRNAs genes, 22 tRNAs and one non-coding region (202 bp) (accession number MF136777). Phylogenetic relationships based on 12 protein-coding gene sequences (PCG) of mitogenomes of a number of trematode and cestode species have shown that T. szidati is the closest genetic relative to nasal bird schistosome T. regenti. The complete mitogenome sequence of T. szidati may serve as a resource for comparative mitogenomics and trematode evolution studies.

Genetics of Atrial Fibrillation: State of the Art in 2017.

Genetic variation is an important determinant of atrial fibrillation (AF) susceptibility. Numerous rare variants in protein-coding sequences of genes have been associated with AF in families and in early-onset cases, and chromosomal loci harbouring common risk variants have been mapped in AF cohorts. Many of these loci are in non-coding regions of the human genome and are thought to contain regulatory sequences that modulate gene expression. Disease genes implicated to date have predominantly encoded cardiac ion channels, with predicted mutation effects on the atrial action potential duration. More recent studies have expanded the spectrum of disease-associated genes to include myocardial structural components and have highlighted an unsuspected role for cardiac transcription factors. These paradigm-shifting discoveries suggest that abnormalities of atrial specification arising during cardiac development might provide a template for AF in later adult life. With the escalating pace of variant discovery, there is an increasing need for mechanistic studies not only to evaluate single variants, but also to determine the collective effects of each person's burden of rare and common genetic variants, co-morbidities and lifestyle factors on the atrial substrate for arrhythmogenesis. Elucidation of an individual's genetic predisposition and modifiable environmental risk factors will facilitate personalised approaches to AF treatment.

The complete mitochondrial genome sequence of Changbai Mountains wild boar (Cetartiodactyla: Suidae)

The Changbai Mountains wild boar (Cetartiodactyla: Suidae) is an omnivorous, gregarious and non-hibernating animal that inhabits in Jilin province of China and is threatened by poaching and deforestation. In this study, we provide the first report of the complete mitogenome of Changbai Mountains wild boar. The total length of mitochondrial genome is 16,720 bp and consists of 22 transfer RNA genes, 13 protein-coding genes, two ribosomal RNA genes, and a major non-coding control region (D-loop region). The D-loop region contains one repeat motif. The full complement of genes in the Changbai Mountains wild boar are located on the H-strand, except for the ND6 gene and eight tRNA genes that are encoded on the L-strand. Phylogenetic analysis based on combining 13 protein-coding gene sequences of several species showed that it is closely related to genus Sus scrofa in north China. The mitochondrial genome of Changbai Mountains wild boar is useful for population genetics and phylogenetic analysis.

A Phylogeny Based on the Concatenated 9 Protein-coding Gene Sequences of Mt Genomes among Three Gibbon Genera

The objective of this study was to reconstruct phylogenetic relationships of gibbons (3 genera, 11 species) deduced from complete sequenced mitochondrial genome sequences. According to conserved (C) / variable (V) sites ratio test, 9 selected protein coding genes were combined into a sequence with 9356 bases long including gaps. A resolved phylogenetic tree was obtained for the mitochondrial genome in the maximum-likelihood and maximum-parsimony analyses. In accordance with all previous morphological and molecular evidence, our results clearly supported monophyly of family Hylobatidae with predominantly strongly supported and each of the three genera with high support based on these coding genes. Among three genera, first Nomascus, next Symphalangus and at last Hylobates diverged, which identical to other previous research based on the whole mitochondrial genome. Our phylogenetic relationships of Nomascus group accord with Chan et al. (2010), only support values of a node of N. gabriellae (90%) were slight lower than Chan's result (97%). Among genus Hylobates, in our result, 6 species are involved and H. pileatus, H. lar, H. klossii, H. agilis, H. moloch and H. muelleri diverged sequentially.