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.     

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.     

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 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.     

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.     

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 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.     

Multilocus phylogeography of the Wedge-billed Woodcreeper Glyphorynchus spirurus (Aves,Furnariidae) in lowland Amazonia: Widespread cryptic diversity and paraphyly reveal a complex diversification pattern

Amazonian rivers function as important barriers to dispersal of Amazonian birds. Studying population genetics of lineages separated by rivers may help us to uncover the dynamics of biological diversification in the Amazon. We reconstructed the phylogeography of the Wedge-billed Woodcreeper, Glyphorynchus spirurus (Furnariidae) in the Amazon basin. Sampling included 134 individuals from 63 sites distributed in eight Amazonian areas of endemism separated by major Amazonian rivers. Nucleotide sequences were generated for five genes: two mtDNA genes (1047 bp for cyt b and 1002 bp for ND2) and three nuclear genes (647 bp from the sex-linked gene ACO, 319 bp from the intron of G3PDH, and 619 bp from intron 2 of MYO). In addition, 37 individuals were randomly selected from the Rondônia and Inambari areas of endemism for genomic fingerprinting, using five ISSR primers. Our results reveal allopatric and well-supported lineages within G. spirurus with high levels of genetic differentiation (p-distances 0.9–6.3%) across opposite banks of major Amazonian rivers. The multilocus phylogenetic reconstructions obtained reveal several incongruences with current subspecies taxonomy. Within currently recognized subspecies, we found high levels of both paraphyly and genetic differentiation, indicating deep divergences and strong isolation consistent with species-level differences. ISSR fingerprinting supports the existence of genetically differentiated populations on opposite sides of the Madeira River. Molecular dating suggests an initial vicariation event isolating populations from the Guiana center of endemism during the Late Miocene/Early Pliocene, while more recent events subdivided Brazilian Shield populations during the Lower Pleistocene.     

Molecular and biochemical characterizations of a novel arthropod endo-β-1,3-glucanase from the Antarctic springtail,Cryptopygus antarcticus,horizontally acquired from bacteria

Collembolan species have been known to have β-1,3-glucanase activity and yet the genes coding such enzymes have not been demonstrated. We report here a novel arthropod endo-β-1,3-glucanase gene CaLam from the Antarctic springtail, Cryptopygus antarcticus. The open reading frame consists of 813 bp encoding 270 amino acids with a putative signal peptide and a typical motif of glycosyl hydrolase family 16 (GHF16), E–I–D–I–T–E. The recombinant protein expressed in E. coli shows the hydrolytic activity toward laminarin (K_m ～ 9.98 mg/mL) with an optimal temperature 50 °C and an optimal pH 6.0. CaLam digests laminarin and laminarioligosaccharides except laminaribiose as an endo-β-1,3-glucanase, releasing glucose, laminaribiose and laminaritriose as the major products. Analyses of molecular phylogeny of CaLam and its protein structure reveal that CaLam is closely related with bacterial β-1,3-glucanases more than with the eukaryotic homologues. Even so, the genomic structure of the CaLam gene consisting of six exons interspersed with approximately 57 to 63 bp introns confirms that it is endogenous in the genome of the Antarctic springtail. These results suggest that CaLam should have been transferred from bacteria to the lineage of the Collembolan species by horizontal gene transfer.     

Analysis of digestion of rice planthopper by Pardosa pseudoannulata based on CO-I gene

In order to systematically study the predatory behavior and digestion regularity of spiders, real-time fluorescence quantification PCR technique was used to detect the number of CO-I genes in Pardosa pseudoannulata after it preyed on rice planthoppers in different temperatures within different periods. At 28 °C, 0, 1, 2, 4, 8, 16, and 24 h after P. pseudoannulata preyed on rich planthopper, DNA was extracted from cephalothorax and abdomen of P. pseudoannulata. Routine PCR and real-time fluorescence PCR techniques were employed for CO-I gene amplification. The results show that: The prey liquid was temporarily stored in the sucking stomach of the spider head within 2 h after prey, and gradually transferred to the midgut of the abdomen with the prolongation of time. After 4 h, CO-I gene residues of rice planthopper in the cephalothorax gradually decreased. The CO-I gene of rice planthopper was basically transferred to the abdomen after 16 h. During 0–1 h, food contained in abdominal midgut and other digestive organs was very small, CO-I gene detection was not obvious. Over time, food entered into the midgut from the sucking stomach for digestion. During 2–4 h, CO-I gene amount increased, at 2–4 h, detected CO-I gene residue reached the peak; but rapidly declined after 8, 16, and 24 h, even it is still detectable. The results at different temperatures reveal that: As the temperature increased from 26 °C to 32 °C, CO-I gene residues of rich planthopper in cephalothorax and abdomen of P. pseudoannulata gradually decreased, which indicated that the digestion rate increased with the increase of temperature with some range. However, when the temperature continued to increase to 34 °C, the digestion rate decreased.     

