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

Diversification of the terrestrial frog genus Anomaloglossus (Anura, Aromobatidae) in the Guiana Shield proceeded from highlands to lowlands, with successive loss and reacquisition of endotrophy

Two main landscapes emerge from the Guiana Shield: the highlands to the west called the Pantepui region and the Amazonian lowlands to the east, both harbouring numerous endemic species. With 32 currently recognized species, the genus Anomaloglossus stands out among Neotropical frogs as one that diversified only within the Guiana Shield both in the highlands and the lowlands. We present a time-calibrated phylogeny obtained by using combined mitogenomic and nuclear DNA, which suggests that the genus originates from Pantepui where extant lineages started diversifying around 21 Ma, and subsequently (ca. 17 Ma) dispersed during the Miocene Climatic Optimum to the lowlands of the eastern Guiana Shield where the ability to produce endotrophic tadpoles evolved. Further diversification within the lowlands in the A. stepheni group notably led to an evolutionary reversal toward exotrophy in one species group during the late Miocene, followed by reacquisition of endotrophy during the Pleistocene. These successive shifts of reproductive mode seem to have accompanied climatic oscillations. Long dry periods might have triggered evolution of exotrophy, whereas wetter climates favoured endotrophic forms, enabling colonization of terrestrial habitats distant from water. Acquisition, loss, and reacquisition of endotrophy makes Anomaloglossus unique among frogs and may largely explain the current species diversity. The micro evolutionary processes involved in these rapid shifts of reproductive mode remain to be revealed.

An Integrated Approach for Synonymization of Rotylenchus rhomboides with R. goodeyi (Nematoda: Hoplolaimidae) Reveals High Intraspecific Mitogenomic Variation

Rotylenchus is a widely distributed, economically important plant-parasitic nematode group whose species-level identification relies largely on limited morphological characters, including character-based tabular keys and molecular data of ribosomal and mitochondrial genes. In this study, a combined morphological and molecular analysis of three populations of Rotylenchus goodeyi from Belgium, Poland, and the Netherlands revealed important character variations of this species, leading to synonymization of R. rhomboides with R. goodeyi and a high nucleotide variation within cox1 gene sequences in these populations. Additional Illumina sequencing of DNA from individuals of the Dutch population revealed two variants of mitogenomes, each approximately 23 Kb in size, differing by approximately 9% and containing 11 protein-coding genes, 2 ribosomal RNA genes, and as many as 29 transfer RNA genes. In addition to the first representative whole-genome shotgun sequence datasets of the genus Rotylenchus, this study also provides the full-length mitogenome and the ribosomal DNA sequences of R. goodeyi.

Meloidogyne graminicola-A Threat to Rice Production: Review Update on Distribution, Biology, Identification, and Management

Simple Summary

New risks to plant health are constantly emerging. Such is the case of the rice root knot nematode Meloidogyne graminicola, adapted to flooded conditions and representing a risk to all types of rice agro-systems. It has been recently detected in Italy and added to the European and Mediterranean Plant Protection Organization (EPPO) Alert List. The presence of this nematode in Europe poses a threat to rice production, as there is a high probability to spread, due to trade activities and climate changes. In view of its importance, an extensive updated review was carried out.

Abstract

Rice (Oryza sativa L.) is one of the main cultivated crops worldwide and represents a staple food for more than half of the world population. Root-knot nematodes (RKNs), Meloidogyne spp., and particularly M. graminicola, are serious pests of rice, being, probably, the most economically important plant-parasitic nematode in this crop. M. graminicola is an obligate sedentary endoparasite adapted to flooded conditions. Until recently, M. graminicola was present mainly in irrigated rice fields in Asia, parts of the Americas, and South Africa. However, in July 2016, it was found in northern Italy in the Piedmont region and in May 2018 in the Lombardy region in the province of Pavia. Following the first detection in the EPPO region, this pest was included in the EPPO Alert List as its wide host range and ability to survive during long periods in environments with low oxygen content, represent a threat for rice production in the European Union. Considering the impact of this nematode on agriculture, a literature review focusing on M. graminicola distribution, biology, identification, and management was conducted.

Genome sequence of the coffee root-knot nematode Meloidogyne exigua

Root-knot nematodes (Meloidogyne spp.) cause serious damages on most crops. Here, we report a high-quality genome sequence of Meloidogyne exigua (population Mex1, Costa Rica), a major pathogen of coffee. Its mitogenome (20,974 bp) was first assembled and annotated. The nuclear genome was then constructed consisting of 206 contigs, with an N50 length of 1.89 Mb and a total assembly length of 42.1 Mb.

