Unzipped genome assemblies of polyploid root-knot nematodes reveal unusual and clade-specific telomeric repeats.

Using long-read sequencing, we assembled and unzipped the polyploid genomes of Meloidogyne incognita, M. javanica and M. arenaria, three of the most devastating plant-parasitic nematodes. We found the canonical nematode telomeric repeat to be missing in these and other Meloidogyne genomes. In addition, we find no evidence for the enzyme telomerase or for orthologs of C. elegans telomere-associated proteins, suggesting alternative lengthening of telomeres. Instead, analyzing our assembled genomes, we identify species-specific composite repeats enriched mostly at one extremity of contigs. These repeats are G-rich, oriented, and transcribed, similarly to canonical telomeric repeats. We confirm them as telomeric using fluorescent in situ hybridization. These repeats are mostly found at one single end of chromosomes in these species. The discovery of unusual and specific complex telomeric repeats opens a plethora of perspectives and highlights the evolutionary diversity of telomeres despite their central roles in senescence, aging, and chromosome integrity.

Ionizing radiation responses appear incidental to desiccation responses in the bdelloid rotifer Adineta vaga.

BACKGROUND: The remarkable resistance to ionizing radiation found in anhydrobiotic organisms, such as some bacteria, tardigrades, and bdelloid rotifers has been hypothesized to be incidental to their desiccation resistance. Both stresses produce reactive oxygen species and cause damage to DNA and other macromolecules. However, this hypothesis has only been investigated in a few species. RESULTS: In this study, we analyzed the transcriptomic response of the bdelloid rotifer Adineta vaga to desiccation and to low- (X-rays) and high- (Fe) LET radiation to highlight the molecular and genetic mechanisms triggered by both stresses. We identified numerous genes encoding antioxidants, but also chaperones, that are constitutively highly expressed, which may contribute to the protection of proteins against oxidative stress during desiccation and ionizing radiation. We also detected a transcriptomic response common to desiccation and ionizing radiation with the over-expression of genes mainly involved in DNA repair and protein modifications but also genes with unknown functions that were bdelloid-specific. A distinct transcriptomic response specific to rehydration was also found, with the over-expression of genes mainly encoding Late Embryogenesis Abundant proteins, specific heat shock proteins, and glucose repressive proteins. CONCLUSIONS: These results suggest that the extreme resistance of bdelloid rotifers to radiation might indeed be a consequence of their capacity to resist complete desiccation. This study paves the way to functional genetic experiments on A. vaga targeting promising candidate proteins playing central roles in radiation and desiccation resistance.

Genome sequence and annotation of Periconia digitata a hopeful biocontrol agent of phytopathogenic oomycetes.

The Periconia fungal genus belongs to the phylum Ascomycota, order Pleosporales, family Periconiaceae. Periconia are found in many habitats, but little is known about their ecology. Several species from this genus produce bioactive molecules. Periconia digitata extracts were shown to be deadly active against the pine wilt nematode. Furthermore, P. digitata was shown to inhibit the plant pathogenic oomycete Phytophthora parasitica. Because P. digitata has great potential as a biocontrol agent and high quality genomic resources are still lacking in the Periconiaceae family, we generated long-read genomic data for P. digitata. Using PacBio Hifi sequencing technology, we obtained a highly-contiguous genome assembled in 13 chromosomes and totaling ca. 39 Mb. In addition, we produced a reference transcriptome, based on 12 different culture conditions, and proteomic data to support the genome annotation. Besides representing a new reference genome within the Periconiaceae, this work will contribute to our better understanding of the Eukaryotic tree of life and opens new possibilities in terms of biotechnological applications.

AvP: A software package for automatic phylogenetic detection of candidate horizontal gene transfers.

