CRISPR/Cas9-generated mutations in a sugar transporter gene reduce cassava susceptibility to bacterial blight.

Bacteria from the genus Xanthomonas are prolific phytopathogens that elicit disease in over 400 plant species. Xanthomonads carry a repertoire of specialized proteins called transcription activator-like (TAL) effectors that promote disease and pathogen virulence by inducing expression of host susceptibility (S) genes. Xanthomonas phaseoli pv. manihotis (Xpm) causes bacterial blight on the staple food crop cassava (Manihot esculenta Crantz). The Xpm effector TAL20 induces ectopic expression of the S gene Manihot esculenta Sugars Will Eventually be Exported Transporter 10a (MeSWEET10a), which encodes a sugar transporter that contributes to cassava bacterial blight susceptibility. We used CRISPR/Cas9 to generate multiple cassava lines with edits to the MeSWEET10a TAL20 effector binding site and/or coding sequence. In several of the regenerated lines, MeSWEET10a expression was no longer induced by Xpm, and in these cases, we observed reduced cassava bacterial blight (CBB) disease symptoms post Xpm infection. Because MeSWEET10a is expressed in cassava flowers, we further characterized the reproductive capability of the MeSWEET10a promoter and coding sequence mutants. Lines were crossed to themselves and to wild-type plants. The results indicated that expression of MeSWEET10a in female, but not male, flowers, is critical to produce viable F1 seed. In the case of promoter mutations that left the coding sequence intact, viable F1 progeny were recovered. Taken together, these results demonstrate that blocking MeSWEET10a induction is a viable strategy for decreasing cassava susceptibility to CBB and that ideal lines will contain promoter mutations that block TAL effector binding while leaving endogenous expression of MeSWEET10a unaltered.

Escalation in the host-pathogen arms race: A host resistance response corresponds to a heightened bacterial virulence response.

The zig-zag model of host-pathogen interaction describes the relative strength of defense response across a spectrum of pathogen-induced plant phenotypes. A stronger defense response results in increased resistance. Here, we investigate the strength of pathogen virulence during disease and place these findings in the context of the zig-zag model. Xanthomonas vasicola pv. holcicola (Xvh) causes sorghum bacterial leaf streak. Despite being widespread, this disease has not been described in detail at the molecular level. We divided diverse sorghum genotypes into three groups based on disease symptoms: water-soaked lesions, red lesions, and resistance. Bacterial growth assays confirmed that these three phenotypes represent a range of resistance and susceptibility. To simultaneously reveal defense and virulence responses across the spectrum of disease phenotypes, we performed dual RNA-seq on Xvh-infected sorghum. Consistent with the zig-zag model, the expression of plant defense-related genes was strongest in the resistance interaction. Surprisingly, bacterial virulence genes related to the type III secretion system (T3SS) and type III effectors (T3Es) were also most highly expressed in the resistance interaction. This expression pattern was observed at multiple time points within the sorghum-Xvh pathosystem. Further, a similar expression pattern was observed in Arabidopsis infected with Pseudomonas syringae for effector-triggered immunity via AvrRps4 but not AvrRpt2. Specific metabolites were able to repress the Xvh virulence response in vitro and in planta suggesting a possible signaling mechanism. Taken together, these findings reveal multiple permutations of the continually evolving host-pathogen arms race from the perspective of host defense and pathogen virulence responses.

Large structural variations in the haplotype-resolved African cassava genome.

