Big Genes, Small Effectors: Pea Aphid Cassette Effector Families Composed From Miniature Exons.

Aphids secrete proteins from their stylets that evidence indicates function similar to pathogen effectors for virulence. Here, we describe two small candidate effector gene families of the pea aphid, Acyrthosiphon pisum, that share highly conserved secretory signal peptide coding regions and divergent non-secretory coding sequences derived from miniature exons. The KQY candidate effector family contains eleven members with additional isoforms, generated by alternative splicing. Pairwise comparisons indicate possible four unique KQY families based on coding regions without the secretory signal region. KQY1a, a representative of the family, is encoded by a 968 bp mRNA and a gene that spans 45.7 kbp of the genome. The locus consists of 37 exons, 33 of which are 15 bp or smaller. Additional KQY members, as well as members of the KHI family, share similar features. Differential expression analyses indicate that the genes are expressed preferentially in salivary glands. Proteomic analysis on salivary glands and saliva revealed 11 KQY members in salivary proteins, and KQY1a was detected in an artificial diet solution after aphid feeding. A single KQY locus and two KHI loci were identified in Myzus persicae, the peach aphid. Of the genes that can be anchored to chromosomes, loci are mostly scattered throughout the genome, except a two-gene region (KQY4/KQY6). We propose that the KQY family expanded in A. pisum through combinatorial assemblies of a common secretory signal cassette and novel coding regions, followed by classical gene duplication and divergence.

Differential Expression of Candidate Salivary Effector Genes in Pea Aphid Biotypes With Distinct Host Plant Specificity.

Effector proteins play crucial roles in determining the outcome of various plant-parasite interactions. Aphids inject salivary effector proteins into plants to facilitate phloem feeding, but some proteins might trigger defense responses in certain plants. The pea aphid, Acyrthosiphon pisum, forms multiple biotypes, and each biotype is specialized to feed on a small number of closely related legume species. Interestingly, all the previously identified biotypes can feed on Vicia faba; hence, it serves as a universal host plant of A. pisum. We hypothesized that the salivary effector proteins have a key role in determining the compatibility between specific host species and A. pisum biotypes and that each biotype produces saliva containing a specific mixture of effector proteins due to differential expression of encoding genes. As the first step to address these hypotheses, we conducted two sets of RNA-seq experiments. RNA-seq analysis of dissected salivary glands (SGs) from reference alfalfa- and pea-specialized A. pisum lines revealed common and line-specific repertoires of candidate salivary effector genes. Based on the results, we created an extended catalogue of A. pisum salivary effector candidates. Next, we used aphid head samples, which contain SGs, to examine biotype-specific expression patterns of candidate salivary genes. RNA-seq analysis of head samples of alfalfa- and pea-specialized biotypes, each represented by three genetically distinct aphid lines reared on either a universal or specific host plant, showed that a majority of the candidate salivary effector genes was expressed in both biotypes at a similar level. Nonetheless, we identified small sets of genes that were differentially regulated in a biotype-specific manner. Little host plant effect (universal vs. specific) was observed on the expression of candidate salivary genes. Analysis of previously obtained genome re-sequenced data of the two biotypes revealed the copy number variations that might explain the differential expression of some candidate salivary genes. In addition, at least four candidate effector genes that were present in the alfalfa biotype but might not be encoded in the pea biotype were identified. This work sets the stage for future functional characterization of candidate genes potentially involved in the determination of plant specificity of pea aphid biotypes.

Fast Evolution and Lineage-Specific Gene Family Expansions of Aphid Salivary Effectors Driven by Interactions with Host-Plants.

Effector proteins play crucial roles in plant-parasite interactions by suppressing plant defenses and hijacking plant physiological responses to facilitate parasite invasion and propagation. Although effector proteins have been characterized in many microbial plant pathogens, their nature and role in adaptation to host plants are largely unknown in insect herbivores. Aphids rely on salivary effector proteins injected into the host plants to promote phloem sap uptake. Therefore, gaining insight into the repertoire and evolution of aphid effectors is key to unveiling the mechanisms responsible for aphid virulence and host plant specialization. With this aim in mind, we assembled catalogues of putative effectors in the legume specialist aphid, Acyrthosiphon pisum, using transcriptomics and proteomics approaches. We identified 3,603 candidate effector genes predicted to be expressed in A. pisum salivary glands (SGs), and 740 of which displayed up-regulated expression in SGs in comparison to the alimentary tract. A search for orthologs in 17 arthropod genomes revealed that SG-up-regulated effector candidates of A. pisum are enriched in aphid-specific genes and tend to evolve faster compared with the whole gene set. We also found that a large fraction of proteins detected in the A. pisum saliva belonged to three gene families, of which certain members show evidence consistent with positive selection. Overall, this comprehensive analysis suggests that the large repertoire of effector candidates in A. pisum constitutes a source of novelties promoting plant adaptation to legumes.

