Complete nucleotide sequence and organization of the mitochondrial genome of eri-silkworm,Samia cynthia ricini (Lepidoptera: Saturniidae)

Samia cynthia ricini is a commercial silk-producing insect that is now reared year-round in Korea, with the expectation of being utilized for diverse purposes. In this report, we present the complete mitochondrial genome (mitogenome) of S. c. ricini. The 15,384-bp long S. cynthia ricini mitogenome was amplified into 26 short fragments using three long overlapping fragments using primers designed from reported lepidopteran mitogenome sequences. The genome comprises 37 genes (13 protein-coding genes, two rRNA genes, and 22 tRNA genes), and one large non-coding region termed the A + T-rich region. The A/T content of the third codon position was 91.7%, which was 18.8% and 21.6% higher than those of first and second codon positions, respectively. The high A/T content in the genome is reflected in codon usage, accounting for 39.5% of A/T-composed codons (TTA, ATT, TTT, and ATA). Unlike a previous report on the start codon for the COI gene, the S. c. ricini COI gene commences with a typical ATT codon. A total of 221 bp of non-coding sequences are dispersed in 17 regions, ranging in size from 1 to 54 bp, which comprise 1.4% of the total genome. One of the non-coding sequence located between tRNA^Gln and ND2 (54 bp) has 77% sequence homology with the 5′-sequence of the neighboring ND2 gene, suggesting partial duplication of the sequence during evolution. The 361-bp long A + T-rich region contains an 18 bp-long poly-T stretch, ATAGA motif, ATTTA element, microsatellite-like A/T sequence, poly-A stretch and one tRNA-like sequence, as typically found in Lepidoptera including Bombycoidea.     

Complete mitochondrial genome sequence of the Arctic gammarid,Onisimus nanseni (Crustacea; Amphipoda): Novel gene structures and unusual control region features

To analyze the mitogenome of the amphipod Onisimus nanseni, we amplified the complete mitogenome of O. nanseni using long-PCR and genome walking techniques. The mitogenome of O. nanseni is circular and contains all the typical mt genes (2 rRNAs, 22 tRNAs, and 13 protein-coding genes). It has two peculiar non-coding regions of 148 bp and 194 bp. The latter can be involved in replication and termination processes. The total length of the pooled protein-coding, rRNA, and tRNA genes is shorter than those of other crustaceans. In addition, the intergenic spacers of the O. nanseni mitogenome are considerably shorter in length than those of other crustaceans. Fourteen adjacent genes overlap, resulting in a compact mitogenomic structure. In the O. nanseni mitogenome, the AT composition is elevated, particularly in the control regions (78.9% AT), as has been demonstrated for two other amphipods. The tRNA order is highly rearranged compared to other arthropod mitogenomes, but the order of protein-coding genes and rRNAs is largely conserved. The gene cluster between the CO1 and CO3 genes is completely conserved among all amphipods compared. This provides insights into the evolution and gene structures of crustacean mitochondrial genomes, particularly in amphipods.     

Complete mitochondrial genome of Coelomactra antiquata (Mollusca: Bivalvia): The first representative from the family Mactridae with novel gene order and unusual tandem repeats

The complete mitochondrial genome plays an important role in the accurate inference of phylogenetic relationships among metazoans. Mactridae, also known as trough shells or duck clams, is an important family of marine bivalve clams in the order Veneroida. Here we present the complete mitochondrial genome sequence of the Xishishe Coelomactra antiquata (Mollusca: Bivalvia), which is the first representative from the family Mactridae. The mitochondrial genome of C. antiquata is of 17,384 bp in length, and encodes 35 genes, including 12 protein-coding, 21 transfer RNA, and 2 ribosomal RNA genes. Compared with the typical gene content of animal mitochondrial genomes, atp8 and tRNAS₂ are missing. Gene order of the mitochondrial genome of C. antiquata is unique compared with others from Veneroida. In the mitochondrial genome of the C. antiquata, a total of 2189 bp of non-coding nucleotides are scattered among 26 non-coding regions. The largest non-coding region contains one section of tandem repeats (99 bp × 11), which is the second largest tandem repeats found in the mitochondrial genomes from Veneroida. The phylogenetic trees based on mitochondrial genomes support the monophyly of Veneridae and Lucinidae, and the relationship at the family level: ((Veneridae + Mactridae) + (Cardiidae + Solecurtidae)) + Lucinidae. The phylogenetic result is consistent with the morphological classification. Meanwhile, bootstrap values are very high (BP = 94–100), suggesting that the evolutionary relationship based on mitochondrial genomes is very reliable.     

