Phylogenetic and genomic analysis of the complete mitochondrial DNA sequence of the spotted asparagus beetle Crioceris duodecimpunctata

We report the complete mitochondrial DNA sequence of the spotted asparagus beetle, Crioceris duodecimpunctata. The genome complement, gene order, and nucleotide composition of this beetle's mitochondrial genome were found to be typical of those reported for other insects. Unusual features of this genome include the substitution of UCU for GCU as the anticodon for tRNA(Ser), an unusual TpsiC loop for the tRNA(Ile) gene, and the identification of a putative ATT start codon for cox1. The utility of complete mitochondrial genome data for phylogenetic inference of the insect orders was tested, and compared to that of cox1 and combined mitochondrial ribosomal DNA sequences. Even though the number of insect orders represented by complete mitochondrial genomes is still limited, several well-established relationships are evident in the phylogenetic analysis of the complete sequences. Monophyly of the orders Diptera, Lepidoptera, and Coleoptera were consistently recovered. Monophyly of the Holometabola was also observed in some (though not all) analyses. The accumulation of complete mitochondrial sequences from a broader array of insect orders holds the promise of clarifying the early diversification of insects.     

First Sequenced Mitochondrial Genome from the Phylum Acanthocephala(Leptorhynchoides thecatus) and Its Phylogenetic Position Within Metazoa

The complete sequence of the mitochondrial genome of Leptorhynchoides thecatus (Acanthocephala) was determined, and a phylogenetic analysis was carried out to determine its placement within Metazoa. The genome is circular, 13,888 bp, and contains at least 36 of the 37 genes typically found in animal mitochondrial genomes. The genes for the large and small ribosomal RNA subunits are shorter than those of most metazoans, and the structures of most of the tRNA genes are atypical. There are two significant noncoding regions (377 and 294 bp), which are the best candidates for a control region; however, these regions do not appear similar to any of the control regions of other animals studied to date. The amino acid and nucleotide sequences of the protein coding genes of L. thecatus and 25 other metazoan taxa were used in both maximum likelihood and maximum parsimony phylogenetic analyses. Results indicate that among taxa with available mitochondrial genome sequences, Platyhelminthes is the closest relative to L. thecatus, which together are the sister taxon of Nematoda; however, long branches and/or base composition bias could be responsible for this result. The monophyly of Ecdysozoa, molting organisms, was not supported by any of the analyses. This study represents the first mitochondrial genome of an acanthocephalan to be sequenced and will allow further studies of systematics, population genetics, and genome evolution.Reviewing Editor: Dr. Rafael Zardoya The entire genome sequence has been deposited with the GenBank Data Libraries under-accession number AY562383.     

The complete mitochondrial genome of Parafronurus youi (Insecta: Ephemeroptera) and phylogenetic position of the Ephemeroptera

The first complete mitochondrial genome of a mayfly, Parafronurus youi (Arthropoda: Insecta: Pterygota: Ephemeroptera: Heptageniidae), was sequenced using a long PCR-based approach. The genome is a circular molecule of 15,481 bp in length, and encodes the set of 38 genes. Among them, 37 genes are found in other conservative insect mitochondrial genomes, and the 38(th) unique gene is trnM-like (trnM2). The duplication-random loss model can be used to explain one of the translocations at least. The A+T content of the control region is 57%, the lowest proportion detected so far in Hexapoda. Based on the nucleotide dataset and the corresponding amino acid dataset of 12 protein-coding genes, Bayesian inference and maximum likelihood analyses yielded stable support for the relationship of the three basal clades of winged insects as Ephemeroptera+(Odonata+Neoptera).     

昆虫基因组及其大小

昆虫基因组大小是由于基因组各种重复序列在扩增、缺失和分化过程中所致的数量差异造成的。这些差异使得昆虫不同类群间、种间和同种的不同种群间表现出基因组大小的不同。目前有59种昆虫已经列入基因组测序计划, 其中6种昆虫（黑腹果蝇Drosophila melanogaster、冈比亚按蚊Anopheles gambiae、家蚕Bombyx mori、意大利蜜蜂Apis mellifera、埃及伊蚊Aedes aegypti和赤拟谷盗Tribolium castaneum）的全基因组序列已经报道。有725种昆虫的基因组大小得到了估计, 大小在0.09~16.93 pg （88~16 558 Mb）之间。本文还介绍了昆虫基因组大小的估计方法, 讨论了昆虫基因组大小的变化及其意义。     

