Chromosomal barcodes for simultaneous tracking of near-isogenic bacterial strains in plant microbiota

DNA-amplicon-based microbiota profiling can estimate species diversity and abundance but cannot resolve genetic differences within individuals of the same species. Here we report the development of modular bacterial tags (MoBacTags) encoding DNA barcodes that enable tracking of near-isogenic bacterial commensals in an array of complex microbiome communities. Chromosomally integrated DNA barcodes are then co-amplified with endogenous marker genes of the community by integrating corresponding primer binding sites into the barcode. We use this approach to assess the contributions of individual bacterial genes to Arabidopsis thaliana root microbiota establishment with synthetic communities that include MoBacTag-labelled strains of Pseudomonas capeferrum. Results show reduced root colonization for certain mutant strains with defects in gluconic-acid-mediated host immunosuppression, which would not be detected with traditional amplicon sequencing. Our work illustrates how MoBacTags can be applied to assess scaling of individual bacterial genetic determinants in the plant microbiota.

Genome-resolved metatranscriptomics reveals conserved root colonization determinants in a synthetic microbiota

The identification of processes activated by specific microbes during microbiota colonization of plant roots has been hampered by technical constraints in metatranscriptomics. These include lack of reference genomes, high representation of host or microbial rRNA sequences in datasets, or difficulty to experimentally validate gene functions. Here, we recolonized germ-free Arabidopsis thaliana with a synthetic, yet representative root microbiota comprising 106 genome-sequenced bacterial and fungal isolates. We used multi-kingdom rRNA depletion, deep RNA-sequencing and read mapping against reference microbial genomes to analyse the in planta metatranscriptome of abundant colonizers. We identified over 3,000 microbial genes that were differentially regulated at the soil-root interface. Translation and energy production processes were consistently activated in planta, and their induction correlated with bacterial strains' abundance in roots. Finally, we used targeted mutagenesis to show that several genes consistently induced by multiple bacteria are required for root colonization in one of the abundant bacterial strains (a genetically tractable Rhodanobacter). Our results indicate that microbiota members activate strain-specific processes but also common gene sets to colonize plant roots.

Cofunctioning of bacterial exometabolites drives root microbiota establishment

Soil-dwelling microbes are the principal inoculum for the root microbiota, but our understanding of microbe-microbe interactions in microbiota establishment remains fragmentary. We tested 39,204 binary interbacterial interactions for inhibitory activities in vitro, allowing us to identify taxonomic signatures in bacterial inhibition profiles. Using genetic and metabolomic approaches, we identified the antimicrobial 2,4-diacetylphloroglucinol (DAPG) and the iron chelator pyoverdine as exometabolites whose combined functions explain most of the inhibitory activity of the strongly antagonistic Pseudomonas brassicacearum R401. Microbiota reconstitution with a core of Arabidopsis thaliana root commensals in the presence of wild-type or mutant strains revealed a root niche-specific cofunction of these exometabolites as root competence determinants and drivers of predictable changes in the root-associated community. In natural environments, both the corresponding biosynthetic operons are enriched in roots, a pattern likely linked to their role as iron sinks, indicating that these cofunctioning exometabolites are adaptive traits contributing to pseudomonad pervasiveness throughout the root microbiota.

Targeted gene deletion with SpCas9 and multiple guide RNAs in Arabidopsis thaliana: four are better than two

BACKGROUND

In plant genome editing, RNA-guided nucleases such as Cas9 from Streptococcus pyogenes (SpCas9) predominantly induce small insertions or deletions at target sites. This can be used for inactivation of protein-coding genes by frame shift mutations. However, in some cases, it may be advantageous to delete larger chromosomal segments. This is achieved by simultaneously inducing double strand breaks upstream and downstream of the segment to be deleted. Experimental approaches for the deletion of larger chromosomal segments have not been systematically evaluated.

RESULTS

We designed three pairs of guide RNAs for deletion of a ~ 2.2 kb chromosomal segment containing the Arabidopsis WRKY30 locus. We tested how the combination of guide RNA pairs and co-expression of the exonuclease TREX2 affect the frequency of wrky30 deletions in editing experiments. Our data demonstrate that compared to one pair of guide RNAs, two pairs increase the frequency of chromosomal deletions. The exonuclease TREX2 enhanced mutation frequency at individual target sites and shifted the mutation profile towards larger deletions. However, TREX2 did not elevate the frequency of chromosomal segment deletions.

CONCLUSIONS

Multiplex editing with at least two pairs of guide RNAs (four guide RNAs in total) elevates the frequency of chromosomal segment deletions at least at the AtWRKY30 locus, and thus simplifies the selection of corresponding mutants. Co-expression of the TREX2 exonuclease can be used as a general strategy to increase editing efficiency in Arabidopsis without obvious negative effects.

