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.

Hybridization and polyploidy enable genomic plasticity without sex in the most devastating plant-parasitic nematodes.

Root-knot nematodes (genus Meloidogyne) exhibit a diversity of reproductive modes ranging from obligatory sexual to fully asexual reproduction. Intriguingly, the most widespread and devastating species to global agriculture are those that reproduce asexually, without meiosis. To disentangle this surprising parasitic success despite the absence of sex and genetic exchanges, we have sequenced and assembled the genomes of three obligatory ameiotic and asexual Meloidogyne. We have compared them to those of relatives able to perform meiosis and sexual reproduction. We show that the genomes of ameiotic asexual Meloidogyne are large, polyploid and made of duplicated regions with a high within-species average nucleotide divergence of ~8%. Phylogenomic analysis of the genes present in these duplicated regions suggests that they originated from multiple hybridization events and are thus homoeologs. We found that up to 22% of homoeologous gene pairs were under positive selection and these genes covered a wide spectrum of predicted functional categories. To biologically assess functional divergence, we compared expression patterns of homoeologous gene pairs across developmental life stages using an RNAseq approach in the most economically important asexually-reproducing nematode. We showed that >60% of homoeologous gene pairs display diverged expression patterns. These results suggest a substantial functional impact of the genome structure. Contrasting with high within-species nuclear genome divergence, mitochondrial genome divergence between the three ameiotic asexuals was very low, signifying that these putative hybrids share a recent common maternal ancestor. Transposable elements (TE) cover a ~1.7 times higher proportion of the genomes of the ameiotic asexual Meloidogyne compared to the sexual relative and might also participate in their plasticity. The intriguing parasitic success of asexually-reproducing Meloidogyne species could be partly explained by their TE-rich composite genomes, resulting from allopolyploidization events, and promoting plasticity and functional divergence between gene copies in the absence of sex and meiosis.

SATQPCR: Website for statistical analysis of real-time quantitative PCR data.

SATQPCR is a web tool providing statistical analysis of real-time quantitative PCR data including all MIQE rules (gene efficiency, selection of reference genes and normalization with them). Our application is a quick tool that provides to the biologist, graphs as well as statistical tables summarizing their results with the chosen methods (t-test or ANOVA with Tukey test). The application is available at http://satqpcr.sophia.inra.fr with a demo dataset. Source code can be found at https://framagit.org/. SUPPLEMENTARY INFORMATION: Tutorials at http://satqpcr.sophia.inra.fr/cgi/help.cgi.

In Vitro Acquisition of Specific Small Interfering RNAs Inhibits the Expression of Some Target Genes in the Plant Ectoparasite Xiphinema index.

Xiphinema index is an important plant parasitic nematode that induces direct damages and specifically transmits the Grapevine fanleaf virus, which is particularly harmful for grapevines. Genomic resources of this nematode species are still limited and no functional gene validation technology is available. RNA interference (RNAi) is a powerful technology to study gene function and here we describe the application of RNAi on several genes in X. index. Soaking the nematodes for 48 h in a suspension containing specific small interfering RNAs resulted in a partial inhibition of the accumulation of some targeted mRNA. However, low reproducible silencing efficiency was observed which could arise from X. index silencing pathway deficiencies. Indeed, essential accustomed proteins for these pathways were not found in the X. index proteome predicted from transcriptomic data. The most reproducible silencing effect was obtained when targeting the piccolo gene potentially involved in endo-exocytosis of synaptic molecules. This represents the first report of gene silencing in a nematode belonging to the Longidoridae family.

Genome-wide expert annotation of the epigenetic machinery of the plant-parasitic nematodes Meloidogyne spp., with a focus on the asexually reproducing species.

