Extensive immune receptor repertoire diversity in disease-resistant rice landraces.

Plants have powerful defense mechanisms and extensive immune receptor repertoires, yet crop monocultures are prone to epidemic diseases. Rice (Oryza sativa) is susceptible to many diseases, such as rice blast caused by Magnaporthe oryzae. Varietal resistance of rice to blast relies on intracellular nucleotide binding, leucine-rich repeat (NLR) receptors that recognize specific pathogen molecules and trigger immune responses. In the Yuanyang terraces in southwest China, rice landraces rarely show severe losses to disease whereas commercial inbred lines show pronounced field susceptibility. Here, we investigate within-landrace NLR sequence diversity of nine rice landraces and eleven modern varieties using complexity reduction techniques. We find that NLRs display high sequence diversity in landraces, consistent with balancing selection, and that balancing selection at NLRs is more pervasive in landraces than modern varieties. Notably, modern varieties lack many ancient NLR haplotypes that are retained in some landraces. Our study emphasizes the value of standing genetic variation that is maintained in farmer landraces as a resource to make modern crops and agroecosystems less prone to disease. The conservation of landraces is, therefore, crucial for ensuring food security in the face of dynamic biotic and abiotic threats.

The unconventional resistance protein PTR recognizes the Magnaporthe oryzae effector AVR-Pita in an allele-specific manner.

Blast disease caused by the fungus Magnaporthe oryzae is one of the most devastating rice diseases. Disease resistance genes such as Pi-ta or Pi-ta2 are critical in protecting rice production from blast. Published work reports that Pi-ta codes for a nucleotide-binding and leucine-rich repeat domain protein (NLR) that recognizes the fungal protease-like effector AVR-Pita by direct binding. However, this model was challenged by the recent discovery that Pi-ta2 resistance, which also relies on AVR-Pita detection, is conferred by the unconventional resistance gene Ptr, which codes for a membrane protein with a cytoplasmic armadillo repeat domain. Here, using NLR Pi-ta and Ptr RNAi knockdown and CRISPR/Cas9 knockout mutant rice lines, we found that AVR-Pita recognition relies solely on Ptr and that the NLR Pi-ta has no role in it, indicating that it is not the Pi-ta resistance gene. Different alleles of Ptr confer different recognition specificities. The A allele of Ptr (PtrA) detects all natural sequence variants of the effector and confers Pi-ta2 resistance, while the B allele of Ptr (PtrB) recognizes a restricted set of AVR-Pita alleles and, thereby, confers Pi-ta resistance. Analysis of the natural diversity in AVR-Pita and of mutant and transgenic strains identified one specific polymorphism in the effector sequence that controls escape from PtrB-mediated resistance. Taken together, our work establishes that the M. oryzae effector AVR-Pita is detected in an allele-specific manner by the unconventional rice resistance protein Ptr and that the NLR Pi-ta has no function in Pi-ta resistance and the recognition of AVR-Pita.

The structural landscape and diversity of Pyricularia oryzae MAX effectors revisited.

Magnaporthe AVRs and ToxB-like (MAX) effectors constitute a family of secreted virulence proteins in the fungus Pyricularia oryzae (syn. Magnaporthe oryzae), which causes blast disease on numerous cereals and grasses. In spite of high sequence divergence, MAX effectors share a common fold characterized by a ß-sandwich core stabilized by a conserved disulfide bond. In this study, we investigated the structural landscape and diversity within the MAX effector repertoire of P. oryzae. Combining experimental protein structure determination and in silico structure modeling we validated the presence of the conserved MAX effector core domain in 77 out of 94 groups of orthologs (OG) identified in a previous population genomic study. Four novel MAX effector structures determined by NMR were in remarkably good agreement with AlphaFold2 (AF2) predictions. Based on the comparison of the AF2-generated 3D models we propose a classification of the MAX effectors superfamily in 20 structural groups that vary in the canonical MAX fold, disulfide bond patterns, and additional secondary structures in N- and C-terminal extensions. About one-third of the MAX family members remain singletons, without strong structural relationship to other MAX effectors. Analysis of the surface properties of the AF2 MAX models also highlights the high variability within the MAX family at the structural level, potentially reflecting the wide diversity of their virulence functions and host targets.

Genome-wide association analysis uncovers rice blast resistance alleles of Ptr and Pia.

