Large tandem duplications affect gene expression, 3D organization, and plant-pathogen response.

Rapid plant genome evolution is crucial to adapt to environmental changes. Chromosomal rearrangements and gene copy number variation (CNV) are two important tools for genome evolution and sources for the creation of new genes. However, their emergence takes many generations. In this study, we show that in Arabidopsis thaliana, a significant loss of ribosomal RNA (rRNA) genes with a past history of a mutation for the chromatin assembly factor 1 (CAF1) complex causes rapid changes in the genome structure. Using long-read sequencing and microscopic approaches, we have identified up to 15 independent large tandem duplications in direct orientation (TDDOs) ranging from 60 kb to 1.44 Mb. Our data suggest that these TDDOs appeared within a few generations, leading to the duplication of hundreds of genes. By subsequently focusing on a line only containing 20% of rRNA gene copies (20rDNA line), we investigated the impact of TDDOs on 3D genome organization, gene expression, and cytosine methylation. We found that duplicated genes often accumulate more transcripts. Among them, several are involved in plant-pathogen response, which could explain why the 20rDNA line is hyper-resistant to both bacterial and nematode infections. Finally, we show that the TDDOs create gene fusions and/or truncations and discuss their potential implications for the evolution of plant genomes.

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.

miR778 mediates gene expression, histone modification, and DNA methylation during cyst nematode parasitism.

Despite the known critical regulatory functions of microRNAs, histone modifications, and DNA methylation in reprograming plant epigenomes in response to pathogen infection, the molecular mechanisms underlying the tight coordination of these components remain poorly understood. Here, we show how Arabidopsis (Arabidopsis thaliana) miR778 coordinately modulates the root transcriptome, histone methylation, and DNA methylation via post-transcriptional regulation of the H3K9 methyltransferases SU(var)3-9 homolog 5 (SUVH5) and SUVH6 upon infection by the beet cyst nematode Heterodera schachtii. miR778 post-transcriptionally silences SUVH5 and SUVH6 upon nematode infection. Manipulation of the expression of miR778 and its two target genes significantly altered plant susceptibility to H. schachtii. RNA-seq analysis revealed a key role of SUVH5 and SUVH6 in reprograming the transcriptome of Arabidopsis roots upon H. schachtii infection. In addition, chromatin immunoprecipitation (ChIP)-seq analysis established SUVH5 and SUVH6 as the main enzymes mediating H3K9me2 deposition in Arabidopsis roots in response to nematode infection. ChIP-seq analysis also showed that these methyltransferases possess distinct DNA binding preferences in that they are targeting transposable elements under noninfected conditions and protein-coding genes in infected plants. Further analyses indicated that H3K9me2 deposition directed by SUVH5 and SUVH6 contributes to gene expression changes both in roots and in nematode feeding sites and preferentially associates with CG DNA methylation. Together, our results uncovered multi-layered epigenetic regulatory mechanisms coordinated by miR778 during Arabidopsis-H. schachtii interactions.

Esophageal Gland RNA-Seq Resource of a Virulent and Avirulent Population of the Soybean Cyst Nematode Heterodera glycines.

