The complex hybrid origins of the root knot nematodes revealed through comparative genomics

Chromosome-level genome assembly of the cereal cyst nematode Heterodera flipjevi

As an economically important plant parasitic nematode (PPN), Heterodera filipjevi causes great damage on wheat, and now it was widely recorded in many countries. While multiple genomes of PPNs have been published, high-quality genome assembly and annotation on H. filipjevi have yet to be performed. This study presents a chromosome-scale genome assembly and annotation for H. filipjevi, utilizing a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing technologies. The genome consists of 9 pseudo-chromosomes that contain 134.19 Mb of sequence, with a scaffold N50 length of 11.88 Mb. In total, 10,036 genes were annotated, representing 75.20% of the total predicted protein-coding genes. Our study provides the first chromosome-scale genome for H. filipjevi, which is also the inaugural high-quality genome of cereal cyst nematodes (CCNs). It provides a valuable genomic resource for further biological research and pest management of cereal cyst nematodes disease.

Identification of Meloidogyne Species on Traditional Chinese Medicine Plants in the Qinling Mountain Area of China and Their Aggressiveness to Different Medicinal Herbs

Root-knot nematodes (Meloidogyne spp.) are plant-parasitic nematodes that cause serious damage worldwide. There are many species of traditional Chinese medicine (TCM) plants, but only a few have been reported to be infected by Meloidogyne species. From 2020 to 2022, a survey was conducted in the Qinling mountain area, which is the main production region of TCM plants in China. Obvious galling symptoms were observed on the root systems of 15 species of TCM plants. Females were collected from diverse diseased TCM plants and subsequently identified at morphological and molecular levels. Among the 20 diseased root samples collected, Meloidogyne hapla populations were identified in 12 samples (60%), and M. incognita populations were identified in eight samples (40%). Among the 15 species of diseased TCM plants, eight species, namely, Scutellaria baicalensis, Leonurus japonicus, Dioscorea zingiberensis, Cornus officinalis, Viola philippica, Achyranthes bidentata, Senecio scandens, and Plantago depressa, were reported to be infected by Meloidogyne species for the first time. The host status of five species of TCM plants for two M. hapla isolates and one M. incognita isolate from TCM plants in this study was then evaluated. Differences in TCM plants' response to nematode infection were apparent when susceptibility was evaluated by the egg counts per gram of fresh root and the reproduction factor of the nematodes. Among the five species of TCM plants tested, Salvia miltiorrhiza and Gynostemma pentaphyllum were the most susceptible, while S. baicalensis and V. philippica were not considered suitable hosts for M. hapla or M. incognita.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

The Plant Parasitic Nematodes Database: A Comprehensive Genomic Data Platform for Plant Parasitic Nematode Research

Plant parasitic nematodes are important phytopathogens that greatly affect the growth of agricultural and forestry plants. Scientists have conducted several studies to prevent and treat the diseases they cause. With the advent of the genomics era, the genome sequencing of plant parasitic nematodes has been considerably accelerated, and a large amount of data has been generated. This study developed the Plant Parasitic Nematodes Database (PPND), a platform to combine these data. The PPND contains genomic, transcriptomic, protein, and functional annotation data, allowing users to conduct BLAST searches and genome browser analyses and download bioinformatics data for in-depth research. PPND will be continuously updated, and new data will be integrated. PPND is anticipated to become a comprehensive genomics data platform for plant parasitic nematode research.

Unzipped chromosome-level genomes reveal allopolyploid nematode origin pattern as unreduced gamete hybridization

The formation and consequences of polyploidization in animals with clonal reproduction remain largely unknown. Clade I root-knot nematodes (RKNs), characterized by parthenogenesis and allopolyploidy, show a widespread geographical distribution and extensive agricultural destruction. Here, we generated 4 unzipped polyploid RKN genomes and identified a putative novel alternative telomeric element. Then we reconstructed 4 chromosome-level assemblies and resolved their genome structures as AAB for triploid and AABB for tetraploid. The phylogeny of subgenomes revealed polyploid RKN origin patterns as hybridization between haploid and unreduced gametes. We also observed extensive chromosomal fusions and homologous gene expression decrease after polyploidization, which might offset the disadvantages of clonal reproduction and increase fitness in polyploid RKNs. Our results reveal a rare pathway of polyploidization in parthenogenic polyploid animals and provide a large number of high-precision genetic resources that could be used for RKN prevention and control.

Counter-attack of biocontrol agents: Environmentally benign Approaches against Root-knot nematodes (Meloidogyne spp.) on Agricultural crops

Root-knot nematodes (Meloidogyne spp.) are obligate sedentary endoparasites, considered severe crop-damaging taxa among all plant-parasitic nematodes globally. Their attacks through parasitic proteins alter the physiology and machinery of the host cells to favour parasitism and reduction in crop yield. Currently, the use of excessive pesticides as a fast remedy to manage this pest is hazardous for both the environment and humans. Keeping this view in mind, there is an urgent need for developing efficient eco-friendly strategies. Bio-control as an eco-friendly is considered the best approach to manage nematodes without disturbing non-target microbes. In bio-control, living agents such as fungi and bacteria are the natural enemies of nematodes and the best substitute for pesticides. Fungi, including nematode-trapping fungi, can sense host signals and produce special trapping devices viz., constricting rings and adhesive knobs/loops, to capture nematodes and kill them. Whereas, endo-parasitic fungi kill nematodes by enzymatic secretions and spore adhesion through their hyphae. Bacteria can also control nematodes by producing antibiotic compounds, competing for nutrients and rhizosphere, production of hydrolytic enzymes viz., chitinases, proteases, lipases, and induction of systemic resistance (ISR) in host plants. Scientists throughout the world are trying to evolve environmentally benign methods that sustain agricultural production and keep nematodes below a threshold level. Whatever methods evolve, in the future the focus should be on important aspects like green approaches for managing nematodes without disturbing human health and the environment.

Biotechnological Tools to Elucidate the Mechanism of Plant and Nematode Interactions

Plant-parasitic nematodes (PPNs) pose a threat to global food security in both the developed and developing worlds. PPNs cause crop losses worth a total of more than USD 150 billion worldwide. The sedentary root-knot nematodes (RKNs) also cause severe damage to various agricultural crops and establish compatible relationships with a broad range of host plants. This review aims to provide a broad overview of the strategies used to identify the morpho-physiological and molecular events that occur during RKN parasitism. It describes the most current developments in the transcriptomic, proteomic, and metabolomic strategies of nematodes, which are important for understanding compatible interactions of plants and nematodes, and several strategies for enhancing plant resistance against RKNs. We will highlight recent rapid advances in molecular strategies, such as gene-silencing technologies, RNA interference (RNAi), and small interfering RNA (siRNA) effector proteins, that are leading to considerable progress in understanding the mechanism of plant-nematode interactions. We also take into account genetic engineering strategies, such as targeted genome editing techniques, the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9) (CRISPR/Cas-9) system, and quantitative trait loci (QTL), to enhance the resistance of plants against nematodes.

Identification, Diversity, and Distribution of Meloidogyne spp. in Vegetable Fields of South Georgia, U.S.A