Adaptive response to olive cultivation in a generalist parasitic nematode (Meloidogyne javanica)

Cultivated plants usually differ from their wild progenitors in several morphological and/or physiological traits. Their microbe communities might also differ because of adaptation to new conditions related to cultivation. To test this hypothesis, we investigated morphological traits in a parthenogenetic root-knot nematode (Meloidogyne javanica) from natural and agricultural environments. Seventeen populations of M. javanica were sampled on cultivated and wild olives in Morocco, then maintained in controlled conditions for a ‘common garden’ experiment. We estimated the genetic variation based on three traits (stylet size, neck width and body width) by a quantitative genetic design (ten families per population and nine individuals per family were measured), and molecular variation was investigated with a mitochondrial marker to identify the genetic lineages of nematode isolates sampled from wild and cultivated olives. Significant morphological differences were detected between individuals from wild vs. cultivated hosts for the three traits, whereas no phylogenetic clustering was observed among isolates collected on those two hosts. Our results thus suggest an adaptive response of the asexual parasite, possibly related to the deep modification of soil nematode communities between natural olive stands and orchards.

The complete mitochondrial genome of Iphisa elegans (Reptilia: Squamata: Gymnophthalmidae)

The complete mitogenome of the lizard Iphisa elegans Gray, 1851 was sequenced using a shotgun approach on an Illumina HiSeq 3000 platform, providing the first mitogenome for Gymnophthalmidae. The genome was 18,622 bp long, with 13 protein-coding genes, two rRNA (12S and 16S), and 22 tRNA, as well as the control region. A maximum likelihood phylogenetic analysis including I. elegans and all other available mitogenomes of Squamata provided a tree in accordance with previous phylogenetic relationships inferred for Squamata.

The mitochondrial genome of Acrobeloides varius (Cephalobomorpha) confirms non-monophyly of Tylenchina (Nematoda)

The infraorder Cephalobomorpha is a diverse and ecologically important nematode group found in almost all terrestrial environments. In a recent nematode classification system based on SSU rDNA, Cephalobomorpha was classified within the suborder Tylenchina with Panagrolaimomorpha, Tylenchomorpha and Drilonematomorpha. However, phylogenetic relationships among species within Tylenchina are not always consistent, and the phylogenetic position of Cephalobomorpha is still uncertain. In this study, in order to examine phylogenetic relationships of Cephalobomorpha with other nematode groups, we determined the complete mitochondrial genome sequence of Acrobeloides varius, the first sequenced representative of Cephalobomorpha, and used this sequence for phylogenetic analyses along with 101 other nematode species. Phylogenetic analyses using amino acid and nucleotide sequence data of 12 protein-coding genes strongly support a sister relationship between the two cephalobomorpha species A. varius and Acrobeles complexus (represented by a partial mt genome sequence). In this mitochondrial genome phylogeny, Cephalobomorpha was sister to all chromadorean species (excluding Plectus acuminatus of Plectida) and separated from Panagrolaimomorpha and Tylenchomorpha, rendering Tylenchina non-monophyletic. Mitochondrial gene order among Tylenchina species is not conserved, and gene clusters shared between A. varius and A. complexus are very limited. Results from phylogenetic analysis and gene order comparison confirms Tylenchina is not monophyletic. To better understand phylogenetic relationships among Tylenchina members, additional mitochondrial genome information is needed from underrepresented taxa representing Panagrolaimomorpha and Cephalobomorpha.

On the Close Relatedness of Two Rice-Parasitic Root-Knot Nematode Species and the Recent Expansion of Meloidogyne graminicola in Southeast Asia