Defense mechanisms against herbivory in Picea: sequence evolution and expression regulation of gene family members in the phenylpropanoid pathway

Background

The insect order Neuroptera encompasses more than 5,700 described species. To date, only three neuropteran mitochondrial genomes have been fully and one partly sequenced. Current knowledge on neuropteran mitochondrial genomes is limited, and new data are strongly required. In the present work, the mitochondrial genome of the ascalaphid owlfly Libelloides macaronius is described and compared with the known neuropterid mitochondrial genomes: Megaloptera, Neuroptera and Raphidioptera. These analyses are further extended to other endopterygotan orders.

Results

The mitochondrial genome of L. macaronius is a circular molecule 15,890 bp long. It includes the entire set of 37 genes usually present in animal mitochondrial genomes. The gene order of this newly sequenced genome is unique among Neuroptera and differs from the ancestral type of insects in the translocation of trnC. The L. macaronius genome shows the lowest A+T content (74.50%) among known neuropterid genomes. Protein-coding genes possess the typical mitochondrial start codons, except for cox1, which has an unusual ACG. Comparisons among endopterygotan mitochondrial genomes showed that A+T content and AT/GC-skews exhibit a broad range of variation among 84 analyzed taxa. Comparative analyses showed that neuropterid mitochondrial protein-coding genes experienced complex evolutionary histories, involving features ranging from codon usage to rate of substitution, that make them potential markers for population genetics/phylogenetics studies at different taxonomic ranks. The 22 tRNAs show variable substitution patterns in Neuropterida, with higher sequence conservation in genes located on the α strand. Inferred secondary structures for neuropterid rrnS and rrnL genes largely agree with those known for other insects. For the first time, a model is provided for domain I of an insect rrnL. The control region in Neuropterida, as in other insects, is fast-evolving genomic region, characterized by AT-rich motifs.

Conclusions

The new genome shares many features with known neuropteran genomes but differs in its low A+T content. Comparative analysis of neuropterid mitochondrial genes showed that they experienced distinct evolutionary patterns. Both tRNA families and ribosomal RNAs show composite substitution pathways. The neuropterid mitochondrial genome is characterized by a complex evolutionary history.     

两种毛蚜线粒体基因组比较分析

【目的】线粒体基因组分析已被应用于昆虫系统发育研究。本研究以蚜科Aphididae重要类群毛蚜亚科物种为代表,测定并比较分析了该类蚜虫的线粒体基因组特征,探讨了基于线粒体基因组信息的蚜虫系统发育关系重建。【方法】以毛蚜亚科三角枫多态毛蚜Periphyllus acerihabitans Zhang和针茅小毛蚜Chaetosiphella stipae Hille Ris Lambers,1947为研究对象,利用长短PCR相结合的方法测定线粒体基因组的序列,分析了基因组的基本特征;基于在线t RNAscan-SE Search Server搜索方法预测了t RNA的二级结构;基于12个物种(本研究获得的2个物种和10个Gen Bank上下载的物种数据)的蛋白编码基因(PCGs)序列,利用最大似然法和贝叶斯法重建了蚜科的系统发育关系。【结果】两种毛蚜均获得了约94%的线粒体基因组数据,P.acerihabitans获得了14 908 bp,控制区为1 205 bp;C.stipae获得了13 893 bp,控制区为609 bp。两种毛蚜同时获得33个基因,包含接近完整的13个蛋白编码基因(PCGs)(nad5不完整),18个tRNA,2个rRNA基因;ka/ks值表明,C.stipae的进化速率更快。从基因组组成、基因排列顺序、核苷酸组成分析、密码子使用情况、t RNA二级结构等特征来分析,两种蚜虫线粒体基因组基本特征相似。系统发育重建结果表明毛蚜亚科、蚜亚科的单系性得到了支持,毛蚜亚科位于蚜科的基部位置。【结论】两种毛蚜线粒体基因组的基本特征相似,符合蚜虫线粒体基因组的一般特征,两种线粒体基因组的长度差异主要来自控制区长度的不同;系统发育重建支持毛蚜亚科与蚜亚科的单系性,毛蚜亚科位于蚜科较为基部的位置。研究结果为蚜虫类系统发育重建提供了参考。     