Horizontal gene transfer (HGT) is the transfer of genes between species outside the transmission from parent to offspring. Due to their impact on the genome and biology of various species, HGTs have gained broader attention, but high-throughput methods to robustly identify them are lacking. One rapid method to identify HGT candidates is to calculate the difference in similarity between the most similar gene in closely related species and the most similar gene in distantly related species. Although metrics on similarity associated with taxonomic information can rapidly detect putative HGTs, these methods are hampered by false positives that are difficult to track. Furthermore, they do not inform on the evolutionary trajectory and events such as duplications. Hence, phylogenetic analysis is necessary to confirm HGT candidates and provide a more comprehensive view of their origin and evolutionary history. However, phylogenetic reconstruction requires several time-consuming manual steps to retrieve the homologous sequences, produce a multiple alignment, construct the phylogeny and analyze the topology to assess whether it supports the HGT hypothesis. Here, we present AvP which automatically performs all these steps and detects candidate HGTs within a phylogenetic framework.

The genome and lifestage-specific transcriptomes of a plant-parasitic nematode and its host reveal susceptibility genes involved in trans-kingdom synthesis of vitamin B5.

Plant-parasitic nematodes are a major threat to crop production in all agricultural systems. The scarcity of classical resistance genes highlights a pressing need to find new ways to develop nematode-resistant germplasm. Here, we sequence and assemble a high-quality phased genome of the model cyst nematode Heterodera schachtii to provide a platform for the first system-wide dual analysis of host and parasite gene expression over time, covering all major parasitism stages. Analysis of the hologenome of the plant-nematode infection site identified metabolic pathways that were incomplete in the parasite but complemented by the host. Using a combination of bioinformatic, genetic, and biochemical approaches, we show that a highly atypical completion of vitamin B5 biosynthesis by the parasitic animal, putatively enabled by a horizontal gene transfer from a bacterium, is required for full pathogenicity. Knockout of either plant-encoded or now nematode-encoded steps in the pathway significantly reduces parasitic success. Our experiments establish a reference for cyst nematodes, further our understanding of the evolution of plant-parasitism by nematodes, and show that congruent differential expression of metabolic pathways in the infection hologenome represents a new way to find nematode susceptibility genes. The approach identifies genome-editing-amenable targets for future development of nematode-resistant crops.

Improvement of eukaryotic protein predictions from soil metagenomes.

During the last decades, metagenomics has highlighted the diversity of microorganisms from environmental or host-associated samples. Most metagenomics public repositories use annotation pipelines tailored for prokaryotes regardless of the taxonomic origin of contigs. Consequently, eukaryotic contigs with intrinsically different gene features, are not optimally annotated. Using a bioinformatics pipeline, we have filtered 7.9 billion contigs from 6,872 soil metagenomes in the JGI's IMG/M database to identify eukaryotic contigs. We have re-annotated genes using eukaryote-tailored methods, yielding 8 million eukaryotic proteins and over 300,000 orphan proteins lacking homology in public databases. Comparing the gene predictions we made with initial JGI ones on the same contigs, we confirmed our pipeline improves eukaryotic proteins completeness and contiguity in soil metagenomes. The improved quality of eukaryotic proteins combined with a more comprehensive assignment method yielded more reliable taxonomic annotation. This dataset of eukaryotic soil proteins with improved completeness, quality and taxonomic annotation reliability is of interest for any scientist aiming at studying the composition, biological functions and gene flux in soil communities involving eukaryotes.

Chromatin Landscape Dynamics in the Early Development of the Plant Parasitic Nematode Meloidogyne incognita.

In model organisms, epigenome dynamics underlies a plethora of biological processes. The role of epigenetic modifications in development and parasitism in nematode pests remains unknown. The root-knot nematode Meloidogyne incognita adapts rapidly to unfavorable conditions, despite its asexual reproduction. However, the mechanisms underlying this remarkable plasticity and their potential impact on gene expression remain unknown. This study provides the first insight into contribution of epigenetic mechanisms to this plasticity, by studying histone modifications in M. incognita. The distribution of five histone modifications revealed the existence of strong epigenetic signatures, similar to those found in the model nematode Caenorhabditis elegans. We investigated their impact on chromatin structure and their distribution relative to transposable elements (TE) loci. We assessed the influence of the chromatin landscape on gene expression at two developmental stages: eggs, and pre-parasitic juveniles. H3K4me3 histone modification was strongly correlated with high levels of expression for protein-coding genes implicated in stage-specific processes during M. incognita development. We provided new insights in the dynamic regulation of parasitism genes kept under histone modifications silencing. In this pioneering study, we establish a comprehensive framework for the importance of epigenetic mechanisms in the regulation of the genome expression and its stability in plant-parasitic nematodes.