A critical step to maximize the usefulness of genome-wide association studies (GWAS) in plant breeding is the identification and validation of candidate genes underlying genetic associations. This is of particular importance in disease resistance breeding where allelic variants of resistance genes often confer resistance to distinct populations, or races, of a pathogen. Here, we perform a genome-wide association analysis of rice blast resistance in 500 genetically diverse rice accessions. To facilitate candidate gene identification, we produce de-novo genome assemblies of ten rice accessions with various rice blast resistance associations. These genome assemblies facilitate the identification and functional validation of novel alleles of the rice blast resistance genes Ptr and Pia. We uncover an allelic series for the unusual Ptr rice blast resistance gene, and additional alleles of the Pia resistance genes RGA4 and RGA5. By linking these associations to three thousand rice genomes we provide a useful tool to inform future rice blast breeding efforts. Our work shows that GWAS in combination with whole-genome sequencing is a powerful tool for gene cloning and to facilitate selection of specific resistance alleles for plant breeding.

Pyricularia oryzae: Lab star and field scourge.

Pyricularia oryzae (syn. Magnaporthe oryzae), is a filamentous ascomycete that causes a major disease called blast on cereal crops, as well as on a wide variety of wild and cultivated grasses. Blast diseases have a tremendous impact worldwide particularly on rice and on wheat, where the disease emerged in South America in the 1980s, before spreading to Asia and Africa. Its economic importance, coupled with its amenability to molecular and genetic manipulation, have inspired extensive research efforts aiming at understanding its biology and evolution. In the past 40 years, this plant-pathogenic fungus has emerged as a major model in molecular plant-microbe interactions. In this review, we focus on the clarification of the taxonomy and genetic structure of the species and its host range determinants. We also discuss recent molecular studies deciphering its lifecycle. TAXONOMY: Kingdom: Fungi, phylum: Ascomycota, sub-phylum: Pezizomycotina, class: Sordariomycetes, order: Magnaporthales, family: Pyriculariaceae, genus: Pyricularia. HOST RANGE: P. oryzae has the ability to infect a wide range of Poaceae. It is structured into different host-specialized lineages that are each associated with a few host plant genera. The fungus is best known to cause tremendous damage to rice crops, but it can also attack other economically important crops such as wheat, maize, barley, and finger millet. DISEASE SYMPTOMS: P. oryzae can cause necrotic lesions or bleaching on all aerial parts of its host plants, including leaf blades, sheaths, and inflorescences (panicles, spikes, and seeds). Characteristic symptoms on leaves are diamond-shaped silver lesions that often have a brown margin and whose appearance is influenced by numerous factors such as the plant genotype and environmental conditions. USEFUL WEBSITES Resources URL Genomic data repositories http://genome.jouy.inra.fr/gemo/ Genomic data repositories http://openriceblast.org/ Genomic data repositories http://openwheatblast.net/ Genome browser for fungi (including P. oryzae) http://fungi.ensembl.org/index.html Comparative genomics database https://mycocosm.jgi.doe.gov/mycocosm/home T-DNA mutant database http://atmt.snu.kr/ T-DNA mutant database http://www.phi-base.org/ SNP and expression data https://fungidb.org/fungidb/app/.

The bZIP transcription factor BIP1 of the rice blast fungus is essential for infection and regulates a specific set of appressorium genes.

The rice blast fungus Magnaporthe oryzae differentiates specialized cells called appressoria that are required for fungal penetration into host leaves. In this study, we identified the novel basic leucine zipper (bZIP) transcription factor BIP1 (B-ZIP Involved in Pathogenesis-1) that is essential for pathogenicity. BIP1 is required for the infection of plant leaves, even if they are wounded, but not for appressorium-mediated penetration of artificial cellophane membranes. This phenotype suggests that BIP1 is not implicated in the differentiation of the penetration peg but is necessary for the initial establishment of the fungus within plant cells. BIP1 expression was restricted to the appressorium by both transcriptional and post-transcriptional control. Genome-wide transcriptome analysis showed that 40 genes were down regulated in a BIP1 deletion mutant. Most of these genes were specifically expressed in the appressorium. They encode proteins with pathogenesis-related functions such as enzymes involved in secondary metabolism including those encoded by the ACE1 gene cluster, small secreted proteins such as SLP2, BAS2, BAS3, and AVR-Pi9 effectors, as well as plant cuticle and cell wall degrading enzymes. Interestingly, this BIP1 network is different from other known infection-related regulatory networks, highlighting the complexity of gene expression control during plant-fungal interactions. Promoters of BIP1-regulated genes shared a GCN4/bZIP-binding DNA motif (TGACTC) binding in vitro to BIP1. Mutation of this motif in the promoter of MGG_08381.7 from the ACE1 gene cluster abolished its appressorium-specific expression, showing that BIP1 behaves as a transcriptional activator. In summary, our findings demonstrate that BIP1 is critical for the expression of early invasion-related genes in appressoria. These genes are likely needed for biotrophic invasion of the first infected host cell, but not for the penetration process itself. Through these mechanisms, the blast fungus strategically anticipates the host plant environment and responses during appressorium-mediated penetration.

