Horizontal gene transfer from a flowering plant to the insular pine Pinus canariensis (Chr. Sm. Ex DC in Buch)

Horizontal gene transfer (HGT) is viewed as very common in the plant mitochondrial (mt) genome, but, to date, only one case of HGT has been found in gymnosperms. Here we report a new case of HGT, in which a mt nad5-1 fragment was transferred from an angiosperm to Pinus canariensis. Quantitative assay and sequence analyses showed that the foreign nad5-1 is located in the mt genome of P. canariensis and is nonfunctional. An extensive survey in the genus Pinus revealed that the angiosperm-derived nad5-1 is restricted to P. canariensis and present across the species'' range. Molecular dating based on chloroplast DNA suggested that the HGT event occurred in the late Miocene after P. canariensis split from its closest relatives, and that the foreign copy became fixed in P. canariensis owing to drift during its colonization of the Canary Islands. The mechanism of this HGT is unclear but it was probably achieved through either direct cell–cell contact or external vectors. Our discovery provides evidence for an important role of HGT in plant mt genome evolution.     

The First Mitochondrial Genome for Caddisfly (Insecta: Trichoptera) with Phylogenetic Implications

The Trichoptera (caddisflies) is a holometabolous insect order with 14,300 described species forming the second most species-rich monophyletic group of animals in freshwater. Hitherto, there is no mitochondrial genome reported of this order. Herein, we describe the complete mitochondrial (mt) genome of a caddisfly species, Eubasilissa regina (McLachlan, 1871). A phylogenomic analysis was carried out based on the mt genomic sequences of 13 mt protein coding genes (PCGs) and two rRNA genes of 24 species belonging to eight holometabolous orders. Both maximum likelihood and Bayesian inference analyses highly support the sister relationship between Trichoptera and Lepidoptera.     

The mitochondrial genome of Baylisascaris procyonis

Background

Baylisascaris procyonis (Nematoda: Ascaridida), an intestinal nematode of raccoons, is emerging as an important helminthic zoonosis due to serious or fatal larval migrans in animals and humans. Despite its significant veterinary and public health impact, the epidemiology, molecular ecology and population genetics of this parasite remain largely unexplored. Mitochondrial (mt) genomes can provide a foundation for investigations in these areas and assist in the diagnosis and control of B. procyonis. In this study, the first complete mt genome sequence of B. procyonis was determined using a polymerase chain reaction (PCR)-based primer-walking strategy.

Methodology/Principal Findings

The circular mt genome (14781 bp) of B. procyonis contained 12 protein-coding, 22 transfer RNA and 2 ribosomal RNA genes congruent with other chromadorean nematodes. Interestingly, the B. procyonis mtDNA featured an extremely long AT-rich region (1375 bp) and a high number of intergenic spacers (17), making it unique compared with other secernentean nematodes characterized to date. Additionally, the entire genome displayed notable levels of AT skew and GC skew. Based on pairwise comparisons and sliding window analysis of mt genes among the available 11 Ascaridida mtDNAs, new primer pairs were designed to amplify specific short fragments of the genes cytb (548 bp fragment) and rrnL (200 bp fragment) in the B. procyonis mtDNA, and tested as possible alternatives to existing mt molecular beacons for Ascaridida. Finally, phylogenetic analysis of mtDNAs provided novel estimates of the interrelationships of Baylisasaris and Ascaridida.

Conclusions/Significance

The complete mt genome sequence of B. procyonis sequenced here should contribute to molecular diagnostic methods, epidemiological investigations and ecological studies of B. procyonis and other related ascaridoids. The information will be important in refining the phylogenetic relationships within the order Ascaridida and enriching the resource of markers for systematic, population genetic and evolutionary biological studies of parasitic nematodes of socio-economic importance.     

Assembly and annotation of full mitochondrial genomes for the corn rootworm species, Diabrotica virgifera virgifera and Diabroticabarberi (Insecta: Coleoptera: Chrysomelidae), using Next Generation Sequence data

Complete mitochondrial genomes for two corn rootworm species, Diabrotica virgifera virgifera (16,747 bp) and Diabrotica barberi (16,632; Insecta: Coleoptera: Chrysomelidae), were assembled from Illumina HiSeq2000 read data. Annotation indicated that the order and orientation of 13 protein coding genes (PCGs), and 22 tRNA and 2 rRNA sequences were in typical of insect mitochondrial genomes. Non-standard nad4 and cox3 stop codons were composed of single T nucleotides and likely completed by adenylation, and atypical TTT start codons was predicted for both D. v. virgifera and D. barberinad1 genes. The D. v. virgifera and D. barberi haplotypes showed 819 variable nucleotide positions within PCG regions (7.36% divergence), which suggest that speciation may have occurred ~ 3.68 million years ago assuming a linear rate of short-term substitution. Phylogenetic analyses of Coleopteran MtD genome show clustering based on family level, and may have the capacity to resolve the evolutionary history within this Order of insects.     

