Comparative analysis of mitochondrial genome data for Necator americanus from two endemic regions reveals substantial genetic variation

Necator americanus is a blood-sucking, intestinal nematode of major human health importance in many tropical and subtropical regions of the world. The aim of the present study was to compare the complete mitochondrial genome sequence from one N. americanus individual from Togo with another from China, in order to estimate the magnitude of genetic variability for different mitochondrial genes and non-coding regions. For the 12 protein genes, this comparison revealed sequence differences at both the nucleotide (3-7%) and amino acid (1-7%) levels. The most conserved of these was the nad4L gene, whereas the nad1 gene was least conserved at both the nucleotide and amino acid levels. Nucleotide differences were also detected in 14 of the 22 transfer RNAs (trns) (1-13%), the AT-rich region ( approximately 8%), non-coding regions (8-25%) and in the small (rrnS) and large (rrnL) subunits of mitochondrial ribosomal RNA (rrn) ( approximately 1%). Comparison of the rrnL sequences among multiple individual worms revealed nine unequivocal nucleotide differences between N. americanus from the two countries. Consistent with previous studies, these findings provide evidence for substantial genetic variation within N. americanus, which may have implications for the transmission and control of hookworm disease.     

Complete Mitochondrial Genome of Haplorchis taichui and Comparative Analysis with Other Trematodes

Mitochondrial genomes have been extensively studied for phylogenetic purposes and to investigate intra- and interspecific genetic variations. In recent years, numerous groups have undertaken sequencing of platyhelminth mitochondrial genomes. Haplorchis taichui (family Heterophyidae) is a trematode that infects humans and animals mainly in Asia, including the Mekong River basin. We sequenced and determined the organization of the complete mitochondrial genome of H. taichui. The mitochondrial genome is 15,130 bp long, containing 12 protein-coding genes, 2 ribosomal RNAs (rRNAs, a small and a large subunit), and 22 transfer RNAs (tRNAs). Like other trematodes, it does not encode the atp8 gene. All genes are transcribed from the same strand. The ATG initiation codon is used for 9 protein-coding genes, and GTG for the remaining 3 (nad1, nad4, and nad5). The mitochondrial genome of H. taichui has a single long non-coding region between trnE and trnG. H. taichui has evolved as being more closely related to Opisthorchiidae than other trematode groups with maximal support in the phylogenetic analysis. Our results could provide a resource for the comparative mitochondrial genome analysis of trematodes, and may yield genetic markers for molecular epidemiological investigations into intestinal flukes.     

The complete mitochondrial genome of Pseudochauhanea macrorchis (Monogenea: Chauhaneidae) revealed a highly repetitive region and a gene rearrangement hot spot in Polyopisthocotylea

The complete mitochondrial genome of Pseudochauhanea macrorchis was determined and compared with other monogenean mitochondrial genomes from GenBank. The circular genome was 15,031 bp in length and encoded 36 genes (12 protein-coding genes, two ribosomal RNAs, and 22 transfer RNAs) typically found in flatworms. Structures of the mitochondrial genome were mostly concordant with that known for Microcotyle sebastis and Polylabris halichoeres, but also contained two noted features-a gene rearrangement hot spot and the highly repetitive region (HRR) in major non-coding region (NCR). The gene rearrangement hot spot located between the cox3 and nad5 genes, including a cluster of tRNA genes, nad6 gene and one major NCR. The HRR seemed to be a unique feature of the polyopisthocotylean mitochondrial genomes. In conclusion, the present study provided new molecular data for future studies of the comparative mitochondrial genomics and also served as a resource of markers for the studies of species populations and monogenean phylogenetics.     

The mitochondrial genome of the thermal dimorphic fungus Penicillium marneffei is more closely related to those of molds than yeasts