Adaptive evolution in virulence effectors of the rice blast fungus Pyricularia oryzae.

Plant pathogens secrete proteins called effectors that target host cellular processes to promote disease. Recently, structural genomics has identified several families of fungal effectors that share a similar three-dimensional structure despite remarkably variable amino-acid sequences and surface properties. To explore the selective forces that underlie the sequence variability of structurally-analogous effectors, we focused on MAX effectors, a structural family of effectors that are major determinants of virulence in the rice blast fungus Pyricularia oryzae. Using structure-informed gene annotation, we identified 58 to 78 MAX effector genes per genome in a set of 120 isolates representing seven host-associated lineages. The expression of MAX effector genes was primarily restricted to the early biotrophic phase of infection and strongly influenced by the host plant. Pangenome analyses of MAX effectors demonstrated extensive presence/absence polymorphism and identified gene loss events possibly involved in host range adaptation. However, gene knock-in experiments did not reveal a strong effect on virulence phenotypes suggesting that other evolutionary mechanisms are the main drivers of MAX effector losses. MAX effectors displayed high levels of standing variation and high rates of non-synonymous substitutions, pointing to widespread positive selection shaping the molecular diversity of MAX effectors. The combination of these analyses with structural data revealed that positive selection acts mostly on residues located in particular structural elements and at specific positions. By providing a comprehensive catalog of amino acid polymorphism, and by identifying the structural determinants of the sequence diversity, our work will inform future studies aimed at elucidating the function and mode of action of MAX effectors.

Does a Similar 3D Structure Mean a Similar Folding Pathway? The Presence of a C-Terminal α-Helical Extension in the 3D Structure of MAX60 Drastically Changes the Folding Pathway Described for Other MAX-Effectors from Magnaporthe oryzae.

Does a similar 3D structure mean a similar folding pathway? This question is particularly meaningful when it concerns proteins sharing a similar 3D structure, but low sequence identity or homology. MAX effectors secreted by the phytopathogenic fungus Magnaporthe oryzae present such characteristics. They share a common 3D structure, a ß-sandwich with the same topology for all the family members, but an extremely low sequence identity/homology. In a previous study, we have investigated the folding of two MAX effectors, AVR-Pia and AVR-Pib, using High-Hydrostatic-Pressure NMR and found that they display a similar folding pathway, with a common folding intermediate. In the present work, we used a similar strategy to investigate the folding conformational landscape of another MAX effector, MAX60, and found a very different folding intermediate. Our analysis strongly supports that the presence of a C-terminal α-helical extension in the 3D structure of MAX60 could be responsible for its different folding pathway.

The RLCK subfamily VII-4 controls pattern-triggered immunity and basal resistance to bacterial and fungal pathogens in rice.

Receptor-like cytoplasmic kinases (RLCKs) mediate the intracellular signaling downstream of pattern-recognition receptors (PRRs). Several RLCKs from subfamily VII of rice (Oryza sativa) have important roles in plant immunity, but the role of RLCK VII-4 in pattern-triggered immune (PTI) signaling and resistance to pathogens has not yet been investigated. Here, we generated by multiplex clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9-mediated genome editing rice sextuple mutant lines where the entire RLCK VII-4 subfamily is inactivated and then analyzed the resulting lines for their response to chitin and flg22 and for their immunity to Xanthomonas oryzae pv. oryzae (Xoo) and Magnaporthe oryzae. Analysis of the rlckvii-4 mutants revealed that they have an impaired reactive oxygen system burst and reduced defense gene expression in response to flg22 and chitin. This indicates that members of the rice RLCK VII-4 subfamily are required for immune signaling downstream of multiple PRRs. Furthermore, we found that the rice RLCK VII-4 subfamily is important for chitin-induced callose deposition and mitogen-activated protein kinase activation and that it is crucial for basal resistance against Xoo and M. oryzae pathogens. This establishes that the RLCK VII-4 subfamily has critical functions in the regulation of multiple PTI pathways in rice and opens the way for deciphering the precise role of its members in the control of rice PTI.