Description of complete mitochondrial genome of the black-veined white,Aporia crataegi (Lepidoptera: Papilionoidea), and comparison to papilionoid species

The black-veined white, Aporia crataegi (Lepidoptera: Papilionoidea) is nearly extinct in South Korea, although substantial numbers of dried specimens are available. One of the common practices used to rescue such endangered species is to launch a re-introduction program after a proper amount of genetic information is analyzed from donor and donee populations. In this study, we sequenced the complete mitochondrial genome (mitogenome) of A. crataegi to accumulate genetic information for subsequent population studies and to further understand the mitogenome evolution in true butterflies, Papilionoidea. The 15,140-bp long A. crataegi mitogenome has typical sets of 37 genes and is the smallest among the true butterfly species, with overall slightly smaller size genes and regions throughout the genome. The A/T content of the genome (81.3%) is the highest in Pieridae, where A. crataegi belongs, but lower than that of the lycaenid species (81.7%–82.7%). Unlike the diversified or modified usage of an anticodon for tRNA^Ser(AGN), the species of Pieridae including A. crataegi all contain GCT that has been hypothesized as being ancestral for Lepidoptera. A total of 111 bp of non-coding sequences are interspersed in 13 regions, ranging in size from 1–49 bp. Among these sequences, relatively longer ones (≥ 16 bp) all have relatively higher sequence identity to other regions of the genome, suggesting partial duplication of the sequences during A. crataegi evolution.     

A physical map of the genome of ethanol fermentative bacterium Zymomonas mobilis ZM4 and localization of genes on the map

A physical map of the Zymomonas mobilis ZM4 genome has been constructed from the results of reciprocal Southern hybridization with PmeI, PacI, and NotI-digested genomic DNA fragments and linking cosmid clones. Restriction enzyme-digested Z. mobilis ZM4 genome was electrophoresed with phage lambda DNA concatemers as a size standard in a Bio-Rad CHEF-DRII pulsed-field gel electrophoresis (PFGE) system. The restriction enzyme PmeI generated 15 fragments (3–625 kb), and PacI produced 19 fragments (7–525 kb). Each size of restriction fragment was calculated by comparison to the size of phage lambda DNA concatemers, and the genome size of Z. mobilis ZM4 was estimated to be 2085.5 kb. The 19 known genes and three rrn operons were localized on the map.     

Mitochondrial tRNA cleavage by tRNA-targeting ribonuclease causes mitochondrial dysfunction observed in mitochondrial disease

Mitochondrial DNA (mtDNA) is a genome possessed by mitochondria. Since reactive oxygen species (ROS) are generated during aerobic respiration in mitochondria, mtDNA is commonly exposed to the risk of DNA damage. Mitochondrial disease is caused by mitochondrial dysfunction, and mutations or deletions on mitochondrial tRNA (mt tRNA) genes are often observed in mtDNA of patients with the disease. Hence, the correlation between mt tRNA activity and mitochondrial dysfunction has been assessed. Then, cybrid cells, which are constructed by the fusion of an enucleated cell harboring altered mtDNA with a ρ⁰ cell, have long been used for the analysis due to difficulty in mtDNA manipulation. Here, we propose a new method that involves mt tRNA cleavage by a bacterial tRNA-specific ribonuclease. The ribonuclease tagged with a mitochondrial-targeting sequence (MTS) was successfully translocated to the mitochondrial matrix. Additionally, mt tRNA cleavage, which resulted in the decrease of cytochrome c oxidase (COX) activity, was observed.     

Next generation molecular diagnosis of mitochondrial disorders