Gnotobiotic Plant Systems for Reconstitution and Functional Studies of the Root Microbiota

Healthy plants host a multi-kingdom community of microbes, which is known as the plant microbiota. Amplicon sequencing technologies for microbial genomic markers were a milestone in revealing the taxonomic composition of the microbiota and its variation associated with a plant host in natural environments. However, this method alone does not allow conclusions to be drawn about functions of these microbial assemblages for the plant. The development of culture collections, which recapitulate natural microbial communities in their diversity, and multiple gnotobiotic plant systems therefore represent a breakthrough in plant-microbiota research such that plants can be inoculated with defined communities to study proposed microbiota functions. These systems provided, for the root microbiota, first insights into mechanisms underlying microbial community establishment and contributions of its microbial members to indirect pathogen protection and mineral nutrition of the host. We argue that the choice of a gnotobiotic system for microbiota reconstitution and subsequent functional analysis depends on the particular plant trait that is influenced by the microbiota. We start by discussing the advantages and limitations of using individual gnotobiotic systems and then describe the general procedures for preparing bacterial cultures from the Arabidopsis thaliana At-R-SPHERE culture collection for inoculation and cocultivation in two gnotobiotic plant growth systems using agar and perlite matrix. Additionally, a protocol for inoculation of plants with opportunistic Pseudomonas pathogens is provided. Lastly, we describe a high-throughput system for visual assessment of roots after inoculation with individual mutants of a transposon library generated from a root-derived bacterial commensal. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Preparation of bacterial cultures from At-R-SPHERE Support Protocol 1: Validation of strains by sequencing hypervariable regions of the 16S rRNA gene Basic Protocol 2: Coinoculation of plants grown on an agar matrix with microbial elicitor and a defined microbial community Alternate Protocol: Inoculation of plants cultivated in a perlite-based growth system Support Protocol 2: Surface sterilization of Arabidopsis thaliana seeds Basic Protocol 3: Inoculation using a Pseudomonas opportunistic pathogen Basic Protocol 4: Assessment of commensal-mediated root phenotypes using phytostrips.

Coordination of microbe-host homeostasis by crosstalk with plant innate immunity

Plants grown in natural soil are colonized by phylogenetically structured communities of microbes known as the microbiota. Individual microbes can activate microbe-associated molecular pattern (MAMP)-triggered immunity (MTI), which limits pathogen proliferation but curtails plant growth, a phenomenon known as the growth-defence trade-off. Here, we report that, in monoassociations, 41% (62 out of 151) of taxonomically diverse root bacterial commensals suppress Arabidopsis thaliana root growth inhibition (RGI) triggered by immune-stimulating MAMPs or damage-associated molecular patterns. Amplicon sequencing of bacterial 16S rRNA genes reveals that immune activation alters the profile of synthetic communities (SynComs) comprising RGI-non-suppressive strains, whereas the presence of RGI-suppressive strains attenuates this effect. Root colonization by SynComs with different complexities and RGI-suppressive activities alters the expression of 174 core host genes, with functions related to root development and nutrient transport. Furthermore, RGI-suppressive SynComs specifically downregulate a subset of immune-related genes. Precolonization of plants with RGI-suppressive SynComs, or mutation of one commensal-downregulated transcription factor, MYB15, renders the plants more susceptible to opportunistic Pseudomonas pathogens. Our results suggest that RGI-non-suppressive and RGI-suppressive root commensals modulate host susceptibility to pathogens by either eliciting or dampening MTI responses, respectively. This interplay buffers the plant immune system against pathogen perturbation and defence-associated growth inhibition, ultimately leading to commensal-host homeostasis.

Disentangling cause and consequence: genetic dissection of the DANGEROUS MIX2 risk locus, and activation of the DM2h NLR in autoimmunity

Nucleotide-binding domain-leucine-rich repeat-type immune receptors (NLRs) protect plants against pathogenic microbes through intracellular detection of effector proteins. However, this comes at a cost, as NLRs can also induce detrimental autoimmunity in genetic interactions with foreign alleles. This may occur when independently evolved genomes are combined in inter- or intraspecific crosses, or when foreign alleles are introduced by mutagenesis or transgenesis. Most autoimmunity-inducing NLRs are encoded within highly variable NLR gene clusters with no known immune functions, which were termed autoimmune risk loci. Whether risk NLRs differ from sensor NLRs operating in natural pathogen resistance and how risk NLRs are activated in autoimmunity is unknown. Here, we analyzed the DANGEROUS MIX2 risk locus, a major autoimmunity hotspot in Arabidopsis thaliana. By gene editing and heterologous expression, we show that a single gene, DM2h, is necessary and sufficient for autoimmune induction in three independent cases of autoimmunity in accession Landsberg erecta. We focus on autoimmunity provoked by an EDS1-yellow fluorescent protein (YFP)NLS fusion protein to characterize DM2h functionally and determine features of EDS1-YFPNLS activating the immune receptor. Our data suggest that risk NLRs function in a manner reminiscent of sensor NLRs, while autoimmunity-inducing properties of EDS1-YFPNLS in this context are unrelated to the protein's functions as an immune regulator. We propose that autoimmunity, at least in some cases, may be caused by spurious, stochastic interactions of foreign alleles with coincidentally matching risk NLRs.