Molecular identification of the first SIFamide receptor

SIFamide is the short name and also the C terminus of the Drosophila neuropeptide AYRKPPFNGSIFamide. SIFamide has been isolated or predicted from various insects and crustaceans, and appears to be extremely well conserved among these arthropods. However, the function of this neuropeptide is still enigmatic. Here, we have identified the Drosophila gene (CG10823) coding for the SIFamide receptor. When expressed in Chinese hamster ovary cells, the receptor is only activated by Drosophila SIFamide (EC(50), 2x10(-8)M) and not by a library of 32 other insect neuropeptides and eight biogenic amines. Database searches revealed SIFamide receptor orthologues in the genomes from the malaria mosquito Anopheles gambiae, the silkworm Bombyx mori, the red flour beetle Tribolium castaneum, and the honey bee Apis mellifera. An alignment of the five insect SIFamide or SIFamide-like receptors showed, again, an impressive sequence conservation (67-77% amino acid sequence identities between the seven-transmembrane areas; 82-87% sequence similarities). The identification of well-conserved SIFamide receptor orthologues in all other insects with a sequenced genome, suggests that the SIFamide/receptor couple must have an essential function in arthropods. This paper is the first report on the identification of a SIFamide receptor.     

The complete sequence of the mitochondrial genome of Daphnia pulex (Cladocera: Crustacea).

The sequence of the mitochondrial DNA (mtDNA) of the branchiopod crustacean Daphnia pulex has been completed. It is 15333bp with an A+T content of 62.3%, and contains the typical complement of 13 protein-coding, 22 transfer RNA (tRNA) and two ribosomal RNA (rRNA) genes. Comparison of this sequence with the sequences of the other eight completely sequenced arthropod mtDNAs showed that gene order and orientation are identical to that of Drosophila but different from Artemia due to the rearrangement of two tRNA genes. Nucleotide composition, codon usage, and amino acid composition are very similar in the crustaceans, but divergent from insects and chelicerates which show a much higher bias towards A+T. However, with few exceptions, the mitochondrial proteins of Daphnia are more similar to those of the dipteran insects (Drosophila and Anopheles) than to those of Artemia, at both the nucleotide and amino acid levels, suggesting that Artemia mtDNA is evolving at an accelerated rate. These results also show that sequence evolution and the evolution of nucleotide composition can be decoupled. Analysis of nucleotide substitution patterns in COII showed that there has been an unbiased acceleration of the overall substitution rate in Artemia. In contrast, the accelerated substitution rate in Apis is due partly to extreme A+T mutation pressure. Secondary structures are proposed for the Daphnia tRNAs and rRNAs. The tRNAs are similar to those of other arthropods but tend to have TPsiC arms that are only 4bp long. The rRNA secondary structures are similar to those proposed for insects except for the absence of a small number of helices in Daphnia. Phylogenetic analysis of second codon positions grouped Daphnia with Artemia, as expected, despite the latter's accelerated divergence rate. In contrast, the unusual pattern of mtDNA divergence in Apis led to a topology in which the holometabolous insects (Anopheles, Drosophila, Apis) appeared to be paraphyletic with respect to the hemimetabolous insect, Locusta, due to the early branching of Apis.     

Complete Mitochondrial Genome of the Free-Living Earwig, Challia fletcheri (Dermaptera: Pygidicranidae) and Phylogeny of Polyneoptera

The insect order Dermaptera, belonging to Polyneoptera, includes ～2,000 extant species, but no dermapteran mitochondrial genome has been sequenced. We sequenced the complete mitochondrial genome of the free-living earwig, Challia fletcheri, compared its genomic features to other available mitochondrial sequences from polyneopterous insects. In addition, the Dermaptera, together with the other known polyneopteran mitochondrial genome sequences (protein coding, ribosomal RNA, and transfer RNA genes), were employed to understand the phylogeny of Polyneoptera, one of the least resolved insect phylogenies, with emphasis on the placement of Dermaptera. The complete mitochondrial genome of C. fletcheri presents the following several unusual features: the longest size in insects is 20,456 bp; it harbors the largest tandem repeat units (TRU) among insects; it displays T- and G-skewness on the major strand and A- and C-skewness on the minor strand, which is a reversal of the general pattern found in most insect mitochondrial genomes, and it possesses a unique gene arrangement characterized by a series of gene translocations and/or inversions. The reversal pattern of skewness is explained in terms of inversion of replication origin. All phylogenetic analyses consistently placed Dermaptera as the sister to Plecoptera, leaving them as the most basal lineage of Polyneoptera or sister to Ephemeroptera, and placed Odonata consistently as the most basal lineage of the Pterygota.     