Root-knot nematode (RKN; Meloidogyne spp.) is the most prevalent plant-parasitic nematode in vegetable fields of Georgia, with an incidence of 67.3%. Because aggressive RKN species are reported in the southeastern United States, molecular-based identification of RKN species was conducted on soil samples taken from a nematode surveillance study in 2018 from 292 RKN-infested vegetable fields in southern Georgia. The RKN-infested soil was potted with tomato cultivar Rutgers, and individual nematode females were isolated from galled roots and subjected to species-specific PCR and mitochondrial haplotype-based RKN species identification. The incidence (%), mean, and maximum relative abundance (second-stage juveniles per 100 cm3 of soil) of the five RKN species identified consisted of M. incognita (91.9, 486, 14,144), M. arenaria (36.0, 707, 14,144), M. floridensis (2.2, 909, 5,264), M. javanica (5.5, 352, 1,488), and M. haplanaria (0.7, 8, 14). A large proportion of fields (29%) had mixed populations of M. incognita and M. arenaria, which may reflect the region's long history of cotton and peanut cultivation. For unknown reasons, mixed populations of M. incognita and M. arenaria were associated with higher population densities. M. incognita is the most important RKN species in vegetable fields, followed by M. arenaria; therefore, pure or mixed populations of these species should be addressed in nematode management programs. Although at a lower incidence, the newly detected species, M. floridensis and M. haplanaria, have the potential to become a major threat since they reproduce on vegetables with Mi-resistant genes.

Population Genomic Analyses Suggest a Hybrid Origin, Cryptic Sexuality, and Decay of Genes Regulating Seed Development for the Putatively Strictly Asexual Kingdonia uniflora (Circaeasteraceae, Ranunculales)

Asexual lineages are perceived to be short-lived on evolutionary timescales. Hence, reports for exceptional cases of putative 'ancient asexuals' usually raise questions about the persistence of such species. So far, there have been few studies to solve the mystery in plants. The monotypic Kingdonia dating to the early Eocene, contains only K. uniflora that has no known definitive evidence for sexual reproduction nor records for having congeneric sexual species, raising the possibility that the species has persisted under strict asexuality for a long period of time. Here, we analyze whole genome polymorphism and divergence in K. uniflora. Our results show that K. uniflora is characterized by high allelic heterozygosity and elevated πN/πS ratio, in line with theoretical expectations under asexual evolution. Allele frequency spectrum analysis reveals the origin of asexuality in K. uniflora occurred prior to lineage differentiation of the species. Although divergence within K. uniflora individuals exceeds that between populations, the topologies of the two haplotype trees, however, fail to match each other, indicating long-term asexuality is unlikely to account for the high allele divergence and K. uniflora may have a recent hybrid origin. Phi-test shows a statistical probability of recombination for the conflicting phylogenetic signals revealed by the split network, suggesting K. uniflora engages in undetected sexual reproduction. Detection of elevated genetic differentiation and premature stop codons (in some populations) in genes regulating seed development indicates mutational degradation of sexuality-specific genes in K. uniflora. This study unfolds the origin and persistence mechanism of a plant lineage that has been known to reproduce asexually and presents the genomic consequences of lack of sexuality.

Rooting Out the Mechanisms of Root-Knot Nematode-Plant Interactions

Root-knot nematodes (RKNs; Meloidogyne spp.) engage in complex parasitic interactions with many different host plants around the world, initiating elaborate feeding sites and disrupting host root architecture. Although RKNs have been the focus of research for many decades, new molecular tools have provided useful insights into the biological mechanisms these pests use to infect and manipulate their hosts. From identifying host defense mechanisms underlying resistance to RKNs to characterizing nematode effectors that alter host cellular functions, the past decade of research has significantly expanded our understanding of RKN-plant interactions, and the increasing number of quality parasite and host genomes promises to enhance future research efforts into RKNs. In this review, we have highlighted recent discoveries, summarized the current understanding within the field, and provided links to new and useful resources for researchers. Our goal is to offer insights and tools to support the study of molecular RKN-plant interactions.

Silencing the conserved small nuclear ribonucleoprotein SmD1 target gene alters susceptibility to root-knot nematodes in plants

Root-knot nematodes (RKNs) are among the most damaging pests of agricultural crops. Meloidogyne is an extremely polyphagous genus of nematodes that can infect thousands of plant species. A few genes for resistance (R-genes) to RKN suitable for use in crop breeding have been identified, but virulent strains and species of RKN have emerged that render these R-genes ineffective. Secretion of RKN effectors targeting plant functions mediates the reprogramming of root cells into specialized feeding cells, the giant cells, essential for RKN development and reproduction. Conserved targets among plant species define the more relevant strategies for controlling nematode infection. The EFFECTOR18 (EFF18) protein from M. incognita interacts with the spliceosomal small nuclear ribonucleoprotein D1 (SmD1) in Arabidopsis (Arabidopsis thaliana), disrupting its function in alternative splicing regulation and modulating the giant cell transcriptome. We show here that EFF18 is a conserved RKN-specific effector that targets this conserved spliceosomal SmD1 protein in Solanaceae. This interaction modulates alternative splicing events produced by tomato (Solanum lycopersicum) in response to M. incognita infection. The alteration of SmD1 expression by virus-induced gene silencing in Solanaceae affects giant cell formation and nematode development. Thus, our work defines a promising conserved SmD1 target gene to develop broad resistance for the control of Meloidogyne spp. in plants.

The causative agents of fascioliasis in animals and humans: Parthenogenetic Fasciola in Asia and other regions

Parthenogenetic Fasciola is the causative agent of fascioliasis in animals and humans and is widely distributed in Asian countries, such as Japan, South Korea, China, Vietnam, Thailand, the Philippines, Myanmar, Bangladesh, Nepal, and India. Parthenogenetic Fasciola geographically originated from central and eastern China, where it exists between the habitats of Fasciola hepatica and Fasciola gigantica; it likely appeared thousands of years ago following hybridization between F. hepatica and F. gigantica. Parthenogenetic Fasciola consists of diploids and triploids that possess nuclear genome of both F. hepatica and F. gigantica and mitochondrial genome of either F. hepatica or F. gigantica. Maternal parents of parthenogenetic Fasciola are either F. hepatica having Fh-C4 haplotype or F. gigantica having Fg-C2 haplotype in mitochondrial NADH dehydrogenase subunit 1 (ND1) nucleotide sequences. Parthenogenetic Fasciola flukes with the Fh-C4 haplotype have spread from China to South Korea and Japan, whereas the flukes with the Fg-C2 haplotype have not only spread to Korea and Japan but also southward to Vietnam, Thailand, the Philippines, Myanmar, Bangladesh, Nepal, and India. Parthenogenetic Fasciola can be distinguished from F. hepatica and F. gigantica using combinational DNA sequence analysis of nuclear phosphoenolpyruvate carboxykinase (pepck) and DNA polymerase delta (pold) along with mitochondrial ND1 markers. The establishment of parthenogenetic Fasciola is expected as follows: parthenogenetic diploids with the Fh-C4 and Fg-C2 haplotypes first appeared based on single or multiple interspecific hybridization events; subsequently, parthenogenetic triploids emerged via backcross events between the maternal parthenogenetic diploid and either paternal bisexual F. hepatica or F. gigantica. Parthenogenetic Fasciola diploids and triploids then survived for thousands of years by clonal parthenogenetic reproduction, and generated descendants with ND1 haplotypes, which were derived from the Fh-C4 and Fg-C2 due to nucleotide substitution. Thus, the emergence of parthenogenetic Fasciola may be due to extremely uncommon and accidental events. Parthenogenetic Fasciola should be treated as a new asexual hybrid species.