We report the complete sequence of the mitochondrial genome of Penicillium marneffei, the first complete mitochondrial DNA sequence of a thermal dimorphic fungus. This 35 kb mitochondrial genome contains the genes encoding ATP synthase subunits 6, 8, and 9 (atp6, atp8, and atp9), cytochrome oxidase subunits I, II, and III (cox1, cox2, and cox3), apocytochrome b (cob), reduced nicotinamide adenine dinucleotide ubiquinone oxireductase subunits (nad1, nad2, nad3, nad4, nad4L, nad5, and nad6), ribosomal protein of the small ribosomal subunit (rps), 28 tRNAs, and small and large ribosomal RNAs. Analysis of gene contents, gene orders, and gene sequences revealed that the mitochondrial genome of P. marneffei is more closely related to those of molds than yeasts.     

Microsatellite mapping of QTL affecting growth, feed consumption, egg production, tonic immobility and body temperature of Japanese quail

Background

Animal mitochondrial genomes are physically separate from the much larger nuclear genomes and have proven useful both for phylogenetic studies and for understanding genome evolution. Within the phylum Arthropoda the subphylum Crustacea includes over 50,000 named species with immense variation in body plans and habitats, yet only 23 complete mitochondrial genomes are available from this subphylum.

Results

I describe here the complete mitochondrial genome of the crustacean Squilla mantis (Crustacea: Malacostraca: Stomatopoda). This 15994-nucleotide genome, the first described from a hoplocarid, contains the standard complement of 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a non-coding AT-rich region that is found in most other metazoans. The gene order is identical to that considered ancestral for hexapods and crustaceans. The 70% AT base composition is within the range described for other arthropods. A single unusual feature of the genome is a 230 nucleotide non-coding region between a serine transfer RNA and the nad1 gene, which has no apparent function. I also compare gene order, nucleotide composition, and codon usage of the S. mantis genome and eight other malacostracan crustaceans. A translocation of the histidine transfer RNA gene is shared by three taxa in the order Decapoda, infraorder Brachyura; Callinectes sapidus, Portunus trituberculatus and Pseudocarcinus gigas. This translocation may be diagnostic for the Brachyura. For all nine taxa nucleotide composition is biased towards AT-richness, as expected for arthropods, and is within the range reported for other arthropods. Codon usage is biased, and much of this bias is probably due to the skew in nucleotide composition towards AT-richness.

Conclusion

The mitochondrial genome of Squilla mantis contains one unusual feature, a 230 base pair non-coding region has so far not been described in any other malacostracan. Comparisons with other Malacostraca show that all nine genomes, like most other mitochondrial genomes, share a bias toward AT-richness and a related bias in codon usage. The nine malacostracans included in this analysis are not representative of the diversity of the class Malacostraca, and additional malacostracan sequences would surely reveal other unusual genomic features that could be useful in understanding mitochondrial evolution in this taxon.     

Ac-SAA-1, an immunodominant 16 kDa surface-associated antigen of infective larvae and adults of Ancylostoma caninum

A cDNA encoding a surface-associated antigen was cloned from an Ancylostoma caninum infective larvae (L(3)) cDNA library by immunoscreening with pooled human immune sera. The sera were obtained from individuals living in an Ancylostoma duodenale hookworm-endemic region of China, who had light intensity infections and high antibody titers against A. caninum L(3). Ancylostoma caninum surface-associated antigen-1 is encoded by an 843 bp mRNA with a predicted open reading frame of 162 amino acids. Recombinant Ancylostoma caninum surface-associated antigen-1 was expressed in Escherichia coli and used to prepare a specific antiserum. A Western blot with anti-Ancylostoma caninum surface-associated antigen-1 specific antiserum showed that native Ancylostoma caninum surface-associated antigen-1 protein is expressed by both L(3) and adult hookworms; RT-PCR confirmed that the mRNA is transcribed in both stages. In adult hookworms, the protein localised to the basal layer of the cuticle and hypodermis of adult worms. Serological analysis determined that recombinant Ancylostoma caninum surface-associated antigen-1 protein is recognised by 61% of human sera from a Necator americanus hookworm endemic area in China, indicating the antigen is immunodominant. Anti-Ancylostoma caninum surface-associated antigen-1 antiserum partially inhibited (46.7%) invasion of hookworm L(3) into dog skin in vitro. Together these results suggest that Ancylostoma caninum surface-associated antigen-1 offers promise as a protective vaccine antigen.     