The soybean cyst nematode Heterodera glycines is the most economically devastating pathogen of soybean in the United States and threatens to become even more damaging through the selection of virulent nematode populations in the field that can overcome natural resistance mechanisms in soybean cultivars. This pathogen, therefore, demands intense transcriptomic/genomic research inquiries into the biology of its parasitic mechanisms. H. glycines delivers effector proteins that are produced in specialized gland cells into the soybean root to enable infection. The study of effector proteins, thus, is particularly promising when exploring novel management options against this pathogen. Here, we announce the availability of a gland cell-specific RNA-seq resource. These data represent an expression snapshot of gland cell activity during early soybean infection of a virulent and an avirulent H. glycines population, providing a unique and highly valuable resource for scientists examining effector biology and nematode virulence.[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.

Phytonematode peptide effectors exploit a host post-translational trafficking mechanism to the ER using a novel translocation signal.

Cyst nematodes induce a multicellular feeding site within roots called a syncytium. It remains unknown how root cells are primed for incorporation into the developing syncytium. Furthermore, it is unclear how CLAVATA3/EMBRYO SURROUNDING REGION (CLE) peptide effectors secreted into the cytoplasm of the initial feeding cell could have an effect on plant cells so distant from where the nematode is feeding as the syncytium expands. Here we describe a novel translocation signal within nematode CLE effectors that is recognized by plant cell secretory machinery to redirect these peptides from the cytoplasm to the apoplast of plant cells. We show that the translocation signal is functionally conserved across CLE effectors identified in nematode species spanning three genera and multiple plant species, operative across plant cell types, and can traffic other unrelated small peptides from the cytoplasm to the apoplast of host cells via a previously unknown post-translational mechanism of endoplasmic reticulum (ER) translocation. Our results uncover a mechanism of effector trafficking that is unprecedented in any plant pathogen to date, andthey illustrate how phytonematodes can deliver effector proteins into host cells and then hijack plant cellular processes for their export back out of the cell to function as external signaling molecules to distant cells.

An Effector from the Cyst Nematode Heterodera schachtii Derepresses Host rRNA Genes by Altering Histone Acetylation.

Cyst nematodes are plant-pathogenic animals that secrete effector proteins into plant root cells to alter host gene expression and reprogram these cells to form specialized feeding sites, known as syncytia. The molecular mechanisms of these effectors are mostly unknown. We determined that the sugar beet cyst nematode (Heterodera schachtii) 32E03 effector protein strongly inhibits the activities of Arabidopsis thaliana histone deacetylases including the HDT1 enzyme, which has a known function in the regulation of rRNA gene expression through chromatin modifications. We determined that plants expressing the 32E03 coding sequence exhibited increased acetylation of histone H3 along the rDNA chromatin. At low 32E03 expression levels, these chromatin changes triggered the derepression of a subset of rRNA genes, which were conducive to H. schachtii parasitism. By contrast, high levels of 32E03 caused profound bidirectional transcription along the rDNA, which triggered rDNA-specific small RNA production leading to RNA-directed DNA methylation and silencing of rDNA, which inhibited nematode development. Our data show that the 32E03 effector alters plant rRNA gene expression by modulating rDNA chromatin in a dose-dependent manner. Thus, the 32E03 effector epigenetically regulates plant gene expression to promote cyst nematode parasitism.

Re-targeting of a plant defense protease by a cyst nematode effector.

Plants mount defense responses during pathogen attacks, and robust host defense suppression by pathogen effector proteins is essential for infection success. 4E02 is an effector of the sugar beet cyst nematode Heterodera schachtii. Arabidopsis thaliana lines expressing the effector-coding sequence showed altered expression levels of defense response genes, as well as higher susceptibility to both the biotroph H. schachtii and the necrotroph Botrytis cinerea, indicating a potential suppression of defenses by 4E02. Yeast two-hybrid analyses showed that 4E02 targets A. thaliana vacuolar papain-like cysteine protease (PLCP) 'Responsive to Dehydration 21A' (RD21A), which has been shown to function in the plant defense response. Activity-based protein profiling analyses documented that the in planta presence of 4E02 does not impede enzymatic activity of RD21A. Instead, 4E02 mediates a re-localization of this protease from the vacuole to the nucleus and cytoplasm, which is likely to prevent the protease from performing its defense function and at the same time, brings it in contact with novel substrates. Yeast two-hybrid analyses showed that RD21A interacts with multiple host proteins including enzymes involved in defense responses as well as carbohydrate metabolism. In support of a role in carbohydrate metabolism of RD21A after its effector-mediated re-localization, we observed cell wall compositional changes in 4E02 expressing A. thaliana lines. Collectively, our study shows that 4E02 removes RD21A from its defense-inducing pathway and repurposes this enzyme by targeting the active protease to different cell compartments.

A new esophageal gland transcriptome reveals signatures of large scale de novo effector birth in the root lesion nematode Pratylenchus penetrans.

BACKGROUND: The root lesion nematode Pratylenchus penetrans is a migratory plant-parasitic nematode responsible for economically important losses in a wide number of crops. Despite the importance of P. penetrans, the molecular mechanisms employed by this nematode to promote virulence remain largely unknown. RESULTS: Here we generated a new and comprehensive esophageal glands-specific transcriptome library for P. penetrans. In-depth analysis of this transcriptome enabled a robust identification of a catalogue of 30 new candidate effector genes, which were experimentally validated in the esophageal glands by in situ hybridization. We further validated the expression of a multifaceted network of candidate effectors during the interaction with different plants. To advance our understanding of the "effectorome" of P. penetrans, we adopted a phylogenetic approach and compared the expanded effector repertoire of P. penetrans to the genome/transcriptome of other nematode species with similar or contrasting parasitism strategies. Our data allowed us to infer plausible evolutionary histories that shaped the effector repertoire of P. penetrans, as well as other close and distant plant-parasitic nematodes. Two remarkable trends were apparent: 1) large scale effector birth in the Pratylenchidae in general and P. penetrans in particular, and 2) large scale effector death in sedentary (endo) plant-parasitic nematodes. CONCLUSIONS: Our study doubles the number of validated Pratylenchus penetrans effectors reported in the literature. The dramatic effector gene gain in P. penetrans could be related to the remarkable ability of this nematode to parasitize a large number of plants. Our data provide valuable insights into nematode parasitism and contribute towards basic understating of the adaptation of P. penetrans and other root lesion nematodes to specific host plants.

