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.     

The complete mitochondrial genome of the endangered Apollo butterfly,Parnassius apollo (Lepidoptera: Papilionidae) and its comparison to other Papilionidae species

The Apollo butterfly, Parnassius apollo is a representative species of the butterfly subfamily Parnassiinae. This charming species is one of the most endangered butterfly species in the world. In this study, we sequenced its complete mitochondrial genome (mitogenome), with the aim of accumulating genetic information for further studies of population genetics and mitogenome evolution in the Papilionidae. The 15,404-bp long mitogenome harbors a typical set of 37 genes and is the largest butterfly mitogenome determined, except for Papilio maraho (16,094 bp). Like many other sequenced lepidopteran species, one tRNA^Trp-like and one tRNA^Leu(UUR)-like sequences were detected in the AT-rich region. A total of 164 bp of non-coding sequences are dispersed in 14 regions throughout the genome. The longest intergenic spacer (68 bp) is located between tRNA^Ser(AGN) and tRNA^Glu, and is the largest spacer at this location among Papilionidae species. This spacer may have resulted from an 8-fold repetition of a TTTCTTCT motif or a 4-fold repetition of a CTTTATTT motif.     

Complete mitochondrial DNA sequence of the ark shell Scapharca broughtonii: An ultra-large metazoan mitochondrial genome

The complete mitochondrial (mt) genome of the ark shell Scapharca broughtonii was determined using long PCR and a genome walking sequencing strategy with genus-specific primers. The S. broughtonii mt genome (GenBank accession number AB729113) contained 12 protein-coding genes (the atp8 gene is missing, as in most bivalves), 2 ribosomal RNA genes, and 42 transfer tRNA genes, in a length of 46,985 nucleotides for the size of mtDNA with only one copy of the heteroplasmic tandem repeat (HTR) unit. Moreover the S. broughtonii mt genome shows size variation; these genomes ranged in size from about 47 kb to about 50 kb because of variation in the number of repeat sequences in the non-coding region. The mt-genome of S. broughtonii is, to date, the longest reported metazoan mtDNA sequence. Sequence duplication in non-coding region and the formation of HTR arrays were two of the factors responsible for the ultra-large size of this mt genome. All the tRNA genes were found within the S. broughtonii mt genome, unlike the other bivalves usually lacking one or more tRNA genes. Twelve additional specimens were used to analyze the patterns of tandem repeat arrays by PCR amplification and agarose electrophoresis. Each of the 12 specimens displayed extensive heteroplasmy and had 8–10 length variants. The motifs of the HTR arrays are about 353–362 bp and the number of repeats ranges from 1 to 11.     

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.     

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.     

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.     

Complete mitochondrial genome of Membranipora grandicella (Bryozoa: Cheilostomatida) determined with next-generation sequencing: The first representative of the suborder Malacostegina

Next-generation sequencing (NGS) has proven a valuable platform for fast and easy obtaining of large numbers of sequences at relatively low cost. In this study we use shot-gun sequencing method on Illumina HiSeq 2000, to obtain enough sequences for the assembly of the bryozoan Membranipora grandicella (Bryozoa: Cheilostomatida) mitochondrial genome, which is the first representative of the suborder Malacostegina. The complete mitochondrial genome is 15,861 bp in length, which is relatively larger than other studied bryozoans. The mitochondrial genome contains 13 protein-coding genes, 2 ribosomal RNAs and 20 transfer RNAs. To investigate the phylogenetic position and the inner relationships of the phylum Bryozoa, phylogenetic trees were constructed with amino acid sequences of 11 PCGs from 30 metazoans. Two superclades of protostomes, namely Lophotrochozoa and Ecdysozoa, are recovered as monophyletic with strong support in both ML and Bayesian analyses. Somewhat to surprise, Bryozoa appears as the sister group of Chaetognatha with moderate or high support. The relationship among five bryozoans is Tubulipora flabellaris + (M. grandicella + (Flustrellidra hispida + (Bugula neritina + Watersipora subtorquata))), which supports for the view that Cheilostomatida is not a natural, monophyletic clade. NGS proved to be a quick and easy method for sequencing a complete mitochondrial genome.     

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.     