Targeted transcriptomics reveals signatures of large-scale independent origins and concerted regulation of effector genes in Radopholus similis.

The burrowing nematode, Radopholus similis, is an economically important plant-parasitic nematode that inflicts damage and yield loss to a wide range of crops. This migratory endoparasite is widely distributed in warmer regions and causes extensive destruction to the root systems of important food crops (e.g., citrus, banana). Despite the economic importance of this nematode, little is known about the repertoire of effectors owned by this species. Here we combined spatially and temporally resolved next-generation sequencing datasets of R. similis to select a list of candidates for the identification of effector genes for this species. We confirmed spatial expression of transcripts of 30 new candidate effectors within the esophageal glands of R. similis by in situ hybridization, revealing a large number of pioneer genes specific to this nematode. We identify a gland promoter motif specifically associated with the subventral glands (named Rs-SUG box), a putative hallmark of spatial and concerted regulation of these effectors. Nematode transcriptome analyses confirmed the expression of these effectors during the interaction with the host, with a large number of pioneer genes being especially abundant. Our data revealed that R. similis holds a diverse and emergent repertoire of effectors, which has been shaped by various evolutionary events, including neofunctionalization, horizontal gene transfer, and possibly by de novo gene birth. In addition, we also report the first GH62 gene so far discovered for any metazoan and putatively acquired by lateral gene transfer from a bacterial donor. Considering the economic damage caused by R. similis, this information provides valuable data to elucidate the mode of parasitism of this nematode.

Chromosome-level genome assembly reveals homologous chromosomes and recombination in asexual rotifer Adineta vaga.

Bdelloid rotifers are notorious as a speciose ancient clade comprising only asexual lineages. Thanks to their ability to repair highly fragmented DNA, most bdelloid species also withstand complete desiccation and ionizing radiation. Producing a well-assembled reference genome is a critical step to developing an understanding of the effects of long-term asexuality and DNA breakage on genome evolution. To this end, we present the first high-quality chromosome-level genome assemblies for the bdelloid Adineta vaga, composed of six pairs of homologous (diploid) chromosomes with a footprint of paleotetraploidy. The observed large-scale losses of heterozygosity are signatures of recombination between homologous chromosomes, either during mitotic DNA double-strand break repair or when resolving programmed DNA breaks during a modified meiosis. Dynamic subtelomeric regions harbor more structural diversity (e.g., chromosome rearrangements, transposable elements, and haplotypic divergence). Our results trigger the reappraisal of potential meiotic processes in bdelloid rotifers and help unravel the factors underlying their long-term asexual evolutionary success.

Movements of transposable elements contribute to the genomic plasticity and species diversification in an asexually reproducing nematode pest.

Despite reproducing without sexual recombination, Meloidogyne incognita is an adaptive and versatile phytoparasitic nematode. This species displays a global distribution, can parasitize a large range of plants, and can overcome plant resistance in a few generations. The mechanisms underlying this adaptability remain poorly known. At the whole-genome level, only a few single nucleotide variations have been observed across different geographical isolates with distinct ranges of compatible hosts. Exploring other factors possibly involved in genomic plasticity is thus important. Transposable elements (TEs), by their repetitive nature and mobility, can passively and actively impact the genome dynamics. This is particularly expected in polyploid hybrid genomes such as the one of M. incognita. Here, we have annotated the TE content of M. incognita, analyzed the statistical properties of this TE landscape, and used whole-genome pool-seq data to estimate the mobility of these TEs across twelve geographical isolates, presenting variations in ranges of compatible host plants. DNA transposons are more abundant than retrotransposons, and the high similarity of TE copies to their consensus sequences suggests they have been at least recently active. We have identified loci in the genome where the frequencies of presence of a TE showed substantial variations across the different isolates. Overall, variations in TE frequencies across isolates followed their phylogenetic divergence, suggesting TEs participate in the species diversification. Compared with the M. incognita reference genome, we detected isolate and lineage-specific de novo insertion of some TEs, including within genic regions or in the upstream regulatory regions. We validated by PCR the insertion of some of these TEs inside genic regions, confirming TE movements have possible functional impacts. Overall, we show DNA transposons can drive genomic plasticity in M. incognita and their role in genome evolution of other parthenogenetic animal deserves further investigation.