Species-specific identification of human hookworms by PCR of the mitochondrial cytochrome oxidase I gene.

Significant differences in the life histories of the human hookworms Ancylostoma duodenale and Necator americanus necessitate their differentiation for epidemiological studies and the design of control programs. Current methods of identification require time-consuming, labor-intensive techniques. A polymerase chain reaction (PCR)-based method that enables rapid species identification is described. The mitochondrial cytochrome oxidase I genes of both species were sequenced, and species-specific primer sets were designed. The primers were used in PCR to amplify 585-bp fragments of the cytochrome oxidase gene from individual hookworm eggs, larvae, and adults. The technique was also able to identify mixed infections containing equal amounts of eggs from each species. The technique is rapid, technically simple, and sensitive and will permit the accurate identification of human hookworms in epidemiological field studies.     

Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing

Little is known about the variations of nematode mitogenomes (mtDNA). Sequencing a complete mtDNA using a PCR approach remains a challenge due to frequent genome reorganizations and low sequence similarities between divergent nematode lineages. Here, a genome skimming approach based on HiSeq sequencing (shotgun) was used to assemble de novo the first complete mtDNA sequence of a root-knot nematode (Meloidogyne graminicola). An AT-rich genome (84.3%) of 20,030 bp was obtained with a mean sequencing depth superior to 300. Thirty-six genes were identified with a semi-automated approach. A comparison with a gene map of the M. javanica mitochondrial genome indicates that the gene order is conserved within this nematode lineage. However, deep genome rearrangements were observed when comparing with other species of the superfamily Hoplolaimoidea. Repeat elements of 111 bp and 94 bp were found in a long non-coding region of 7.5 kb, as similarly reported in M. javanica and M. hapla. This study points out the power of next generation sequencing to produce complete mitochondrial genomes, even without a reference sequence, and possibly opening new avenues for species/race identification, phylogenetics and population genetics of nematodes.     

The complete mitochondrial genome of the gall-forming fly, Fergusonina taylori Nelson and Yeates (Diptera: Fergusoninidae)

The monogeneric family Fergusoninidae consists of gall-forming flies that, together with Fergusobia (Tylenchida: Neotylenchidae) nematodes, form the only known mutualistic association between insects and nematodes. In this study, the entire 16,000 bp mitochondrial genome of Fergusonina taylori Nelson and Yeates was sequenced. The circular genome contains one encoding region including 27 genes and one non-coding A+T-rich region. The arrangement of the protein-coding, ribosomal RNA (rRNA) and transfer RNA (tRNA) genes was the same as that found in the ancestral insect. Nucleotide composition is highly A+T biased. All of the protein initiation codons are ATN, except for nad1 which begins with TTT. All 22 tRNA anticodons of F. taylori match those observed in Drosophila yakuba, and all form the typical cloverleaf structure except for tRNA-Ser((AGN)) which lacks a dihydrouridine (DHU) arm. Secondary structural features of the rRNA genes of Fergusonina are similar to those proposed for other insects, with minor modifications. The mitochondrial genome of Fergusonina presented here may prove valuable for resolving the sister group to the Fergusoninidae, and expands the available mtDNA data sources for acalyptrates overall.     

The complete mitochondrial genome of Anisakis simplex (Ascaridida: Nematoda) and phylogenetic implications