A role for Arabidopsis growth-regulating factors 1 and 3 in growth-stress antagonism.

Growth-regulating factors (GRFs) belong to a small family of transcription factors that are highly conserved in plants. GRFs regulate many developmental processes and plant responses to biotic and abiotic stimuli. Despite the importance of GRFs, a detailed mechanistic understanding of their regulatory functions is still lacking. In this study, we used ChIP sequencing (ChIP-seq) to identify genome-wide binding sites of Arabidopsis GRF1 and GRF3, and correspondingly their direct downstream target genes. RNA-sequencing (RNA-seq) analysis revealed that GRF1 and GRF3 regulate the expression of a significant number of the identified direct targets. The target genes unveiled broad regulatory functions of GRF1 and GRF3 in plant growth and development, phytohormone biosynthesis and signaling, and the cell cycle. Our analyses also revealed that clock core genes and genes with stress- and defense-related functions are most predominant among the GRF1- and GRF3-bound targets, providing insights into a possible role for these transcription factors in mediating growth-defense antagonism and integrating environmental stimuli into developmental programs. Additionally, GRF1 and GRF3 target molecular nodes of growth-defense antagonism and modulate the levels of defense- and development-related hormones in opposite directions. Taken together, our results point to GRF1 and GRF3 as potential key determinants of plant fitness under stress conditions.

The genome of the soybean cyst nematode (Heterodera glycines) reveals complex patterns of duplications involved in the evolution of parasitism genes.

BACKGROUND: Heterodera glycines, commonly referred to as the soybean cyst nematode (SCN), is an obligatory and sedentary plant parasite that causes over a billion-dollar yield loss to soybean production annually. Although there are genetic determinants that render soybean plants resistant to certain nematode genotypes, resistant soybean cultivars are increasingly ineffective because their multi-year usage has selected for virulent H. glycines populations. The parasitic success of H. glycines relies on the comprehensive re-engineering of an infection site into a syncytium, as well as the long-term suppression of host defense 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 of effector acquisition, diversification, and selection need to be understood before effective control strategies can be developed, but the lack of an annotated genome has been a major roadblock. RESULTS: Here, we use PacBio long-read technology to assemble a H. glycines genome of 738 contigs into 123 Mb with annotations for 29,769 genes. The genome contains significant numbers of repeats (34%), tandem duplicates (18.7 Mb), and horizontal gene transfer events (151 genes). A large number of putative effectors (431 genes) were identified in the genome, many of which were found in transposons. CONCLUSIONS: This advance provides a glimpse into the host and parasite interplay by revealing a diversity of mechanisms that give rise to virulence genes in the soybean cyst nematode, including: tandem duplications containing over a fifth of the total gene count, virulence genes hitchhiking in transposons, and 107 horizontal gene transfers not reported in other plant parasitic nematodes thus far. Through extensive characterization of the H. glycines genome, we provide new insights into H. glycines biology and shed light onto the mystery underlying complex host-parasite interactions. This genome sequence is an important prerequisite to enable work towards generating new resistance or control measures against H. glycines.

Homeostasis in the soybean miRNA396-GRF network is essential for productive soybean cyst nematode infections.

Heterodera glycines, the soybean cyst nematode, penetrates soybean roots and migrates to the vascular cylinder where it forms a feeding site called the syncytium. MiRNA396 (miR396) targets growth-regulating factor (GRF) genes, and the miR396-GRF1/3 module is a master regulator of syncytium development in model cyst nematode H. schachtii infection of Arabidopsis. Here, we investigated whether this regulatory system operates similarly in soybean roots and is likewise important for H. glycines infection. We found that a network involving nine MIR396 and 23 GRF genes is important for normal development of soybean roots and that GRF function is specified in the root apical meristem by miR396. All MIR396 genes are down-regulated in the syncytium during its formation phase while, specifically, 11 different GRF genes are up-regulated. The switch to the syncytium maintenance phase coincides with up-regulation of MIR396 and down-regulation of the 11 GRF genes specifically via post-transcriptional regulation by miR396. Furthermore, interference with the miR396-GRF6/8-13/15-17/19 regulatory network, through either overexpression or knockdown experiments, does not affect the number of H. glycines juveniles that enter the vascular cylinder to initiate syncytia, but specifically inhibits efficient H. glycines development to adult females. Therefore, homeostasis in the miR396-GRF6/8-13/15-17/19 regulatory network is essential for productive H. glycines infections.