Root-knot nematodes (Meloidogyne spp.) a threat to agriculture in Mexico: biology, current control strategies, and perspectives

Root-knot nematodes (RKN) are sedentary parasites of the roots of plants and are considered some of the most damaging pests in agriculture. Since RKN target the root vascular system, they provoke host nutrient deprivation and defective water transport, causing above-ground symptoms of growth stunting, wilting, chlorosis, and reduced crop yields. In Mexico RKN infestations are primarily dealt with by treating with synthetic chemically based nematicides that are preferred by farmers over available bioproducts. However, due to environmental and human health concerns chemical control is increasingly restricted. Biological control of RKNs can help reduce the use of chemical nematicides as it is achieved with antagonistic organisms, mainly bacteria, fungi, other nematodes, or consortia of diverse microorganisms, which control nematodes directly by predation and parasitism at different stages: eggs, juveniles, or adults; or indirectly by the action of toxic diffusible inhibitory metabolites. The need to increase agricultural production and reduce negative environmental impact creates an opportunity for optimizing biological control agents to suppress nematode populations, but this endeavour remains challenging as researchers around the world try to understand diverse control mechanisms, nematode and microbe life cycles, ecology, metabolite production, predatory behaviours, molecular and biochemical interactions, in order to generate attractive products with the approval of local regulatory bodies. Here, we provide a brief review of the biology of the genus Meloidogyne, biological control strategies, and a comparison between chemical and bioproducts in the Mexican market, and guidelines emitted by national agencies to ensure safety and effectiveness of new developments.

Identification of individual root-knot nematodes using low coverage long-read sequencing

Root-knot nematodes (RKN; genus Meloidogyne) are polyphagous plant pathogens of great economic importance to agriculturalists globally. These species are small, diverse, and can be challenging for accurate taxonomic identification. Many of the most important crop pests confound analysis with simple genetic marker loci as they are polyploids of likely hybrid origin. Here we take a low-coverage, long-read genome sequencing approach to characterisation of individual root-knot nematodes. We demonstrate library preparation for Oxford Nanopore Technologies Flongle sequencing of low input DNA from individual juveniles and immature females, multiplexing up to twelve samples per flow cell. Taxonomic identification with Kraken 2 (a k-mer-based taxonomic assignment tool) is shown to reliably identify individual nematodes to species level, even within the very closely related Meloidogyne incognita group. Our approach forms a robust, low-cost, and scalable method for accurate RKN species diagnostics.

A chromosomal assembly of the soybean cyst nematode genome

The soybean cyst nematode (Heterodera glycines) is a sedentary plant parasite that exceeds billion USD annually in yield losses. This problem is exacerbated by H. glycines populations overcoming the limited sources of natural resistance in soybean and by the lack of effective and safe alternative treatments. Although there are genetic determinants that render soybeans resistant to nematode genotypes, resistant soybeans are increasingly ineffective because their multiyear usage has selected for virulent H. glycines populations. Successful H. glycines infection relies on the comprehensive re-engineering of soybean root cells into a syncytium, as well as the long-term suppression of host defences to ensure syncytial viability. At the forefront of these complex molecular interactions are effectors, the proteins secreted by H. glycines into host root tissues. The mechanisms that control genomic effector acquisition, diversification, and selection are important insights needed for the development of essential novel control strategies. As a foundation to obtain this understanding, we created a nine-scaffold, 158 Mb pseudomolecule assembly of the H. glycines genome using PacBio, Chicago, and Hi-C sequencing. A Mikado consensus gene prediction produced an annotation of 22,465 genes using short- and long-read expression data. To evaluate assembly and annotation quality, we cross-examined synteny among H. glycines assemblies, and compared BUSCO across related species. To describe the predicted proteins involved in H. glycines' secretory pathway, we contrasted expression between preparasitic and parasitic stages with functional gene information. Here, we present the results from our assembly and annotation of the H. glycines genome and contribute this resource to the scientific community.

Chromosome Fusion Affects Genetic Diversity and Evolutionary Turnover of Functional Loci but Consistently Depends on Chromosome Size

Major changes in chromosome number and structure are linked to a series of evolutionary phenomena, including intrinsic barriers to gene flow or suppression of recombination due to chromosomal rearrangements. However, chromosome rearrangements can also affect the fundamental dynamics of molecular evolution within populations by changing relationships between linked loci and altering rates of recombination. Here, we build chromosome-level assembly Eueides isabella and, together with a recent chromosome-level assembly of Dryas iulia, examine the evolutionary consequences of multiple chromosome fusions in Heliconius butterflies. These assemblies pinpoint fusion points on 10 of the 20 autosomal chromosomes and reveal striking differences in the characteristics of fused and unfused chromosomes. The ten smallest autosomes in D. iulia and E. isabella, which have each fused to a longer chromosome in Heliconius, have higher repeat and GC content, and longer introns than predicted by their chromosome length. When fused, these characteristics change to become more in line with chromosome length. The fusions also led to reduced diversity, which likely reflects increased background selection and selection against introgression between diverging populations, following a reduction in per-base recombination rate. We further show that chromosome size and fusion impact turnover rates of functional loci at a macroevolutionary scale. Together these results provide further evidence that chromosome fusion in Heliconius likely had dramatic effects on population level processes shaping rates of neutral and adaptive divergence. These effects may have impacted patterns of diversification in Heliconius, a classic example of an adaptive radiation.

Genus-level evolutionary relationships of FAR proteins reflect the diversity of lifestyles of free-living and parasitic nematodes

Background

Nematodes are a widespread and diverse group comprising free-living and parasitic species, some of which have major detrimental effects on crops, animals, and human health. Genomic comparisons of nematodes may help reveal the genetic bases for the evolution of parasitic lifestyles. Fatty acid and retinol-binding proteins (FARs) are thought to be unique to nematodes and play essential roles in their development, reproduction, infection, and possibly parasitism through promoting the uptake, transport, and distribution of lipid and retinol. However, the evolution of FAR family proteins across the phylum Nematoda remains elusive.

Results

We report here the evolutionary relationship of the FAR gene family across nematodes. No FAR was found in Trichocephalida species and Romanomermis culicivorax from Clade I, and FAR could be found in species from Clades III, IV, and V. FAR proteins are conserved in Clade III species and separated into three clusters. Tandem duplications and high divergence events lead to variable richness and low homology of FARs in Steinernema of Clade IVa, Strongyloides of Clade IVb, and intestinal parasitic nematodes from Clades Vc and Ve. Moreover, different richness and sequence variations of FARs in pine wood, root-knot, stem, and cyst nematodes might be determined by reproduction mode or parasitism. However, murine lungworm Angiostrongylus and bovine lungworm Dictyocaulus viviparus from Clade Vd have only 3–4 orthologs of FAR. RNA-seq data showed that far genes, especially far-1 and far-2, were highly expressed in most nematodes. Angiostrongylus cantonensis FAR-1 and FAR-3 have low sequence homology and distinct ligand-binding properties, leading to differences in the cavity volume of proteins. These data indicate that FAR proteins diverged early and experienced low selective pressure to form genus-level diversity. The far genes are present in endophyte or root-colonized bacteria of Streptomyces, Kitasatospora sp., Bacillus subtilis, and Lysobacter, suggesting that bacterial far genes might be derived from plant-parasitic nematodes by horizontal gene transfer.

Conclusions