双翅目昆虫线粒体基因组结构特点及其测序通用引物的设计和应用

研究双翅目昆虫线粒体基因组的结构特点, 并设计其测序的通用引物, 为今后双翅目昆虫线粒体基因组的研究提供参考和依据。利用比较基因组学和生物信息学方法, 分析了已经完全测序的26个双翅目昆虫线粒体基因组的结构特点、 碱基组成和保守区, 并据此设计了双翅目昆虫基因组测序的通用引物。结果表明： 双翅目昆虫线粒体基因组长14 503~19 517 bp, 其结构保守, 含有37个编码基因, 包括13个蛋白质编码基因, 22个tRNA编码基因和2个rRNA编码基因, 此外还包含一段长度差异很大的非编码区(AT富含区)。基因组内基因排列次序稳定, 除个别基因外, 其余都与黑腹果蝇Drosophila melanogaster基因排列次序一致。基因组的碱基组成不均衡, AT含量在72.59%~85.15%之间, 碱基使用存在偏向性, 偏好使用AC碱基。全基因组的核苷酸和氨基酸序列保守, 共鉴定了11个保守区。在保守区内共设计了26对双翅目线粒体基因组测序通用引物, 扩增的目标片段都在1 200 bp以内。将该套通用引物用于葱蝇Delia antiqua线粒体全基因组测序, 结果证明其高效、 合用。     

The mitochondrial genomes of the human hookworms, Ancylostoma duodenale and Necator americanus (Nematoda: Secernentea).

The complete mitochondrial genome sequences were determined for two species of human hookworms, Ancylostoma duodenale (13,721 bp) and Necator americanus (13,604 bp). The circular hookworm genomes are amongst the smallest reported to date for any metazoan organism. Their relatively small size relates mainly to a reduced length in the AT-rich region. Both hookworm genomes encode 12 protein, two ribosomal RNA and 22 transfer RNA genes, but lack the ATP synthetase subunit 8 gene, which is consistent with three other species of Secernentea studied to date. All genes are transcribed in the same direction and have a nucleotide composition high in A and T, but low in G and C. The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. For both hookworm species, genes were arranged in the same order as for Caenorhabditis elegans, except for the presence of a non-coding region between genes nad3 and nad5. In A. duodenale, this non-coding region is predicted to form a stem-and-loop structure which is not present in N. americanus. The mitochondrial genome structure for both hookworms differs from Ascaris suum only in the location of the AT-rich region, whereas there are substantial differences when compared with Onchocerca volvulus, including four gene or gene-block translocations and the positions of some transfer RNA genes and the AT-rich region. Based on genome organisation and amino acid sequence identity, A. duodenale and N. americanus were more closely related to C. elegans than to A. suum or O. volvulus (all secernentean nematodes), consistent with a previous phylogenetic study using ribosomal DNA sequence data. Determination of the complete mitochondrial genome sequences for two human hookworms (the first members of the order Strongylida ever sequenced) provides a foundation for studying the systematics, population genetics and ecology of these and other nematodes of socio-economic importance.     

Identification of the Drosophila and Tribolium receptors for the recently discovered insect RYamide neuropeptides