Suppression or Activation of Immune Responses by Predicted Secreted Proteins of the Soybean Rust Pathogen Phakopsora pachyrhizi.

Rust fungi, such as the soybean rust pathogen Phakopsora pachyrhizi, are major threats to crop production. They form specialized haustoria that are hyphal structures intimately associated with host-plant cell membranes. These haustoria have roles in acquiring nutrients and secreting effector proteins that manipulate host immune systems. Functional characterization of effector proteins of rust fungi is important for understanding mechanisms that underlie their virulence and pathogenicity. Hundreds of candidate effector proteins have been predicted for rust pathogens, but it is not clear how to prioritize these effector candidates for further characterization. There is a need for high-throughput approaches for screening effector candidates to obtain experimental evidence for effector-like functions, such as the manipulation of host immune systems. We have focused on identifying effector candidates with immune-related functions in the soybean rust fungus P. pachyrhizi. To facilitate the screening of many P. pachyrhizi effector candidates (named PpECs), we used heterologous expression systems, including the bacterial type III secretion system, Agrobacterium infiltration, a plant virus, and a yeast strain, to establish an experimental pipeline for identifying PpECs with immune-related functions and establishing their subcellular localizations. Several PpECs were identified that could suppress or activate immune responses in nonhost Nicotiana benthamiana, N. tabacum, Arabidopsis, tomato, or pepper plants.

A Small Cysteine-Rich Protein from the Asian Soybean Rust Fungus, Phakopsora pachyrhizi, Suppresses Plant Immunity.

The Asian soybean rust fungus, Phakopsora pachyrhizi, is an obligate biotrophic pathogen causing severe soybean disease epidemics. Molecular mechanisms by which P. pachyrhizi and other rust fungi interact with their host plants are poorly understood. The genomes of all rust fungi encode many small, secreted cysteine-rich proteins (SSCRP). While these proteins are thought to function within the host, their roles are completely unknown. Here, we present the characterization of P. pachyrhizi effector candidate 23 (PpEC23), a SSCRP that we show to suppress plant immunity. Furthermore, we show that PpEC23 interacts with soybean transcription factor GmSPL12l and that soybean plants in which GmSPL12l is silenced have constitutively active immunity, thereby identifying GmSPL12l as a negative regulator of soybean defenses. Collectively, our data present evidence for a virulence function of a rust SSCRP and suggest that PpEC23 is able to suppress soybean immune responses and physically interact with soybean transcription factor GmSPL12l, a negative immune regulator.

The novel cyst nematode effector protein 30D08 targets host nuclear functions to alter gene expression in feeding sites.

Cyst nematodes deliver effector proteins into host cells to manipulate cellular processes and establish a metabolically hyperactive feeding site. The novel 30D08 effector protein is produced in the dorsal gland of parasitic juveniles, but its function has remained unknown. We demonstrate that expression of 30D08 contributes to nematode parasitism, the protein is packaged into secretory granules and it is targeted to the plant nucleus where it interacts with SMU2 (homolog of suppressor of mec-8 and unc-52 2), an auxiliary spliceosomal protein. We show that SMU2 is expressed in feeding sites and an smu2 mutant is less susceptible to nematode infection. In Arabidopsis expressing 30D08 under the SMU2 promoter, several genes were found to be alternatively spliced and the most abundant functional classes represented among differentially expressed genes were involved in RNA processing, transcription and binding, as well as in development, and hormone and secondary metabolism, representing key cellular processes known to be important for feeding site formation. In conclusion, we demonstrated that the 30D08 effector is secreted from the nematode and targeted to the plant nucleus where its interaction with a host auxiliary spliceosomal protein may alter the pre-mRNA splicing and expression of a subset of genes important for feeding site formation.

Cooperative Regulatory Functions of miR858 and MYB83 during Cyst Nematode Parasitism.