Two ancestral genes shaped the Xanthomonas campestris TAL effector gene repertoire.

Xanthomonas transcription activator-like effectors (TALEs) are injected inside plant cells to promote host susceptibility by enhancing transcription of host susceptibility genes. TALE-encoding (tal) genes were thought to be absent from Brassicaceae-infecting Xanthomonas campestris (Xc) genomes based on four reference genomic sequences. We discovered tal genes in 26 of 49 Xc strains isolated worldwide and used a combination of single molecule real time (SMRT) and tal amplicon sequencing to yield a near-complete description of the TALEs found in Xc (Xc TALome). The 53 sequenced tal genes encode 21 distinct DNA binding domains that sort into seven major DNA binding specificities. In silico analysis of the Brassica rapa promoterome identified a repertoire of predicted TALE targets, five of which were experimentally validated using quantitative reverse transcription polymerase chain reaction. The Xc TALome shows multiple signs of DNA rearrangements that probably drove its evolution from two ancestral tal genes. We discovered that Tal12a and Tal15a of Xcc strain Xca5 contribute together in the development of disease symptoms on susceptible B. oleracea var. botrytis cv Clovis. This large and polymorphic repertoire of TALEs opens novel perspectives for elucidating TALE-mediated susceptibility of Brassicaceae to black rot disease and for understanding the molecular processes underlying TALE evolution.

Optimization of Agroinfiltration in Pisum sativum Provides a New Tool for Studying the Salivary Protein Functions in the Pea Aphid Complex.

Aphids are piercing-sucking insect pests and feed on phloem sap. During feeding, aphids inject a battery of salivary proteins into host plant. Some of these proteins function like effectors of microbial pathogens and influence the outcome of plant-aphid interactions. The pea aphid (Acyrthosiphon pisum) is the model aphid and encompasses multiple biotypes each specialized to one or a few legume species, providing an opportunity to investigate the underlying mechanisms of the compatibility between plants and aphid biotypes. We aim to identify the aphid factors that determine the compatibility with host plants, hence involved in the host plant specialization process, and hypothesize that salivary proteins are one of those factors. Agrobacterium-mediated transient gene expression is a powerful tool to perform functional analyses of effector (salivary) proteins in plants. However, the tool was not established for the legume species that A. pisum feeds on. Thus, we decided to optimize the method for legume plants to facilitate the functional analyses of A. pisum salivary proteins. We screened a range of cultivars of pea (Pisum sativum) and alfalfa (Medicago sativa). None of the M. sativa cultivars was suitable for agroinfiltration under the tested conditions; however, we established a protocol for efficient transient gene expression in two cultivars of P. sativum, ZP1109 and ZP1130, using A. tumefaciens AGL-1 strain and the pEAQ-HT-DEST1 vector. We confirmed that the genes are expressed from 3 to 10 days post-infiltration and that aphid lines of the pea adapted biotype fed and reproduced on these two cultivars while lines of alfalfa and clover biotypes did not. Thus, the pea biotype recognizes these two cultivars as typical pea plants. By using a combination of ZP1109 and an A. pisum line, we defined an agroinfiltration procedure to examine the effect of in planta expression of selected salivary proteins on A. pisum fitness and demonstrated that transient expression of one candidate salivary gene increased the fecundity of the aphids. This result confirms that the agroinfiltration can be used to perform functional analyses of salivary proteins in P. sativum and consequently to study the molecular mechanisms underlying host specialization in the pea aphid complex.

Corrigendum: Optimization of Agroinfiltration in Pisum sativum Provides a New Tool for Studying the Salivary Protein Functions in the Pea Aphid Complex.

[This corrects the article on p. 1171 in vol. 7, PMID: 27555856.].

Genomics and transcriptomics of Xanthomonas campestris species challenge the concept of core type III effectome.