Data from these comparative analyses have provided insights into genus-level diversity of FAR proteins in the phylum Nematoda. FAR diversification provides a glimpse into the complicated evolution history across free-living and parasitic nematodes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01111-3.

Nematode Genome Announcement: Draft Genome of Meloidogyne chitwoodi, an Economically Important Pest of Potato in the Pacific Northwest

Meloidogyne chitwoodi is one of the most devastating pests of potato in the U.S. Pacific Northwest (PNW). Nematode-infected tubers develop external as well as internal defects, making the potato tubers unmarketable, and resulting in economic losses. Draft genome assemblies of three M. chitwoodi genotypes-race 1, race 2 and race 1 pathotype Roza-were generated using Illumina and PacBio Sequel RS II sequencing. The final assemblies consist of 30, 39, and 38 polished contigs for race 1, race 2 and race 1 pathotype Roza, respectively, with average N₅₀ of 2.37 Mb and average assembled genome size of approximately 47.41 Mb. On average, 10,508 genes were annotated for each genome. Benchmarking universal single-copy ortholog (BUSCO) analysis indicated that 69.80% of the BUSCOs were complete whereas 68.80, 0.93, and 12.67% were single copy, duplicated, and fragmented, respectively. These highly contiguous genomes will enrich resources to study potato-nematode interactions and enhance breeding efforts to develop nematode-resistant potato varieties for the PNW.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

How Do Pathogens Evolve Novel Virulence Activities?

This article is part of the Top 10 Unanswered Questions in MPMI invited review series.We consider the state of knowledge on pathogen evolution of novel virulence activities, broadly defined as anything that increases pathogen fitness with the consequence of causing disease in either the qualitative or quantitative senses, including adaptation of pathogens to host immunity and physiology, host species, genotypes, or tissues, or the environment. The evolution of novel virulence activities as an adaptive trait is based on the selection exerted by hosts on variants that have been generated de novo or arrived from elsewhere. In addition, the biotic and abiotic environment a pathogen experiences beyond the host may influence pathogen virulence activities. We consider host-pathogen evolution, host range expansion, and external factors that can mediate pathogen evolution. We then discuss the mechanisms by which pathogens generate and recombine the genetic variation that leads to novel virulence activities, including DNA point mutation, transposable element activity, gene duplication and neofunctionalization, and genetic exchange. In summary, if there is an (epi)genetic mechanism that can create variation in the genome, it will be used by pathogens to evolve virulence factors. Our knowledge of virulence evolution has been biased by pathogen evolution in response to major gene resistance, leaving other virulence activities underexplored. Understanding the key driving forces that give rise to novel virulence activities and the integration of evolutionary concepts and methods with mechanistic research on plant-microbe interactions can help inform crop protection.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genomic Features of Parthenogenetic Animals

Evolution without sex is predicted to impact genomes in numerous ways. Case studies of individual parthenogenetic animals have reported peculiar genomic features that were suggested to be caused by their mode of reproduction, including high heterozygosity, a high abundance of horizontally acquired genes, a low transposable element load, or the presence of palindromes. We systematically characterized these genomic features in published genomes of 26 parthenogenetic animals representing at least 18 independent transitions to asexuality. Surprisingly, not a single feature was systematically replicated across a majority of these transitions, suggesting that previously reported patterns were lineage-specific rather than illustrating the general consequences of parthenogenesis. We found that only parthenogens of hybrid origin were characterized by high heterozygosity levels. Parthenogens that were not of hybrid origin appeared to be largely homozygous, independent of the cellular mechanism underlying parthenogenesis. Overall, despite the importance of recombination rate variation for the evolution of sexual animal genomes, the genome-wide absence of recombination does not appear to have had the dramatic effects which are expected from classical theoretical models. The reasons for this are probably a combination of lineage-specific patterns, the impact of the origin of parthenogenesis, and a survivorship bias of parthenogenetic lineages.

The Meloidogyne incognita Nuclear Effector MiEFF1 Interacts With Arabidopsis Cytosolic Glyceraldehyde-3-Phosphate Dehydrogenases to Promote Parasitism

Root-knot nematodes are obligate endoparasites that maintain a biotrophic relationship with their hosts over a period of several weeks. They induce the differentiation of root cells into specialized multinucleate hypertrophied feeding cells known as giant cells. Nematode effectors synthesized in the esophageal glands and injected into the plant tissue through the syringe-like stylet play a key role in giant cell ontogenesis. The Meloidogyne incognita MiEFF1 is one of the rare effectors of phytopathogenic nematodes to have been located in vivo in feeding cells. This effector specifically targets the giant cell nuclei. We investigated the Arabidopsis functions modulated by this effector, by using a yeast two-hybrid approach to identify its host targets. We characterized a universal stress protein (USP) and cytosolic glyceraldehyde-3-phosphate dehydrogenases (GAPCs) as the targets of MiEFF1. We validated the interaction of MiEFF1 with these host targets in the plant cell nucleus, by bimolecular fluorescence complementation (BiFC). A functional analysis with Arabidopsis GUS reporter lines and knockout mutant lines showed that GAPCs were induced in giant cells and that their non-metabolic functions were required for root-knot nematode infection. These susceptibility factors are potentially interesting targets for the development of new root-knot nematode control strategies.

Speciation and adaptive evolution reshape antioxidant enzymatic system diversity across the phylum Nematoda

Background

Nematodes have evolved to survive in diverse ecological niches and can be a serious burden on agricultural economy, veterinary medicine, and public health. Antioxidant enzymes in parasitic nematodes play a critical role in defending against host oxidative stress. However, the features of the evolution of antioxidant enzymes in the phylum Nematoda remain elusive.

Results

Here, we systematically investigated the evolution and gene expression of antioxidant enzymes in the genomes of 59 nematodes and transcriptomes of 20 nematodes. Catalase has been independently lost in several orders, suggesting that it is unnecessary for some nematodes. Unlike in mammals, phospholipid hydroperoxide glutathione peroxidase is widely distributed in nematodes, among which it has evolved independently. We found that superoxide dismutase (SOD) has been present throughout nematode evolutionary process, and the extracellular isoform (SOD3) is diverged from the corresponding enzyme in mammals and has undergone duplication and differentiation in several nematodes. Moreover, the evolution of intracellular and extracellular SOD isoforms in filaria strongly indicates that extracellular SOD3 originated from intracellular SOD1 and underwent rapid evolution to form the diversity of extracellular SOD3. We identify a novel putative metal-independent extracellular SOD presenting independently in Steinernema and Strongyloididae lineage that featured a high expression level in Strongyloides larvae. Sequence divergence of SOD3 between parasitic nematodes and their closest free-living nematode, the specifically high expression in the parasitic female stage, and presence in excretory-secretory proteome of Strongyloides suggest that SOD3 may be related with parasitism.

Conclusions

This study advances our understanding of the complex evolution of antioxidant enzymes across Nematoda and provides targets for controlling parasitic nematode diseases.

Comparative Genomics Reveals Novel Target Genes towards Specific Control of Plant-Parasitic Nematodes