Meloidogyne graminicola is a facultative meiotic parthenogenetic root-knot nematode (RKN) that seriously threatens agriculture worldwide. We have little understanding of its origin, genomic structure, and intraspecific diversity. Such information would offer better knowledge of how this nematode successfully damages rice in many different environments. Previous studies on nuclear ribosomal DNA (nrDNA) suggested a close phylogenetic relationship between M. graminicola and Meloidogyne oryzae, despite their different modes of reproduction and geographical distribution. In order to clarify the evolutionary history of these two species and explore their molecular intraspecific diversity, we sequenced the genome of 12 M. graminicola isolates, representing populations of worldwide origins, and two South American isolates of M. oryzae. k-mer analysis of their nuclear genome and the detection of divergent homologous genomic sequences indicate that both species show a high proportion of heterozygous sites (ca. 1–2%), which had never been previously reported in facultative meiotic parthenogenetic RKNs. These analyses also point to a distinct ploidy level in each species, compatible with a diploid M. graminicola and a triploid M. oryzae. Phylogenetic analyses of mitochondrial genomes and three nuclear genomic sequences confirm close relationships between these two species, with M. graminicola being a putative parent of M. oryzae. In addition, comparative mitogenomics of those 12 M. graminicola isolates with a Chinese published isolate reveal only 15 polymorphisms that are phylogenetically non-informative. Eight mitotypes are distinguished, the most common one being shared by distant populations from Asia and America. This low intraspecific diversity, coupled with a lack of phylogeographic signal, suggests a recent worldwide expansion of M. graminicola.

Nematode Genome Announcement: A Draft Genome for Rice Root-Knot Nematode, Meloidogyne graminicola

The rice root-knot nematode Meloidogyne graminicola has emerged as a devastating pest of rice in South-East Asian countries. Here we present a draft genome sequence for M. graminicola , assembled using data from short and long insert libraries sequenced on Illumina GAIIx sequencing platform.

The mitochondrial genome of the oribatid mite Paraleius leontonychus: new insights into tRNA evolution and phylogenetic relationships in acariform mites

Bilaterian mitochondrial (mt) genomes are circular molecules that typically contain 37 genes. To date, only a single complete mitogenome sequence is available for the species-rich sarcoptiform mite order Oribatida. We sequenced the mitogenome of Paraleius leontonychus, another species of this suborder. It is 14,186 bp long and contains 35 genes, including only 20 tRNAs, lacking tRNA^Gly and tRNA^Tyr. Re-annotation of the mitogenome of Steganacarus magnus increased the number of mt tRNAs for this species to 12. As typical for acariform mites, many tRNAs are highly truncated in both oribatid species. The total number of tRNAs and the number of tRNAs with a complete cloverleaf-like structure in P. leontonychus, however, clearly exceeds the numbers previously reported for Sarcoptiformes. This indicates, contrary to what has been previously assumed, that reduction of tRNAs is not a general characteristic for sarcoptiform mites. Compared to other Sarcoptiformes, the two oribatid species have the least rearranged mt genome with respect to the pattern observed in Limulus polyphemus, a basal arachnid species. Phylogenetic analysis of the newly sequenced mt genome and previously published data on other acariform mites confirms paraphyly of the Oribatida and an origin of the Astigmata within the Oribatida.

A Hypersensitivity-Like Response to Meloidogyne graminicola in Rice (Oryza sativa)

Meloidogyne graminicola is a major plant-parasitic nematode affecting rice cultivation in Asia. Resistance to this nematode was found in the African rice genotypes Oryza glaberrima and O. longistaminata; however, due to interspecific hybrid sterility, the introgression of resistance genes in the widely consumed O. sativa varieties remains challenging. Recently, resistance was found in O. sativa and, here, we report for the first time the histological and genetic characterization of the resistance to M. graminicola in Zhonghua 11, an O. sativa variety. Bright-light microscopy and fluorescence observations of the root tissue of this variety revealed that the root cells surrounding the nematode displayed a hypersensitivity-like reaction with necrotic cells at early stages of infection when nematodes are migrating in the root's mesoderm. An accumulation of presumably phenolic compounds in the nematodes' neighboring root cells was also observed. In addition, at a later stage of infection, not only were few feeding sites observed but also the giant cells were underdeveloped, underlining an incompatible interaction. Furthermore, we generated a hybrid O. sativa population by crossing Zhonghua 11 with the susceptible O. sativa variety IR64 in order to describe the genetic background of this resistance. Our data suggested that the resistance to M. graminicola infection was qualitative rather than quantitative and, therefore, major resistance genes must be involved in this infection process. The full characterization of the defense mechanism and the preliminary study of the genetic inheritance of novel sources of resistance to Meloidogyne spp. in rice constitute a major step toward their use in crop breeding.

Reproductive, pathogenic and genotypic characterisation of five Meloidogyne graminicola populations from the Philippines on susceptible and resistant rice varieties