Cassava (Manihot esculenta Crantz, 2n = 36) is a global food security crop. It has a highly heterozygous genome, high genetic load, and genotype-dependent asynchronous flowering. It is typically propagated by stem cuttings and any genetic variation between haplotypes, including large structural variations, is preserved by such clonal propagation. Traditional genome assembly approaches generate a collapsed haplotype representation of the genome. In highly heterozygous plants, this results in artifacts and an oversimplification of heterozygous regions. We used a combination of Pacific Biosciences (PacBio), Illumina, and Hi-C to resolve each haplotype of the genome of a farmer-preferred cassava line, TME7 (Oko-iyawo). PacBio reads were assembled using the FALCON suite. Phase switch errors were corrected using FALCON-Phase and Hi-C read data. The ultralong-range information from Hi-C sequencing was also used for scaffolding. Comparison of the two phases revealed >5000 large haplotype-specific structural variants affecting over 8 Mb, including insertions and deletions spanning thousands of base pairs. The potential of these variants to affect allele-specific expression was further explored. RNA-sequencing data from 11 different tissue types were mapped against the scaffolded haploid assembly and gene expression data are incorporated into our existing easy-to-use web-based interface to facilitate use by the broader plant science community. These two assemblies provide an excellent means to study the effects of heterozygosity, haplotype-specific structural variation, gene hemizygosity, and allele-specific gene expression contributing to important agricultural traits and further our understanding of the genetics and domestication of cassava.

Current status and impending progress for cassava structural genomics.

KEY MESSAGE: We demystify recent advances in genome assemblies for the heterozygous staple crop cassava (Manihot esculenta), and highlight key cassava genomic resources. Cassava, Manihot esculenta Crantz, is a crop of societal and agricultural importance in tropical regions around the world. Genomics provides a platform for accelerated improvement of cassava's nutritional and agronomic traits, as well as for illuminating aspects of cassava's history including its path towards domestication. The highly heterozygous nature of the cassava genome is widely recognized. However, the full extent and context of this heterozygosity has been difficult to reveal because of technological limitations within genome sequencing. Only recently, with several new long-read sequencing technologies coming online, has the genomics community been able to tackle some similarly difficult genomes. In light of these recent advances, we provide this review to document the current status of the cassava genome and genomic resources and provide a perspective on what to look forward to in the coming years.

Identification and genomic characterization of major effect bacterial blight resistance locus (BB-13) in Upland cotton (Gossypium hirsutum L.).

KEY MESSAGE: Identification and genomic characterization of major resistance locus against cotton bacterial blight (CBB) using GWAS and linkage mapping to enable genomics-based development of durable CBB resistance and gene discovery in cotton. Cotton bacterial leaf blight (CBB), caused by Xanthomonas citri subsp. malvacearum (Xcm), has periodically been a damaging disease in the USA. Identification and deployment of genetic resistance in cotton cultivars is the most economical and efficient means of reducing crop losses due to CBB. In the current study, genome-wide association study (GWAS) of CBB resistance using an elite diversity panel of 380 accessions, genotyped with the cotton single nucleotide polymorphism (SNP) 63 K array, and phenotyped with race-18 of CBB, localized the CBB resistance to a 2.01-Mb region in the long arm of chromosome D02. Molecular genetic mapping using an F6 recombinant inbred line (RIL) population showed the CBB resistance in cultivar Arkot 8102 was controlled by a single locus (BB-13). The BB-13 locus was mapped within the 0.95-cM interval near the telomeric region in the long arm of chromosome D02. Flanking SNP markers, i04890Gh and i04907Gh of the BB-13 locus, identified from the combined linkage analysis and GWAS, targeted it to a 371-Kb genomic region. Candidate gene analysis identified thirty putative gene sequences in the targeted genomic region. Nine of these putative genes and two NBS-LRR genes adjacent to the targeted region were putatively involved in plant disease resistance and are possible candidate genes for BB-13 locus. Genetic mapping and genomic targeting of the BB13 locus in the current study will help in cloning the CBB-resistant gene and establishing the molecular genetic architecture of the BB-13 locus towards developing durable resistance to CBB in cotton.

Genomics-enabled analysis of the emergent disease cotton bacterial blight.