Mitochondrial and Nuclear Ribosomal DNA Evidence Supports the Existence of a New Trichuris Species in the Endangered Fran?ois’ Leaf-Monkey

The whipworm of humans, Trichuris trichiura, is responsible for a neglected tropical disease (NTD) of major importance in tropical and subtropical countries of the world. Whipworms also infect animal hosts, including pigs, dogs and non-human primates, cause clinical disease (trichuriasis) similar to that of humans. Although Trichuris species are usually considered to be host specific, it is not clear whether non-human primates are infected with T. trichiura or other species. In the present study, we sequenced the complete mitochondrial (mt) genome as well as the first and second internal transcribed spacers (ITS-1 and ITS-2) of Trichuris from the François’ leaf-monkey (langur), and compared them with homologous sequences from human- and pig-derived Trichuris. In addition, sequence comparison of a conserved mt ribosomal gene among multiple individual whipworms revealed substantial nucleotide differences among these three host species but limited sequence variation within each of them. The molecular data indicate that the monkey-derived whipworm is a separate species from that of humans. Future work should focus on detailed population genetic and morphological studies (by electron microscopy) of whipworms from various non-humans primates and humans.     

Multiple primer PCR for the identification of anisakid nematodes from Taiwan Strait

There were six major larval anisakid species found in commercial marine fishes caught in the Minnan fishing ground in the Taiwan Strait: Anisakis physeteris, Anisakis pegreffii, Raphidascaris trichiuri, Contracaecum aduncum, Contracaecum muraenesoxi, Contracaecum sp. For rapid identification of the parasite species above, a single and a multiple primer PCR (multiplex PCR) method, using specific primers based on aligned sequences of the internal transcribed spacer ITS-1, 5.8S, and ITS-2 of nuclear ribosomal DNA, were jointly used for the rapid identification of these anisakid larvae. The primers yielded distinct PCR products for each of the anisakid nematodes, providing rapid and accurate tools for identifying anisakid nematodes with distinct geographical distribution.     

Genotyping Echinococcus granulosus from dogs from Western Iran

Cystic echinococcosis is a zoonotic infection caused by the dog tapeworm, Echinococcus granulosus. In the present study, adults of E. granulosus (n = 20) were collected from 71 dogs from Western Iran and were genetically characterized using DNA sequencing of the partial mitochondrial cytochrome c oxidase subunit 1 (cox1) and NADH dehydrogenase 1 (nad1). Consensus sequences were obtained for cox1 (366) and nad1 (471) genes. Phylogenetic analysis of concatenated nad1 and cox1 nucleotide sequence data was performed using Bayesian Inference approach. Overall, the dog isolates indicated nine different sequences in cox1 and seven in nad1 genes. Three genotypes (G1 [75%], G2 [10%] and G3 [15%]) were identified from the isolates. The G2 sequences indicated 100% homology with reference G2 sequence in both cox1 (Genbank accession number M84662) and nad1 (AJ237633) genes. G3 sequences showed 100% homology with G3 reference sequence in nad1 (AJ237633), but displayed two different cox1 profiles, each having 99% homology with reference G3 sequence (M84663). In the phylogenetic tree all of the isolates were grouped into a distinct cluster corresponding to the G1–G3 complex with relevant reference sequences. The presence of G1 genotype (sheep strain) of E. granulosus sensu stricto as dominant genotype in dogs is emphasized. To the best of our knowledge, this study established the first record of E. granulosus sensu stricto, G2 genotype in Iran.     

Transmission of lungworms of harbour porpoises and harbour seals: Molecular tools determine potential vertebrate intermediate hosts