Mitochondrial disorders are by far the most genetically heterogeneous group of diseases, involving two genomes, the 16.6 kb mitochondrial genome and ~ 1500 genes encoded in the nuclear genome. For maternally inherited mitochondrial DNA disorders, a complete molecular diagnosis requires several different methods for the detection and quantification of mtDNA point mutations and large deletions. For mitochondrial disorders caused by autosomal recessive, dominant, and X-linked nuclear genes, the diagnosis has relied on clinical, biochemical, and molecular studies to point to a group of candidate genes followed by stepwise Sanger sequencing of the candidate genes one-by-one. The development of Next Generation Sequencing (NGS) has revolutionized the diagnostic approach. Using massively parallel sequencing (MPS) analysis of the entire mitochondrial genome, mtDNA point mutations and deletions can be detected and quantified in one single step. The NGS approach also allows simultaneous analyses of a group of genes or the whole exome, thus, the mutations in causative gene(s) can be identified in one-step. New approaches make genetic analyses much faster and more efficient. Huge amounts of sequencing data produced by the new technologies brought new challenges to bioinformatics, analytical pipelines, and interpretation of numerous novel variants. This article reviews the clinical utility of next generation sequencing for the molecular diagnoses of complex dual genome mitochondrial disorders.     

Complete mitochondrial genome sequence of <Emphasis Type="Italic">Cucullaea labiata</Emphasis> (Arcoida: Cucullaeidae) and phylogenetic implications

The complete mitochondrial genome of Cucullaea labiata (Arcoida: Cucullaeidae) was firstly determined in this study in order to better understand the phylogenetic relationship between Cucullaeidae and Arcidae. The C. labiata mitochondrial genome was 25,845 bp in size and contained 12 protein-coding genes, 2 rRNA and 22 tRNA genes. The number and the location of the tRNA genes were different from three Arcidae species (Scapharca broughtonii, Scapharca kagoshimensis and Tegillarca granosa). Gene arrangement also differed dramatically. The length of the non-coding regions was 10,559 bp, in which the largest one (6057 bp) included eight point nine copies of a 659 bp repeat motif. The number of repeated sequences was different in different individuals, similar to the findings from the mitochondrial genome of S. broughtonii and Placopecten magellanicus. One intron was found in cox1 gene both in CL_98 and in CL_99 individuals of C. labiata. The reason why mitochondrial introns are retained so scarcely in bivalve taxa needs further research. Phylogenetic analyses based on 12 concatenated amino acid sequences of protein-coding genes supported Cucullaeidae was the sister group of Arcidae.     

The complete mitogenome of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) and comparison with other lepidopteran insects

The complete mitochondrial genome (mitogenome) of Bombyx mori strain Dazao (Lepidoptera: Bombycidae) was determined to be 15,653 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a A + T-rich region. It has the typical gene organization and order of mitogenomes from lepidopteran insects. The AT skew of this mitogenome was slightly positive and the nucleotide composition was also biased toward A + T nucleotides (81.31%). All PCGs were initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which was initiated by CGA. The cox1 and cox2 genes had incomplete stop codons consisting of just a T. All the tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The A + T-rich region of the mitogenome was 495 bp in length and consisted of several features common to the lepidopteras. Phylogenetic analysis showed that the B. mori Dazao was close to Bombycidae.     

The complete mitochondrial genome of the mantid shrimp Oratosquilla oratoria (Crustacea: Malacostraca: Stomatopoda): Novel non-coding regions features and phylogenetic implications of the Stomatopoda

The complete mitochondrial (mt) genome sequence of Oratosquilla oratoria (Crustacea: Malacostraca: Stomatopoda) was determined; a circular molecule of 15,783 bp in length. The gene content and arrangement are consistent with the pancrustacean ground pattern. The mt control region of O. oratoria is characterized by no GA-block near the 3′ end and different position of [TA(A)]n-blocks compared with other reported Stomatopoda species. The sequence of the second hairpin structure is relative conserved which suggests this region may be a synapomorphic character for the Stomatopoda. In addition, a relative large intergenic spacer (101 bp) with higher A + T content than that in control region was identified between the tRNA^Glu and tRNA^Phe genes. Phylogenetic analyses based on the current dataset of complete mt genomes strongly support the Stomatopoda is closely related to Euphausiacea. They in turn cluster with Penaeoidea and Caridea clades while other decapods form a separate group, which rejects the monophyly of Decapoda. This challenges the suitability of Stomatopoda as an outgroup of Decapoda in phylogenetic analyses. The basal position of Stomatopoda within Eumalacostraca according to the morphological characters is also questioned.     