Mitochondrial distribution patterns in canine oocytes as related to the reproductive cycle stage

This study investigates the mitochondrial (mt) distribution in canine ovarian oocytes examined at recovery time, as related to the reproductive cycle stage, and in oviductal oocytes. Ovarian Germinal Vesicle (GV) stage oocytes were recovered from bitches in anestrous (A, n = 2), follicular phase (F, n = 4), ovulation (O, n = 2), early luteal (EL, n = 7) and mid/late luteal phase (MLL, n = 2). Oviductal GV, metaphase I (MI) or MII stage oocytes were recovered from six bitches between 56 and 110 h after ovulation. Mitochondria were revealed by using MitoTracker Orange CMTM Ros and confocal microscopy. In ovarian oocytes, three mt distribution patterns were found: (I) small aggregates diffused throughout the cytoplasm; (II) diffused tubular networks; (III) pericortical tubular networks. Significantly higher rates of oocytes showing heterogeneous mt patterns (II + III) were obtained from bitches in F (75%) and in O (96%) compared with bitches in A (31%; F vs. A: P < 0.05; O vs. A: P < 0.001), in EL (61%; O vs. EL: P < 0.01), or in MLL (0%; F vs. MLL: P < 0.05; O vs. MLL: P < 0.001). Fluorescence intensity did not vary according to mt distribution pattern except that it was lower in oocytes recovered in EL phase and showing small mt aggregations (P < 0.001). The majority of ovulated MII stage oocytes (79%) showed diffused tubular mt network. We conclude that mt distribution pattern of canine ovarian immature oocytes changes in relation to reproductive cycle stage and that patterns observed in oocytes recovered from bitches in periovulatory phases are heterogeneous and similar to those of in vivo matured oocytes.     

Limited phylogenetic distribution of a long tandem-repeat cluster in the mitochondrial control region in Bubo (Aves,Strigidae) and cluster variation in Blakiston’s fish owl (Bubo blakistoni)

To investigate the phylogenetic position of Blakiston’s fish owl (Bubo blakistoni), we sequenced the mitochondrial (mt) DNA control region and cytochrome b (cyt b) for nine Bubo species. Maximum-likelihood analyses of combined control region and cyt b sequences, and cyt b sequences alone, showed that species formerly placed in genus Ketupa comprise a monophyletic group. Unexpectedly, we discovered a long cluster of 20–25 tandem repeat units 77 or 78 bp long in the third control region domain in four of the nine Bubo species for which the control region was sequenced (B. blakistoni, B. flavipes, and B. ketupu in the Ketupa clade; B. lacteus), leading to overall control region lengths of 3.0–3.8 kpb estimated from agarose gel electrophoresis. The control region in B. lacteus is the longest (3.8 kbp) reported to date in vertebrates. Sequencing of eight repeat units at each end of the cluster in 20 B. blakistoni individuals detected several types of repeat units 77 or 78 bp long, and six patterns in the order of unit types. The occurrence of a repeat cluster in all three species examined in the Ketupa clade suggests their common ancestor also had a cluster, whereas a maximum parsimony tree showed repeat-unit types grouping by species, rather than by paralog groups, suggesting independent origins of the clusters. We reconcile these results with a turnover model, in which the range in cluster-length variation and unit types at the 5′ end are hypothetically functionally constrained by the protein-binding function of the control region, but otherwise there is a continual turnover of units in evolutionary time, with new unit types arising through mutations, proliferating by duplication of single and double repeat blocks, and being lost through deletion. Estimated free energies for reconstructed secondary structures of single and especially pairs of repeat units were higher than for homologous single-unit blocks in species lacking a repeat cluster, supporting slipped-strand mispairing as the mechanism of cluster turnover.     

Cloning and sequencing of three new putative toxin genes from Clostridium bifermentans CH18