Cotton bacterial blight (CBB), an important disease of (Gossypium hirsutum) in the early 20th century, had been controlled by resistant germplasm for over half a century. Recently, CBB re-emerged as an agronomic problem in the United States. Here, we report analysis of cotton variety planting statistics that indicate a steady increase in the percentage of susceptible cotton varieties grown each year since 2009. Phylogenetic analysis revealed that strains from the current outbreak cluster with race 18 Xanthomonas citri pv. malvacearum (Xcm) strains. Illumina based draft genomes were generated for thirteen Xcm isolates and analyzed along with 4 previously published Xcm genomes. These genomes encode 24 conserved and nine variable type three effectors. Strains in the race 18 clade contain 3 to 5 more effectors than other Xcm strains. SMRT sequencing of two geographically and temporally diverse strains of Xcm yielded circular chromosomes and accompanying plasmids. These genomes encode eight and thirteen distinct transcription activator-like effector genes. RNA-sequencing revealed 52 genes induced within two cotton cultivars by both tested Xcm strains. This gene list includes a homeologous pair of genes, with homology to the known susceptibility gene, MLO. In contrast, the two strains of Xcm induce different clade III SWEET sugar transporters. Subsequent genome wide analysis revealed patterns in the overall expression of homeologous gene pairs in cotton after inoculation by Xcm. These data reveal important insights into the Xcm-G. hirsutum disease complex and strategies for future development of resistant cultivars.

Simultaneous CRISPR/Cas9-mediated editing of cassava eIF4E isoforms nCBP-1 and nCBP-2 reduces cassava brown streak disease symptom severity and incidence.

Cassava brown streak disease (CBSD) is a major constraint on cassava yields in East and Central Africa and threatens production in West Africa. CBSD is caused by two species of positive-sense RNA viruses belonging to the family Potyviridae, genus Ipomovirus: Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). Diseases caused by the family Potyviridae require the interaction of viral genome-linked protein (VPg) and host eukaryotic translation initiation factor 4E (eIF4E) isoforms. Cassava encodes five eIF4E proteins: eIF4E, eIF(iso)4E-1, eIF(iso)4E-2, novel cap-binding protein-1 (nCBP-1), and nCBP-2. Protein-protein interaction experiments consistently found that VPg proteins associate with cassava nCBPs. CRISPR/Cas9-mediated genome editing was employed to generate ncbp-1, ncbp-2, and ncbp-1/ncbp-2 mutants in cassava cultivar 60444. Challenge with CBSV showed that ncbp-1/ncbp-2 mutants displayed delayed and attenuated CBSD aerial symptoms, as well as reduced severity and incidence of storage root necrosis. Suppressed disease symptoms were correlated with reduced virus titre in storage roots relative to wild-type controls. Our results demonstrate the ability to modify multiple genes simultaneously in cassava to achieve tolerance to CBSD. Future studies will investigate the contribution of remaining eIF4E isoforms on CBSD and translate this knowledge into an optimized strategy for protecting cassava from disease.

Correction: A Genetic Screen Identifies a Requirement for Cysteine-Rich-Receptor-Like Kinases in Rice NH1 (OsNPR1)-Mediated Immunity.

[This corrects the article DOI: 10.1371/journal.pgen.1006049.].

A Genetic Screen Identifies a Requirement for Cysteine-Rich-Receptor-Like Kinases in Rice NH1 (OsNPR1)-Mediated Immunity.

Systemic acquired resistance, mediated by the Arabidopsis NPR1 gene and the rice NH1 gene, confers broad-spectrum immunity to diverse pathogens. NPR1 and NH1 interact with TGA transcription factors to activate downstream defense genes. Despite the importance of this defense response, the signaling components downstream of NPR1/NH1 and TGA proteins are poorly defined. Here we report the identification of a rice mutant, snim1, which suppresses NH1-mediated immunity and demonstrate that two genes encoding previously uncharacterized cysteine-rich-receptor-like kinases (CRK6 and CRK10), complement the snim1 mutant phenotype. Silencing of CRK6 and CRK10 genes individually in the parental genetic background recreates the snim1 phenotype. We identified a rice mutant in the Kitaake genetic background with a frameshift mutation in crk10; this mutant also displays a compromised immune response highlighting the important role of crk10. We also show that elevated levels of NH1 expression lead to enhanced CRK10 expression and that the rice TGA2.1 protein binds to the CRK10 promoter. These experiments demonstrate a requirement for CRKs in NH1-mediated immunity and establish a molecular link between NH1 and induction of CRK10 expression.

