Ksenia V. Krasileva

Adaptive Evolution Has Targeted the C-Terminal Domain of the RXLR Effectors of Plant Pathogenic Oomycetes

Oomycete plant pathogens deliver effector proteins inside host cells to modulate plant defense circuitry and to enable parasitic colonization. These effectors are defined by a conserved motif, termed RXLR (for Arg, any amino acid, Leu, Arg), that is located downstream of the signal peptide and that has been implicated in host translocation. Because the phenotypes of RXLR effectors extend to plant cells, their genes are expected to be the direct target of the evolutionary forces that drive the antagonistic interplay between pathogen and host. We used the draft genome sequences of three oomycete plant pathogens, Phytophthora sojae, Phytophthora ramorum, and Hyaloperonospora parasitica, to generate genome-wide catalogs of RXLR effector genes and determine the extent to which these genes are under positive selection. These analyses revealed that the RXLR sequence is overrepresented and positionally constrained in the secretome of Phytophthora relative to other eukaryotes. The three examined plant pathogenic oomycetes carry complex and diverse sets of RXLR effector genes that have undergone relatively rapid birth and death evolution. We obtained robust evidence of positive selection in more than two-thirds of the examined paralog families of RXLR effectors. Positive selection has acted for the most part on the C-terminal region, consistent with the view that RXLR effectors are modular, with the N terminus involved in secretion and host translocation and the C-terminal domain dedicated to modulating host defenses inside plant cells.

Activation of an Arabidopsis Resistance Protein Is Specified by the in Planta Association of Its Leucine-Rich Repeat Domain with the Cognate Oomycete Effector

The Arabidopsis disease resistance protein RPP1 recognizes the ATR1 effector protein from Hyaloperonospora arabidopsidis. This works shows that the molecular basis of recognition is mediated by the in planta association of the LRR domain of the RPP1 protein with the ATR1 effector protein. The in planta association of specific alleles of ATR1 leads to the activation of plant immune responses. Activation of plant immunity relies on recognition of pathogen effectors by several classes of plant resistance proteins. To discover the underlying molecular mechanisms of effector recognition by the Arabidopsis thaliana RECOGNITION OF PERONOSPORA PARASITICA1 (RPP1) resistance protein, we adopted an Agrobacterium tumefaciens–mediated transient protein expression system in tobacco (Nicotiana tabacum), which allowed us to perform coimmunoprecipitation experiments and mutational analyses. Herein, we demonstrate that RPP1 associates with its cognate effector ARABIDOPSIS THALIANA RECOGNIZED1 (ATR1) in a recognition-specific manner and that this association is a prerequisite step in the induction of the hypersensitive cell death response of host tissue. The leucine-rich repeat (LRR) domain of RPP1 mediates the interaction with ATR1, while the Toll/Interleukin1 Receptor (TIR) domain facilitates the induction of the hypersensitive cell death response. Additionally, we demonstrate that mutations in the TIR and nucleotide binding site domains, which exhibit loss of function for the induction of the hypersensitive response, are still able to associate with the effector in planta. Thus, our data suggest molecular epistasis between signaling activity of the TIR domain and the recognition function of the LRR and allow us to propose a model for ATR1 recognition by RPP1.

Comparative genomics reveals diversity among xanthomonads infecting tomato and pepper

BackgroundBacterial spot of tomato and pepper is caused by four Xanthomonas species and is a major plant disease in warm humid climates. The four species are distinct from each other based on physiological and molecular characteristics. The genome sequence of strain 85-10, a member of one of the species, Xanthomonas euvesicatoria (Xcv) has been previously reported. To determine the relationship of the four species at the genome level and to investigate the molecular basis of their virulence and differing host ranges, draft genomic sequences of members of the other three species were determined and compared to strain 85-10.ResultsWe sequenced the genomes of X. vesicatoria (Xv) strain 1111 (ATCC 35937), X. perforans (Xp) strain 91-118 and X. gardneri (Xg) strain 101 (ATCC 19865). The genomes were compared with each other and with the previously sequenced Xcv strain 85-10. In addition, the molecular features were predicted that may be required for pathogenicity including the type III secretion apparatus, type III effectors, other secretion systems, quorum sensing systems, adhesins, extracellular polysaccharide, and lipopolysaccharide determinants. Several novel type III effectors from Xg strain 101 and Xv strain 1111 genomes were computationally identified and their translocation was validated using a reporter gene assay. A homolog to Ax21, the elicitor of XA21-mediated resistance in rice, and a functional Ax21 sulfation system were identified in Xcv. Genes encoding proteins with functions mediated by type II and type IV secretion systems have also been compared, including enzymes involved in cell wall deconstruction, as contributors to pathogenicity.ConclusionsComparative genomic analyses revealed considerable diversity among bacterial spot pathogens, providing new insights into differences and similarities that may explain the diverse nature of these strains. Genes specific to pepper pathogens, such as the O-antigen of the lipopolysaccharide cluster, and genes unique to individual strains, such as novel type III effectors and bacteriocin genes, have been identified providing new clues for our understanding of pathogen virulence, aggressiveness, and host preference. These analyses will aid in efforts towards breeding for broad and durable resistance in economically important tomato and pepper cultivars.