Molecular engineering of plant immune receptors for tailored crop disease resistance.

The specific recognition of pathogen effectors by intracellular nucleotide-binding and leucine-rich repeat domain receptors (NLRs) is an important component of plant immunity. Creating NLRs with new bespoke recognition specificities is a major goal in molecular plant pathology as it promises to provide unlimited resources for the resistance of crops against diseases. Recent breakthrough discoveries on the structure and molecular activity of NLRs begin to enable their knowledge-guided molecular engineering. First, studies succeeded to extend or change effector recognition specificities by modifying, in a structure-guided manner, the NLR domains that directly bind effectors. By modifying the LRR domain of the singleton NLR Sr35 or the unconventional decoy domains of the helper NLRs RGA5 or Pik-1, receptors that detected other or additional effectors were created.

Rapid evolution of an RNA virus to escape recognition by a rice nucleotide-binding and leucine-rich repeat domain immune receptor.

Viral diseases are a major limitation for crop production, and their control is crucial for sustainable food supply. We investigated by a combination of functional genetics and experimental evolution the resistance of rice to the rice yellow mottle virus (RYMV), which is among the most devastating rice pathogens in Africa, and the mechanisms underlying the extremely fast adaptation of the virus to its host. We found that the RYMV3 gene that protects rice against the virus codes for a nucleotide-binding and leucine-rich repeat domain immune receptor (NLRs) from the Mla-like clade of NLRs. RYMV3 detects the virus by forming a recognition complex with the viral coat protein (CP). The virus escapes efficiently from detection by mutations in its CP, some of which interfere with the formation of the recognition complex. This study establishes that NLRs also confer in monocotyledonous plants immunity to viruses, and reveals an unexpected functional diversity for NLRs of the Mla clade that were previously only known as fungal disease resistance proteins. In addition, it provides precise insight into the mechanisms by which viruses adapt to plant immunity and gives important knowledge for the development of sustainable resistance against viral diseases of cereals.

Correction: Maintenance of divergent lineages of the Rice Blast Fungus Pyricularia oryzae through niche separation, loss of sex and post-mating genetic incompatibilities.

[This corrects the article DOI: 10.1371/journal.ppat.1010687.].

Maintenance of divergent lineages of the Rice Blast Fungus Pyricularia oryzae through niche separation, loss of sex and post-mating genetic incompatibilities.

Many species of fungal plant pathogens coexist as multiple lineages on the same host, but the factors underlying the origin and maintenance of population structure remain largely unknown. The rice blast fungus Pyricularia oryzae is a widespread model plant pathogen displaying population subdivision. However, most studies of natural variation in P. oryzae have been limited in genomic or geographic resolution, and host adaptation is the only factor that has been investigated extensively as a contributor to population subdivision. In an effort to complement previous studies, we analyzed genetic and phenotypic diversity in isolates of the rice blast fungus covering a broad geographical range. Using single-nucleotide polymorphism genotyping data for 886 isolates sampled from 152 sites in 51 countries, we showed that population subdivision of P. oryzae in one recombining and three clonal lineages with broad distributions persisted with deeper sampling. We also extended previous findings by showing further population subdivision of the recombining lineage into one international and three Asian clusters, and by providing evidence that the three clonal lineages of P. oryzae were found in areas with different prevailing environmental conditions, indicating niche separation. Pathogenicity tests and bioinformatic analyses using an extended set of isolates and rice varieties indicated that partial specialization to rice subgroups contributed to niche separation between lineages, and differences in repertoires of putative virulence effectors were consistent with differences in host range. Experimental crosses revealed that female sterility and early post-mating genetic incompatibilities acted as strong additional barriers to gene flow between clonal lineages. Our results demonstrate that the spread of a fungal pathogen across heterogeneous habitats and divergent populations of a crop species can lead to niche separation and reproductive isolation between distinct, widely distributed, lineages.