BACKGROUND: The renewed interest in epigenetics has led to the understanding that both the environment and individual lifestyle can directly interact with the epigenome to influence its dynamics. Epigenetic phenomena are mediated by DNA methylation, stable chromatin modifications and non-coding RNA-associated gene silencing involving specific proteins called epigenetic factors. Multiple organisms, ranging from plants to yeast and mammals, have been used as model systems to study epigenetics. The interactions between parasites and their hosts are models of choice to study these mechanisms because the selective pressures are strong and the evolution is fast. The asexually reproducing root-knot nematodes (RKN) offer different advantages to study the processes and mechanisms involved in epigenetic regulation. RKN genomes sequencing and annotation have identified numerous genes, however, which of those are involved in the adaption to an environment and potentially relevant to the evolution of plant-parasitism is yet to be discovered. RESULTS: Here, we used a functional comparative annotation strategy combining orthology data, mining of curated genomics as well as protein domain databases and phylogenetic reconstructions. Overall, we show that (i) neither RKN, nor the model nematode Caenorhabditis elegans possess any DNA methyltransferases (DNMT) (ii) RKN do not possess the complete machinery for DNA methylation on the 6th position of adenine (6mA) (iii) histone (de)acetylation and (de)methylation pathways are conserved between C. elegans and RKN, and the corresponding genes are amplified in asexually reproducing RKN (iv) some specific non-coding RNA families found in plant-parasitic nematodes are dissimilar from those in C. elegans. In the asexually reproducing RKN Meloidogyne incognita, expression data from various developmental stages supported the putative role of these proteins in epigenetic regulations. CONCLUSIONS: Our results refine previous predictions on the epigenetic machinery of model species and constitute the most comprehensive description of epigenetic factors relevant to the plant-parasitic lifestyle and/or asexual mode of reproduction of RKN. Providing an atlas of epigenetic factors in RKN is an informative resource that will enable researchers to explore their potential role in adaptation of these parasites to their environment.

VapC10 toxin of the legume symbiont Sinorhizobium meliloti targets tRNASer and controls intracellular lifestyle.

The soil bacterium Sinorhizobium meliloti can establish a nitrogen-fixing symbiosis with the model legume Medicago truncatula. The rhizobia induce the formation of a specialized root organ called nodule, where they differentiate into bacteroids and reduce atmospheric nitrogen into ammonia. Little is known on the mechanisms involved in nodule senescence onset and in bacteroid survival inside the infected plant cells. Although toxin-antitoxin (TA) systems have been shown to promote intracellular survival within host cells in human pathogenic bacteria, their role in symbiotic bacteria was rarely investigated. S. meliloti encodes several TA systems, mainly of the VapBC family. Here we present the functional characterization, through a multidisciplinary approach, of the VapBC10 TA system of S. meliloti. Following a mapping by overexpression of an RNase in Escherichia coli (MORE) RNA-seq analysis, we demonstrated that the VapC10 toxin is an RNase that cleaves the anticodon loop of two tRNASer. Thereafter, a bioinformatics approach was used to predict VapC10 targets in bacteroids. This analysis suggests that toxin activation triggers a specific proteome reprogramming that could limit nitrogen fixation capability and viability of bacteroids. Accordingly, a vapC10 mutant induces a delayed senescence in nodules, associated to an enhanced bacteroid survival. VapBC10 TA system could contribute to S. meliloti adaptation to symbiotic lifestyle, in response to plant nitrogen status.

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.

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.

BABA-induced pathogen resistance: a multi-omics analysis of the tomato response reveals a hyper-receptive status involving ethylene.

Prior exposure to microbial-associated molecular patterns or specific chemical compounds can promote plants into a primed state with stronger defence responses. ß-aminobutyric acid (BABA) is an endogenous stress metabolite that induces resistance protecting various plants towards diverse stresses. In this study, by integrating BABA-induced changes in selected metabolites with transcriptome and proteome data, we generated a global map of the molecular processes operating in BABA-induced resistance (BABA-IR) in tomato. BABA significantly restricts the growth of the pathogens Oidium neolycopersici and Phytophthora parasitica but not Botrytis cinerea. A cluster analysis of the upregulated processes showed that BABA acts mainly as a stress factor in tomato. The main factor distinguishing BABA-IR from other stress conditions was the extensive induction of signaling and perception machinery playing a key role in effective resistance against pathogens. Interestingly, the signalling processes and immune response activated during BABA-IR in tomato differed from those in Arabidopsis with substantial enrichment of genes associated with jasmonic acid (JA) and ethylene (ET) signalling and no change in Asp levels. Our results revealed key differences between the effect of BABA on tomato and other model plants studied until now. Surprisingly, salicylic acid (SA) is not involved in BABA downstream signalization whereas ET and JA play a crucial role.

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.

A hierarchical classification of polysaccharide lyases for glycogenomics.