Three new open reading frames were found downstream from cbm71, a toxin gene from Clostridium bifermentans malaysia (Cbm) strain CH18. The first one (91 bp downstream) called cbm72, is 1857 bp long and encodes a 71 727-Da protein (Cbm72) with a sequence similar to that of Bacillus thuringiensis delta-endotoxins. This protein shows no significant toxicity to mosquito larvae. The two others, cbm17.1 (462 bp) and cbm17.2 (459 bp), are copies of the same gene encoding Cbm P18 and P16 polypeptides and located 426 bp and 1022 bp downstream from cbm72, respectively. They encode 17 189-Da and 17 451-Da proteins with sequences 44.6% similar to that of Aspergillus fumigatus hemolysin; however, they were not hemolytic in the conditions tested.     

Expression of allograft inflammatory factor-1 (AIF-1) in response to bacterial challenge and tissue injury in the pearl oyster,Pinctada martensii

Allograft inflammatory factor-1 (AIF-1), an interferon (IFN)-γ-inducible calcium-binding cytokine, is associated with the inflammatory response and defense. We cloned and analyzed the expression pattern of the AIF-1 gene of the pearl oyster Pinctada martensii, hereafter designated PmAIF-1. The full-length PmAIF-1 cDNA is 946 bp in length and consists of a 5′-untranslated region (UTR) of 120 bp, a 3′-UTR of 376 bp, and an open reading frame (ORF) of 450 bp encoding a polypeptide of 149 amino acids with an estimated molecular mass of 17 kDa. Sequence analysis reveals that PmAIF-1 contains two EF hand Ca⁺²-binding motifs like those in previously characterized AIF-1s while alignment with known AIF-1 protein sequences reveals higher similarity to invertebrate orthologs than to those of vertebrates.Quantitative PCR analysis reveals that PmAIF-1 is constitutively expressed, with the highest expression detected in hemocytes, and the expression level of PmAIF-1 mRNA was significantly up-regulated in hemocytes, gill, digestive gland under bacterial challenge and tissue injury. After challenged by gram-negative bacteria Vibrio alginolyticus and Vibrio parahaemolyticus, gram-positive bacteria Bacillus subtilis, the expression level of this gene in hemocytes were all up-regulated and reached the maximum point at 12 h (5.80 folds, P < 0.01), 6 h (5.02 folds, P < 0.01) and 12 h (5.49 folds, P < 0.01), respectively. Under shell damage and mantle injury, PmAIF-1 mRNA increased gradually in the first 3 h and reached a peak of expression at 6 h post-injury. These findings suggest that PmAIF-1 is an acute-response protein involved in the innate immune responses of pearl oysters, and provide general information about the mechanisms of innate immune defense against bacterial infection in pearl oysters.     

Generation and analysis of expressed sequence tags from adductor muscle of Japanese scallop Mizuhopecten yessoensis

A normalized cDNA library was constructed from the adductor muscle of M. yessoensis and acquired 4595 high quality expressed sequence tags (ESTs). After clustering and assembly of the ESTs, 3061 unigenes containing 654 contigs and 2407 singletons were identified. The contig length ranged from 266 bp to 2364 bp and the average length of these contigs was 544 bp. Blastx nonredundant protein database analysis showed that 1522 unigenes had significant homology to known genes (E value ≤ 10^? 5). By comparing to Clusters of Orthologous Groups (COG) categories, 460 unigenes were annotated (E value ≤ 10^? 10). Using Kyoto Encyclopedia of Genes and Genomes (KEGG), 345 of 3061 unigenes were assigned into 103 pathways (E value ≤ 10^? 5). For InterProScan searches, 1237 unigenes were annotated containing 727 different types of protein domains. 941 of the 1237 unigenes were annotated for Gene Ontology (GO) classification using Uniprot2GO associations in any category (biological, cellular, and molecular). By sequences comparability and analysis of Blastx NCBI nonredundant protein database and KEGG, 66 unigenes were identified that may be involved in genetic information processing based on the known knowledge. The study provides a material basis as useful information for the genomic analysis of shellfish.     

Cloning of Wap65 in sea bass (Dicentrarchus labrax) and sea bream (Sparus aurata) and expression in sea bass tissues