Phylogenetic structure of the Thomomys bottae–umbrinus complex in North America

The phylogeography of the Thomomys bottae–umbrinus complex in the United States and Mexico was assessed with sequences of the mitochondrial cytochrome b gene. These sequences were obtained from 225 individuals representing 108 locations over the range, including 56 sequences from GenBank. 110 (500 bp) sequences were used for Bayesian inference and neighbor-joining analyses, and 34 (1140 bp) specimens from the main clades obtained from the Bayesian inference were used in maximum-parsimony and maximum-likelihood analyses. The different analyses indicate significant variation within the species complex that averages 13% among major groups of genetic differences among Thomomys bottae–umbrinus. The overall pattern of geographic variation is not concordant with the current taxonomy. To the contrary, eight monophyletic groups are supported by all analyses and can be considered phylogenetic species. Overall divergence among these groups appears influenced by historical biogeographic events active during the Pliocene and Pleistocene.     

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.     

Copper,zinc superoxide dismutase of Curcuma aromatica is a kinetically stable protein

This study details on cloning and characterization of Cu,Zn superoxide dismutase (Ca–Cu,Zn SOD) from a medicinally important plant species Curcuma aromatica. Ca–Cu,Zn SOD was 692 bp with an open reading frame of 459 bp. Expression of the gene in Escherichia coli cells followed by purification yielded the enzyme with K_m of 0.047 ± 0.008 μM and V_max of 1250 ± 24 units/mg of protein. The enzyme functioned (i) across a temperature range of −10 to +80 °C with temperature optima at 20 °C; and (ii) at pH range of 6–9 with optimum activity at pH 7.8. Ca–Cu,Zn SOD retained 50% of the maximum activity after autoclaving, and was stable at a wide storage pH ranging from 3 to 10. The enzyme tolerated varying concentrations of denaturating agent, reductants, inhibitors, trypsin, was fairly resistant to inactivation at 80 °C for 180 min (k_d, 6.54 ± 0.17 × 10⁻³ min⁻¹; t_1/2, 106.07 ± 2.68 min), and had midpoint of thermal transition (T_m) of 70.45 °C. The results suggested Ca–Cu,Zn SOD to be a kinetically stable protein that could be used for various industrial applications.     

Gene cloning and expression of a detergent stable alkaline protease from Aspergillus clavatus ES1

The genes, cDNA alpES1 and alpES1, encoding Aspergillus clavatus ES1 alkaline protease were amplified from complementary DNA (cDNA) and genomic DNA, respectively, cloned in pCR^®II-TOPO plasmid and then sequenced. Sequence analysis of the cDNA alpES1 gene revealed an open reading frame (ORF) of 1212 bp encoding a pre–pro-protein of 403 amino acid residues consisting of a 21-aa signal peptide, a 100-aa pro-peptide and a 282-aa mature protein with a calculated molecular weight of 28.5 kDa. Compared to the cDNA alpES1 gene, the alpES1 gene contained three introns, which had 53, 57 and 54 bp, respectively. The cDNA alpES1 gene was then sub-cloned in pET-30b(+) and expressed in Escherichia coli BL21 (λDE3). The purified recombinant protease had a molecular weight of about 32 kDa estimated by SDS-PAGE. Kinetic parameters, K_m and k_cat values of the recombinant AlpES1 for casein, were 0.23 mM and 12.38 min⁻¹, respectively. The catalytic efficiency (k_cat/K_m) was 53.82 min⁻¹ mM⁻¹.     

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.