Allele exchange at the EPSPS locus confers glyphosate tolerance in cassava.

Effective weed control can protect yields of cassava (Manihot esculenta) storage roots. Farmers could benefit from using herbicide with a tolerant cultivar. We applied traditional transgenesis and gene editing to generate robust glyphosate tolerance in cassava. By comparing promoters regulating expression of transformed 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) genes with various paired amino acid substitutions, we found that strong constitutive expression is required to achieve glyphosate tolerance during in vitro selection and in whole cassava plants. Using strategies that exploit homologous recombination (HR) and nonhomologous end-joining (NHEJ) DNA repair pathways, we precisely introduced the best-performing allele into the cassava genome, simultaneously creating a promoter swap and dual amino acid substitutions at the endogenous EPSPS locus. Primary EPSPS-edited plants were phenotypically normal, tolerant to high doses of glyphosate, with some free of detectable T-DNA integrations. Our methods demonstrate an editing strategy for creating glyphosate tolerance in crop plants and demonstrate the potential of gene editing for further improvement of cassava.

CG gene body DNA methylation changes and evolution of duplicated genes in cassava.

DNA methylation is important for the regulation of gene expression and the silencing of transposons in plants. Here we present genome-wide methylation patterns at single-base pair resolution for cassava (Manihot esculenta, cultivar TME 7), a crop with a substantial impact in the agriculture of subtropical and tropical regions. On average, DNA methylation levels were higher in all three DNA sequence contexts (CG, CHG, and CHH, where H equals A, T, or C) than those of the most well-studied model plant Arabidopsis thaliana. As in other plants, DNA methylation was found both on transposons and in the transcribed regions (bodies) of many genes. Consistent with these patterns, at least one cassava gene copy of all of the known components of Arabidopsis DNA methylation pathways was identified. Methylation of LTR transposons (GYPSY and COPIA) was found to be unusually high compared with other types of transposons, suggesting that the control of the activity of these two types of transposons may be especially important. Analysis of duplicated gene pairs resulting from whole-genome duplication showed that gene body DNA methylation and gene expression levels have coevolved over short evolutionary time scales, reinforcing the positive relationship between gene body methylation and high levels of gene expression. Duplicated genes with the most divergent gene body methylation and expression patterns were found to have distinct biological functions and may have been under natural or human selection for cassava traits.

Foundational and Translational Research Opportunities to Improve Plant Health.

Reader Comments | Submit a Comment The white paper reports the deliberations of a workshop focused on biotic challenges to plant health held in Washington, D.C. in September 2016. Ensuring health of food plants is critical to maintaining the quality and productivity of crops and for sustenance of the rapidly growing human population. There is a close linkage between food security and societal stability; however, global food security is threatened by the vulnerability of our agricultural systems to numerous pests, pathogens, weeds, and environmental stresses. These threats are aggravated by climate change, the globalization of agriculture, and an over-reliance on nonsustainable inputs. New analytical and computational technologies are providing unprecedented resolution at a variety of molecular, cellular, organismal, and population scales for crop plants as well as pathogens, pests, beneficial microbes, and weeds. It is now possible to both characterize useful or deleterious variation as well as precisely manipulate it. Data-driven, informed decisions based on knowledge of the variation of biotic challenges and of natural and synthetic variation in crop plants will enable deployment of durable interventions throughout the world. These should be integral, dynamic components of agricultural strategies for sustainable agriculture.

Gene expression atlas for the food security crop cassava.

Cassava (Manihot esculenta) feeds c. 800 million people world-wide. Although this crop displays high productivity under drought and poor soil conditions, it is susceptible to disease, postharvest deterioration and the roots contain low nutritional content. Here, we provide molecular identities for 11 cassava tissue/organ types through RNA-sequencing and develop an open access, web-based interface for further interrogation of the data. Through this dataset, we consider the physiology of cassava. Specifically, we focus on identification of the transcriptional signatures that define the massive, underground storage roots used as a food source and the favored target tissue for transgene integration and genome editing, friable embryogenic callus (FEC). Further, we identify promoters able to drive strong expression in multiple tissue/organs. The information gained from this study is of value for both conventional and biotechnological improvement programs.