Plant-parasitic nematodes cause extensive annual yield losses to worldwide agricultural production. Most cultivated plants have no known resistance against nematodes and the few bearing a resistance gene can be overcome by certain species. Chemical methods that have been deployed to control nematodes have largely been banned from use due to their poor specificity and high toxicity. Hence, there is an urgent need for the development of cleaner and more specific control methods. Recent advances in nematode genomics, including in phytoparasitic species, provide an unprecedented opportunity to identify genes and functions specific to these pests. Using phylogenomics, we compared 61 nematode genomes, including 16 for plant-parasitic species and identified more than 24,000 protein families specific to these parasites. In the genome of Meloidogyne incognita, one of the most devastating plant parasites, we found ca. 10,000 proteins with orthologs restricted only to phytoparasitic species and no further homology in protein databases. Among these phytoparasite-specific proteins, ca. 1000 shared the same properties as known secreted effectors involved in essential parasitic functions. Of these, 68 were novel and showed strong expression during the endophytic phase of the nematode life cycle, based on both RNA-seq and RT-qPCR analyses. Besides effector candidates, transcription-related and neuro-perception functions were enriched in phytoparasite-specific proteins, revealing interesting targets for nematode control methods. This phylogenomics analysis constitutes a unique resource for the further understanding of the genetic basis of nematode adaptation to phytoparasitism and for the development of more efficient control methods.

Genome assembly and annotation of Meloidogyne enterolobii, an emerging parthenogenetic root-knot nematode

Root-knot nematodes (genus Meloidogyne) are plant parasites causing huge economic loss in the agricultural industry and affecting severely numerous developing countries. Control methods against these plant pests are sparse, the preferred one being the deployment of plant cultivars bearing resistance genes against Meloidogyne species. However, M. enterolobii is not controlled by the resistance genes deployed in the crop plants cultivated in Europe. The recent identification of this species in Europe is thus a major concern. Here, we sequenced the genome of M. enterolobii using short and long-read technologies. The genome assembly spans 240 Mbp with contig N50 size of 143 kbp, enabling high-quality annotations of 59,773 coding genes, 4,068 non-coding genes, and 10,944 transposable elements (spanning 8.7% of the genome). We validated the genome size by flow cytometry and the structure, quality and completeness by bioinformatics metrics. This ensemble of resources will fuel future projects aiming at pinpointing the genome singularities, the origin, diversity, and adaptive potential of this emerging plant pest.

Measurement(s)	genome • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing assay • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Meloidogyne enterolobii

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12410363

Genome structure and content of the rice root-knot nematode (Meloidogyne graminicola)

Discovered in the 1960s, Meloidogyne graminicola is a root-knot nematode species considered as a major threat to rice production. Yet, its origin, genomic structure, and intraspecific diversity are poorly understood. So far, such studies have been limited by the unavailability of a sufficiently complete and well-assembled genome. In this study, using a combination of Oxford Nanopore Technologies and Illumina sequencing data, we generated a highly contiguous reference genome (283 scaffolds with an N50 length of 294 kb, totaling 41.5 Mb). The completeness scores of our assembly are among the highest currently published for Meloidogyne genomes. We predicted 10,284 protein-coding genes spanning 75.5% of the genome. Among them, 67 are identified as possibly originating from horizontal gene transfers (mostly from bacteria), which supposedly contribute to nematode infection, nutrient processing, and plant defense manipulation. Besides, we detected 575 canonical transposable elements (TEs) belonging to seven orders and spanning 2.61% of the genome. These TEs might promote genomic plasticity putatively related to the evolution of M. graminicola parasitism. This high-quality genome assembly constitutes a major improvement regarding previously available versions and represents a valuable molecular resource for future phylogenomic studies of Meloidogyne species. In particular, this will foster comparative genomic studies to trace back the evolutionary history of M. graminicola and its closest relatives.

Genomic evidence for a hybrid origin of the yeast opportunistic pathogen Candida albicans

Background

Opportunistic yeast pathogens of the genus Candida are an important medical problem. Candida albicans, the most prevalent Candida species, is a natural commensal of humans that can adopt a pathogenic behavior. This species is highly heterozygous and cannot undergo meiosis, adopting instead a parasexual cycle that increases genetic variability and potentially leads to advantages under stress conditions. However, the origin of C. albicans heterozygosity is unknown, and we hypothesize that it could result from ancestral hybridization. We tested this idea by analyzing available genomes of C. albicans isolates and comparing them to those of hybrid and non-hybrid strains of other Candida species.

Results

Our results show compelling evidence that C. albicans is an evolved hybrid. The genomic patterns observed in C. albicans are similar to those of other hybrids such as Candida orthopsilosis MCO456 and Candida inconspicua, suggesting that it also descends from a hybrid of two divergent lineages. Our analysis indicates that most of the divergence between haplotypes in C. albicans heterozygous blocks was already present in a putative heterozygous ancestor, with an estimated 2.8% divergence between homeologous chromosomes. The levels and patterns of ancestral heterozygosity found cannot be fully explained under the paradigm of vertical evolution and are not consistent with continuous gene flux arising from lineage-specific events of admixture.

Conclusions

Although the inferred level of sequence divergence between the putative parental lineages (2.8%) is not clearly beyond current species boundaries in Saccharomycotina, we show here that all analyzed C. albicans strains derive from a single hybrid ancestor and diverged by extensive loss of heterozygosity. This finding has important implications for our understanding of C. albicans evolution, including the loss of the sexual cycle, the origin of the association with humans, and the evolution of virulence traits.

Attempt to Silence Genes of the RNAi Pathways of the Root-Knot Nematode, Meloidogyne incognita Results in Diverse Responses Including Increase and No Change in Expression of Some Genes

Control of plant-parasitic nematodes (PPNs) via host-induced gene silencing (HIGS) involves rational selection of genes and detailed assessment of effects of a possible knockdown on the nematode. Some genes by nature may be very important for the survival of the nematode that knockdown may be resisted. Possible silencing and effects of 20 such genes involved in the RNA interference (RNAi) pathways of Meloidogyne incognita were investigated in this study using long double-stranded RNAs (dsRNAs) as triggers. Two of the genes, ego-1 and mes-2, could not be knocked down. Expression of six genes (xpo-1, pash-1, xpo-2, rha-1, ekl-4, and csr-1) were significantly upregulated after RNAi treatment whereas for 12 of the genes, significant knockdown was achieved and with the exception of mes-2 and mes-6, RNAi was accompanied by defective phenotypes in treated nematodes including various degrees of paralysis and abnormal behaviors and movement such as curling, extreme wavy movements, and twitching. These abnormalities resulted in up to 75% reduction in infectivity of a tomato host, the most affected being the J2s previously treated with dsRNA of the gfl-1 gene. For 10 of the genes, effects of silencing in the J2s persisted as the adult females isolated from galls were under-developed, elongated, and transparent compared to the normal saccate, white adult females. Following RNAi of ego-1, smg-2, smg-6, and eri-1, reduced expression and/or the immediate visible effects on the J2s were not permanent as the nematodes infected and developed normally in tomato hosts. Equally intriguing was the results of RNAi of the mes-2 gene where the insignificant change in gene expression and behavior of treated J2s did not mean the nematodes were not affected as they were less effective in infecting host plants. Attempt to silence drsh-1, mut-7, drh-3, rha-1, pash-1, and vig-1 through HIGS led to reduction in nematode infestation by up to 89%. Our results show that genes may respond to RNAi knockdown differently so an exhaustive assessment of target genes as targets for nematode control via RNAi is imperative.

A Single, Shared Triploidy in Three Species of Parasitic Nematodes