Efficient Genome-Wide Detection and Cataloging of EMS-Induced Mutations Using Exome Capture and Next-Generation Sequencing

The authors combined high-throughput global sequencing targeted to genes and custom-designed bioinformatics tools to catalog and characterize chemically induced mutations in rice and wheat. They demonstrate that this approach can be used to develop large-scale induced mutation resources with relatively small investments and is applicable to resource-limited and polyploid organisms. Chemical mutagenesis efficiently generates phenotypic variation in otherwise homogeneous genetic backgrounds, enabling functional analysis of genes. Advances in mutation detection have brought the utility of induced mutant populations on par with those produced by insertional mutagenesis, but systematic cataloguing of mutations would further increase their utility. We examined the suitability of multiplexed global exome capture and sequencing coupled with custom-developed bioinformatics tools to identify mutations in well-characterized mutant populations of rice (Oryza sativa) and wheat (Triticum aestivum). In rice, we identified ∼18,000 induced mutations from 72 independent M2 individuals. Functional evaluation indicated the recovery of potentially deleterious mutations for >2600 genes. We further observed that specific sequence and cytosine methylation patterns surrounding the targeted guanine residues strongly affect their probability to be alkylated by ethyl methanesulfonate. Application of these methods to six independent M2 lines of tetraploid wheat demonstrated that our bioinformatics pipeline is applicable to polyploids. In conclusion, we provide a method for developing large-scale induced mutation resources with relatively small investments that is applicable to resource-poor organisms. Furthermore, our results demonstrate that large libraries of sequenced mutations can be readily generated, providing enhanced opportunities to study gene function and assess the effect of sequence and chromatin context on mutations.

Separating homeologs by phasing in the tetraploid wheat transcriptome

BackgroundThe high level of identity among duplicated homoeologous genomes in tetraploid pasta wheat presents substantial challenges for de novo transcriptome assembly. To solve this problem, we develop a specialized bioinformatics workflow that optimizes transcriptome assembly and separation of merged homoeologs. To evaluate our strategy, we sequence and assemble the transcriptome of one of the diploid ancestors of pasta wheat, and compare both assemblies with a benchmark set of 13,472 full-length, non-redundant bread wheat cDNAs.ResultsA total of 489 million 100 bp paired-end reads from tetraploid wheat assemble in 140,118 contigs, including 96% of the benchmark cDNAs. We used a comparative genomics approach to annotate 66,633 open reading frames. The multiple k-mer assembly strategy increases the proportion of cDNAs assembled full-length in a single contig by 22% relative to the best single k-mer size. Homoeologs are separated using a post-assembly pipeline that includes polymorphism identification, phasing of SNPs, read sorting, and re-assembly of phased reads. Using a reference set of genes, we determine that 98.7% of SNPs analyzed are correctly separated by phasing.ConclusionsOur study shows that de novo transcriptome assembly of tetraploid wheat benefit from multiple k-mer assembly strategies more than diploid wheat. Our results also demonstrate that phasing approaches originally designed for heterozygous diploid organisms can be used to separate the close homoeologous genomes of tetraploid wheat. The predicted tetraploid wheat proteome and gene models provide a valuable tool for the wheat research community and for those interested in comparative genomic studies.

An improved assembly and annotation of the allohexaploid wheat genome identifies complete families of agronomic genes and provides genomic evidence for chromosomal translocations

Advances in genome sequencing and assembly technologies are generating many high-quality genome sequences, but assemblies of large, repeat-rich polyploid genomes, such as that of bread wheat, remain fragmented and incomplete. We have generated a new wheat whole-genome shotgun sequence assembly using a combination of optimized data types and an assembly algorithm designed to deal with large and complex genomes. The new assembly represents >78% of the genome with a scaffold N50 of 88.8 kb that has a high fidelity to the input data. Our new annotation combines strand-specific Illumina RNA-seq and Pacific Biosciences (PacBio) full-length cDNAs to identify 104,091 high-confidence protein-coding genes and 10,156 noncoding RNA genes. We confirmed three known and identified one novel genome rearrangements. Our approach enables the rapid and scalable assembly of wheat genomes, the identification of structural variants, and the definition of complete gene models, all powerful resources for trait analysis and breeding of this key global crop.

High-throughput genomic sequencing of cassava bacterial blight strains identifies conserved effectors to target for durable resistance