Five populations of Meloidogyne graminicola isolated from different rice-growing areas in the Philippines were characterised. The populations showed little phenotypic variability of second-stage juveniles and female perineal pattern. Differences in reproduction among M. graminicola populations were not observed on mature resistant Oryza glaberrima varieties ‘TOG5674’, ‘TOG5675’, ‘RAM131’ and ‘CG14’, or on susceptible O. sativa varieties ‘IR64’ and ‘UPLRi-5’. In all infected rice varieties, plant growth and yield-contributing traits showed no differences among the populations. A search on M. graminicola populations from the Philippines for single-nucleotide polymorphism on the sequences of Internal Transcribed Spacer (ITS) of rDNA genes and mtDNA indicated only few points of heteroplasmy. Nematode reproduction and disease induction of the five M. graminicola populations in the Philippines exerted the same level of aggressiveness and virulence. The absence of resistance-breaking populations of M. graminicola is important for the maintenance of durability of resistance to this important rice pathogen.

Mitochondrial genome diversity in dagger and needle nematodes (Nematoda: Longidoridae)

Dagger and needle nematodes included in the family Longidoridae (viz. Longidorus, Paralongidorus, and Xiphinema) are highly polyphagous plant-parasitic nematodes in wild and cultivated plants and some of them are plant-virus vectors (nepovirus). The mitochondrial (mt) genomes of the dagger and needle nematodes, Xiphinema rivesi, Xiphinema pachtaicum, Longidorus vineacola and Paralongidorus litoralis were sequenced in this study. The four circular mt genomes have an estimated size of 12.6, 12.5, 13.5 and 12.7 kb, respectively. Up to date, the mt genome of X. pachtaicum is the smallest genome found in Nematoda. The four mt genomes contain 12 protein-coding genes (viz. cox1-3, nad1-6, nad4L, atp6 and cob) and two ribosomal RNA genes (rrnL and rrnS), but the atp8 gene was not detected. These mt genomes showed a gene arrangement very different within the Longidoridae species sequenced, with the exception of very closely related species (X. americanum and X. rivesi). The sizes of non-coding regions in the Longidoridae nematodes were very small and were present in a few places in the mt genome. Phylogenetic analysis of all coding genes showed a closer relationship between Longidorus and Paralongidorus and different phylogenetic possibilities for the three Xiphinema species.

Meloidogyne graminicola: a major threat to rice agriculture

TAXONOMY

Superkingdom Eukaryota; Kingdom Metazoa; Phylum Nematoda; Class Chromadorea; Order Tylenchida; Suborder Tylenchina; Infraorder Tylenchomorpha; Superfamily Tylenchoidea; Family Meloidogynidae; Subfamily Meloidogyninae; Genus Meloidogyne.

BIOLOGY

Microscopic non-segmented roundworm. Plant pathogen; obligate sedentary endoparasitic root-knot nematode. Reproduction: facultative meiotic parthenogenetic species in which amphimixis can occur at a low frequency (c. 0.5%); relatively fast life cycle completed in 19-27 days on rice depending on the temperature range.

HOST RANGE

Reported to infect over 100 plant species, including cereals and grass plants, as well as dicotyledonous plants. Main host: rice (Oryza sativa).

SYMPTOMS

Characteristic hook-shaped galls (root swellings), mainly formed at the root tips of infected plants. Alteration of the root vascular system causes disruption of water and nutrient transport, stunting, chlorosis and loss of vigour, resulting in poor growth and reproduction of the plants with substantial yield losses in crops.

DISEASE CONTROL

Nematicides, chemical priming, constant immersion of rice in irrigated fields, crop rotation with resistant or non-host plants, use of nematode-free planting material. Some sources of resistance to Meloidogyne graminicola have been identified in African rice species (O. glaberrima and O. longistaminata), as well as in a few Asian rice cultivars.

AGRONOMIC IMPORTANCE

Major threat to rice agriculture, particularly in Asia. Adapted to flooded conditions, Meloidogyne graminicola causes problems in all types of rice agrosystems.

The Mitochondrial Genomes of the Zoonotic Canine Filarial Parasites Dirofilaria (Nochtiella) repens and Candidatus Dirofilaria (Nochtiella) hongkongensis Provide Evidence for Presence of Cryptic Species

Background

Cutaneous dirofilariosis is a canine mosquito-borne zoonosis that can cause larva migrans disease in humans. Dirofilaria repens is considered an emerging pathogen occurring with high prevalence in Mediterranean areas and many parts of tropical Asia. In Hong Kong, a second species, Candidatus Dirofilaria hongkongensis, has been reported. The present study aimed to compare mitochondrial genomes from these parasites and to obtain population genetic information.