The complete mitochondrial genome of Sarcoptes scabiei var. nyctereutis from the Japanese raccoon dog: Prediction and detection of two transfer RNAs (tRNA-A and tRNA-Y)

Sarcoptes scabiei (Acari: Sarcoptidae) causes a common contagious skin disease that affects many mammals. Here, the complete mitochondrial genome of a mite, S. scabiei var. nyctereutis, from Japanese wild raccoon dogs was analyzed. The 13,837 bp circular genome contained 13 protein-coding genes, two rRNA genes, and 22 tRNA genes. For the first time, two tRNAs (alanine and tyrosine), that were thought to be absent in scabies mites from other animals, were predicted to have short, non-cloverleaf structures by in silico annotation and detected by RT-PCR, sequencing, and northern analysis. The mitochondrial genome structure of S. scabiei is similar to that of Psoroptes cuniculi and Dermatophagoides farinae. While small and unusual tRNA genes seem to be common among acariform mites, further experimental evidence for their presence is needed. Furthermore, through an analysis of the cox1 gene, we have provided new evidence to confirm the transmission of this mite between different animal hosts.     

Species trees for spiny lizards (Genus Sceloporus): Identifying points of concordance and conflict between nuclear and mitochondrial data

Spiny lizards (genus Sceloporus) represent one of the most diverse and species rich clades of squamate reptiles in continental North America. Sceloporus contains 90+ species, which are partitioned into 21 species groups containing anywhere from one to 15 species. Despite substantial progress towards elucidating the phylogeographic patterns for many species of Sceloporus, efforts to resolve the phylogenetic relationships among the major species groups remain limited. In this study, the phylogenetic relationships of 53 species of Sceloporus, representing all 21 species groups, are estimated based on four nuclear genes (BDNF, PNN, R35, RAG-1; >3.3 kb) and contrasted with a new mitochondrial DNA genealogy based on six genes (12S, ND1, ND4, and the histidine, serine, and leucine tRNA genes; >2.5 kb). Species trees estimated from the nuclear loci using data concatenation or a coalescent-based inference method result in concordant topologies, but the coalescent approach provides lower resolution and support. When comparing nuclear versus mtDNA-based topologies for Sceloporus species groups, conflicting relationships outnumber concordant relationships. Incongruence is not restricted to weak or unresolved nodes as might be expected under a scenario of rapid diversification, but extends to conflicts involving strongly support clades. The points of concordance and conflict between the nuclear and mtDNA data are discussed, and arguments for preferring the species trees estimated from the multilocus nuclear data are presented.     

Shigella flexneri type-specific antigen V: cloning,sequencing and characterization of the glucosyl transferase gene of temperate bacteriophage SfV

With lysogeny by bacteriophage SfV, Shigella flexneri serotype Y is converted to serotype 5a. The glucosyl transferase gene (gtr) from bacteriophage SfV of S. flexneri, involved in serotype-specific conversion, was cloned and characterized. The DNA sequence of a 3.7 kb EcoRI–BamHI fragment of bacteriophage SfV which includes the gtr gene was determined. This gene, encoding a polypeptide of 417 aa with 47.67 kDa molecular mass, caused partial serotype conversion of S. flexneri from serotype Y to type V antigen as demonstrated by Western blotting and the sensitivity of the hybrid strain to phage Sf6. The deduced protein of the partially sequenced open reading frame upstream of the gtr showed similarity to various glycosyl transferases of other bacteria. Orf3, separated from the gtr by a non-coding region and transcribed convergently, codes for a 167 aa (18.8 kDa) protein found to have homology with tail fibre genes of phage lambda and P2.     