Genome Expression Dynamics Reveal the Parasitism Regulatory Landscape of the Root-Knot Nematode Meloidogyne incognita and a Promoter Motif Associated with Effector Genes.

Root-knot nematodes (genus Meloidogyne) are the major contributor to crop losses caused by nematodes. These nematodes secrete effector proteins into the plant, derived from two sets of pharyngeal gland cells, to manipulate host physiology and immunity. Successful completion of the life cycle, involving successive molts from egg to adult, covers morphologically and functionally distinct stages and will require precise control of gene expression, including effector genes. The details of how root-knot nematodes regulate transcription remain sparse. Here, we report a life stage-specific transcriptome of Meloidogyne incognita. Combined with an available annotated genome, we explore the spatio-temporal regulation of gene expression. We reveal gene expression clusters and predicted functions that accompany the major developmental transitions. Focusing on effectors, we identify a putative cis-regulatory motif associated with expression in the dorsal glands, providing an insight into effector regulation. We combine the presence of this motif with several other criteria to predict a novel set of putative dorsal gland effectors. Finally, we show this motif, and thereby its utility, is broadly conserved across the Meloidogyne genus, and we name it Mel-DOG. Taken together, we provide the first genome-wide analysis of spatio-temporal gene expression in a root-knot nematode and identify a new set of candidate effector genes that will guide future functional analyses.

Dissecting protein domain variability in the core RNA interference machinery of five insect orders.

RNA interference (RNAi)-mediated gene silencing can be used to control specific insect pest populations. Unfortunately, the variable efficiency in the knockdown levels of target genes has narrowed the applicability of this technology to a few species. Here, we examine the current state of knowledge regarding the miRNA (micro RNA) and siRNA (small interfering RNA) pathways in insects and investigate the structural variability at key protein domains of the RNAi machinery. Our goal was to correlate domain variability with mechanisms affecting the gene silencing efficiency. To this end, the protein domains of 168 insect species, encompassing the orders Coleoptera, Diptera, Hemiptera, Hymenoptera, and Lepidoptera, were analysed using our pipeline, which takes advantage of meticulous structure-based sequence alignments. We used phylogenetic inference and the evolutionary rate coefficient (K) to outline the variability across domain regions and surfaces. Our results show that four domains, namely dsrm, Helicase, PAZ and Ribonuclease III, are the main contributors of protein variability in the RNAi machinery across different insect orders. We discuss the potential roles of these domains in regulating RNAi-mediated gene silencing and the role of loop regions in fine-tuning RNAi efficiency. Additionally, we identified several order-specific singularities which indicate that lepidopterans have evolved differently from other insect orders, possibly due to constant coevolution with plants and viruses. In conclusion, our results highlight several variability hotspots that deserve further investigation in order to improve the application of RNAi technology in the control of insect pests.

Recent Advances in Population Genomics of Plant-Parasitic Nematodes.

Plant-parasitic nematodes are a costly burden of crop production. Ubiquitous in nature, phytoparasitic nematodes are associated with nearly every important agricultural crop and represent a significant constraint on global food security. Population genetics is a key discipline in plant nematology to understand aspects of the life strategies of these parasites, in particular their modes of reproduction, geographic origins, evolutionary histories, and dispersion abilities. Advances in high-throughput sequencing technologies have enabled a recent but active effort in genomic analyses of plant-parasitic nematodes. Such genomic approaches applied to multiple populations are providing new insights into the molecular and evolutionary processes that underpin the establishment of these nematodes and into a better understanding of the genetic and mechanistic basis of their pathogenicity and adaptation to their host plants. In this review, we attempt to update information about genome resources and genotyping techniques useful for nematologists who are thinking about initiating population genomics or genome sequencing projects. This review is intended also to foster the development of population genomics in plant-parasitic nematodes through highlighting recent publications that illustrate the potential for this approach to identify novel molecular markers or genes of interest and improve our knowledge of the genome variability, pathogenicity, and evolutionary potential of plant-parasitic nematodes.