Carbohydrate-active enzymes face huge substrate diversity in a highly selective manner using only a limited number of available folds. They are therefore subjected to multiple divergent and convergent evolutionary events. This and their frequent modularity render their functional annotation in genomes difficult in a number of cases. In the present paper, a classification of polysaccharide lyases (the enzymes that cleave polysaccharides using an elimination instead of a hydrolytic mechanism) is shown thoroughly for the first time. Based on the analysis of a large panel of experimentally characterized polysaccharide lyases, we examined the correlation of various enzyme properties with the three levels of the classification: fold, family and subfamily. The resulting hierarchical classification, which should help annotate relevant genes in genomic efforts, is available and constantly updated at the Carbohydrate-Active Enzymes Database (http://www.cazy.org).

The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics.

The Carbohydrate-Active Enzyme (CAZy) database is a knowledge-based resource specialized in the enzymes that build and breakdown complex carbohydrates and glycoconjugates. As of September 2008, the database describes the present knowledge on 113 glycoside hydrolase, 91 glycosyltransferase, 19 polysaccharide lyase, 15 carbohydrate esterase and 52 carbohydrate-binding module families. These families are created based on experimentally characterized proteins and are populated by sequences from public databases with significant similarity. Protein biochemical information is continuously curated based on the available literature and structural information. Over 6400 proteins have assigned EC numbers and 700 proteins have a PDB structure. The classification (i) reflects the structural features of these enzymes better than their sole substrate specificity, (ii) helps to reveal the evolutionary relationships between these enzymes and (iii) provides a convenient framework to understand mechanistic properties. This resource has been available for over 10 years to the scientific community, contributing to information dissemination and providing a transversal nomenclature to glycobiologists. More recently, this resource has been used to improve the quality of functional predictions of a number genome projects by providing expert annotation. The CAZy resource resides at URL: http://www.cazy.org/.

Complete genome sequence of the complex carbohydrate-degrading marine bacterium, Saccharophagus degradans strain 2-40 T.

The marine bacterium Saccharophagus degradans strain 2-40 (Sde 2-40) is emerging as a vanguard of a recently discovered group of marine and estuarine bacteria that recycles complex polysaccharides. We report its complete genome sequence, analysis of which identifies an unusually large number of enzymes that degrade >10 complex polysaccharides. Not only is this an extraordinary range of catabolic capability, many of the enzymes exhibit unusual architecture including novel combinations of catalytic and substrate-binding modules. We hypothesize that many of these features are adaptations that facilitate depolymerization of complex polysaccharides in the marine environment. This is the first sequenced genome of a marine bacterium that can degrade plant cell walls, an important component of the carbon cycle that is not well-characterized in the marine environment.

Overlapping genes produce proteins with unusual sequence properties and offer insight into de novo protein creation.

It is widely assumed that new proteins are created by duplication, fusion, or fission of existing coding sequences. Another mechanism of protein birth is provided by overlapping genes. They are created de novo by mutations within a coding sequence that lead to the expression of a novel protein in another reading frame, a process called "overprinting." To investigate this mechanism, we have analyzed the sequences of the protein products of manually curated overlapping genes from 43 genera of unspliced RNA viruses infecting eukaryotes. Overlapping proteins have a sequence composition globally biased toward disorder-promoting amino acids and are predicted to contain significantly more structural disorder than nonoverlapping proteins. By analyzing the phylogenetic distribution of overlapping proteins, we were able to confirm that 17 of these had been created de novo and to study them individually. Most proteins created de novo are orphans (i.e., restricted to one species or genus). Almost all are accessory proteins that play a role in viral pathogenicity or spread, rather than proteins central to viral replication or structure. Most proteins created de novo are predicted to be fully disordered and have a highly unusual sequence composition. This suggests that some viral overlapping reading frames encoding hypothetical proteins with highly biased composition, often discarded as noncoding, might in fact encode proteins. Some proteins created de novo are predicted to be ordered, however, and whenever a three-dimensional structure of such a protein has been solved, it corresponds to a fold previously unobserved, suggesting that the study of these proteins could enhance our knowledge of protein space.

Characterization of Two Satellite DNA Families in the Genome of the Oomycete Plant Pathogen Phytophthora parasitica.

Satellite DNA is a class of repetitive sequences that are organized in long arrays of tandemly repeated units in most eukaryotes. Long considered as selfish DNA, satellite sequences are now proposed to contribute to genome integrity. Despite their potential impact on the architecture and evolution of the genome, satellite DNAs have not been investigated in oomycetes due to the paucity of genomic data and the difficulty of assembling highly conserved satellite arrays. Yet gaining knowledge on the structure and evolution of genomes of oomycete pathogens is crucial to understanding the mechanisms underlying adaptation to their environment and to proposing efficient disease control strategies. A de novo assembly of the genome of Phytophthora parasitica, an important oomycete plant pathogen, led to the identification of several families of tandemly repeated sequences varying in size, copy number, and sequence conservation. Among them, two abundant families, designated as PpSat1 and PpSat2, displayed typical features of satellite DNA and were collectively designated as PpSat. These two satellite families differ by their length, sequence, organization, genomic environment, and evolutionary dynamics. PpSat1, but not PpSat2, presented homologs among oomycetes. This observation, as well as the characterization of transcripts of PpSat families, suggested that these satellite DNA families likely play a conserved role within this important group of pathogens.

Transcriptomic and Ultrastructural Signatures of K+-Induced Aggregation in Phytophthora parasitica Zoospores.

Most pathogenic oomycetes of the genus Phytophthora spread in water films as flagellated zoospores. Zoospores perceive and produce signals attracting other zoospores, resulting in autoaggregation in vitro or biofilm formation on plant surface. The mechanisms underlying intercellular communication and consequent attraction, adhesion and aggregation are largely unknown. In Phytophthora parasitica, the perception of a K+ gradient induces coordinated motion and aggregation. To define cellular and molecular events associated with oomycete aggregation, we combined transcriptomic and ultrastructural analyses. Results indicate involvement of electroception in K+ sensing. They establish that the transcriptome repertoire required for swimming and aggregation is already fully functional at zoospore release. At the time points analyzed, aggregates are mainly constituted of zoospores. They produce vesicular and fibrillary material discharged at cell-to-cell contacts. Consistently, the signature of transcriptome dynamics during transition to aggregates is an upregulation of genes potentially related to vesicular trafficking. Moreover, transcriptomic and functional analyses show a strong enhancement of carbonic anhydrase activity, indicating that pH homeostasis may contribute to aggregation by acting on both zoospore movement and adhesion. This study poses the molecular and cellular bases of aggregative behavior within oomycetes and expands the current knowledge of ion perception-mediated dissemination of propagules in the rhizosphere.

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.

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.

FOLy: an integrated database for the classification and functional annotation of fungal oxidoreductases potentially involved in the degradation of lignin and related aromatic compounds.

The breakdown of lignin by fungi is a key step during carbon recycling in terrestrial ecosystems. This process is of great interest for green and white biotechnological applications. Given the importance of these enzymatic processes, we have classified the enzymes potentially involved in lignin catabolism into sequence-based families and integrated them in a newly developed database, designated Fungal Oxidative Lignin enzymes (FOLy). Families were defined after sequence similarity searches starting from protein sequences and validated by the convergence of results with biochemical experiments reported in the literature. The resulting database was applied as a tool for the functional annotation of genomes from different fungi, namely (i) the Basidiomycota Coprinopsis cinerea, Phanerochaete chrysosporium and Ustilago maydis and (ii) the Ascomycota Aspergillus nidulans and Trichoderma reesei. Genomic comparison of the oxidoreductases of these fungi revealed significant differences in the putative enzyme arsenals. Two Ascomycota fungal genomes were annotated and new candidate genes were identified that could be useful for lignin degradation and (or) melanin synthesis, and their function investigated experimentally. This database efforts aims at providing the means to get new insights for the understanding and biotechnological exploitation of the lignin degradation. A WWW server giving access to the routinely updated FOLy classifications of enzymes potentially involved in lignin degradation can be found at http://foly.esil.univ-mrs.fr.

Alienness: Rapid Detection of Candidate Horizontal Gene Transfers across the Tree of Life.

Horizontal gene transfer (HGT) is the transmission of genes between organisms by other means than parental to offspring inheritance. While it is prevalent in prokaryotes, HGT is less frequent in eukaryotes and particularly in Metazoa. Here, we propose Alienness, a taxonomy-aware web application available at http://alienness.sophia.inra.fr. Alienness parses BLAST results against public libraries to rapidly identify candidate HGT in any genome of interest. Alienness takes as input the result of a BLAST of a whole proteome of interest against any National Center for Biotechnology Information (NCBI) protein library. The user defines recipient (e.g., Metazoa) and donor (e.g., bacteria, fungi) branches of interest in the NCBI taxonomy. Based on the best BLAST E-values of candidate donor and recipient taxa, Alienness calculates an Alien Index (AI) for each query protein. An AI > 0 indicates a better hit to candidate donor than recipient taxa and a possible HGT. Higher AI represent higher gap of E-values between candidate donor and recipient and a more likely HGT. We confirmed the accuracy of Alienness on phylogenetically confirmed HGT of non-metazoan origin in plant-parasitic nematodes. Alienness scans whole proteomes to rapidly identify possible HGT in any species of interest and thus fosters exploration of HGT more easily and largely across the tree of life.