Methods and Findings

Complete mitochondrial genomes were obtained by PCR and Sanger sequencing or ILLUMINA sequencing for four worms. Cytochrome oxidase subunit 1 sequences identified three as D. repens (all from Europe) and one as C. D. hongkongensis (from India). Mitochondrial genomes have the same organization as in other spirurid nematodes but a higher preference for thymine in the coding strand. Phylogenetic analysis was in contradiction to current taxonomy of the Onchocercidae but in agreement with a recent multi-locus phylogenetic analysis using both mitochondrial and nuclear markers. D. repens and C. D. hongkongensis sequences clustered together and were the common sister group to Dirofilaria immitis. Analysis of a 2.5 kb mitochondrial genome fragment from macrofilaria or canine blood samples from Europe (42), Thailand (2), India (1) and Vietnam (1) revealed only small genetic differences in the D. repens samples including all European and the Vietnam sample. The Indian C. D. hongkongensis and the two Thai samples formed separate clusters and differences were comparatively large.

Conclusion

Genetic differences between Dirofilaria spp. causing cutaneous disease can be considerable whereas D. repens itself was genetically quite homogenous. C. D. hongkongensis was identified for the first time from the Indian subcontinent. The full mitochondrial genome sequence strengthens the hypothesis that it represents an independent species and the Thai samples might represent another cryptic species, Candidatus Dirofilaria sp. ‘Thailand II’, or a quite divergent population of C. D. hongkongensis.

The mitochondrial genomes of the zoonotic parasite species Dirofilaria repens and Candidatus Dirofilaria hongkongensis were characterized and compared to the genomes of other filariae. The resulting phylogeny is largely in agreement with recent molecular data. C. D. hongkongensis was placed as a sister group to D. repens and both as a common sister to D. immitis. The large genetic difference between D. repens and C. D. hongkongensis further supports the hypothesis that both are distinct valid species. Two canine samples from Thailand that were diagnosed as D. repens are either from a C. D. hongkongensis population that is quite divergent from the Indian population or might represent another currently unrecognized species in the genus.

Shotgun sequencing of the mitochondrial genome of the Aldabra giant tortoise (Aldabrachelys gigantea)

The complete mitochondrial genome of the Aldabra giant tortoise [Aldabrachelys gigantea (Schweigger, 1812): Reptilia, Testudines, Testudinidae] was sequenced using a shotgun approach on an Illumina HiSeq 2500 platform (Illumina Inc., San Diego, CA). This genome was 16 467 bp long and presents the typical organization found in vertebrates. The mean coverage of sequencing was 116×. A phylogenetic analysis of the Testudinidae confirms the placement of Aldabrachelys in an Indian Ocean group (including Madagascar). This mitogenome constitutes a reference for ancient DNA analyses of the extinct Madagascan lineages of Aldabrachelys.

The complete mitochondrial genome of Anomaloglossus baeobatrachus (Amphibia: Anura: Aromobatidae)

The complete mitogenome of the rocket frog Anomaloglossus baeobatrachus was sequenced using a shotgun approach on an Illumina HiSeq 2500 platform (Illumina Inc., San Diego, CA), providing the first mitogenome for this genus. The genome was 17,572 bp long and presents the typical organization found in other neobatrachian anurans. A phylogenetic analysis including A. baeobatrachus and all other available mitogenomes of Hyloidea provided relationships in accordance with previous phylogenetic studies.

Gnathostoma spinigerum Mitochondrial Genome Sequence: a Novel Gene Arrangement and its Phylogenetic Position within the Class Chromadorea

Human gnathostomiasis is an emerging food-borne parasitic disease caused by nematodes in the genus Gnathostoma. In spite of their significance as pathogens, these parasites remain poorly understood at the molecular level. In the present study, we sequenced the mitochondrial (mt) genome of G. spinigerum, which infects a range of definitive hosts including dogs, cats, tigers, leopards and humans. The mt genome of G. spinigerum is 14,079 bp in size and shows substantial changes in gene order compared to other nematodes studied to date. Phylogenetic analyses of mt genome sequences by Bayesian inference (BI) revealed that the infraorder Gnathostomatomorpha (represented by G. spinigerum) is closely related to the infraorder Ascaridomorpha. G. spinigerum is the first species from the infraorder Gnathostomatomorpha for which a complete mt genome has been sequenced. The new data will help understand the evolution, population genetics and systematics of this medically important group of parasites.