We determined the nucleotide sequence of the complete mitochondrial genome of the nematode species Anisakis simplex. The genome is circular, 13,916 bp in size and conforms to the general characteristics of nematode mitochondrial DNAs. The gene arrangement of A. simplex is the same as that of Ascaris suum and almost identical to those of rhabditid species with a minor exception concerning the relative position of the AT-rich and non-coding regions and radically different from those of spirurid species. Along with comparisons of gene arrangement, phylogenetic analyses (maximum parsimony, neighbour joining and maximum likelihood methods) based on concatenated amino acid sequences of 12 protein-coding genes from 13 nematode species provided strong support for the sister-group relationship between Ascaridida and Rhabditida. The Shimodaira-Hasegawa and Templeton's tests both rejected the alternative hypothesis of a closer relationship between Ascaridida and Spirurida. These results contradicted the traditional view of nematode classification and a recent molecular phylogenetic study of 18S rDNA data that assigned Ascaridida and Spirurida as being a sister-group. Mapping of gene arrangement across the phylogenetic tree lead to the assumption that the conserved gene arrangement found in Ascaridida-Rhabditida members might have been acquired after the most recent common ancestor of ascaridid/rhabditid members branched off from the basal stock of the rhabditid lineage.     

Genes and genomes of Necator americanus and related hookworms

The human hookworms (Necator americanus and Ancylostoma duodenale) infect over one billion people. The phylogenetic relationships of the human hookworms suggest independent acquisition of the human host. The hookworms probably have a haploid chromosome number n = 6, and an XX-XO sex determination mechanism is likely to be used. Genetic and molecular research on hookworms is in its infancy, but several important genes and gene products have already been identified. Of note are cathepsin genes, a family of secreted proteins known as Ancylostoma activation-associated proteins and a family of anticoagulants. The inception of an expressed sequence tag program on the human hookworm, N. americanus, promises to yield many new genes with novel functions in the biology of these important parasites.     

Analysis of three leafminers' complete mitochondrial genomes

Liriomyza trifolii (Burgess), Liriomyza huidobrensis (Blanchard), and Liriomyza bryoniae (Kaltenbach), are three closely related and economically important leafminer pests in the world. This study examined the complete mitochondrial genomes of L. trifolii, L. huidobrensis and L. bryoniae, which were 16141 bp, 16236 bp and 16183 bp in length, respectively. All of them displayed 37 typical animal mitochondrial genes and an A + T-rich region. The genomes were highly compact with only 60–68 bp of non-coding intergenic spacer. However, considerable differences in the A + T-rich region were detected among the three species. Results of this study also showed the two ribosomal RNA genes of the three species had very limited variable sites and thus should not provide much information in the study of population genetics of these species. Data generated from three leafminers' complete mitochondrial genomes should provide valuable information in studying phylogeny of Diptera, and developing genetic markers for species identification in leafminers.     

The complete mitochondrial genome and phylogenetic analysis of Nyctereutes procyonoides

The complete mitochondrial genome sequence of the raccoon dog (Nyctereutes procyonoides) was determined by using the long and accurate polymerase chain reaction. The entire mitochondrial genome sequence is 16,713 bp in length contains two ribosomal RNA genes, 13 protein-coding genes, 22 transfer RNA genes and 1 control region. Most mitochondrial genes are encoded on the H strand, except for the ND6 gene and 8 tRNA genes. The base compositions of mitochondrial genomes present clearly A–T skew. All the transfer RNA genes can be folded into the typical cloverleaf-shaped structure except tRNA-Ser (AGY), which lacks the dihydrouridine arm. Protein-coding genes mainly initiate with ATG and terminate with TAA. Some reading frame intervals and overlaps are found in the mitochondrial genome. The control region can be divided into three domains: the extended termination associated sequences (ETASs) domain, the central conserved domain and the conserved sequence blocks (CSBs) domain. Three conserved sequence blocks (CSBs) and one extended termination associated sequences (ETAS-1) is found in the control region. The phylogenetic analysis based on the concatenated data set of 14 genes in the mitochondrial genome of Canidae shows that the raccoon dog has close phylogenetic position with the red fox (Vulpes vulpes) and they constitute a clade which has an equil evolutionary position with the clade formed by the genera Canis and Cuon.     

The mitochondrial genome of Strongyloides stercoralis (Nematoda) - idiosyncratic gene order and evolutionary implications