Quantitative, Image-Based Phenotyping Methods Provide Insight into Spatial and Temporal Dimensions of Plant Disease.

Plant disease symptoms exhibit complex spatial and temporal patterns that are challenging to quantify. Image-based phenotyping approaches enable multidimensional characterization of host-microbe interactions and are well suited to capture spatial and temporal data that are key to understanding disease progression. We applied image-based methods to investigate cassava bacterial blight, which is caused by the pathogen Xanthomonas axonopodis pv. manihotis (Xam). We generated Xam strains in which individual predicted type III effector (T3E) genes were mutated and applied multiple imaging approaches to investigate the role of these proteins in bacterial virulence. Specifically, we quantified bacterial populations, water-soaking disease symptoms, and pathogen spread from the site of inoculation over time for strains with mutations in avrBs2, xopX, and xopK as compared to wild-type Xam ∆avrBs2 and ∆xopX both showed reduced growth in planta and delayed spread through the vasculature system of cassava. ∆avrBs2 exhibited reduced water-soaking symptoms at the site of inoculation. In contrast, ∆xopK exhibited enhanced induction of disease symptoms at the site of inoculation but reduced spread through the vasculature. Our results highlight the importance of adopting a multipronged approach to plant disease phenotyping to more fully understand the roles of T3Es in virulence. Finally, we demonstrate that the approaches used in this study can be extended to many host-microbe systems and increase the dimensions of phenotype that can be explored.

Multilocus sequence analysis of xanthomonads causing bacterial spot of tomato and pepper plants reveals strains generated by recombination among species and recent global spread of Xanthomonas gardneri.

Four Xanthomonas species are known to cause bacterial spot of tomato and pepper, but the global distribution and genetic diversity of these species are not well understood. A collection of bacterial spot-causing strains from the Americas, Africa, Southeast Asia, and New Zealand were characterized for genetic diversity and phylogenetic relationships using multilocus sequence analysis of six housekeeping genes. By examining strains from different continents, we found unexpected phylogeographic patterns, including the global distribution of a single multilocus haplotype of X. gardneri, possible regional differentiation in X. vesicatoria, and high species diversity on tomato in Africa. In addition, we found evidence of multiple recombination events between X. euvesicatoria and X. perforans. Our results indicate that there have been shifts in the species composition of bacterial spot pathogen populations due to the global spread of dominant genotypes and that recombination between species has generated genetic diversity in these populations.

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.

Xanthomonas axonopodis virulence is promoted by a transcription activator-like effector-mediated induction of a SWEET sugar transporter in cassava.

The gene-for-gene concept has historically been applied to describe a specific resistance interaction wherein single genes from the host and the pathogen dictate the outcome. These interactions have been observed across the plant kingdom and all known plant microbial pathogens. In recent years, this concept has been extended to susceptibility phenotypes in the context of transcription activator-like (TAL) effectors that target SWEET sugar transporters. However, because this interaction has only been observed in rice, it was not clear whether the gene-for-gene susceptibility was unique to that system. Here, we show, through a combined systematic analysis of the TAL effector complement of Xanthomonas axonopodis pv. manihotis and RNA sequencing to identify targets in cassava, that TAL20Xam668 specifically induces the sugar transporter MeSWEET10a to promote virulence. Designer TAL effectors (dTALE) complement TAL20Xam668 mutant phenotypes, demonstrating that MeSWEET10a is a susceptibility gene in cassava. Sucrose uptake-deficient X. axonopodis pv. manihotis bacteria do not lose virulence, indicating that sucrose may be cleaved extracellularly and taken up as hexoses into X. axonopodis pv. manihotis. Together, our data suggest that pathogen hijacking of plant nutrients is not unique to rice blight but also plays a role in bacterial blight of the dicot cassava.