Intraspecific variability of the facultative meiotic parthenogenetic root-knot nematode (Meloidogyne graminicola) from rice fields in Vietnam

Twenty years ago, the facultative meiotic parthenogenetic root-knot nematode (RKN), Meloidogyne graminicola, was recognised as an important rice pathogen in South Vietnam. Although this country is one of the most important rice exporters worldwide, a comprehensive picture of the occurrence of M. graminicola in Vietnamese rice fields is still not available. Therefore a nematode survey was carried out with the aim of better understanding the geographical distribution, and the pathogenic and genetic variability of the RKN in Vietnam. From the fields surveyed in a range of ecosystems, 21 RKN populations were recovered from infected rice roots. A diagnostic SCAR marker was developed showing that all Vietnamese populations belong to M. graminicola. Furthermore, sequencing of the Internal Transcribed Spacer (ITS) of the rDNA genes confirmed this identification. These populations were then characterised using morphometrics and pathogenicity tests (host plant range diversity, reproduction and virulence diversity) revealing intraspecific variability. We showed that morphometric traits are mainly genetically heritable characters with significant differences among the studied populations. Finally, a distinctive trait signature was found for the populations isolated from the upland rice cultures. All together, our study reveals the prevalence of M. graminicola populations in Vietnamese rice. Further investigations need to be developed to explore the population dynamics and evolutionary history of this species in South East Asia.

PCR-free shotgun sequencing of the stone loach mitochondrial genome (Barbatula barbatula)

The complete mitochondrial genome of the stone loach Barbatula barbatula (Linnaeus, 1758) (Actinopterygii: Cypriniformes: Nemacheilidae) has been sequenced using a genome-skimming approach on an Illumina Hiseq 2500 platform. The mitochondrial genome of B. barbatula was determined to be 16,630 bp long and presents an organization typical of vertebrate mitogenomes. The mean coverage was 82× with a minimum coverage of 33× for the control region and 52× for the remaining part of the genome. A phylogenetic analysis of the Nemacheilidae family shows the monophyly of the Barbatula genus with strong support.

Soup to Tree: The Phylogeny of Beetles Inferred by Mitochondrial Metagenomics of a Bornean Rainforest Sample

In spite of the growth of molecular ecology, systematics and next-generation sequencing, the discovery and analysis of diversity is not currently integrated with building the tree-of-life. Tropical arthropod ecologists are well placed to accelerate this process if all specimens obtained through mass-trapping, many of which will be new species, could be incorporated routinely into phylogeny reconstruction. Here we test a shotgun sequencing approach, whereby mitochondrial genomes are assembled from complex ecological mixtures through mitochondrial metagenomics, and demonstrate how the approach overcomes many of the taxonomic impediments to the study of biodiversity. DNA from approximately 500 beetle specimens, originating from a single rainforest canopy fogging sample from Borneo, was pooled and shotgun sequenced, followed by de novo assembly of complete and partial mitogenomes for 175 species. The phylogenetic tree obtained from this local sample was highly similar to that from existing mitogenomes selected for global coverage of major lineages of Coleoptera. When all sequences were combined only minor topological changes were induced against this reference set, indicating an increasingly stable estimate of coleopteran phylogeny, while the ecological sample expanded the tip-level representation of several lineages. Robust trees generated from ecological samples now enable an evolutionary framework for ecology. Meanwhile, the inclusion of uncharacterized samples in the tree-of-life rapidly expands taxon and biogeographic representation of lineages without morphological identification. Mitogenomes from shotgun sequencing of unsorted environmental samples and their associated metadata, placed robustly into the phylogenetic tree, constitute novel DNA “superbarcodes” for testing hypotheses regarding global patterns of diversity.

Comparative and evolutionary analyses of Meloidogyne spp. Based on mitochondrial genome sequences