BACKGROUND: The bacterial species Xanthomonas campestris infects a wide range of Brassicaceae. Specific pathovars of this species cause black rot (pv. campestris), bacterial blight of stock (pv. incanae) or bacterial leaf spot (pv. raphani). RESULTS: In this study, we extended the genomic coverage of the species by sequencing and annotating the genomes of strains from pathovar incanae (CFBP 1606R and CFBP 2527R), pathovar raphani (CFBP 5828R) and a pathovar formerly named barbareae (CFBP 5825R). While comparative analyses identified a large core ORFeome at the species level, the core type III effectome was limited to only three putative type III effectors (XopP, XopF1 and XopAL1). In Xanthomonas, these effector proteins are injected inside the plant cells by the type III secretion system and contribute collectively to virulence. A deep and strand-specific RNA sequencing strategy was adopted in order to experimentally refine genome annotation for strain CFBP 5828R. This approach also allowed the experimental definition of novel ORFs and non-coding RNA transcripts. Using a constitutively active allele of hrpG, a master regulator of the type III secretion system, a HrpG-dependent regulon of 141 genes co-regulated with the type III secretion system was identified. Importantly, all these genes but seven are positively regulated by HrpG and 56 of those encode components of the Hrp type III secretion system and putative effector proteins. CONCLUSIONS: This dataset is an important resource to mine for novel type III effector proteins as well as for bacterial genes which could contribute to pathogenicity of X. campestris.

The Decoy Substrate of a Pathogen Effector and a Pseudokinase Specify Pathogen-Induced Modified-Self Recognition and Immunity in Plants.

In plants, host response to pathogenic microbes is driven both by microbial perception and detection of modified-self. The Xanthomonas campestris effector protein AvrAC/XopAC uridylylates the Arabidopsis BIK1 kinase to dampen basal resistance and thereby promotes bacterial virulence. Here we show that PBL2, a paralog of BIK1, is similarly uridylylated by AvrAC. However, in contrast to BIK1, PBL2 uridylylation is specifically required for host recognition of AvrAC to trigger immunity, but not AvrAC virulence. PBL2 thus acts as a decoy and enables AvrAC detection. AvrAC recognition also requires the RKS1 pseudokinase of the ZRK family and the NOD-like receptor ZAR1, which is known to recognize the Pseudomonas syringae effector HopZ1a. ZAR1 forms a stable complex with RKS1, which specifically recruits PBL2 when the latter is uridylylated by AvrAC, triggering ZAR1-mediated immunity. The results illustrate how decoy substrates and pseudokinases can specify and expand the capacity of the plant immune system.

Genomics of adaptation to host-plants in herbivorous insects.

Herbivorous insects represent the most species-rich lineages of metazoans. The high rate of diversification in herbivorous insects is thought to result from their specialization to distinct host-plants, which creates conditions favorable for the build-up of reproductive isolation and speciation. These conditions rely on constraints against the optimal use of a wide range of plant species, as each must constitute a viable food resource, oviposition site and mating site for an insect. Utilization of plants involves many essential traits of herbivorous insects, as they locate and select their hosts, overcome their defenses and acquire nutrients while avoiding intoxication. Although advances in understanding insect-plant molecular interactions have been limited by the complexity of insect traits involved in host use and the lack of genomic resources and functional tools, recent studies at the molecular level, combined with large-scale genomics studies at population and species levels, are revealing the genetic underpinning of plant specialization and adaptive divergence in non-model insect herbivores. Here, we review the recent advances in the genomics of plant adaptation in hemipterans and lepidopterans, two major insect orders, each of which includes a large number of crop pests. We focus on how genomics and post-genomics have improved our understanding of the mechanisms involved in insect-plant interactions by reviewing recent molecular discoveries in sensing, feeding, digesting and detoxifying strategies. We also present the outcomes of large-scale genomics approaches aimed at identifying loci potentially involved in plant adaptation in these insects.

Genome sequence of Xanthomonas fuscans subsp. fuscans strain 4834-R reveals that flagellar motility is not a general feature of xanthomonads.

BACKGROUND: Xanthomonads are plant-associated bacteria responsible for diseases on economically important crops. Xanthomonas fuscans subsp. fuscans (Xff) is one of the causal agents of common bacterial blight of bean. In this study, the complete genome sequence of strain Xff 4834-R was determined and compared to other Xanthomonas genome sequences. RESULTS: Comparative genomics analyses revealed core characteristics shared between Xff 4834-R and other xanthomonads including chemotaxis elements, two-component systems, TonB-dependent transporters, secretion systems (from T1SS to T6SS) and multiple effectors. For instance a repertoire of 29 Type 3 Effectors (T3Es) with two Transcription Activator-Like Effectors was predicted. Mobile elements were associated with major modifications in the genome structure and gene content in comparison to other Xanthomonas genomes. Notably, a deletion of 33 kbp affects flagellum biosynthesis in Xff 4834-R. The presence of a complete flagellar cluster was assessed in a collection of more than 300 strains representing different species and pathovars of Xanthomonas. Five percent of the tested strains presented a deletion in the flagellar cluster and were non-motile. Moreover, half of the Xff strains isolated from the same epidemic than 4834-R was non-motile and this ratio was conserved in the strains colonizing the next bean seed generations. CONCLUSIONS: This work describes the first genome of a Xanthomonas strain pathogenic on bean and reports the existence of non-motile xanthomonads belonging to different species and pathovars. Isolation of such Xff variants from a natural epidemic may suggest that flagellar motility is not a key function for in planta fitness.