The complete mitochondrial genome sequence of the parasitic nematode Strongyloides stercoralis was determined, and its organisation and structure compared with other nematodes for which complete mitochondrial sequence data were available. The mitochondrial genome of S. stercoralis is 13,758 bp in size and contains 36 genes (all transcribed in the clockwise direction) but lacks the atp8 gene. This genome has a high T content (55.9%) and a low C content (8.3%). Corresponding to this T content, there are 16 (poly-T) tracts of >/=12 Ts distributed across the genome. In protein-coding genes, the T bias is greatest (76.4%) at the third codon position compared with the first and second codon positions. Also, the C content is higher at the first (9.3%) and second (13.4%) codon positions than at the third (2%) position. These nucleotide biases have a significant effect on predicted codon usage patterns and, hence, on amino acid compositions of the mitochondrial proteins. Interestingly, six of the 12 protein-coding genes are predicted to employ a unique initiation codon (TTT), which has not yet been reported for any other animal mitochondrial genome. The secondary structures predicted for the 22 transfer RNA (trn) genes and the two ribosomal RNA (rrn) genes are similar to those of other nematodes. In contrast, the gene arrangement in the mitochondrial genome of S. stercoralis is different from all other nematodes studied to date, revealing only a limited number of shared gene boundaries (atp6-nad2 and cox2-rrnL). Evolutionary analyses of mitochondrial nucleotide and amino acid sequence data sets for S. stercoralis and seven other nematodes demonstrate that the mitochondrial genome provides a rich source of phylogenetically informative characters. In conclusion, the S. stercoralis mitochondrial genome, with its unique gene order and characteristics, should provide a resource for comparative mitochondrial genomics and systematics studies of parasitic nematodes.     

The complete mitochondrial genome of the black coral Chrysopathes formosa (Cnidaria:Anthozoa:Antipatharia) supports classification of antipatharians within the subclass Hexacorallia

Black corals comprise a globally distributed shallow- and deep-water taxon whose phylogenetic position within the Anthozoa has been debated. We sequenced the complete mitochondrial genome of the antipatharian Chrysopathes formosa to further evaluate its phylogenetic relationships. The circular mitochondrial genome (18,398 bp) consists of 13 energy pathway protein-coding genes and two ribosomal RNAs, but only two transfer RNA genes (trnM and trnW), as well as a group I intron within the nad5 gene that contains the only copies of nad1 and nad3. No novel genes were found in the antipatharian mitochondrial genome. Gene order and genome content are most similar to those of the sea anemone Metridium senile (subclass Hexacorallia), with differences being the relative location of three contiguous genes (cox2-nad4-nad6) and absence (from the antipatharian) of a group I intron within the cox1 gene. Phylogenetic analyses of multiple protein-coding genes support classifying the Antipatharia within the subclass Hexacorallia and not the subclass Ceriantipatharia; however, the sister-taxon relationships of black corals within Hexacorallia remain inconclusive.     

Complete sequence of the mitochondrial genome of Taenia saginata: comparison with T. solium and T. asiatica

The complete sequence of the Taenia saginata mitochondrial genome was determined, and its organization and structure were compared to other human-tropic Taenia tapeworms for which complete mitochondrial sequence data were available. The mitochondrial genome was 13,670 bp long, contained 12 protein-coding genes, two ribosomal RNAs (rRNAs, a small and a large subunit), and 22 transfer RNAs (tRNAs). It did not encode the atp8 gene. Overlapping regions were found between nad4L and nad4, nad1 and trnN, and cox1 and trnT. The ATG initiation codon was used for 10 protein-coding genes, and the GTG initiation codon was used for the remaining 2 genes (nad4 and atp6). The size of the protein-coding genes of the three human Taenia tapeworms did not vary, except for Taenia solium nad1 (891 aa) and nad4 (1212 aa) and Taenia asiatica cox2 (576 aa). The tRNA genes were 57-75 bp long, and the predicted secondary structures of 18 of these genes had typical clover-leaf shapes with paired dihydrouridine (DHU) arms. The genes in all human Taenia tapeworms for the two mitochondrial rRNA subunits rrnL and rrnS are separated by trnC. The putative T. saginata rrnL and rrnS are 972 and 732 bp long, respectively. The non-coding regions of the mt genome of T. saginata consisted of 2 regions: a short non-coding region (SNR, 66 nucleotides) and a long non-coding region (LNR, 159 nucleotides). The overall sequence difference in the full mitochondrial genome between T. saginata and T. asiatica was 4.6%, while T. solium differed by 11%. In conclusion, the complete sequence of the T. saginata mitochondrial genome will serve as a resource for comparative mitochondrial genomics and systematic studies of the parasitic cestodes.     