Harbour porpoises (Phocoena phocoena) and harbour seals (Phoca vitulina) from German waters are infected by six species of lungworms (Metastrongyloidea). These nematodes parasitise the respiratory tract, are pathogenic and often cause secondary bacterial infections. In spite of their clinical and epidemiological significance, the life cycle and biology of lungworms in the marine environment is still largely unknown. Regions of ribosomal DNA (ITS-2) of all lungworms parasitising harbour porpoises and harbour seals in German waters were sequenced to characterise and compare the different species. The phylogenetic relationship among the lungworm species was analysed by means of their ITS-2 nucleotide sequences and the species-specific traits of the ITS-2 were used to screen wild fish as possible intermediate hosts for larval lungworms. Molecular markers were developed to identify larval nematodes via in-situ hybridisation of tissues of harbour porpoise and harbour seal prey fish. Potential wild intermediate fish hosts from the North Sea were dissected and found to harbour larval nematodes. Histological examination and in-situ hybridisation of tissue samples from these fish showed lungworm larvae within the intestinal wall. Based on larval ITS-2 nucleotide sequences, larval nematodes were identified as Pseudalius inflexus and Parafilaroides gymnurus. Turbot (Psetta maxima) bred and raised in captivity were experimentally infected with live L1s of Otostrongylus circumlitus and ensheathed larvae were recovered from the gastrointestinal tract of turbot and identified using molecular tools. Our results show that fish intermediate hosts play a role in the transmission of metastrongyloid nematodes of harbour porpoises and harbour seals.     

Mitochondrial gene sequences alone or combined with ITS region sequences provide firm molecular criteria for the classification of Lecanicillium species

The recent revision of Verticillium sect. Prostrata led to the introduction of the genus Lecanicillium, which comprises the majority of the entomopathogenic strains. Sixty-five strains previously classified as Verticillium lecanii or Verticillium sp. from different geographical regions and hosts were examined and their phylogenetic relationships were determined using sequences from three mitochondrial (mt) genes [the small rRNA subunit (rns), the NADH dehydrogenase subunits 1 (nad1) and 3 (nad3)] and the ITS region. In general, single gene phylogenetic trees differentiated and placed the strains examined in well-supported (by BS analysis) groups of L. lecanii, L. longisporum, L. muscarium, and L. nodulosum, although in some cases a few uncertainties still remained. nad1 was the most informative single gene in phylogenetic analyses and was also found to contain group I introns with putative open reading frames (ORFs) encoding for GIY–YIG endonucleases. The combined use of mt gene sequences resolved taxonomic uncertainties arisen from ITS analysis and, alone or in combination with ITS sequences, helped in placing uncharacterised Verticillium lecanii and Verticillium sp. firmly into Lecanicillium species. Combined gene data from all the mt genes and all the mt genes and the ITS region together, were very similar. Furthermore, a relaxed correlation with host specificity—at least for Homoptera—was indicated for the rns and the combined mt gene sequences. Thus, the usefulness of mt gene sequences as a convenient molecular tool in phylogenetic studies of entomopathogenic fungi was demonstrated.     

Complete Mitochondrial Genome of Echinostoma hortense (Digenea: Echinostomatidae)

Echinostoma hortense (Digenea: Echinostomatidae) is one of the intestinal flukes with medical importance in humans. However, the mitochondrial (mt) genome of this fluke has not been known yet. The present study has determined the complete mt genome sequences of E. hortense and assessed the phylogenetic relationships with other digenean species for which the complete mt genome sequences are available in GenBank using concatenated amino acid sequences inferred from 12 protein-coding genes. The mt genome of E. hortense contained 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 non-coding region. The length of the mt genome of E. hortense was 14,994 bp, which was somewhat smaller than those of other trematode species. Phylogenetic analyses based on concatenated nucleotide sequence datasets for all 12 protein-coding genes using maximum parsimony (MP) method showed that E. hortense and Hypoderaeum conoideum gathered together, and they were closer to each other than to Fasciolidae and other echinostomatid trematodes. The availability of the complete mt genome sequences of E. hortense provides important genetic markers for diagnostics, population genetics, and evolutionary studies of digeneans.     

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

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

The Complete Mitochondrial Genome of Gossypium hirsutum and Evolutionary Analysis of Higher Plant Mitochondrial Genomes

Background

Mitochondria are the main manufacturers of cellular ATP in eukaryotes. The plant mitochondrial genome contains large number of foreign DNA and repeated sequences undergone frequently intramolecular recombination. Upland Cotton (Gossypium hirsutum L.) is one of the main natural fiber crops and also an important oil-producing plant in the world. Sequencing of the cotton mitochondrial (mt) genome could be helpful for the evolution research of plant mt genomes.

Methodology/Principal Findings