The root-knot nematodes of the genus Meloidogyne are important and damaging parasites capable of infecting most flowering plants. Within this genus, several species of the Meloidogyne incognita group show evidence of paleopolyploidy in their genomes. We used our software tool POInT, the Polyploidy Orthology Inference Tool, to phylogenetically model the gene losses that followed that polyploidy. These models, and simulations based on them, show that three of these species (M. incognita, M. arenaria and M. javanica) descend from a single common hybridization event that yielded triplicated genomes with three distinguishable subgenomes. While one of the three subgenomes shows elevated gene loss rates relative to the other two, this subgenome does not show elevated sequence divergence. In all three species, ancestral loci where two of the three gene copies have been lost are less likely to have orthologs in Caenorhabditis elegans that are lethal when knocked down than are ancestral loci with surviving duplicate copies.

Genome sequence of the root-knot nematode Meloidogyne luci

Root-knot nematodes from the genus Meloidogyne are polyphagous plant endoparasites and agricultural pests of global importance. Here, we report the high-quality genome sequence of Meloidogyne luci population SI-Smartno V13. The resulting genome assembly of M. luci SI-Smartno V13 consists of 327 contigs, with an N50 contig length of 1,711,905 bp and a total assembly length of 209.16 Mb.

Signatures of the Evolution of Parthenogenesis and Cryptobiosis in the Genomes of Panagrolaimid Nematodes

Most animal species reproduce sexually and fully parthenogenetic lineages are usually short lived in evolution. Still, parthenogenesis may be advantageous as it avoids the cost of sex and permits colonization by single individuals. Panagrolaimid nematodes have colonized environments ranging from arid deserts to Arctic and Antarctic biomes. Many are obligatory meiotic parthenogens, and most have cryptobiotic abilities, being able to survive repeated cycles of complete desiccation and freezing. To identify systems that may contribute to these striking abilities, we sequenced and compared the genomes and transcriptomes of parthenogenetic and outcrossing panagrolaimid species, including cryptobionts and non-cryptobionts. The parthenogens are triploids, most likely originating through hybridization. Adaptation to cryptobiosis shaped the genomes of panagrolaimid nematodes and is associated with the expansion of gene families and signatures of selection on genes involved in cryptobiosis. All panagrolaimids have acquired genes through horizontal gene transfer, some of which are likely to contribute to cryptobiosis.

The genome of the migratory nematode, Radopholus similis, reveals signatures of close association to the sedentary cyst nematodes

Radopholus similis, commonly known as the burrowing nematode, is an important pest of myriad crops and ornamentals including banana (Musa spp.) and Citrus spp. In order to characterize the potential role of putative effectors encoded by R. similis genes we compared predicted proteins from a draft R. similis genome with other plant-parasitic nematodes in order to define the suite of excreted/secreted proteins that enable it to function as a parasite and to ascertain the phylogenetic position of R. similis in the Tylenchida order. Identification and analysis of candidate genes encoding for key plant cell-wall degrading enzymes including GH5 cellulases, PL3 pectate lyases and GH28 polygalactouranase revealed a pattern of occurrence similar to other PPNs, although with closest phylogenetic associations to the sedentary cyst nematodes. We also observed the absence of a suite of effectors essential for feeding site formation in the cyst nematodes. Clustering of various orthologous genes shared by R. similis with other nematodes showed higher overlap with the cyst nematodes than with the root-knot or other migratory endoparasitic nematodes. The data presented here support the hypothesis that R. similis is evolutionarily closer to the cyst nematodes, however, differences in the effector repertoire delineate ancient divergence of parasitism, probably as a consequence of niche specialization. These similarities and differences further underscore distinct evolutionary relationships during the evolution of parasitism in this group of nematodes.

Parthenogenesis and developmental constraints

The absence of a paternal contribution in an unfertilized ovum presents two developmental constraints against the evolution of parthenogenesis. We discuss the constraint caused by the absence of a centrosome and the one caused by the missing set of chromosomes and how they have been broken in specific taxa. They are examples of only a few well-underpinned examples of developmental constraints acting at macro-evolutionary scales in animals. Breaking of the constraint of the missing chromosomes is the best understood and generally involves rare occasions of drastic changes of meiosis. These drastic changes can be best explained by having been induced, or at least facilitated, by sudden cytological events (e.g., repeated rounds of hybridization, endosymbiont infections, and contagious infections). Once the genetic and developmental machinery is in place for regular or obligate parthenogenesis, shifts to other types of parthenogenesis can apparently rather easily evolve, for example, from facultative to obligate parthenogenesis, or from pseudoarrhenotoky to haplodiploidy. We argue that the combination of the two developmental constraints forms a near-absolute barrier against the gradual evolution from sporadic to obligate or regular facultative parthenogenesis, which can probably explain why the occurrence of the highly advantageous mode of regular facultative parthenogenesis is so rare and entirely absent in vertebrates.

Rice susceptibility to root-knot nematodes is enhanced by the Meloidogyne incognita MSP18 effector gene

MAIN CONCLUSION

This study revealed novel insights into the function of MSP18 effector during root-knot nematode parasitism in rice roots. MSP18 may modulate host immunity and enhance plant susceptibility to Meloidogyne spp. Rice (Oryza sativa) production is seriously impacted by root-knot nematodes (RKN), including Meloidogyne graminicola, Meloidogyne incognita, and Meloidogyne javanica, in upland and irrigated culture systems. Successful plant infection by RKN is likely achieved by releasing into the host cells some effector proteins to suppress the activation of immune responses. Here, we conducted a series of functional analyses to assess the role of the Meloidogyne-secreted protein (MSP) 18 from M. incognita (Mi-MSP18) during rice infection by RKN. Developmental expression profiles of M. javanica and M. graminicola showed that the MSP18 gene is up-regulated throughout nematode parasitic stages in rice. Reproduction of M. javanica and M. graminicola is enhanced in rice plants overexpressing Mi-MSP18, indicating that the Mi-MSP18 protein facilitates RKN parasitism. Transient expression assays in onion cells suggested that Mi-MSP18 is localized to the cytoplasm of the host cells. In tobacco, Mi-MSP18 suppressed the cell death induced by the INF1 elicitin, suggesting that Mi-MSP18 can interfere with the plant defense pathways. The data obtained in this study highlight Mi-MSP18 as a novel RKN effector able to enhance plant susceptibility and modulate host immunity.

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.

Unravelling the Genetic Diversity among Cassava Bemisia tabaci Whiteflies Using NextRAD Sequencing

Bemisia tabaci threatens production of cassava in Africa through vectoring viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). B. tabaci sampled from cassava in eight countries in Africa were genotyped using NextRAD sequencing, and their phylogeny and population genetics were investigated using the resultant single nucleotide polymorphism (SNP) markers. SNP marker data and short sequences of mitochondrial DNA cytochrome oxidase I (mtCOI) obtained from the same insect were compared. Eight genetically distinct groups were identified based on mtCOI, whereas phylogenetic analysis using SNPs identified six major groups, which were further confirmed by PCA and multidimensional analyses. STRUCTURE analysis identified four ancestral B. tabaci populations that have contributed alleles to the six SNP-based groups. Significant gene flows were detected between several of the six SNP-based groups. Evidence of gene flow was strongest for SNP-based groups occurring in central Africa. Comparison of the mtCOI and SNP identities of sampled insects provided a strong indication that hybrid populations are emerging in parts of Africa recently affected by the severe CMD pandemic. This study reveals that mtCOI is not an effective marker at distinguishing cassava-colonizing B. tabaci haplogroups, and that more robust SNP-based multilocus markers should be developed. Significant gene flows between populations could lead to the emergence of haplogroups that might alter the dynamics of cassava virus spread and disease severity in Africa. Continuous monitoring of genetic compositions of whitefly populations should be an essential component in efforts to combat cassava viruses in Africa.