Cassava bacterial blight (CBB), incited by Xanthomonas axonopodis pv. manihotis (Xam), is the most important bacterial disease of cassava, a staple food source for millions of people in developing countries. Here we present a widely applicable strategy for elucidating the virulence components of a pathogen population. We report Illumina-based draft genomes for 65 Xam strains and deduce the phylogenetic relatedness of Xam across the areas where cassava is grown. Using an extensive database of effector proteins from animal and plant pathogens, we identify the effector repertoire for each sequenced strain and use a comparative sequence analysis to deduce the least polymorphic of the conserved effectors. These highly conserved effectors have been maintained over 11 countries, three continents, and 70 y of evolution and as such represent ideal targets for developing resistance strategies.

Uncovering hidden variation in polyploid wheat

Significance Pasta and bread wheat are polyploid species that carry multiple copies of each gene. Therefore, loss-of-function mutations in one gene copy are frequently masked by functional copies on other genomes. We sequenced the protein coding regions of 2,735 mutant lines and developed a public database including more than 10 million mutations. Researchers and breeders can search this database online, identify mutations in the different copies of their target gene, and request seeds to study gene function or improve wheat varieties. Mutations are being used to improve the nutritional value of wheat, increase the size of the wheat grains, and generate additional variability in flowering genes to improve wheat adaptation to new and changing environments. Comprehensive reverse genetic resources, which have been key to understanding gene function in diploid model organisms, are missing in many polyploid crops. Young polyploid species such as wheat, which was domesticated less than 10,000 y ago, have high levels of sequence identity among subgenomes that mask the effects of recessive alleles. Such redundancy reduces the probability of selection of favorable mutations during natural or human selection, but also allows wheat to tolerate high densities of induced mutations. Here we exploited this property to sequence and catalog more than 10 million mutations in the protein-coding regions of 2,735 mutant lines of tetraploid and hexaploid wheat. We detected, on average, 2,705 and 5,351 mutations per tetraploid and hexaploid line, respectively, which resulted in 35–40 mutations per kb in each population. With these mutation densities, we identified an average of 23–24 missense and truncation alleles per gene, with at least one truncation or deleterious missense mutation in more than 90% of the captured wheat genes per population. This public collection of mutant seed stocks and sequence data enables rapid identification of mutations in the different copies of the wheat genes, which can be combined to uncover previously hidden variation. Polyploidy is a central phenomenon in plant evolution, and many crop species have undergone recent genome duplication events. Therefore, the general strategy and methods developed herein can benefit other polyploid crops.

Sequence Divergent RXLR Effectors Share a Structural Fold Conserved across Plant Pathogenic Oomycete Species

The availability of genome sequences for some of the most devastating eukaryotic plant pathogens has led a revolution in our understanding of how these parasites cause disease, and how their hosts respond to invasion [1]–[7]. One of the most significant discoveries from the genome sequences of plant pathogenic oomycetes is the plethora of putative translocated effector proteins these organisms encode. Many effector genes display signatures of rapid evolution and tend to reside in dynamic regions of the pathogen genomes. Once inside the host, effector proteins modulate cellular processes, mainly suppressing plant immunity [8]–[12]. Effectors can also be recognized directly or indirectly by the plant immune system through the action of disease resistance (R) proteins [13], [14].

Comparative analysis of plant immune receptor architectures uncovers host proteins likely targeted by pathogens.

BackgroundPlants deploy immune receptors to detect pathogen-derived molecules and initiate defense responses. Intracellular plant immune receptors called nucleotide-binding leucine-rich repeat (NLR) proteins contain a central nucleotide-binding (NB) domain followed by a series of leucine-rich repeats (LRRs), and are key initiators of plant defense responses. However, recent studies demonstrated that NLRs with non-canonical domain architectures play an important role in plant immunity. These composite immune receptors are thought to arise from fusions between NLRs and additional domains that serve as “baits” for the pathogen-derived effector proteins, thus enabling pathogen recognition. Several names have been proposed to describe these proteins, including “integrated decoys” and “integrated sensors”. We adopt and argue for “integrated domains” or NLR-IDs, which describes the product of the fusion without assigning a universal mode of action.ResultsWe have scanned available plant genome sequences for the full spectrum of NLR-IDs to evaluate the diversity of integrations of potential sensor/decoy domains across flowering plants, including 19 crop species. We manually curated wheat and brassicas and experimentally validated a subset of NLR-IDs in wild and cultivated wheat varieties. We have examined NLR fusions that occur in multiple plant families and identified that some domains show re-occurring integration across lineages. Domains fused to NLRs overlap with previously identified pathogen targets confirming that they act as baits for the pathogen. While some of the integrated domains have been previously implicated in disease resistance, others provide new targets for engineering durable resistance to plant pathogens.ConclusionsWe have built a robust reproducible pipeline for detecting variable domain architectures in plant immune receptors across species. We hypothesize that NLR-IDs that we revealed provide clues to the host proteins targeted by pathogens, and that this information can be deployed to discover new sources of disease resistance.