We utilized 454 technology for sequencing and combined with Fosmid library of the Gossypium hirsutum mt genome screening and positive clones sequencing and conducted a series of evolutionary analysis on Cycas taitungensis and 24 angiosperms mt genomes. After data assembling and contigs joining, the complete mitochondrial genome sequence of G. hirsutum was obtained. The completed G.hirsutum mt genome is 621,884 bp in length, and contained 68 genes, including 35 protein genes, four rRNA genes and 29 tRNA genes. Five gene clusters are found conserved in all plant mt genomes; one and four clusters are specifically conserved in monocots and dicots, respectively. Homologous sequences are distributed along the plant mt genomes and species closely related share the most homologous sequences. For species that have both mt and chloroplast genome sequences available, we checked the location of cp-like migration and found several fragments closely linked with mitochondrial genes.

Conclusion

The G. hirsutum mt genome possesses most of the common characters of higher plant mt genomes. The existence of syntenic gene clusters, as well as the conservation of some intergenic sequences and genic content among the plant mt genomes suggest that evolution of mt genomes is consistent with plant taxonomy but independent among different species.     

New record of Ascaridia nymphii (Secernentea: Ascaridiidae) from macaw parrot,Ara chloroptera,in China

Present study was performed to identify the species of ascarids from macaw parrot, Ara chloroptera, in China. Total 6 ascarids (3 males and 3 females) were collected in the feces of 3 macaws at Guangzhou Zoo in Guangdong Province, China. Their morphological characteristics with dimensions were observed under a light microscope, and their genetic characters were analyzed with the partial 18S rDNA, ITS rDNA and nad4 gene sequences, respectively. Results showed that all worms have no interlabia but male worms have two alate spicules, well-developed precloacal sucker and a tail with ventrolateral caudal alae and 11 pairs of papillae. The partial 18S rDNA, ITS rDNA and nad4 sequences were 831 bp, 1015 bp and 394 bp in length, respectively. They showed the highest similarity of 99.8% (18S rDNA) with Ascaridia nymphii, 93.8% identities (ITS rDNA) with A. columbae and 98.5% to 99.5% identities (nad4) with Ascaridia sp. from infected parrot. All Ascaridia nematodes from the macaws were clustered into one clade and formed monophyletic group of Ascaridia with A. columbae and A. galli in two phylogenetic trees. It is observed that the combining morphological and sequencing data from three loci, the present Ascaridia species was identified as Ascaridia nymphii, which is the first record of A. nymphii from macaw parrot in China.     

The First Mitochondrial Genome for the Fishfly Subfamily Chauliodinae and Implications for the Higher Phylogeny of Megaloptera

Megaloptera are a basal holometabolous insect order with larvae exclusively predacious and aquatic. The evolutionary history of Megaloptera attracts great interest because of its antiquity and important systematic status in Holometabola. However, due to the difficulties identifying morphological apomorphies for the group, controversial hypotheses on the monophyly and higher phylogeny of Megaloptera have been proposed. Herein, we describe the complete mitochondrial (mt) genome of a fishfly species, Neochauliodes punctatolosus Liu & Yang, 2006, representing the first mt genome of the subfamily Chauliodinae. A phylogenomic analysis was carried out based on the mt genomic sequences of 13 mt protein-coding genes (PCGs) and two rRNA genes of nine Neuropterida species, comprising all three orders of Neuropterida and all families and subfamilies of Megaloptera. Both maximum likelihood and Bayesian inference analyses highly support the monophyly of Megaloptera, which was recovered as the sister of Neuroptera. Within Megaloptera, the sister relationship between Corydalinae and Chauliodinae was corroborated. The divergence time estimation suggests that stem lineage of Neuropterida and Coleoptera separated in the Early Permian. The interordinal divergence within Neuropterida might have occurred in the Late Permian.     

The complete mitochondrial genome of Macrobrachium nipponense