Comparative Genomics of Apomictic Root-Knot Nematodes: Hybridization, Ploidy, and Dynamic Genome Change

The root-knot nematodes (genus Meloidogyne) are important plant parasites causing substantial agricultural losses. The Meloidogyne incognita group (MIG) of species, most of which are obligatory apomicts (mitotic parthenogens), are extremely polyphagous and important problems for global agriculture. While understanding the genomic basis for their variable success on different crops could benefit future agriculture, analyses of their genomes are challenging due to complex evolutionary histories that may incorporate hybridization, ploidy changes, and chromosomal fragmentation. Here, we sequence 19 genomes, representing five species of key root-knot nematodes collected from different geographic origins. We show that a hybrid origin that predated speciation within the MIG has resulted in each species possessing two divergent genomic copies. Additionally, the apomictic MIG species are hypotriploids, with a proportion of one genome present in a second copy. The hypotriploid proportion varies among species. The evolutionary history of the MIG genomes is revealed to be very dynamic, with noncrossover recombination both homogenizing the genomic copies, and acting as a mechanism for generating divergence between species. Interestingly, the automictic MIG species M. floridensis differs from the apomict species in that it has become homozygous throughout much of its genome.

Plant-Parasitic Nematodes and Food Security in Sub-Saharan Africa

Sub-Saharan Africa (SSA) is a region beset with challenges, not least its ability to feed itself. Low agricultural productivity, exploding populations, and escalating urbanization have led to declining per capita food availability. In order to reverse this trend, crop production systems must intensify, which brings with it an elevated threat from pests and diseases, including plant-parasitic nematodes. A holistic systems approach to pest management recognizes disciplinary integration. However, a critical under-representation of nematology expertise is a pivotal shortcoming, especially given the magnitude of the threat nematodes pose under more intensified systems. With more volatile climates, efficient use of water by healthy root systems is especially crucial. Within SSA, smallholder farming systems dominate the agricultural landscape, where a limited understanding of nematode problems prevails. This review provides a synopsis of current nematode challenges facing SSA and presents the opportunities to overcome current shortcomings, including a means to increase nematology capacity.

High-Quality Genome Sequence of the Root-Knot Nematode Meloidogyne arenaria Genotype A2-O

Root-knot nematodes (Meloidogyne spp.) cause serious damage to many crops globally. We report the high-quality genome sequence of Meloidogyne arenaria genotype A2-O. The genome assembly of M. arenaria A2-O is composed of 2,224 contigs with an N₅₀ contig length of 204,551 bp and a total assembly length of 284.05 Mb.

Host-mediated RNAi of a Notch-like receptor gene in Meloidogyne incognita induces nematode resistance

GLP-1 (abnormal germline proliferation) is a Notch-like receptor protein that plays an essential role in pharyngeal development. In this study, an orthologue of Caenorhabditis elegans glp-1 was identified in Meloidogyne incognita. A computational analysis revealed that the orthologue contained almost all the domains present in the C. elegans gene: specifically, the LIN-12/Notch repeat, the ankyrin repeat, a transmembrane domain and different ligand-binding motifs were present in orthologue, but the epidermal growth factor-like motif was not observed. An expression analysis showed differential expression of glp-1 throughout the life cycle of M. incognita, with relatively higher expression in the egg stage. To evaluate the silencing efficacy of Mi-glp-1, transgenic Arabidopsis plants carrying double-stranded RNA constructs of glp-1 were generated, and infection of these plants with M. incognita resulted in a 47-50% reduction in the numbers of galls, females and egg masses. Females obtained from the transgenic RNAi lines exhibited 40-60% reductions in the transcript levels of the targeted glp-1 gene compared with females isolated from the control plants. Second-generation juveniles (J2s), which were descendants of the infected females from the transgenic lines, showed aberrant phenotypes. These J2s exhibited a significant decrease in the overall distance from the stylet to the metacorpus region, and this effect was accompanied by disruption around the metacorporeal bulb of the pharynx. The present study suggests a role for this gene in organ (pharynx) development during embryogenesis in M. incognita and its potential use as a target in the management of nematode infestations in plants.

Comparative genomics of bdelloid rotifers: Insights from desiccating and nondesiccating species

Bdelloid rotifers are a class of microscopic invertebrates that have existed for millions of years apparently without sex or meiosis. They inhabit a variety of temporary and permanent freshwater habitats globally, and many species are remarkably tolerant of desiccation. Bdelloids offer an opportunity to better understand the evolution of sex and recombination, but previous work has emphasised desiccation as the cause of several unusual genomic features in this group. Here, we present high-quality whole-genome sequences of 3 bdelloid species: Rotaria macrura and R. magnacalcarata, which are both desiccation intolerant, and Adineta ricciae, which is desiccation tolerant. In combination with the published assembly of A. vaga, which is also desiccation tolerant, we apply a comparative genomics approach to evaluate the potential effects of desiccation tolerance and asexuality on genome evolution in bdelloids. We find that ancestral tetraploidy is conserved among all 4 bdelloid species, but homologous divergence in obligately aquatic Rotaria genomes is unexpectedly low. This finding is contrary to current models regarding the role of desiccation in shaping bdelloid genomes. In addition, we find that homologous regions in A. ricciae are largely collinear and do not form palindromic repeats as observed in the published A. vaga assembly. Consequently, several features interpreted as genomic evidence for long-term ameiotic evolution are not general to all bdelloid species, even within the same genus. Finally, we substantiate previous findings of high levels of horizontally transferred nonmetazoan genes in both desiccating and nondesiccating bdelloid species and show that this unusual feature is not shared by other animal phyla, even those with desiccation-tolerant representatives. These comparisons call into question the proposed role of desiccation in mediating horizontal genetic transfer.

A Genome Resequencing-Based Genetic Map Reveals the Recombination Landscape of an Outbred Parasitic Nematode in the Presence of Polyploidy and Polyandry

The parasitic nematode Haemonchus contortus is an economically and clinically important pathogen of small ruminants, and a model system for understanding the mechanisms and evolution of traits such as anthelmintic resistance. Anthelmintic resistance is widespread and is a major threat to the sustainability of livestock agriculture globally; however, little is known about the genome architecture and parameters such as recombination that will ultimately influence the rate at which resistance may evolve and spread. Here, we performed a genetic cross between two divergent strains of H. contortus, and subsequently used whole-genome resequencing of a female worm and her brood to identify the distribution of genome-wide variation that characterizes these strains. Using a novel bioinformatic approach to identify variants that segregate as expected in a pseudotestcross, we characterized linkage groups and estimated genetic distances between markers to generate a chromosome-scale F1 genetic map. We exploited this map to reveal the recombination landscape, the first for any helminth species, demonstrating extensive variation in recombination rate within and between chromosomes. Analyses of these data also revealed the extent of polyandry, whereby at least eight males were found to have contributed to the genetic variation of the progeny analyzed. Triploid offspring were also identified, which we hypothesize are the result of nondisjunction during female meiosis or polyspermy. These results expand our knowledge of the genetics of parasitic helminths and the unusual life-history of H. contortus, and enhance ongoing efforts to understand the genetic basis of resistance to the drugs used to control these worms and for related species that infect livestock and humans throughout the world. This study also demonstrates the feasibility of using whole-genome resequencing data to directly construct a genetic map in a single generation cross from a noninbred nonmodel organism with a complex lifecycle.

