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.

Small RNA-Based Antiviral Defense in the Phytopathogenic Fungus Colletotrichum higginsianum.

Even though the fungal kingdom contains more than 3 million species, little is known about the biological roles of RNA silencing in fungi. The Colletotrichum genus comprises fungal species that are pathogenic for a wide range of crop species worldwide. To investigate the role of RNA silencing in the ascomycete fungus Colletotrichum higginsianum, knock-out mutants affecting genes for three RNA-dependent RNA polymerase (RDR), two Dicer-like (DCL), and two Argonaute (AGO) proteins were generated by targeted gene replacement. No effects were observed on vegetative growth for any mutant strain when grown on complex or minimal media. However, Δdcl1, Δdcl1Δdcl2 double mutant, and Δago1 strains showed severe defects in conidiation and conidia morphology. Total RNA transcripts and small RNA populations were analyzed in parental and mutant strains. The greatest effects on both RNA populations was observed in the Δdcl1, Δdcl1Δdcl2, and Δago1 strains, in which a previously uncharacterized dsRNA mycovirus [termed Colletotrichum higginsianum non-segmented dsRNA virus 1 (ChNRV1)] was derepressed. Phylogenetic analyses clearly showed a close relationship between ChNRV1 and members of the segmented Partitiviridae family, despite the non-segmented nature of the genome. Immunoprecipitation of small RNAs associated with AGO1 showed abundant loading of 5'U-containing viral siRNA. C. higginsianum parental and Δdcl1 mutant strains cured of ChNRV1 revealed that the conidiation and spore morphology defects were primarily caused by ChNRV1. Based on these results, RNA silencing involving ChDCL1 and ChAGO1 in C. higginsianum is proposed to function as an antiviral mechanism.

Roles and programming of Arabidopsis ARGONAUTE proteins during Turnip mosaic virus infection.

In eukaryotes, ARGONAUTE proteins (AGOs) associate with microRNAs (miRNAs), short interfering RNAs (siRNAs), and other classes of small RNAs to regulate target RNA or target loci. Viral infection in plants induces a potent and highly specific antiviral RNA silencing response characterized by the formation of virus-derived siRNAs. Arabidopsis thaliana has ten AGO genes of which AGO1, AGO2, and AGO7 have been shown to play roles in antiviral defense. A genetic analysis was used to identify and characterize the roles of AGO proteins in antiviral defense against Turnip mosaic virus (TuMV) in Arabidopsis. AGO1, AGO2 and AGO10 promoted anti-TuMV defense in a modular way in various organs, with AGO2 providing a prominent antiviral role in leaves. AGO5, AGO7 and AGO10 had minor effects in leaves. AGO1 and AGO10 had overlapping antiviral functions in inflorescence tissues after systemic movement of the virus, although the roles of AGO1 and AGO10 accounted for only a minor amount of the overall antiviral activity. By combining AGO protein immunoprecipitation with high-throughput sequencing of associated small RNAs, AGO2, AGO10, and to a lesser extent AGO1 were shown to associate with siRNAs derived from silencing suppressor (HC-Pro)-deficient TuMV-AS9, but not with siRNAs derived from wild-type TuMV. Co-immunoprecipitation and small RNA sequencing revealed that viral siRNAs broadly associated with wild-type HC-Pro during TuMV infection. These results support the hypothesis that suppression of antiviral silencing during TuMV infection, at least in part, occurs through sequestration of virus-derived siRNAs away from antiviral AGO proteins by HC-Pro. These findings indicate that distinct AGO proteins function as antiviral modules, and provide a molecular explanation for the silencing suppressor activity of HC-Pro.

Functional analysis of three Arabidopsis ARGONAUTES using slicer-defective mutants.

In RNA-directed silencing pathways, ternary complexes result from small RNA-guided ARGONAUTE (AGO) associating with target transcripts. Target transcripts are often silenced through direct cleavage (slicing), destabilization through slicer-independent turnover mechanisms, and translational repression. Here, wild-type and active-site defective forms of several Arabidopsis thaliana AGO proteins involved in posttranscriptional silencing were used to examine several AGO functions, including small RNA binding, interaction with target RNA, slicing or destabilization of target RNA, secondary small interfering RNA formation, and antiviral activity. Complementation analyses in ago mutant plants revealed that the catalytic residues of AGO1, AGO2, and AGO7 are required to restore the defects of Arabidopsis ago1-25, ago2-1, and zip-1 (AGO7-defective) mutants, respectively. AGO2 had slicer activity in transient assays but could not trigger secondary small interfering RNA biogenesis, and catalytically active AGO2 was necessary for local and systemic antiviral activity against Turnip mosaic virus. Slicer-defective AGOs associated with miRNAs and stabilized AGO-miRNA-target RNA ternary complexes in individual target coimmunoprecipitation assays. In genome-wide AGO-miRNA-target RNA coimmunoprecipitation experiments, slicer-defective AGO1-miRNA associated with target RNA more effectively than did wild-type AGO1-miRNA. These data not only reveal functional roles for AGO1, AGO2, and AGO7 slicer activity, but also indicate an approach to capture ternary complexes more efficiently for genome-wide analyses.