The complete mitochondrial (mt) genome sequence plays an important role in the accurate determination of phylogenetic relationships among metazoans. Herein, we determined the complete mt genome sequence, structure and organization of Macrobrachium nipponense (M. nipponense) (GenBank ID: NC_015073.1) and compared it to that of Macrobrachium lanchesteri (M. lanchesteri) and Macrobrachium rosenbergii (M. rosenbergii). The 15,806 base pair (bp) M. nipponense mt genome, which is comprised of 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and 2 ribosomal RNAs (rRNAs), is slightly larger than that of M. lanchesteri (15,694 bp, GenBank ID: NC_012217.1) and M. rosenbergii (15,772 bp, GenBank ID: NC_006880.1). The M. nipponense genome contains a high AT content (66.0%), which is a common feature among metazoan mt genomes. Compared with M. lanchesteri and M. rosenbergii, we found a peculiar non-coding region of 950 bp with a microsatellite-like (TA)₆ element and many hairpin structures. The 13 PCGs are comprised of a total of 3707 codons, excluding incomplete termination codons, and the most frequently used amino acid is Leu (16.0%). The predicted start codons in the M. nipponense mt genome include ATG, ATC and ATA. Seven PCGs use TAA as a stop codon, whereas two use TAG, three use T and only one uses TA. Twenty-three of the genes are encoded on the L strand, and ND1, ND4, ND5, ND4L, 12S rRNA, 16S rRNA, tRNA^His, tRNA^Pro, tRNA^Phe, tRNA^Val, tRNA^Gln, tRNA^Cys, tRNA^Tyr and a tRNA^Leu are encoded on the H strand. The two rRNAs of M. nipponense and M. rosenbergii are encoded on the H strand, whereas the M. lanchesteri rRNAs are encoded on the L stand.     

The mitochondrial genome sequence analysis of Ophidascaris baylisi from the Burmese python (Python molurus bivittatus)

Ophidascaris species are parasitic roundworms that inhabit the python gut, resulting in severe granulomatous lesions or even death. However, the classification and nomenclature of these roundworms are still controversial. Our study aims to identify a snake roundworm from the Burmese python (Python molurus bivittatus) and analyze the mitochondrial genome. We identified this roundworm as Ophidascaris baylisi based on the morphology and cytochrome c oxidase subunit I (cox1) sequence. Ophidascaris baylisi complete mitochondrial genome was 14,784 bp in length, consisting of two non-coding regions and 36 mitochondrial genes (12 protein-coding genes, 22 tRNA genes, and two rRNA genes). The protein-coding genes used TTG, ATG, ATT, or TTA as start codons and TAG, TAA, or T as stop codons. All tRNA genes showed a TV-loop structure, except trnS1^AGN and trnS2^UCN revealed a D-loop structure. The mitochondrial large ribosomal subunit 16S (rrnL) and small ribosomal subunit 12S (rrnS) were 956 bp and 700 bp long, respectively. Phylogenetic analysis based on O. baylisi mitochondrial protein-coding genes demonstrated that O. baylisi clustered with the family Ascarididae members and was most closely related to Ophidascaris wangi. These results may enhance the nematode mitochondrial genome database and provide valuable molecular markers for further research on the taxonomy, phylogeny, and genetic relationships of Ophidascaris nematodes.     

Mitochondrial genome sequence diversity of Indian Plasmodium falciparum isolates

We have analysed the whole mitochondrial (mt) genome sequences (each ~6 kilo nucleotide base pairs in length) of four field isolates of the malaria parasite Plasmodium falciparum collected from different locations in India. Comparative genomic analyses of mt genome sequences revealed three novel India-specific single nucleotide polymorphisms. In general, high mt genome diversity was found in Indian P. falciparum, at a level comparable to African isolates. A population phylogenetic tree placed the presently sequenced Indian P. falciparum with the global isolates, while a previously sequenced Indian isolate was an outlier. Although this preliminary study is limited to a few numbers of isolates, the data have provided fundamental evidence of the mt genome diversity and evolutionary relationships of Indian P. falciparum with that of global isolates.     

The complete nucleotide sequence of the mitochondrial genome of <Emphasis Type="Italic">Phthonandria atrilineata</Emphasis> (Lepidoptera: Geometridae)

Using long-polymerase chain reaction (Long-PCR) method, we determined the complete nucleotide sequence of the mitochondrial genome (mitogenome) of Phthonandria atrilineata. The complete mtDNA from P. atrilineata was 15,499 base pairs in length and contained 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The P. atrilineata genes were in the same order and orientation as the completely sequenced mitogenomes of other lepidopteran species. The nucleotide composition of P. atrilineata mitogenome was biased toward A + T nucleotides (81.02%), and the 13 PCGs show different A + T contents that range from 73.25% (cox1) to 92.12% (atp8). Phthonandria had the canonical set of 22 tRNA genes, that fold in the typical cloverleaf structure described for metazoan mt tRNAs, with the unique exception of trnS(AGN). The phylogenetic relationships were reconstructed with the concatenated sequences of the 13 PCGs of the mitochondrial genome, which confirmed that P. atrilineata is most closely related to the superfamily Bombycoidea.