Highly efficient multiplex editing: one-shot generation of 8× Nicotiana benthamiana and 12× Arabidopsis mutants

Genome editing by RNA-guided nucleases, such as SpCas9, has been used in numerous different plant species. However, to what extent multiple independent loci can be targeted simultaneously by multiplexing has not been well documented. Here, we developed a toolkit, based on a highly intron-optimized zCas9i gene, which allows assembly of nuclease constructs expressing up to 32 single guide RNAs (sgRNAs). We used this toolkit to explore the limits of multiplexing in two major model species, and report on the isolation of transgene-free octuple (8×) Nicotiana benthamiana and duodecuple (12×) Arabidopsis thaliana mutant lines in a single generation (T1 and T2 , respectively). We developed novel counter-selection markers for N. benthamiana, most importantly Sl-FAST2, comparable to the well-established Arabidopsis seed fluorescence marker, and FCY-UPP, based on the production of toxic 5-fluorouracil in the presence of a precursor. Targeting eight genes with an array of nine different sgRNAs and relying on FCY-UPP for selection of non-transgenic T1 , we identified N. benthamiana mutant lines with astonishingly high efficiencies: All analyzed plants carried mutations in all genes (approximately 112/116 target sites edited). Furthermore, we targeted 12 genes by an array of 24 sgRNAs in A. thaliana. Efficiency was significantly lower in A. thaliana, and our results indicate Cas9 availability is the limiting factor in such higher-order multiplexing applications. We identified a duodecuple mutant line by a combination of phenotypic screening and amplicon sequencing. The resources and results presented provide new perspectives for how multiplexing can be used to generate complex genotypes or to functionally interrogate groups of candidate genes.

An EDS1-SAG101 Complex Is Essential for TNL-Mediated Immunity in Nicotiana benthamiana

Heterodimeric complexes containing the lipase-like protein ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1) are regarded as central regulators of plant innate immunity. In this context, a complex of EDS1 with PHYTOALEXIN DEFICIENT4 (PAD4) is required for basal resistance and signaling downstream of immune receptors containing an N-terminal Toll-interleukin-1 receptor-like domain (TNLs) in Arabidopsis (Arabidopsis thaliana). Here we analyze EDS1 functions in the model Solanaceous plant Nicotiana benthamiana (Nb). Stable Nb mutants deficient in EDS1 complexes are not impaired in basal resistance, a finding which contradicts a general role for EDS1 in immunity. In Nb, PAD4 demonstrated no detectable immune functions, but TNL-mediated resistance responses required EDS1 complexes incorporating a SENESCENCE ASSOCIATED GENE101 (SAG101) isoform. Intriguingly, SAG101 is restricted to those genomes also encoding TNL receptors, and we propose it may be required for TNL-mediated immune signaling in most plants, except the Brassicaceae. Transient complementation in Nb was used for accelerated mutational analyses while avoiding complex biotic interactions. We identify a large surface essential for EDS1-SAG101 immune functions that extends from the N-terminal lipase domains to the C-terminal EDS1-PAD4 domains and might mediate interaction partner recruitment. Furthermore, this work demonstrates the value of genetic resources in Nb, which will facilitate elucidation of EDS1 functions.

Optimized Cas9 expression systems for highly efficient Arabidopsis genome editing facilitate isolation of complex alleles in a single generation

Genetic resources for the model plant Arabidopsis comprise mutant lines defective in almost any single gene in reference accession Columbia. However, gene redundancy and/or close linkage often render it extremely laborious or even impossible to isolate a desired line lacking a specific function or set of genes from segregating populations. Therefore, we here evaluated strategies and efficiencies for the inactivation of multiple genes by Cas9-based nucleases and multiplexing. In first attempts, we succeeded in isolating a mutant line carrying a 70 kb deletion, which occurred at a frequency of ~ 1.6% in the T₂ generation, through PCR-based screening of numerous individuals. However, we failed to isolate a line lacking Lhcb1 genes, which are present in five copies organized at two loci in the Arabidopsis genome. To improve efficiency of our Cas9-based nuclease system, regulatory sequences controlling Cas9 expression levels and timing were systematically compared. Indeed, use of DD45 and RPS5a promoters improved efficiency of our genome editing system by approximately 25-30-fold in comparison to the previous ubiquitin promoter. Using an optimized genome editing system with RPS5a promoter-driven Cas9, putatively quintuple mutant lines lacking detectable amounts of Lhcb1 protein represented approximately 30% of T₁ transformants. These results show how improved genome editing systems facilitate the isolation of complex mutant alleles, previously considered impossible to generate, at high frequency even in a single (T₁) generation.