MicroRNAs (miRNAs) recently have been established as key regulators of transcriptome reprogramming that define cell function and identity. Nevertheless, the molecular functions of the greatest number of miRNA genes remain to be determined. Here, we report cooperative regulatory functions of miR858 and its MYB83 transcription factor target gene in transcriptome reprogramming during Heterodera cyst nematode parasitism of Arabidopsis (Arabidopsis thaliana). Gene expression analyses and promoter-GUS fusion assays documented a role of miR858 in posttranscriptional regulation of MYB83 in the Heterodera schachtii-induced feeding sites, the syncytia. Constitutive overexpression of miR858 interfered with H. schachtii parasitism of Arabidopsis, leading to reduced susceptibility, while reduced miR858 abundance enhanced plant susceptibility. Similarly, MYB83 expression increases were conducive to nematode infection because overexpression of a noncleavable coding sequence of MYB83 significantly increased plant susceptibility, whereas a myb83 mutation rendered the plants less susceptible. In addition, RNA-seq analysis revealed that genes involved in hormone signaling pathways, defense response, glucosinolate biosynthesis, cell wall modification, sugar transport, and transcriptional control are the key etiological factors by which MYB83 facilitates nematode parasitism of Arabidopsis. Furthermore, we discovered that miR858-mediated silencing of MYB83 is tightly regulated through a feedback loop that might contribute to fine-tuning the expression of more than a thousand of MYB83-regulated genes in the H. schachtii-induced syncytium. Together, our results suggest a role of the miR858-MYB83 regulatory system in finely balancing gene expression patterns during H. schachtii parasitism of Arabidopsis to ensure optimal cellular function.

The cyst nematode effector protein 10A07 targets and recruits host posttranslational machinery to mediate its nuclear trafficking and to promote parasitism in Arabidopsis.

Plant-parasitic cyst nematodes synthesize and secrete effector proteins that are essential for parasitism. One such protein is the 10A07 effector from the sugar beet cyst nematode, Heterodera schachtii, which is exclusively expressed in the nematode dorsal gland cell during all nematode parasitic stages. Overexpression of H. schachtii 10A07 in Arabidopsis thaliana produced a hypersusceptible phenotype in response to H. schachtii infection along with developmental changes reminiscent of auxin effects. The 10A07 protein physically associates with a plant kinase and the IAA16 transcription factor in the cytoplasm and nucleus, respectively. The interacting plant kinase (IPK) phosphorylates 10A07 at Ser-144 and Ser-231 and mediates its trafficking from the cytoplasm to the nucleus. Translocation to the nucleus is phosphorylation dependent since substitution of Ser-144 and Ser-231 by alanine resulted in exclusive cytoplasmic accumulation of 10A07. IPK and IAA16 are highly upregulated in the nematode-induced syncytium (feeding cells), and deliberate manipulations of their expression significantly alter plant susceptibility to H. schachtii in an additive fashion. An inactive variant of IPK functioned antagonistically to the wild-type IPK and caused a dominant-negative phenotype of reduced plant susceptibility. Thus, exploitation of host processes to the advantage of the parasites is one mechanism by which cyst nematodes promote parasitism of host plants.

A soybean cyst nematode resistance gene points to a new mechanism of plant resistance to pathogens.

Soybean (Glycine max (L.) Merr.) is an important crop that provides a sustainable source of protein and oil worldwide. Soybean cyst nematode (Heterodera glycines Ichinohe) is a microscopic roundworm that feeds on the roots of soybean and is a major constraint to soybean production. This nematode causes more than US$1 billion in yield losses annually in the United States alone, making it the most economically important pathogen on soybean. Although planting of resistant cultivars forms the core management strategy for this pathogen, nothing is known about the nature of resistance. Moreover, the increase in virulent populations of this parasite on most known resistance sources necessitates the development of novel approaches for control. Here we report the map-based cloning of a gene at the Rhg4 (for resistance to Heterodera glycines 4) locus, a major quantitative trait locus contributing to resistance to this pathogen. Mutation analysis, gene silencing and transgenic complementation confirm that the gene confers resistance. The gene encodes a serine hydroxymethyltransferase, an enzyme that is ubiquitous in nature and structurally conserved across kingdoms. The enzyme is responsible for interconversion of serine and glycine and is essential for cellular one-carbon metabolism. Alleles of Rhg4 conferring resistance or susceptibility differ by two genetic polymorphisms that alter a key regulatory property of the enzyme. Our discovery reveals an unprecedented plant resistance mechanism against a pathogen. The mechanistic knowledge of the resistance gene can be readily exploited to improve nematode resistance of soybean, an increasingly important global crop.