Author Correction: Genome assembly and annotation of Meloidogyne enterolobii, an emerging parthenogenetic root-knot nematode.

A Correction to this paper has been published: https://doi.org/10.1038/s41597-020-00747-0.

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 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.

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.

Characterization of raffinose metabolism genes uncovers a wild Arachis galactinol synthase conferring tolerance to abiotic stresses.

Raffinose family oligosaccharides (RFOs) are implicated in plant regulatory mechanisms of abiotic stresses tolerance and, despite their antinutritional proprieties in grain legumes, little information is available about the enzymes involved in RFO metabolism in Fabaceae species. In the present study, the systematic survey of legume proteins belonging to five key enzymes involved in the metabolism of RFOs (galactinol synthase, raffinose synthase, stachyose synthase, alpha-galactosidase, and beta-fructofuranosidase) identified 28 coding-genes in Arachis duranensis and 31 in A. ipaënsis. Their phylogenetic relationships, gene structures, protein domains, and chromosome distribution patterns were also determined. Based on the expression profiling of these genes under water deficit treatments, a galactinol synthase candidate gene (AdGolS3) was identified in A. duranensis. Transgenic Arabidopsis plants overexpressing AdGolS3 exhibited increased levels of raffinose and reduced stress symptoms under drought, osmotic, and salt stresses. Metabolite and expression profiling suggested that AdGolS3 overexpression was associated with fewer metabolic perturbations under drought stress, together with better protection against oxidative damage. Overall, this study enabled the identification of a promising GolS candidate gene for metabolic engineering of sugars to improve abiotic stress tolerance in crops, whilst also contributing to the understanding of RFO metabolism in legume species.

Evolutionarily conserved plant genes responsive to root-knot nematodes identified by comparative genomics.

Root-knot nematodes (RKNs, genus Meloidogyne) affect a large number of crops causing severe yield losses worldwide, more specifically in tropical and sub-tropical regions. Several plant species display high resistance levels to Meloidogyne, but a general view of the plant immune molecular responses underlying resistance to RKNs is still lacking. Combining comparative genomics with differential gene expression analysis may allow the identification of widely conserved plant genes involved in RKN resistance. To identify genes that are evolutionary conserved across plant species, we used OrthoFinder to compared the predicted proteome of 22 plant species, including important crops, spanning 214 Myr of plant evolution. Overall, we identified 35,238 protein orthogroups, of which 6,132 were evolutionarily conserved and universal to all the 22 plant species (PLAnts Common Orthogroups-PLACO). To identify host genes responsive to RKN infection, we analyzed the RNA-seq transcriptome data from RKN-resistant genotypes of a peanut wild relative (Arachis stenosperma), coffee (Coffea arabica L.), soybean (Glycine max L.), and African rice (Oryza glaberrima Steud.) challenged by Meloidogyne spp. using EdgeR and DESeq tools, and we found 2,597 (O. glaberrima), 743 (C. arabica), 665 (A. stenosperma), and 653 (G. max) differentially expressed genes (DEGs) during the resistance response to the nematode. DEGs' classification into the previously characterized 35,238 protein orthogroups allowed identifying 17 orthogroups containing at least one DEG of each resistant Arachis, coffee, soybean, and rice genotype analyzed. Orthogroups contain 364 DEGs related to signaling, secondary metabolite production, cell wall-related functions, peptide transport, transcription regulation, and plant defense, thus revealing evolutionarily conserved RKN-responsive genes. Interestingly, the 17 DEGs-containing orthogroups (belonging to the PLACO) were also universal to the 22 plant species studied, suggesting that these core genes may be involved in ancestrally conserved immune responses triggered by RKN infection. The comparative genomic approach that we used here represents a promising predictive tool for the identification of other core plant defense-related genes of broad interest that are involved in different plant-pathogen interactions.