Insight into the structure and molecular mode of action of plant paired NLR immune receptors.

The specific recognition of pathogen effectors by intracellular nucleotide-binding domain and leucine-rich repeat receptors (NLRs) is an important component of plant immunity. NLRs have a conserved modular architecture and can be subdivided according to their signaling domain that is mostly a coiled-coil (CC) or a Toll/Interleukin1 receptor (TIR) domain into CNLs and TNLs. Single NLR proteins are often sufficient for both effector recognition and immune activation. However, sometimes, they act in pairs, where two different NLRs are required for disease resistance. Functional studies have revealed that in these cases one NLR of the pair acts as a sensor (sNLR) and one as a helper (hNLR). The genes corresponding to such resistance protein pairs with one-to-one functional co-dependence are clustered, generally with a head-to-head orientation and shared promoter sequences. sNLRs in such functional NLR pairs have additional, non-canonical and highly diverse domains integrated in their conserved modular architecture, which are thought to act as decoys to trap effectors. Recent structure-function studies on the Arabidopsis thaliana TNL pair RRS1/RPS4 and on the rice CNL pairs RGA4/RGA5 and Pik-1/Pik-2 are unraveling how such protein pairs function together. Focusing on these model NLR pairs and other recent examples, this review highlights the distinctive features of NLR pairs and their various fascinating mode of action in pathogen effector perception. We also discuss how these findings on NLR pairs pave the way toward improved plant disease resistance.

The activity of the RGA5 sensor NLR from rice requires binding of its integrated HMA domain to effectors but not HMA domain self-interaction.

The rice nucleotide-binding (NB) and leucine-rich repeat (LRR) domain immune receptors (NLRs) RGA4 and RGA5 form a helper NLR/sensor NLR (hNLR/sNLR) pair that specifically recognizes the effectors AVR-Pia and AVR1-CO39 from the blast fungus Magnaporthe oryzae. While RGA4 contains only canonical NLR domains, RGA5 has an additional unconventional heavy metal-associated (HMA) domain integrated after its LRR domain. This RGA5HMA domain binds the effectors and is crucial for their recognition. Investigation of the three-dimensional structure of the AVR1-CO39/RGA5HMA complex by X-ray crystallography identified a candidate surface for effector binding in the HMA domain and showed that the HMA domain self-interacts in the absence of effector through the same surface. Here, we investigated the relevance of this HMA homodimerization for RGA5 function and the role of the RGA5HMA effector-binding and self-interaction surface in effector recognition. By analysing structure-informed point mutations in the RGA5HMA -binding surface in protein interaction studies and in Nicotiana benthamiana cell death assays, we found that HMA self-interaction does not contribute to RGA5 function. However, the effector-binding surface of RGA5HMA identified by X-ray crystallography is crucial for both in vitro and in vivo effector binding as well as effector recognition. These results support the current hypothesis that noncanonical integrated domains of NLRs act primarily as effector traps and deepen our understanding of the sNLRs' function within NLR pairs.

1H, 13C, 15 N backbone and side-chain NMR assignments for three MAX effectors from Magnaporthe oryzae.

Effectors are small and very diverse proteins secreted by fungi and translocated in plant cells during infection. Among them, MAX effectors (for Magnaporthe Avrs and ToxB) were identified as a family of effectors that share an identical fold topology despite having highly divergent sequences. They are mostly secreted by ascomycetes from the Magnaporthe genus, a fungus that causes the rice blast, a plant disease leading to huge crop losses. As rice is the first source of calories in many countries, especially in Asia and Africa, this constitutes a threat for world food security. Hence, a better understanding of these effectors, including structural and functional characterization, constitutes a strategic milestone in the fight against phytopathogen fungi and may give clues for the development of resistant varieties of rice. We report here the near complete 1H, 15 N and 13C NMR resonance assignment of three new putative MAX effectors (MAX47, MAX60 and MAX67). Secondary structure determination using TALOS-N and CSI.3 demonstrates a high content of ß-strands in all the three proteins, in agreement with the canonic ß-sandwich structure of MAX effectors. This preliminary study provides foundations for further structural characterization, that will help in turn to improve sequence predictions of other MAX effectors through data mining.