Molecular characterization of a cloned dolphin mitochondrial genome

Summary DNA clones have been isolated that span the complete mitochondrial (mt) genome of the dolphin,Cephalorhynchus commersonii. Hybridization experiments with purified primate mtDNA probes have established that there is close resemblance in the general organization of the dolphin mt genome and the terrestrial mammalian mt genomes. Sequences covering 2381 bp of the dolphin mt genome from the major noncoding region, three tRNA genes, and parts of the genes encoding cytochrome b, NADH dehydrogenase subunit 3 (ND3), and 16S rRNA have been compared with corresponding regions from other mammalian genomes. There is a general tendency throughout the sequenced regions for greater similarity between dolphin and bovine mt genomes than between dolphin and rodent or human mt genomes.     

Characterisation of the complete mitochondrial genome of Helice wuana (Grapsoidea: Varunidae) and comparison with other Brachyuran crabs

The mitochondrial genome (mitogenome) provides important information for phylogenetic analysis and understanding evolutionary origins. Herein, we sequenced, annotated, and characterised the mitogenome of the crab Helice wuana to better understand its molecular evolution and phylogeny. The 16,359 bp mitogenome includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one control region. The genome composition is highly A + T biased 68.42%, and exhibits a negative AT–skew (? 0.036) and GC–skew (? 0.269) among Brachyura crabs. Gene rearrangements were detected, as was tandem duplication followed by random loss, which explains the translocation of mitochondrial genes. Phylogenetic analysis showed that H. wuana and H. tientsinensis clustered on one branch with high nodal support values. These results confirm that the placement of H. wuana within the Varunidae family of Thoracotrematan crabs. This study will provided a better understanding for gene rearrangements and crab evolution in the further.     

BAC library for the amphipod crustacean,Parhyale hawaiensis

Bacterial artificial chromosomes (BACs) are capable of propagating large fragments of DNA and have become an invaluable tool for studying genome biology. To fill a phylogenetic gap in available genomic sequence and to complement ongoing molecular and genetic studies, we have generated a BAC library for the marine amphipod crustacean, Parhyale hawaiensis. The library was generated from genomic DNA isolated from whole adult animals and comprises 129,024 individual clones. The median insert size is ～ 140 kb and the genomic coverage is approximately five genome equivalents. We have screened the Parhyale BAC library for developmentally relevant genes and characterized the genomic organization of four genes: a hedgehog ortholog, and three Pax3/7 paralogs. Preliminary analysis suggests that introns are larger and more prevalent in Parhyale than in other arthropods whose genomes have been sequenced, which may partly account for the large genome size in Parhyale.     

A unique tRNA gene family and a novel,highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS^− AGY and trnQ genes, a trnM isomerism, but it lacks trnS^− CUN. We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP → trnS^− AGY) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, “Oxidored_q1_N”, which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.     

Complete mtDNA of Meretrix lusoria (Bivalvia: Veneridae) reveals the presence of an atp8 gene, length variation and heteroplasmy in the control region

The complete nucleotide sequence of the mitochondrial genome of the clam Meretrix lusoria (Bivalvia: Veneridae) was determined. It comprises 20,268 base pairs (bp) and contains 13 protein-coding genes, including ATPase subunit 8 (atp8), two ribosomal RNAs, 22 transfer RNAs, and a non-coding control region. The atp8 encodes a protein of 39 amino acids. All genes are encoded on the same strand. A putative control region (CR or D-loop) was identified in the major non-coding region (NCR) between the tRNA^Gly and tRNA^Gln. A 1087 bp tandem repeat fragment was identified that comprises nearly 11 copies of a 101 bp motif and accounts for approximately 41% of the NCR. The 101 bp tandem repeat motif of the NCR can be folded into a stem–loop secondary structure. Samples of eight individuals from Hainan and Fujian provinces were collected and their NCR regions were successfully amplified and sequenced. The data revealed a highly polymorphic VNTR (variable number of tandem repeats) associated with high levels of heteroplasmy in the D-loop region. The size of the CR ranged from 1942 to 3354 bp depending upon the copy number of the repeat sequence.