Further analysis of barley MORC1 using a highly efficient RNA-guided Cas9 gene-editing system

Microrchidia (MORC) proteins comprise a family of proteins that have been identified in prokaryotes and eukaryotes. They are defined by two hallmark domains: a GHKL-type ATPase and an S5-fold. In plants, MORC proteins were first discovered in a genetic screen for Arabidopsis thaliana mutants compromised for resistance to a viral pathogen. Subsequent studies expanded their role in plant immunity and revealed their involvement in gene silencing and genome stabilization. Little is known about the role of MORC proteins of cereals, especially because knockout (KO) mutants were not available and assessment of loss of function relied only on RNAi strategies, which were arguable, given that MORC proteins in itself are influencing gene silencing. Here, we used a Streptococcus pyogenes Cas9 (SpCas9)-mediated KO strategy to functionally study HvMORC1, one of the current seven MORC members of barley. Using a novel barley RNA Pol III-dependent U3 small nuclear RNA (snRNA) promoter to drive expression of the synthetic single guide RNA (sgRNA), we achieved a very high mutation frequency in HvMORC1. High frequencies of mutations were detectable by target sequencing in the callus, the T0 generation (77%) and T1 generation (70%-100%), which constitutes an important improvement of the gene-editing technology in cereals. Corroborating and extending earlier findings, SpCas9-edited hvmorc1-KO barley, in clear contrast to Arabidopsis atmorc1 mutants, had a distinct phenotype of increased disease resistance to fungal pathogens, while morc1 mutants of either plant showed de-repressed expression of transposable elements (TEs), substantiating that plant MORC proteins contribute to genome stabilization in monocotyledonous and dicotyledonous plants.

Peripheral infrastructure vectors and an extended set of plant parts for the Modular Cloning system

Standardized DNA assembly strategies facilitate the generation of multigene constructs from collections of building blocks in plant synthetic biology. A common syntax for hierarchical DNA assembly following the Golden Gate principle employing Type IIs restriction endonucleases was recently developed, and underlies the Modular Cloning and GoldenBraid systems. In these systems, transcriptional units and/or multigene constructs are assembled from libraries of standardized building blocks, also referred to as phytobricks, in several hierarchical levels and by iterative Golden Gate reactions. Here, a toolkit containing further modules for the novel DNA assembly standards was developed. Intended for use with Modular Cloning, most modules are also compatible with GoldenBraid. Firstly, a collection of approximately 80 additional phytobricks is provided, comprising e.g. modules for inducible expression systems, promoters or epitope tags. Furthermore, DNA modules were developed for connecting Modular Cloning and Gateway cloning, either for toggling between systems or for standardized Gateway destination vector assembly. Finally, first instances of a "peripheral infrastructure" around Modular Cloning are presented: While available toolkits are designed for the assembly of plant transformation constructs, vectors were created to also use coding sequence-containing phytobricks directly in yeast two hybrid interaction or bacterial infection assays. The presented material will further enhance versatility of hierarchical DNA assembly strategies.

Generation of chromosomal deletions in dicotyledonous plants employing a user-friendly genome editing toolkit

Genome editing facilitated by Cas9-based RNA-guided nucleases (RGNs) is becoming an increasingly important and popular technique for reverse genetics in both model and non-model species. So far, RGNs were mainly applied for the induction of point mutations, and one major challenge consists in the detection of genome-edited individuals from a mutagenized population. Also, point mutations are not appropriate for functional dissection of non-coding DNA. Here, the multiplexing capacity of a newly developed genome editing toolkit was exploited for the induction of inheritable chromosomal deletions at six different loci in Nicotiana benthamiana and Arabidopsis. In both species, the preferential formation of small deletions was observed, suggesting reduced efficiency with increasing deletion size. Importantly, small deletions (<100 bp) were detected at high frequencies in N. benthamiana T0 and Arabidopsis T2 populations. Thus, targeting of small deletions by paired nucleases represents a simple approach for the generation of mutant alleles segregating as size polymorphisms in subsequent generations. Phenotypically selected deletions of up to 120 kb occurred at low frequencies in Arabidopsis, suggesting larger population sizes for the discovery of valuable alleles from addressing gene clusters or non-coding DNA for deletion by programmable nucleases.