Molecular taxonomy and evolution of nematodes have been recently the focus of several studies. Mitochondrial sequences were proposed as an alternative for precise identification of Meloidogyne species, to study intraspecific variability and to follow maternal lineages. We characterized the mitochondrial genomes (mtDNAs) of the root knot nematodes M. floridensis, M. hapla and M. incognita. These were AT rich (81–83%) and highly compact, encoding 12 proteins, 2 rRNAs, and 22 tRNAs. Comparisons with published mtDNAs of M. chitwoodi, M. incognita (another strain) and M. graminicola revealed that they share protein and rRNA gene order but differ in the order of tRNAs. The mtDNAs of M. floridensis and M. incognita were strikingly similar (97–100% identity for all coding regions). In contrast, M. floridensis, M. chitwoodi, M. hapla and M. graminicola showed 65–84% nucleotide identity for coding regions. Variable mitochondrial sequences are potentially useful for evolutionary and taxonomic studies. We developed a molecular taxonomic marker by sequencing a highly-variable ~2 kb mitochondrial region, nad5-cox1, from 36 populations of root-knot nematodes to elucidate relationships within the genus Meloidogyne. Isolates of five species formed monophyletic groups and showed little intraspecific variability. We also present a thorough analysis of the mitochondrial region cox2-rrnS. Phylogenies based on either mitochondrial region had good discrimination power but could not discriminate between M. arenaria, M. incognita and M. floridensis.

Mitochondrial genome plasticity among species of the nematode genus Meloidogyne (Nematoda: Tylenchina)

The mitochondrial (mt) genomes of the plant-parasitic root-knot nematodes Meloidogyne arenaria, Meloidogyne enterolobii and Meloidogyne javanica were sequenced and compared with those of three other root-knot nematode species in order to explore the mt genome plasticity within Meloidogyne. The mt genomes of M. arenaria, M. enterolobii and M. javanica are circular, with an estimated size of 18.8, 18.9 and 19.6 kb, respectively. Compared to other nematodes these mt genomes are larger, due to the presence of large non-coding regions. The mt genome architecture within the genus Meloidogyne varied in the position of trn genes and in the position, length and nucleotide composition of non-coding regions. These variations were observed independent of the species' natural environments or reproductive modes. M. enterolobii showed three main non-coding regions whereas Meloidogyne chitwoodi, Meloidogyne incognita, M. javanica and M. arenaria had two non-coding regions, and Meloidogyne graminicola had a unique large non-coding region interrupted by two trn genes. trn genes were positioned in different regions of the mt genomes in M. chitwoodi, M. enterolobii and M. graminicola, whereas the trn gene order was identical between M. arenaria, M. incognita and M. javanica. Importantly, M. graminicola had extra copies of trnV and trnS2. High divergence levels between the two copies of each trn might indicate duplication events followed by random loss and mutations in the anticodon. Tree-based methods based on amino acid sequences of 12 mt protein-coding genes support the monophyly for the tropical and mitotic parthenogenetic species, M. arenaria, M. enterolobii, M. incognita and M. javanica and for a clade that includes the meiotic parthenogenetic species, M. chitwoodi and M. graminicola. A comparison of the mt genome architecture in plant-parasitic nematodes and phylogenetic analyses support that Pratylenchus is the most recent ancestor of root-knot nematodes.

Sequencing of the mitochondrial genome of the avocado lace bug Pseudacysta perseae (Heteroptera, Tingidae) using a genome skimming approach

Lace bugs (Tingidae) are a family of phytophagous heteropterans, some of which are important agricultural and forestry pests. They currently comprise around 2500 species distributed worldwide, for which only one mitochondrial genome has been described so far. We sequenced the complete mitochondrial genome and the nuclear ribosomal gene segment of the avocado lace bug Pseudacysta perseae using a genome skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including the one of the sycamore lace bug Corythucha ciliata, were retrieved to allow for comparisons and phylogenetic analyses. P. perseae mitochondrial genome was determined to be 15,850 bp long, and presented the typical organisation of insect mitogenomes. The phylogenetic analysis placed P. perseae as a sister to C. ciliata but did not confirm the monophyly of Miroidae including Tingidae. Our results contradicted widely accepted phylogenetic hypothesis, which highlights the limits of analyses based on mitochondrial data only. Shotgun sequencing approaches should provide substantial improvements in harmonizing mitochondrial and nuclear databases.

Shotgun assembly of the assassin bug Brontostoma colossus mitochondrial genome (Heteroptera, Reduviidae)