Improving cassava bacterial blight resistance by editing the epigenome.

Pathogens rely on expression of host susceptibility (S) genes to promote infection and disease. As DNA methylation is an epigenetic modification that affects gene expression, blocking access to S genes through targeted methylation could increase disease resistance. Xanthomonas phaseoli pv. manihotis, the causal agent of cassava bacterial blight (CBB), uses transcription activator-like20 (TAL20) to induce expression of the S gene MeSWEET10a. In this work, we direct methylation to the TAL20 effector binding element within the MeSWEET10a promoter using a synthetic zinc-finger DNA binding domain fused to a component of the RNA-directed DNA methylation pathway. We demonstrate that this methylation prevents TAL20 binding, blocks transcriptional activation of MeSWEET10a in vivo and that these plants display decreased CBB symptoms while maintaining normal growth and development. This work therefore presents an epigenome editing approach useful for crop improvement.

Mutations in DNA polymerase δ subunit 1 co-segregate with CMD2-type resistance to Cassava Mosaic Geminiviruses.

Cassava mosaic disease (CMD) suppresses cassava yields across the tropics. The dominant CMD2 locus confers resistance to cassava mosaic geminiviruses. It has been reported that CMD2-type landraces lose resistance after regeneration through de novo morphogenesis. As full genome bisulfite sequencing failed to uncover an epigenetic mechanism for this loss of resistance, whole genome sequencing and genetic variant analysis was performed and the CMD2 locus was fine-mapped to a 190 kilobase interval. Collectively, these data indicate that CMD2-type resistance is caused by a nonsynonymous, single nucleotide polymorphism in DNA polymerase δ subunit 1 (MePOLD1) located within this region. Virus-induced gene silencing of MePOLD1 in a CMD-susceptible cassava variety produced a recovery phenotype typical of CMD2-type resistance. Analysis of other CMD2-type cassava varieties identified additional candidate resistance alleles within MePOLD1. Genetic variation of MePOLD1, therefore, could represent an important genetic resource for resistance breeding and/or genome editing, and elucidating mechanisms of resistance to geminiviruses.

Global position analysis of the Pseudomonas aeruginosa quorum-sensing transcription factor LasR.

In Pseudomonas aeruginosa quorum sensing (QS), the transcriptional regulator LasR controls the expression of more than 300 genes. Several of these genes are activated indirectly via a second, subordinate QS regulator, RhlR. Conserved sequence elements upstream of individual other genes have been shown to bind LasR in vitro. To comprehensively identify all regions that are bound by LasR in vivo, we employed chromatin immunoprecipitation in conjunction with microarray analysis. We identified 35 putative promoter regions that direct the expression of up to 74 genes. In vitro DNA binding studies allowed us to distinguish between cooperative and non-cooperative LasR binding sites, and allowed us to build consensus sequences according to the mode of binding. Five promoter regions were not previously recognized as QS-controlled. Two of the associated transcript units encode proteins involved in the cold-shock response and in Psl exopolysaccharide synthesis respectively. The LasR regulon includes seven genes encoding transcriptional regulators, while secreted factors and secretion machinery are the most over-represented functional categories overall. This supports the notion that the core function of LasR is to co-ordinate the production of extracellular factors, although many of its effects on global gene expression are likely mediated indirectly by regulatory genes under its control.

Hiding in plain sight: New virus genomes discovered via a systematic analysis of fungal public transcriptomes.

The distribution and diversity of RNA viruses in fungi is incompletely understood due to the often cryptic nature of mycoviral infections and the focused study of primarily pathogenic and/or economically important fungi. As most viruses that are known to infect fungi possess either single-stranded or double-stranded RNA genomes, transcriptomic data provides the opportunity to query for viruses in diverse fungal samples without any a priori knowledge of virus infection. Here we describe a systematic survey of all transcriptomic datasets from fungi belonging to the subphylum Pezizomycotina. Using a simple but effective computational pipeline that uses reads discarded during normal RNA-seq analyses, followed by identification of a viral RNA-dependent RNA polymerase (RdRP) motif in de novo assembled contigs, 59 viruses from 44 different fungi were identified. Among the viruses identified, 88% were determined to be new species and 68% are, to our knowledge, the first virus described from the fungal species. Comprehensive analyses of both nucleotide and inferred protein sequences characterize the phylogenetic relationships between these viruses and the known set of mycoviral sequences and support the classification of up to four new families and two new genera. Thus the results provide a deeper understanding of the scope of mycoviral diversity while also increasing the distribution of fungal hosts. Further, this study demonstrates the suitability of analyzing RNA-seq data to facilitate rapid discovery of new viruses.

Mutagenesis of the carboxy terminal protease CtpA decreases desiccation tolerance in Rhizobium leguminosarum.