The complementary DNA encoding WAP65 protein was cloned from the liver of two fish species sea bass (Dicentrarchus labrax) and sea bream (Sparus aurata). Full-length cDNA sequences were obtained from reverse transcribed total RNA, followed by 5′ and 3′ rapid amplification of cDNA end (RACE) experiments. The full-length cDNA sequence of D. labrax is 1709 bp and the coding sequence is flanked by a 67 bp 5′-UTR and a 358 bp 3′-UTR. The full-length cDNA sequence of S. aurata is 1599 bp, and the coding sequence is flanked by a 48 bp 5′-UTR and a 273 bp 3′-UTR. The deduced amino acid putative primary sequences are composed of 427 and 425 amino acid residues for D. labrax and S. aurata, respectively. They display high homologies with previously described fish WAP65 and other hemopexin-like proteins from rabbit (Oryctolagus cuniculus). Expression of Wap65 has proved to be a natural physiological adaptive answer of teleost fish to warm temperature acclimation. In all fish species studied to date, Wap65 was found expressed mainly by the liver, although other tissues seem able to express Wap65 in response to a warm temperature acclimation, in a specie specific manner. Here, we investigate the tissue specific expression of Wap65 in D. labrax and S. aurata in response to a warm temperature acclimation, by RT-PCR analysis.     

ITS2 sequence variations among members of Diachasmimorpha longicaudata complex (Hymenoptera: Braconidae) in Thailand

The parasitoid Diachasmimorpha longicaudata complex in Thailand contains at least 3 cryptic species informally designated as species D. longicaudata A, B and BB. DNA sequence data of nuclear ITS2 (second internal transcribed spacer) were used to characterize members of this D. longicaudata complex. The polymerase chain reaction (PCR) amplicon of ITS2 region of D. longicaudata B (≈ 650 bp) clearly differentiated this species from A and BB (amplicon of ≈ 590 bp). Sequence alignment of individual parasitoids revealed that low intraspecies differences ranged from 0.457 to 3.991%, but interspecies differences ranged from 7.566 to 12.989%. Phylogenetic trees constructed using Neighbor-Joining (NJ) and Maximum Parsimony (MP) methods, taking the parasitoid Psyttalia concolor complex as an outgroup, revealed that D. longicaudata A, B and BB formed a monophyletic group, with species A and BB being more closely related than species B. ITS2 characterization of D. longicaudata complex has revealed an interesting divergence of the three cryptic sibling species in Thailand.     

Characterization and PCR optimization of the thermostable family B DNA polymerase from Thermococcus guaymasensis

The gene encoding Thermococcus guaymasensis DNA polymerase (Tgu DNA polymerase) was cloned and sequenced. The 2328 bp Tgu DNA polymerase gene encoded a 775 amino acid residue protein. Alignment of the entire amino acid sequence revealed a high degree of sequence homology between Tgu DNA polymerase and other archaeal family B DNA polymerases. The Tgu DNA polymerase gene was expressed under the control of the T7lac promoter on pET-22b(+) in Escherichia coli BL21-CodonPlus(DE3)-RIL. The expressed enzyme was then purified by heat treatment followed by two steps of chromatography. The optimum pH and temperature were 7.5 and 80 °C, respectively. The optimal buffer for PCR with Tgu DNA polymerase consisted of 50 mM Tris–HCl (pH 8.2), 4 mM MgCl₂, 50 mM KCl, and 0.02% Triton X-100. Tgu DNA polymerase revealed 4-fold higher fidelity (3.17 × 10^?6) than Taq DNA polymerase (12.13 × 10^?6) and a faster amplification rate than Taq and Pfu DNA polymerases. Tgu DNA polymerase had an extension rate of 30 bases/s and a processivity of 150 nucleotides (nt). Thus, Tgu DNA polymerase has some faster elongation rate and a higher processivity than Pfu DNA polymerase. Use of different ratios of Taq and Tgu DNA polymerases determined that a ratio of 4:1 efficiently facilitated long PCR (approximately 15 kb) and a 3-fold lower error rate (4.44 × 10^?6) than Taq DNA polymerase.