New features of Asian Crassostrea oyster mitochondrial genomes: A novel alloacceptor tRNA gene recruitment and two novel ORFs

A feasible way to perform evolutionary analyses is to compare characters divergent enough to observe significant differences, but sufficiently similar to exclude saturation of the differences that occurred. Thus, comparisons of invertebrate mitochondrial (mt) genomes at low taxonomic levels can be extremely helpful in investigating patterns of variation and evolutionary dynamics of genomes, as intermediate stages of the process may be identified. Fortunately, in this study, we newly sequenced the mt genome of the eighth member of Asian Crassostrea oysters which can provide necessary intermediate characters for us to believe that the variation of Crassostrea mt genomes is considerably greater than previously acknowledged. Several new features of Asian Crassostrea oyster mitochondrial genomes were revealed, and our results are particularly significant as they 1) suggest a novel model of alloacceptor tRNA gene recruitment, namely "vertical" tRNA gene recruitment, which can be successfully used to explain the origination of the unusually additional trnK and trnQ genes (annotated as trnK(2) and trnQ(2) respectively) in the mt genomes of the five Asian oysters, and we speculate that this recruitment progress may be a common phenomenon in the evolution of the tRNA multigene family; 2) reveal the existence of two additional, lineage-specific, mtDNA-encoded genes that may originate from duplication of nad2 followed by rapid evolutionary change. Each of these two genes encodes a unique amino terminal signal peptide, thus each might possess an unknown function; and 3) identify for the first time the atp8 gene in oysters. The present study thus gives further credence to the comparison of congeneric bivalves as a meaningful strategy to investigate mt genomic evolutionary trends in genome organization, tRNA multigene family, and gene loss and/or duplication that are difficult to undertake at higher taxonomic levels. In particular, our study provides new evidence for the identification and characterization of ORFs in the "non-coding region" of animal mt genomes.     

Complete Mitochondrial Genomes of Carcinoscorpius rotundicauda and Tachypleus tridentatus (Xiphosura,Arthropoda) and Implications for Chelicerate Phylogenetic Studies

Horseshoe crabs (order Xiphosura) are often referred to as an ancient order of marine chelicerates and have been considered as keystone taxa for the understanding of chelicerate evolution. However, the mitochondrial genome of this order is only available from a single species, Limulus polyphemus. In the present study, we analyzed the complete mitochondrial genomes from two Asian horseshoe crabs, Carcinoscorpius rotundicauda and Tachypleus tridentatus to offer novel data for the evolutionary relationship within Xiphosura and their position in the chelicerate phylogeny. The mitochondrial genomes of C. rotundicauda (15,033 bp) and T. tridentatus (15,006 bp) encode 13 protein-coding genes, two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. Overall sequences and genome structure of two Asian species were highly similar to that of Limulus polyphemus, though clear differences among three were found in the stem-loop structure of the putative control region. In the phylogenetic analysis with complete mitochondrial genomes of 43 chelicerate species, C. rotundicauda and T. tridentatus were recovered as a monophyly, while L. polyphemus solely formed an independent clade. Xiphosuran species were placed at the basal root of the tree, and major other chelicerate taxa were clustered in a single monophyly, clearly confirming that horseshoe crabs composed an ancestral taxon among chelicerates. By contrast, the phylogenetic tree without the information of Asian horseshoe crabs did not support monophyletic clustering of other chelicerates. In conclusion, our analyses may provide more robust and reliable perspective on the study of evolutionary history for chelicerates than earlier analyses with a single Atlantic species.