To better understand the role of proteases in Rhizobium leguminosarum biovar viciae, a gene with homology to the carboxy-terminal protease (CtpA), which belongs to a novel group of serine proteases, was studied. The ctpA gene was cloned and mutated using allelic exchange and a gusA reporter gene was used to study ctpA expression. Mutational analysis shows that ctpA is critical for the viability of R. leguminosarum when cells are grown on complex semi-solid media but is dispensable when cells are grown in complex liquid media and that this is likely due to an increase in susceptibility to desiccation on semi-solid media. The ctpA mutant also displayed an increased sensitivity to detergents, indicating an alteration in the permeability of the cell envelope. This is the first characterization of a ctpA gene within the Rhizobiaceae and the first report of a ctpA mutant that exhibits an increased sensitivity to desiccation.

PlantCV v2: Image analysis software for high-throughput plant phenotyping.

Systems for collecting image data in conjunction with computer vision techniques are a powerful tool for increasing the temporal resolution at which plant phenotypes can be measured non-destructively. Computational tools that are flexible and extendable are needed to address the diversity of plant phenotyping problems. We previously described the Plant Computer Vision (PlantCV) software package, which is an image processing toolkit for plant phenotyping analysis. The goal of the PlantCV project is to develop a set of modular, reusable, and repurposable tools for plant image analysis that are open-source and community-developed. Here we present the details and rationale for major developments in the second major release of PlantCV. In addition to overall improvements in the organization of the PlantCV project, new functionality includes a set of new image processing and normalization tools, support for analyzing images that include multiple plants, leaf segmentation, landmark identification tools for morphometrics, and modules for machine learning.

A Versatile Phenotyping System and Analytics Platform Reveals Diverse Temporal Responses to Water Availability in Setaria.

Phenotyping has become the rate-limiting step in using large-scale genomic data to understand and improve agricultural crops. Here, the Bellwether Phenotyping Platform for controlled-environment plant growth and automated multimodal phenotyping is described. The system has capacity for 1140 plants, which pass daily through stations to record fluorescence, near-infrared, and visible images. Plant Computer Vision (PlantCV) was developed as open-source, hardware platform-independent software for quantitative image analysis. In a 4-week experiment, wild Setaria viridis and domesticated Setaria italica had fundamentally different temporal responses to water availability. While both lines produced similar levels of biomass under limited water conditions, Setaria viridis maintained the same water-use efficiency under water replete conditions, while Setaria italica shifted to less efficient growth. Overall, the Bellwether Phenotyping Platform and PlantCV software detected significant effects of genotype and environment on height, biomass, water-use efficiency, color, plant architecture, and tissue water status traits. All ∼ 79,000 images acquired during the course of the experiment are publicly available.

Preparation of Multiplexed Small RNA Libraries From Plants.

High-throughput sequencing is a powerful tool for exploring small RNA populations in plants. The ever-increasing output from an Illumina Sequencing System allows for multiplexing multiple samples while still obtaining sufficient data for small RNA discovery and characterization. Here we describe a protocol for generating multiplexed small RNA libraries for sequencing up to 12 samples in one lane of an Illumina HiSeq System single-end, 50 base pair run. RNA ligases are used to add the 3' and 5' adaptors to purified small RNAs; ligation products that lack a small RNA molecule (adaptor-adaptor products) are intentionally depleted. After cDNA synthesis, a linear PCR step amplifies the DNA fragments. The 3' PCR primers used here include unique 6-nucleotide sequences to allow for multiplexing up to 12 samples.

Phytophthora have distinct endogenous small RNA populations that include short interfering and microRNAs.

In eukaryotes, RNA silencing pathways utilize 20-30-nucleotide small RNAs to regulate gene expression, specify and maintain chromatin structure, and repress viruses and mobile genetic elements. RNA silencing was likely present in the common ancestor of modern eukaryotes, but most research has focused on plant and animal RNA silencing systems. Phytophthora species belong to a phylogenetically distinct group of economically important plant pathogens that cause billions of dollars in yield losses annually as well as ecologically devastating outbreaks. We analyzed the small RNA-generating components of the genomes of P. infestans, P. sojae and P. ramorum using bioinformatics, genetic, phylogenetic and high-throughput sequencing-based methods. Each species produces two distinct populations of small RNAs that are predominantly 21- or 25-nucleotides long. The 25-nucleotide small RNAs were primarily derived from loci encoding transposable elements and we propose that these small RNAs define a pathway of short-interfering RNAs that silence repetitive genetic elements. The 21-nucleotide small RNAs were primarily derived from inverted repeats, including a novel microRNA family that is conserved among the three species, and several gene families, including Crinkler effectors and type III fibronectins. The Phytophthora microRNA is predicted to target a family of amino acid/auxin permeases, and we propose that 21-nucleotide small RNAs function at the post-transcriptional level. The functional significance of microRNA-guided regulation of amino acid/auxin permeases and the association of 21-nucleotide small RNAs with Crinkler effectors remains unclear, but this work provides a framework for testing the role of small RNAs in Phytophthora biology and pathogenesis in future work.