One year ago, we discovered a new family of insect RYamide neuropeptides, which has the C-terminal consensus sequence FFXXXRYamide, and which is widely occurring in most insects, including the fruitfly Drosophila melanogaster and the red flour beetle Tribolium castaneum (F. Hauser et al., J. Proteome Res. 9 (2010) 5296–5310). Here, we identify a Drosophila G-protein-coupled receptor (GPCR) coded for by gene CG5811 and its Tribolium GPCR ortholog as insect RYamide receptors. The Drosophila RYamide receptor is equally well activated (EC₅₀, 1 × 10⁻⁹ M) by the two Drosophila RYamide neuropeptides: RYamide-1 (PVFFVASRYamide) and RYamide-2 (NEHFFLGSRYamide), both contained in a preprohormone coded for by gene CG40733. The Tribolium receptor shows a somewhat higher affinity to Tribolium RYamide-2 (ADAFFLGPRYamide; EC₅₀, 5 × 10⁻⁹ M) than to Tribolium RYamide-1 (VQNLATFKTMMRYamide; EC₅₀, 7 × 10⁻⁸ M), which might be due to the fact that the last peptide does not completely follow the RYamide consensus sequence rule. There are other neuropeptides in insects that have similar C-terminal sequences (RWamide or RFamide), such as the FMRFamides, sulfakinins, myosuppressins, neuropeptides F, and the various short neuropeptides F. Amazingly, these neuropeptides show no cross-reactivity to the Tribolium RYamide receptor, while the Drosophila RYamide receptor is only very slightly activated by high concentrations (>10⁻⁶ M) of neuropeptide F and short neuropeptide F-1, showing that the two RYamide receptors are quite specific for activation by insect RYamides, and that the sequence FFXXXRYamide is needed for effective insect RYamide receptor activation. Phylogenetic tree analyses and other amino acid sequence comparisons show that the insect RYamide receptors are not closely related to any other known insect or invertebrate/vertebrate receptors, including mammalian neuropeptide Y and insect neuropeptide F and short neuropeptide F receptors. Gene expression data published in Flybase (www.flybase.org) show that the Drosophila CG5811 gene is significantly expressed in the hindgut of adult flies, suggesting a role of insect RYamides in digestion or water reabsorption.     

The mitochondrial genome of Ruspolia dubia (Orthoptera: Conocephalidae) contains a short A+T-rich region of 70 bp in length.

The complete sequence (14 971 bp) of the Ruspolia dubia mitochondrial genome was determined and annotated. The genome contains the gene content, base composition, and codon usage typical of metazoan mitochondrial genomes. All 37 genes are conserved in the positions observed most frequently in insect mitochondrial genome structures. The secondary structures of both small subunit and large subunit rRNA were predicted. The most unusual features found were the initiation codon (TTA) of COI and a short A+T-rich region of 70 bp in length. In addition, a short, highly conserved polythymidine stretch that was previously described in Orthoptera and Diptera was also present in the A+T-rich region.     

Horizontal Escape of the Novel Tc1-Like Lepidopteran Transposon TCp3.2 into Cydia pomonella Granulovirus

We characterized an insertion mutant of the baculovirus Cydia pomonella granulovirus (CpGV), which contained a transposable element of 3.2 kb. This transposon, termed TCp3.2, has unusually long inverted terminal repeats (ITRs) of 756 bp and encodes a defective gene for a putative transposase. Amino acid sequence comparison of the defective transposase gene revealed a distant relationship to a putative transposon in Caenorhabditis elegans which also shares some similarity of the ITRs. Maximum parsimony analysis of the predicted amino acid sequences of Tc1- and mariner-like transposases available from the GenBank data base grouped TCp3.2 within the superfamily of Tc1-like transposons. DNA hybridization indicated that TCp3.2 originated from the genome of Cydia pomonella, which is the natural host of CpGV, and is present in less than 10 copies in the C. pomonella genome. The transposon TCp3.2 most likely was inserted into the viral genome during infection of host larvae. TCp3.2 and the recently characterized Tc1-like transposon TC14.7 (Jehle et al. 1995), which was also found in a CpGV mutant, represent a new family of transposons found in baculovirus genomes. The occasional horizontal escape of different types of host transposons into baculovirus genomes evokes the question about the possible role of baculoviruses as an interspecies vector in the horizontal transmission of insect transposons. Received: 27 February 1997 / Accepted: 16 May 1997     

Sequence and comparative analysis of the maize NB mitochondrial genome