The complete mitochondrial genome of the assassin bug Brontostoma colossus (Distant, 1902) (Heteroptera: Reduviidae) has been sequenced using a genome-skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including five assassin bug species, were retrieved to allow for comparisons and phylogenetic analyses. The mitochondrial genome of B. colossus was determined to be 16,625 bp long, and consists of 13 protein-coding genes (PCGs), 23 transfer-RNA genes (tRNAs), two ribosomal-RNA genes (rRNAs), and one control region. The nucleotide composition is biased toward adenine and thymine (A+T=73.4%). Overall, architecture, nucleotide composition and genome asymmetry are similar among all available assassin bug mitogenomes. All PCGs have usual start-codons (Met and Ile). Three T and two TA incomplete termination codons were identified adjacent to tRNAs, which was consistent with the punctuation model for primary transcripts processing followed by 3' polyadenylation of mature mRNA. All tRNAs exhibit the classic clover-leaf secondary structure except for tRNASer(AGN) in which the DHU arm forms a simple loop. Two notable features are present in the B. colossus mitogenome: (i) a 131 bp duplicated unit including the complete tRNAArg gene, resulting in 23 potentially functional tRNAs in total, and (ii) a 857 bp duplicated region comprising 277 bp of the srRNA gene and 580 bp of the control region. A phylogenetic analysis based on 55 true bug mitogenomes confirmed that B. colossus belongs to Reduviidae, but contradicted a widely accepted hypothesis. This highlights the limits of phylogenetic analyses based on mitochondrial data only.

The complete mitochondrial genome of Meloidogyne graminicola (Tylenchina): a unique gene arrangement and its phylogenetic implications

Meloidogyne graminicola is one of the most economically important plant parasitic-nematodes (PPNs). In the present study, we determined the complete mitochondrial (mt) DNA genome sequence of this plant pathogen. Compared with other PPNs genera, this genome (19,589 bp) is only slightly smaller than that of Pratylenchus vulnus (21,656 bp). The nucleotide composition of the whole mtDNA sequence of M. graminicola is significantly biased toward A and T, with T being the most favored nucleotide and C being the least favored. The A+T content of the entire genome is 83.51%. The mt genome of M. graminicola contains 36 genes (lacking atp8) that are transcribed in the same direction. The gene arrangement of the mt genome of M. graminicola is unique. A total of 21 out of 22 tRNAs possess a DHU loop only, while tRNASer(AGN) lacks a DHU loop. The two large noncoding regions (2,031 bp and 5,063 bp) are disrupted by tRNASer(UCN). Phylogenetic analysis based on concatenated amino acid sequences of 12 protein-coding genes support the monophylies of the three orders Rhabditida, Mermithida and Trichinellida, the suborder Rhabditina and the three infraorders Spiruromorpha, Oxyuridomorpha and Ascaridomorpha, but do not support the monophylies of the two suborders Spirurina and Tylenchina, and the three infraorders Rhabditomorpha, Panagrolaimomorpha and Tylenchomorpha. The four Tylenchomorpha species including M. graminicola, P. vulnus, H. glycines and R. similis from the superfamily Tylenchoidea are placed within a well-supported monophyletic clade, but far from the other two Tylenchomorpha species B. xylophilus and B. mucronatus of Aphelenchoidea. In the clade of Tylenchoidea, M. graminicola is sister to P. vulnus, and H. glycines is sister to R. similis, which suggests root-knot nematodes has a closer relationship to Pratylenchidae nematodes than to cyst nematodes.

The complete mitochondrial genome of Aphelenchoides besseyi (Nematoda: Aphelenchoididae), the first sequenced representative of the subfamily Aphelenchoidinae

The complete mitochondrial genome (mitogenome) ofAphelenchoides besseyiis 16 216 bp in size and has the typical organisation of nematode mitogenomes of Chromadorea, including 12 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and the AT-rich non-coding region. The nucleotide composition of the mitogenome ofA. besseyiis AT-biased (80.0%) and the AT skew is −0.289. The most common start codon forA. besseyiis ATT. Thenad3andnad4Lgenes have an incomplete stop codon consisting of just a T and the other PCGs stop with the full stop codons. All the tRNA genes display a non-typical cloverleaf structure of mitochondrial tRNA. The AT-rich non-coding region contains ten tandem repeat units with four different regions. Phylogenetic analysis based on concatenated amino acid sequences of 12 protein-coding genes showed that three Tylenchomorpha species, includingA. besseyi,Bursaphelenchus mucronatusandB. xylophilusfrom the superfamily Aphelenchoidea, are placed within a well-supported monophyletic clade, but far from the other six Tylenchomorpha speciesMeloidogyne chitwoodi,M. graminicola,M. incognita,Pratylenchus vulnus,Heterodera glycinesandRadopholus similisof Tylenchoidea. This phylogeny suggests thatAphelenchoideshas a close relative relationship withBursaphelenchusand that the Tylenchomorpha is not monophyletic.