Arabidopsis miR827 mediates post-transcriptional gene silencing of its ubiquitin E3 ligase target gene in the syncytium of the cyst nematode Heterodera schachtii to enhance susceptibility.

MicroRNAs (miRNAs) are a major class of small non-coding RNAs with emerging functions in biotic and abiotic interactions. Here, we report on a new functional role of Arabidopsis miR827 and its NITROGEN LIMITATION ADAPTATION (NLA) target gene in mediating plant susceptibility to the beet cyst nematode Heterodera schachtii. Cyst nematodes are sedentary endoparasites that induce the formation of multinucleated feeding structures termed syncytia in the roots of host plants. Using promoter:GUS fusion assays we established that miR827 was activated in the initial feeding cells and this activation was maintained in the syncytium during all sedentary stages of nematode development. Meanwhile, the NLA target gene, which encodes an ubiquitin E3 ligase enzyme, was post-transcriptionally silenced in the syncytium to permanently suppress its activity during all nematode parasitic stages. Overexpression of miR827 in Arabidopsis resulted in hyper-susceptibility to H. schachtii. In contrast, inactivation of miR827 activity through target mimicry or by overexpression a miR827-resistant cDNA of NLA produced the opposite phenotype of reduced plant susceptibility to H. schachtii. Gene expression analysis of several pathogenesis-related genes together with Agrobacterium-mediated transient expression in Nicotiana benthamiana provided strong evidence that miR827-mediated downregulation of NLA to suppress basal defense pathways. In addition, using yeast two-hybrid screens we identified several candidates of NLA-interacting proteins that are involved in a wide range of biological processes and molecular functions, including three pathogenesis-related proteins. Taken together, we conclude that nematode-activated miR827 in the syncytium is necessary to suppress immune responses in order to establish infection and cause disease.

Horizontal gene transfer of acetyltransferases, invertases and chorismate mutases from different bacteria to diverse recipients.

BACKGROUND: Hoplolaimina plant-parasitic nematodes (PPN) are a lineage of animals with many documented cases of horizontal gene transfer (HGT). In a recent study, we reported on three likely HGT candidate genes in the soybean cyst nematode Heterodera glycines, all of which encode secreted candidate effectors with putative functions in the host plant. Hg-GLAND1 is a putative GCN5-related N-acetyltransferase (GNAT), Hg-GLAND13 is a putative invertase (INV), and Hg-GLAND16 is a putative chorismate mutase (CM), and blastp searches of the non-redundant database resulted in highest similarity to bacterial sequences. Here, we searched nematode and non-nematode sequence databases to identify all the nematodes possible that contain these three genes, and to formulate hypotheses about when they most likely appeared in the phylum Nematoda. We then performed phylogenetic analyses combined with model selection tests of alternative models of sequence evolution to determine whether these genes were horizontally acquired from bacteria. RESULTS: Mining of nematode sequence databases determined that GNATs appeared in Hoplolaimina PPN late in evolution, while both INVs and CMs appeared before the radiation of the Hoplolaimina suborder. Also, Hoplolaimina GNATs, INVs and CMs formed well-supported clusters with different rhizosphere bacteria in the phylogenetic trees, and the model selection tests greatly supported models of HGT over descent via common ancestry. Surprisingly, the phylogenetic trees also revealed additional, well-supported clusters of bacterial GNATs, INVs and CMs with diverse eukaryotes and archaea. There were at least eleven and eight well-supported clusters of GNATs and INVs, respectively, from different bacteria with diverse eukaryotes and archaea. Though less frequent, CMs from different bacteria formed supported clusters with multiple different eukaryotes. Moreover, almost all individual clusters containing bacteria and eukaryotes or archaea contained species that inhabit very similar niches. CONCLUSIONS: GNATs were horizontally acquired late in Hoplolaimina PPN evolution from bacteria most similar to the saprophytic and plant-pathogenic actinomycetes. INVs and CMs were horizontally acquired from bacteria most similar to rhizobacteria and Burkholderia soil bacteria, respectively, before the radiation of Hoplolaimina. Also, these three gene groups appear to have been frequent subjects of HGT from different bacteria to numerous, diverse lineages of eukaryotes and archaea, which suggests that these genes may confer important evolutionary advantages to many taxa. In the case of Hoplolaimina PPN, this advantage likely was an improved ability to parasitize plants.