The NB mitochondrial genome found in most fertile varieties of commercial maize (Zea mays subsp. mays) was sequenced. The 569,630-bp genome maps as a circle containing 58 identified genes encoding 33 known proteins, 3 ribosomal RNAs, and 21 tRNAs that recognize 14 amino acids. Among the 22 group II introns identified, 7 are trans-spliced. There are 121 open reading frames (ORFs) of at least 300 bp, only 3 of which exist in the mitochondrial genome of rice (Oryza sativa). In total, the identified mitochondrial genes, pseudogenes, ORFs, and cis-spliced introns extend over 127,555 bp (22.39%) of the genome. Integrated plastid DNA accounts for an additional 25,281 bp (4.44%) of the mitochondrial DNA, and phylogenetic analyses raise the possibility that copy correction with DNA from the plastid is an ongoing process. Although the genome contains six pairs of large repeats that cover 17.35% of the genome, small repeats (20-500 bp) account for only 5.59%, and transposable element sequences are extremely rare. MultiPip alignments show that maize mitochondrial DNA has little sequence similarity with other plant mitochondrial genomes, including that of rice, outside of the known functional genes. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly three-fourths of the maize NB mitochondrial genome is still of unknown origin and function.     

The complete mitochondrial genome of the cockroach <Emphasis Type="Italic">Eupolyphaga sinensis</Emphasis> (Blattaria: Polyphagidae) and the phylogenetic relationships within the Dictyoptera

We present the complete mitochondrial DNA sequence of Eupolyphaga sinensis. This closed circular molecule is 15553 bp long and consists of 37 genes that encode for 13 inner membrane proteins, 2 ribosomal RNAs and 22 transfer RNAs. The genome shares the gene order and orientation with previously known Blattaria mitochondrial genomes. All tRNAs could be folded into the typical cloverleaf secondary structure, but the tRNASer (AGN) appears to be missing the DHU arm. The A + T-rich region is 857 bp long and longer than other cockroaches. Based on the concatenated amino acid sequences of all protein coding genes of E. sinensis in conjunction with those 23 other arthropod sequences, we reconstruct the phylogenetic tree. Phylogenetic analyses shows that Blataria (including Isoptera) and the Mantodea are sister groups. Furthermore the relationship of the three basal clades of winged insects are different from the three previous hypotheses ((Ephemeroptera + Odonata) +Neoptera, Ephemeroptera + (Odonata + Neoptera), Odonata + (Ephemeroptera +Neoptera)). The Ephemeroptera (Parafronurus youi) clusters with the Plecoptera (Pteronarcys princes).     

The partial mitochondrial genome of the Cephalothrix rufifrons (Nemertea, Palaeonemertea): characterization and implications for the phylogenetic position of Nemertea

A continuous 10.1kb fragment of the Cephalothrix rufifrons (Nemertea, Palaeonemertea) mitochondrial genome was sequenced and characterized to further assess organization of protostome mitochondrial genomes and evaluate the phylogenetic potential of gene arrangement and amino acid characters. The genome is A-T rich (72%), and this biased base composition is partly reflected in codon usage. Inferred tRNA secondary structures are typical of those reported for other metazoan mitochondrial DNAs. The arrangement of the 26 genes contained in the fragment exhibits marked similarity to those of many protostome taxa, most notably molluscs with highly conserved arrangements and a phoronid. Separate and simultaneous phylogenetic analyses of inferred amino acid sequences and gene adjacencies place the nemertean within the protostomes among coelomate lophotrochozoan taxa, but do not find a well-supported sister taxon link.     

Overlapping of the VP2-VP3 gene and the VP1 gene in the SV40 genome.

The nucleotide sequence of the SV40 Hind E fragment has been determined mainly by the partial chemical degradation procedure of Maxam and Gilbert (1977). The sequence of the strand with the same polarity as the late messenger RNA shows only one open reading frame for translation. Considering that VP3 corresponds to the carbosyl terminal part of VP2, and considering various evidence which indicates that the SV40 Hind E segment is part of the amino acid sequence of VP2-VP3. It continues clockwise in Hind K, where it terminates with a UAA signal. The latter is located 110 nucleotides beyond the initiation signal for the major structural protein VP1 (Fiers et al., 1975; Van de Voorde et al., 1976). Hence this small overlapping region of the genome codes for the synthesis of three different proteins in two different reading frames. The deduced amino acid sequence covers a major part of the vp3 poly peptide, and the amino acid composition is in good agreement with published values (Greenaway and Levine, 1973).