xopAC-triggered immunity against Xanthomonas depends on Arabidopsis receptor-like cytoplasmic kinase genes PBL2 and RIPK.

Xanthomonas campestris pv. campestris (Xcc) colonizes the vascular system of Brassicaceae and ultimately causes black rot. In susceptible Arabidopsis plants, XopAC type III effector inhibits by uridylylation positive regulators of the PAMP-triggered immunity such as the receptor-like cytoplasmic kinases (RLCK) BIK1 and PBL1. In the resistant ecotype Col-0, xopAC is a major avirulence gene of Xcc. In this study, we show that both the RLCK interaction domain and the uridylyl transferase domain of XopAC are required for avirulence. Furthermore, xopAC can also confer avirulence to both the vascular pathogen Ralstonia solanacearum and the mesophyll-colonizing pathogen Pseudomonas syringae indicating that xopAC-specified effector-triggered immunity is not specific to the vascular system. In planta, XopAC-YFP fusions are localized at the plasma membrane suggesting that XopAC might interact with membrane-localized proteins. Eight RLCK of subfamily VII predicted to be localized at the plasma membrane and interacting with XopAC in yeast two-hybrid assays have been isolated. Within this subfamily, PBL2 and RIPK RLCK genes but not BIK1 are important for xopAC-specified effector-triggered immunity and Arabidopsis resistance to Xcc.

Natural genetic variation of Xanthomonas campestris pv. campestris pathogenicity on arabidopsis revealed by association and reverse genetics.

ABSTRACT The pathogenic bacterium Xanthomonas campestris pv. campestris, the causal agent of black rot of Brassicaceae, manipulates the physiology and the innate immunity of its hosts. Association genetic and reverse-genetic analyses of a world panel of 45 X. campestris pv. campestris strains were used to gain understanding of the genetic basis of the bacterium's pathogenicity to Arabidopsis thaliana. We found that the compositions of the minimal predicted type III secretome varied extensively, with 18 to 28 proteins per strain. There were clear differences in aggressiveness of those X. campestris pv. campestris strains on two Arabidopsis natural accessions. We identified 3 effector genes (xopAC, xopJ5, and xopAL2) and 67 amplified fragment length polymorphism (AFLP) markers that were associated with variations in disease symptoms. The nature and distribution of the AFLP markers remain to be determined, but we observed a low linkage disequilibrium level between predicted effectors and other significant markers, suggesting that additional genetic factors make a meaningful contribution to pathogenicity. Mutagenesis of type III effectors in X. campestris pv. campestris confirmed that xopAC functions as both a virulence and an avirulence gene in Arabidopsis and that xopAM functions as a second avirulence gene on plants of the Col-0 ecotype. However, we did not detect the effect of any other effector in the X. campestris pv. campestris 8004 strain, likely due to other genetic background effects. These results highlight the complex genetic basis of pathogenicity at the pathovar level and encourage us to challenge the agronomical relevance of some virulence determinants identified solely in model strains. IMPORTANCE The identification and understanding of the genetic determinants of bacterial virulence are essential to be able to design efficient protection strategies for infected plants. The recent availability of genomic resources for a limited number of pathogen isolates and host genotypes has strongly biased our research toward genotype-specific approaches. Indeed, these do not consider the natural variation in both pathogens and hosts, so their applied relevance should be challenged. In our study, we exploited the genetic diversity of Xanthomonas campestris pv. campestris, the causal agent of black rot on Brassicaceae (e.g., cabbage), to mine for pathogenicity determinants. This work evidenced the contribution of known and unknown loci to pathogenicity relevant at the pathovar level and identified these virulence determinants as prime targets for breeding resistance to X. campestris pv. campestris in Brassicaceae.

Genome Sequence of Xanthomonas campestris pv. campestris Strain Xca5.

An annotated high-quality draft genome sequence for Xanthomonas campestris pv. campestris race 1 strain Xca5 (originally described as X. campestris pv. armoraciae), the causal agent of black rot on Brassicaceae plants, has been determined. This genome sequence is a valuable resource for comparative genomics within the campestris pathovar.