Characterization of Putative Effectors from the Cereal Cyst Nematode Heterodera avenae

Few molecular details of effectors of Heterodera avenae parasitism are known. We performed a high-throughput sequencing analysis of the H. avenae transcriptome at five developmental stages. A total of 82,549 unigenes were ultimately obtained, and 747 transcripts showed best hits to genes putatively encoding carbohydrate-active enzymes in plant-parasitic nematodes that play an important role in the invasion process. A total of 1,480 unigenes were homologous to known phytonematode effectors, and 63 putative novel effectors were identified in the H. avenae transcriptomes. Twenty-three unigenes were analyzed by qRT-PCR and confirmed to be highly expressed during at least one developmental stage. For in situ hybridization, 17 of the 22 tested putative effectors were specifically expressed and located in the subventral gland cells, and five putative novel effectors were specifically expressed in the dorsal gland. Furthermore, 115 transcripts were found to have putative lethal RNA interference (RNAi) phenotypes. Three target genes with lethal RNAi phenotypes and two of the four tested putative effectors were associated with a decrease in the number of cysts through in vitro RNAi technology. These transcriptomic data lay a foundation for further studies of interactions of H. avenae with cereal and H. avenae parasitic control.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class

BACKGROUND

Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes.

RESULTS

The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder.

CONCLUSIONS

Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.

On the origin of asexual species by means of hybridization and drift

Several species of asexuals appear to have existed for millions of years. This is despite the prevalent view that natural selection is weakened without gene exchange, which should cause these organisms to rapidly go extinct. In theory, one can identify evolutionary long-lived asexuals from their allelic sequence divergence, also known as the 'Meselson effect', which leads to elevated within-individual diversity. Yet evidence that this phenomenon exists is mixed. Furthermore, several confounding factors can lead to similar outcomes, including the formation of asexual species by hybridization. Disentangling these factors has proved to be tricky, but Ament-Velásquez et al. (2016) have provided an elegant solution in this issue of Molecular Ecology. They studied transcriptomes and mitochondrial DNA from the Lineus genus of nemertean worms, which contains both sexual and asexual types, to first show that the asexual L. pseudolactus is a hybrid between a sexual and an asexual species. After isolating out diversity arising from this hybridization, they find subsequent evidence for the Meselson effect. This study sets a new standard for differentiating between the complex causes and consequences of asexuality.

The Transcriptomes of Xiphinema index and Longidorus elongatus Suggest Independent Acquisition of Some Plant Parasitism Genes by Horizontal Gene Transfer in Early-Branching Nematodes

Nematodes have evolved the ability to parasitize plants on at least four independent occasions, with plant parasites present in Clades 1, 2, 10 and 12 of the phylum. In the case of Clades 10 and 12, horizontal gene transfer of plant cell wall degrading enzymes from bacteria and fungi has been implicated in the evolution of plant parasitism. We have used ribonucleic acid sequencing (RNAseq) to generate reference transcriptomes for two economically important nematode species, Xiphinema index and Longidorus elongatus, representative of two genera within the early-branching Clade 2 of the phylum Nematoda. We used a transcriptome-wide analysis to identify putative horizontal gene transfer events. This represents the first in-depth transcriptome analysis from any plant-parasitic nematode of this clade. For each species, we assembled ~30 million Illumina reads into a reference transcriptome. We identified 62 and 104 transcripts, from X. index and L. elongatus, respectively, that were putatively acquired via horizontal gene transfer. By cross-referencing horizontal gene transfer prediction with a phylum-wide analysis of Pfam domains, we identified Clade 2-specific events. Of these, a GH12 cellulase from X. index was analysed phylogenetically and biochemically, revealing a likely bacterial origin and canonical enzymatic function. Horizontal gene transfer was previously shown to be a phenomenon that has contributed to the evolution of plant parasitism among nematodes. Our findings underline the importance and the extensiveness of this phenomenon in the evolution of plant-parasitic life styles in this speciose and widespread animal phylum.

Differences in the genetic control of early egg development and reproduction between C. elegans and its parthenogenetic relative D. coronatus

BACKGROUND

The free-living nematode Diploscapter coronatus is the closest known relative of Caenorhabditis elegans with parthenogenetic reproduction. It shows several developmental idiosyncracies, for example concerning the mode of reproduction, embryonic axis formation and early cleavage pattern (Lahl et al. in Int J Dev Biol 50:393-397, 2006). Our recent genome analysis (Hiraki et al. in BMC Genomics 18:478, 2017) provides a solid foundation to better understand the molecular basis of developmental idiosyncrasies in this species in an evolutionary context by comparison with selected other nematodes. Our genomic data also yielded indications for the view that D. coronatus is a product of interspecies hybridization.

RESULTS

In a genomic comparison between D. coronatus, C. elegans, other representatives of the genus Caenorhabditis and the more distantly related Pristionchus pacificus and Panagrellus redivivus, certain genes required for central developmental processes in C. elegans like control of meiosis and establishment of embryonic polarity were found to be restricted to the genus Caenorhabditis. The mRNA content of early D. coronatus embryos was sequenced and compared with similar stages in C. elegans and Ascaris suum. We identified 350 gene families transcribed in the early embryo of D. coronatus but not in the other two nematodes. Looking at individual genes transcribed early in D. coronatus but not in C. elegans and A. suum, we found that orthologs of most of these are present in the genomes of the latter species as well, suggesting heterochronic shifts with respect to expression behavior. Considerable genomic heterozygosity and allelic divergence lend further support to the view that D. coronatus may be the result of an interspecies hybridization. Expression analysis of early acting single-copy genes yields no indication for silencing of one parental genome.

CONCLUSIONS

Our comparative cellular and molecular studies support the view that the genus Caenorhabditis differs considerably from the other studied nematodes in its control of development and reproduction. The easy-to-culture parthenogenetic D. coronatus, with its high-quality draft genome and only a single chromosome when haploid, offers many new starting points on the cellular, molecular and genomic level to explore alternative routes of nematode development and reproduction.

Genome Architecture and Evolution of a Unichromosomal Asexual Nematode

Asexual reproduction in animals, though rare, is the main or exclusive mode of reproduction in some long-lived lineages. The longevity of asexual clades may be correlated with the maintenance of heterozygosity by mechanisms that rearrange genomes and reduce recombination. Asexual species thus provide an opportunity to gain insight into the relationship between molecular changes, genome architecture, and cellular processes. Here we report the genome sequence of the parthenogenetic nematode Diploscapter pachys with only one chromosome pair. We show that this unichromosomal architecture is shared by a long-lived clade of asexual nematodes closely related to the genetic model organism Caenorhabditis elegans. Analysis of the genome assembly reveals that the unitary chromosome arose through fusion of six ancestral chromosomes, with extensive rearrangement among neighboring regions. Typical nematode telomeres and telomeric protection-encoding genes are lacking. Most regions show significant heterozygosity; homozygosity is largely concentrated to one region and attributed to gene conversion. Cell-biological and molecular evidence is consistent with the absence of key features of meiosis I, including synapsis and recombination. We propose that D. pachys preserves heterozygosity and produces diploid embryos without fertilization through a truncated meiosis. As a prelude to functional studies, we demonstrate that D. pachys is amenable to experimental manipulation by RNA interference.