Combining High-Pressure NMR and Geometrical Sampling to Obtain a Full Topological Description of Protein Folding Landscapes: Application to the Folding of Two MAX Effectors from Magnaporthe oryzae.

Despite advances in experimental and computational methods, the mechanisms by which an unstructured polypeptide chain regains its unique three-dimensional structure remains one of the main puzzling questions in biology. Single-molecule techniques, ultra-fast perturbation and detection approaches and improvement in all-atom and coarse-grained simulation methods have greatly deepened our understanding of protein folding and the effects of environmental factors on folding landscape. However, a major challenge remains the detailed characterization of the protein folding landscape. Here, we used high hydrostatic pressure 2D NMR spectroscopy to obtain high-resolution experimental structural information in a site-specific manner across the polypeptide sequence and along the folding reaction coordinate. We used this residue-specific information to constrain Cyana3 calculations, in order to obtain a topological description of the entire folding landscape. This approach was used to describe the conformers populating the folding landscape of two small globular proteins, AVR-Pia and AVR-Pib, that belong to the structurally conserved but sequence-unrelated MAX effectors superfamily. Comparing the two folding landscapes, we found that, in spite of their divergent sequences, the folding pathway of these two proteins involves a similar, inescapable, folding intermediate, even if, statistically, the routes used are different.

New recognition specificity in a plant immune receptor by molecular engineering of its integrated domain.

Plant nucleotide-binding and leucine-rich repeat domain proteins (NLRs) are immune sensors that recognize pathogen effectors. Here, we show that molecular engineering of the integrated decoy domain (ID) of an NLR can extend its recognition spectrum to a new effector. We relied for this on detailed knowledge on the recognition of the Magnaporthe oryzae effectors AVR-PikD, AVR-Pia, and AVR1-CO39 by, respectively, the rice NLRs Pikp-1 and RGA5. Both receptors detect their effectors through physical binding to their HMA (Heavy Metal-Associated) IDs. By introducing into RGA5_HMA the AVR-PikD binding residues of Pikp-1_HMA, we create a high-affinity binding surface for this effector. RGA5 variants carrying this engineered binding surface perceive the new ligand, AVR-PikD, and still recognize AVR-Pia and AVR1-CO39 in the model plant N. benthamiana. However, they do not confer extended disease resistance specificity against M. oryzae in transgenic rice plants. Altogether, our study provides a proof of concept for the design of new effector recognition specificities in NLRs through molecular engineering of IDs.

A novel robust and high-throughput method to measure cell death in Nicotiana benthamiana leaves by fluorescence imaging.

Assessing immune responses and cell death in Nicotiana benthamiana leaf agro-infiltration assays is a powerful and widely used experimental approach in molecular plant pathology. Here, we describe a reliable high-throughput protocol to quantify strong, macroscopically visible cell death responses in N. benthamiana agro-infiltration assays. The method relies on measuring the reduction of leaf autofluorescence in the red spectrum upon cell death induction and provides quantitative data suitable for straightforward statistical analysis. Two different well-established model nucleotide-binding and leucine-rich repeat domain proteins (NLRs) were used to ensure the genericity of the approach. Its accuracy and versatility were compared to visual scoring of the cell death response and standard methods commonly used to characterize NLR activities in N. benthamiana. A discussion of the advantages and limitations of our method compared to other protocols demonstrates its robustness and versatility and provides an effective means to select the best-suited protocol for a defined experiment.

Precision Breeding Made Real with CRISPR: Illustration through Genetic Resistance to Pathogens.

Since its discovery as a bacterial adaptive immune system and its development for genome editing in eukaryotes, the CRISPR technology has revolutionized plant research and precision crop breeding. The CRISPR toolbox holds great promise in the production of crops with genetic disease resistance to increase agriculture resilience and reduce chemical crop protection with a strong impact on the environment and public health. In this review, we provide an extensive overview on recent breakthroughs in CRISPR technology, including the newly developed prime editing system that allows precision gene editing in plants. We present how each CRISPR tool can be selected for optimal use in accordance with its specific strengths and limitations, and illustrate how the CRISPR toolbox can foster the development of genetically pathogen-resistant crops for sustainable agriculture.