Search PubMedSearch

SEARCH · Search PubMed

Results for “plant Immunity”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

High-throughput glycan array screening reveals rhamnogalacturonan-I as a ligand for Arabidopsis leucine-rich repeat receptor kinases involved in plant immunity.

The plant cell wall not only serves as a physical barrier against pathogens but, when damaged, also functions as a source of cell wall-derived molecules that play crucial roles in plant immunity as damage-associated molecular patterns. While oligogalacturonides from homogalacturonan are well-studied damage-associated molecular patterns, the immune-signaling potential of other cell wall components remains largely unexplored. Conventional genetic and biochemical approaches aimed at identifying ligand-receptor pairs in plant immunity have been limited by the vast diversity of potential ligand molecules and functional redundancy of putative receptors. In this study, we developed a high-throughput screening pipeline that simultaneously examines multiple interactions between plant cell wall-derived glycans and >350 extracellular domains of receptor kinases and receptor-like proteins in Arabidopsis, resulting in the screening of >40 000 interactions. We discovered a group of leucine-rich repeat receptor kinases named ARMs (AWARENESS of RG-I MAINTENANCES) that interact with rhamnogalacturonan-I (RG-I), a major component of pectin. RG-I treatment induced pattern-triggered immunity responses with distinct kinetics compared to oligogalacturonide responses. We identified RG-I oligosaccharide structures required for interaction with ARM receptors and immune activation and found that ARM receptors function redundantly in plant immunity. Collectively, our work provides a powerful platform for discovering glycan-receptor pairs in plants, facilitating a more comprehensive understanding of cell wall surveillance mechanisms in plant immunity.

Arabidopsis

A CsWRKY46-CsPBL9-CsARI1 tripartite regulatory module coordinates H2O2 production and callose deposition in citrus fruit immunity.

Plant immunity against pathogens involves multiple immune responses and intricate regulatory networks. However, how immune networks are deployed in fruit remains poorly understood. Here, we show that citrus fruit immune responses, including hydrogen peroxide (H2O2) production and callose deposition, are multiply regulated by transcriptional activation, phosphorylation, and ubiquitination. Citrus sinensis genes encoding nicotinamide adenine dinucleotide phosphate (NADPH) oxidase CsRBOHG and callose synthase CsCalS5, responsible for H2O2 production and callose deposition, respectively, are transcriptionally activated by CsWRKY46. Phosphorylation-enhanced activity of CsRBOHG by CsPBL9 enhances immunity. RING1-IBR-RING2 (RBR)-type E3 ligase CsARI1, acting as an immune brake, ubiquitinates CsRBOHG and CsCalS5 for degradation. Interestingly, CsARI1 also shows a moonlight function wherein it interacts with CsPBL9 in a non-ubiquitination manner, disrupting CsPBL9's interaction with CsRBOHG. This CsARI1-CsPBL9 interaction is stimulated by H2O2 as feedback. Moreover, H2O2 contributes to callose deposition, indicating an interplay between two immune responses. Our study reveals a tripartite regulatory hub orchestrating self-linked immunity in citrus fruit.

CP: plants

Pan-analysis of intra- and inter-species diversity reveals a group of highly variable immune receptor genes in rice.

Plant immune receptors and their natural variations play a central role in combating disease-causing pathogens. These immune receptors include intracellular nucleotide-binding leucine-rich repeat (LRR) receptors (NLRs) and cell-surface pattern recognition receptors (PRRs) that can be further classified as receptor-like proteins (RLPs) and receptor-like kinases (RLKs). Although the NLRome has been characterized, the repertoire and extent of diversity of PRRome remain undetermined in rice. In this study, we examined the diversity of immune receptor genes using high-quality genomes of 309 rice accessions from 8 species within the genus Oryza. A total of 376 310 immune receptor genes were identified, including 149 592 NLR-coding genes and 226 718 PRR coding genes. Shannon entropy analysis revealed a set of immune receptors that display significant intra-species and inter-species diversity in rice. In general, RLPs are more variable than RLKs, while NLRs and LRR-RLPs are more variable than LRR-RLKs. Additionally, NLR and PRR genes exhibit contrasting shoot/root expression patterns, with NLRs generally skewed towards root expression. Furthermore, we found that the size of the LRR-RLK gene families correlates with local annual precipitation, suggesting a stronger selection pressure on LRR-RLK genes in rice accessions grown under wet conditions than dry conditions. In sum, this pan-genomic analysis not only reveals the extensive diversity of the immune receptor repertoires in rice but also provides potential target genes for improving disease resistance in rice.

Oryza

The bZIP54 (GBF2)-SARD1 module regulates salicylic acid-mediated resistance to Pst DC3000 in Arabidopsis.

Salicylic acid (SA)-mediated defense responses are crucial for plant immunity, yet transcription factors (TFs) that coordinate SA biosynthesis with immune activation remain incompletely characterized. Here, a basic leucine zipper (bZIP) TF, bZIP54, was identified as a positive regulator in response to Pseudomonas syringae pv. tomato (Pst) DC3000. Consistent with this finding, bZIP54 regulated SA accumulation and a suite of SA-related defense genes following Pst DC3000 infection. Mechanistically, bZIP54 directly bound to a G-box-like motif in the SARD1 promoter, activating its expression-an interaction that was further enhanced by SA. Genetic analysis demonstrated that SARD1 operates downstream of bZIP54 to confer resistance to Pst DC3000. Additionally, bZIP54 also contributed to defense against the fungal pathogen Sclerotinia sclerotiorum, indicating a broader role in plant immunity. Together, these findings revealed a bZIP54-SARD1 regulatory module, thus providing insights into the transcriptional networks governing disease resistance in Arabidopsis.

Arabidopsis

Microbial partnerships and molecular mechanisms in plant stress physiology for climate-resilient and sustainable farming.

Plant-microbial partnerships and their underlying molecular mechanisms are indispensable, natural drivers of improved nutrient acquisition and stress tolerance in the face of climate-driven environmental challenges. Modern multi-omics tools, when coupled with artificial intelligence and synthetic biology, enable the precise design of targeted bioinoculants and synthetic microbial consortia. Translating these advanced microbiome-based strategies into scalable, field-level agricultural applications provides a sustainable path toward securing global food production while maintaining soil health. Global climate change imposes multifaceted abiotic and biotic stresses on crops, disrupting physiological and molecular processes and threatening agricultural productivity. Plant-associated microbes represent an underexplored yet powerful ally in enhancing crop resilience. This review presents current knowledge of plant-microbe interactions and the molecular mechanisms governing plant stress physiology, with an emphasis on climate-resilient and sustainable farming. Hence, ever-changing environmental cues pose a significant burden on agricultural productivity, and plant-associated microbial communities modulate a cascade of physiological and molecular responses, including production of phytohormones, signaling, regulation of reactive oxygen species homeostasis, and activation of plant immune responses to help plants withstand stress and enhance productivity. Moreover, root exudates, phytohormones, and quorum sensing mediate the central communication networks, facilitating plant-microbe cross talk. Additionally, the advances in OMICs approaches aid in disentangling the molecular underpinnings of these interactions by providing mechanistic insights and potential candidate gene targets for crop improvement and stress resilience. In the post-genomic era, integrating artificial intelligence and big data analysis to optimize microbiome-based strategies for sustainable agriculture is a new frontier for disentangling plant-microbe symbiosis to improve soil health, enhance crop yields, and improve stress tolerance. Thus, by integrating the ecological, physiological, and molecular perspectives, this review highlights the transformative potential of harnessing plant-microbe symbiosis for climate-resilient and sustainable agriculture.

Stress, Physiological

Integrated Functional Characterization of Hemileia vastatrix Effector Candidates Reveals Coordinated Immune Suppression, Sequential Deployment and Compartment-Specific Targeting.

Coffee leaf rust, caused by the obligate biotrophic fungus Hemileia vastatrix, remains the most destructive disease of coffee worldwide. Although genomic and transcriptomic studies have identified a large number of candidate effectors, experimental evidence supporting their biological roles during infection remains limited. Here, we integrated functional assays, temporal expression profiling during coffee infection and subcellular localization analyses to investigate the biological properties of 44 H. vastatrix effector candidates (HvECs). Using the Pseudomonas fluorescens EtHAn effector delivery system in Nicotiana benthamiana, 15 HvECs consistently suppressed pattern-triggered immunity (PTI), indicating that immune suppression is a widespread property among the H. vastatrix effector repertoire, as assessed in this heterologous system. Five HvECs also attenuated AvrB-triggered effector-triggered immunity (ETI), and three suppressed both PTI and ETI, suggesting that a subset of HvECs targets conserved regulatory nodes shared by these interconnected immune pathways. Temporal expression profiling revealed sequential deployment of HvECs throughout infection, with distinct subsets predominating during pre-biotrophic development, host penetration or biotrophic colonization, consistent with stage-specific functions during fungal pathogenesis. Subcellular localization analyses further showed that HvECs preferentially accumulated in the nucleus or chloroplasts, compartments known as central hubs of plant immune regulation. This study provides the most comprehensive functional characterization of H. vastatrix effector candidates to date, establishes a biologically informed framework for prioritizing candidates for future identification of avirulence determinants recognized by SH resistance genes, and advances our understanding of how the coffee rust fungus orchestrates immune suppression across time and cellular space during pathogenesis.

Nicotiana

Genomic analysis of regulatory mechanisms governing EPS66A biosynthesis in Streptomyces changanensis HL-66.

Streptomyces changanensis HL-66 produces the α-(1,4)/(1,6)-glucan exopolysaccharide EPS66A, a potent plant immune elicitor with promising applications in plant protection. However, its low native fermentation yield limits large-scale application. To investigate the biosynthetic potential and regulatory mechanisms underlying EPS66A production, the whole genome of HL-66 was sequenced and analyzed. The HL-66 genome is 6.82 Mb in size, with a GC content of 74%, and encodes 6081 predicted functional genes. Among these, 1390 genes were annotated to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, 4187 were assigned to Gene Ontology (GO) terms, and 143 were classified into Clusters of Orthologous Groups (COG) categories. antiSMASH analysis identified 22 secondary metabolite biosynthetic gene clusters, including multiple polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) clusters. Functional analyses revealed that the glycosyltransferase gene (GTy) and the global regulatory gene (bldD) are involved in EPS66A biosynthesis. bldD is involved in morphological development and EPS66A production, whereas GTy specifically regulates EPS66A production without affecting growth or development. In both in vivo and potted-plant experiments, EPS66A (200 μg/mL) significantly reduced the severity of tobacco mosaic virus, apple anthracnose leaf spot, walnut bacterial leaf spot, and jujube anthracnose, achieving control efficacies of 90.21%, 87.95%, 77.41%, and 68.55%, respectively, and outperforming a commercial chitosan oligosaccharide control. These findings provide new insights into the genetic architecture and regulatory mechanisms of EPS66A biosynthesis and support its development as a polysaccharide-based green pesticide.

Streptomyces

Myosin XI-mediated BIK1 recruitment to nanodomains facilitates FLS2-BIK1 complex formation during innate immunity in Arabidopsis.

Plants rely on immune receptor complexes at the cell surface to perceive microbial molecules and transduce these signals into the cell to regulate immunity. Various immune receptors and associated proteins are often dynamically distributed in specific nanodomains on the plasma membrane (PM). However, the exact molecular mechanism and functional relevance of this nanodomain targeting in plant immunity regulation remain largely unknown. By utilizing high spatiotemporal resolution imaging and single-particle tracking analysis, we show that myosin XIK interacts with remorin to recruit and stabilize PM-associated kinase BOTRYTIS-INDUCED KINASE 1 (BIK1) within immune receptor FLAGELLIN SENSING 2 (FLS2)-containing nanodomains. This recruitment facilitates FLS2/BIK1 complex formation, leading to the full activation of BIK1-dependent defense responses upon ligand perception. Collectively, our findings provide compelling evidence that myosin XI functions as a molecular scaffold to enable a spatially confined complex assembly within nanodomains. This ensures the presence of a sufficient quantity of preformed immune receptor complex for efficient signaling transduction from the cell surface.

Arabidopsis

Leaf Rust in Rye: From Pathogen Biology to Host Defense and Resistance Breeding.

Leaf rust (LR), caused by Puccinia recondita f. sp. secalis (Prs), is considered one of the most dangerous rye (Secale cereale L.) diseases, causing yield losses exceeding 35%. This review summarizes all currently available data about this disease: pathogen characteristics (including its life cycle, natural variation, and disease symptoms), resistance resources, and the background of the plant immune response at the genome, transcriptome, and metabolome levels. The research conducted so far has allowed for the identification of dozens of genes that play a significant role in the rye immune response to Prs infection. Among them, genes encoding NBS-LRR proteins (including SECCE1Rv1G0014220, the most likely Pr3 candidate), glycosyltransferase, β-1,3-glucanase, 1-deoxy-D-xylulose 5-phosphate synthase, β-1,3-glucanase, UDP-glycosyltransferase, pathogenesis-related protein 1, ammonium transporter, and cytochrome P450 enzymes are candidates for seedling and all-stage resistance, whereas ScLr_ABC25 currently represents the most promising candidate associated with adult-plant resistance. Among the metabolites differentially accumulated in response to Prs, those related to phenylpropanoids, diterpenoids, and thiamine branches seem to play the most important role in the immune response. Finally, we suggest how the knowledge acquired so far about the rye-Prs interaction can be used in modern breeding programs aimed at obtaining cultivars with enhanced resistance to LR, such as through the use of functional gene markers and/or metabolic biomarker-assisted selection and, in the more distant future, by developing and applying new genomic techniques for precise editing of resistance and susceptibility genes, engineering synthetic immune receptors and decoys, and pan-genomic exploration for identification of rare or lineage-specific resistance alleles. [Formula: see text] Copyright © 2026 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Plant Diseases

CRISPR/Cas9-driven double modification of grapevine MLO6-7 imparts powdery mildew resistance, while editing of NPR3 augments powdery and downy mildew tolerance.

The implementation of genome editing strategies in grapevine is the easiest way to improve sustainability and resilience while preserving the original genotype. Among others, the Mildew Locus-O (MLO) genes have already been reported as good candidates to develop powdery mildew-immune plants. A never-explored grapevine target is NPR3, a negative regulator of the systemic acquired resistance. We report the exploitation of a cisgenic approach with the Cre-lox recombinase technology to generate grapevine-edited plants with the potential to be transgene-free while preserving their original genetic background. The characterization of three edited lines for each target demonstrated immunity development against Erysiphe necator in MLO6-7-edited plants. Concomitantly, a significant improvement of resilience, associated with increased leaf thickness and specific biochemical responses, was observed in defective NPR3 lines against E. necator and Plasmopara viticola. Transcriptomic analysis revealed that both MLO6-7 and NPR3 defective lines modulated their gene expression profiles, pointing to distinct though partially overlapping responses. Furthermore, targeted metabolite analysis highlighted an overaccumulation of stilbenes coupled with an improved oxidative scavenging potential in both editing targets, likely protecting the MLO6-7 mutants from detrimental pleiotropic effects. Finally, the Cre-loxP approach allowed the recovery of one MLO6-7 edited plant with the complete removal of transgene. Taken together, our achievements provide a comprehensive understanding of the molecular and biochemical adjustments occurring in double MLO-defective grape plants. In parallel, the potential of NPR3 mutants for multiple purposes has been demonstrated, raising new questions on its wide role in orchestrating biotic stress responses.

Vitis

The genome of Lespedeza potaninii reveals biased subgenome evolution and drought adaptation.

Lespedeza potaninii, a xerophytic subshrub belonging to the legume family, is native to the Tengger Desert and is highly adapted to drought. It has important ecological value due to its drought adaptability, but the underlying molecular mechanisms remain largely unknown. Here, we report a 1.24 Gb chromosome-scale assembly of the L. potaninii genome (contig N50 = 15.75 Mb). Our results indicate that L. potaninii underwent an allopolyploid event with 2 subgenomes, A and B, presenting asymmetric evolution and B subgenome dominance. We estimate that the 2 diploid progenitors of L. potaninii diverged around 3.6 million years ago (MYA) and merged around 1.0 MYA. We revealed that the expansion of hub genes associated with drought responses, such as the binding partner 1 of accelerated cell death 11 (ACD11) (BPA1), facilitated environmental adaptations of L. potaninii to desert habitats. We found a novel function of the BPA1 family in abiotic stress tolerance in addition to the known role in regulating the plant immune response, which could improve drought tolerance by positively regulating reactive oxygen species homeostasis in plants. We revealed that bZIP transcription factors could bind to the BPA1 promoter and activate its transcription. Our work fills the genomic data gap in the Lespedeza genus and the tribe Desmodieae, which should provide theoretical support both in the study of drought tolerance and in the molecular breeding of legume crops.

Genome, Plant

Revealing Functional Traits of Insect Pest Suppressive Rhizobacterial Strains Through Comparative Genomics.

Root inoculation with rhizobacteria is an emerging strategy to enhance plant resistance to aphid herbivory, yet the microbial functional traits underpinning these responses remain poorly characterised. Here, we present a comparative genomic analysis of five rhizobacteria (Acidovorax radicis N35, Bacillus subtilis B171, Bacillus velezensis FZB42, Rhizobium radiobacter F4 and Pseudomonas simiae WCS417r) that suppress aphids when inoculated onto barley. As expected, functional variation largely reflected phylogenetic relatedness; however, candidate traits implicated in modulation of plant immune defences were conserved across all strains, including biosynthesis of 2,3-butanediol, riboflavin and salicylic acid. Additional shared functions, linked to plant defence signalling, included phytoene and squalene biosynthesis (absent in P. simiae) and N-acyl homoserine lactone quorum sensing (absent in Bacillus spp.). Strain-specific traits were also identified, including surfactin production in Bacillus spp. and hydrogen cyanide biosynthesis in A. radicis and P. simiae. Comparison with a broader collection of rhizobacteria revealed that many putative plant-beneficial functions identified were widely conserved, including among closely related phytopathogens. This extensive functional overlap suggests aphid suppression cannot be explained solely by presence or absence of broad functional traits, but rather by specific trait combinations, regulatory differences, or context-dependent expression. This highlights the need for genome-informed approaches for bioinoculant discovery.

Animals

A cooperative regulatory module between TAGL2 and JMJC1 activates specific defense genes against root-knot nematodes in tomato.

Plant-parasitic nematodes (PPNs) threaten global food security. Although epigenetic modifications are crucial for plant immunity, how histone modifiers contribute to root-knot nematodes (RKNs, Meloidogyne incognita) resistance remains unclear. Here, using genetic, molecular and biochemical approaches, we investigated the epigenetic and transcriptional mechanisms underlying RKN resistance mediated by the histone demethylase (HDM) JMJC1 and the MADS-box transcription factor TAGL2 in tomato (Solanum lycopersicum). We identified JMJC1 as an RKN-induced positive defense regulator targeting H3K9me3 and H3K27me3 histone marks. JMJC1 physically interacts with TAGL2, which also positively regulates RKN resistance. Transcriptomic analysis indicated that TAGL2 regulates multiple layers of the plant defense network, transcriptionally activating representative genes from distinct pathways (including PUB10, bHLH98, CCaMK, and SAUR3), which we validated as positive regulators of RKN resistance via virus-induced gene silencing (VIGS). At the chromatin level, TAGL2 and JMJC1 co-regulate these loci, associating with localized H3K9me3 and H3K27me3 reduction. Furthermore, TAGL2 directly activates JMJC1 transcription, establishing a positive feedback loop that amplifies immune signaling. Our findings reveal a cooperative model wherein a HDM and a transcription factor coordinate at specific loci to fine-tune multiple defense layers at both epigenetic and transcriptional levels, providing insights for breeding durable nematode-resistant plants.

Solanum lycopersicum

Rapidly evolving aphid gall effector proteins exhibit saposin-like folds.

Many insects manipulate plants by injecting effector proteins. In one extreme example of this molecular "hijacking", Hormaphis cornu aphids inject bicycle proteins into Hamamelis virginiana (Witch Hazel), contributing to the development of novel organs called galls. Bicycle proteins share no amino acid sequence similarity with proteins of known function. Here, we report the crystal structures of two divergent bicycle proteins. Both proteins contain saposin-like folds: one with multiple disulfide bonds exhibits a helix swap; the other has no disulfide bonds and possesses two tandem domains. To explore the structural evolution of bicycle proteins, we predicted bicycle protein structures with Alphafold2 (AF2). While AF2 did not recover the two experimental structures using existing databases, it succeeded after we provided multiple sequence alignments (MSAs) containing protein sequences encoded in new genome sequences from closely related aphid species. Using this customized approach at scale, we generated 2400 high-confidence predictions for bicycle proteins from seven aphid species. This dataset revealed that bicycle proteins without cysteines are outliers in fold space and appear to have evolved from ancestral proteins with disulfide-bonded saposin-like folds. While all bicycle proteins contain predicted saposin-like folds, they display a vast diversity of structural and physicochemical properties. While this diversity thwarts prediction of conserved functions encoded in structure, it suggests that bicycle proteins have evolved to target diverse plant processes and/or to evade plant immune surveillance.

AlphaFold predictions

Adaptation of the Cyst Nematode Globodera pallida to the Colinear Potato Resistant QTLs GpaVvrn and GpaVspl Involved Distinct Genomic Regions and Absence of Cross-Virulence.

The use of alternative methods to control cyst nematode populations has accelerated since the ban of chemical nematicides in Europe. The resistant QTL GpaVvrn, derived from the wild species Solanum vernei, is widely present in resistant European potato cultivars and provides strong protection against Globodera pallida populations although a risk of resistance breakdown has already been demonstrated in both experimental evolution studies and field populations. The wild relative S. sparsipilum, harbouring the resistant QTL GpaVspl, would be an interesting alternative source of resistance to control virulent G. pallida. The goal of the present study was to understand the genomics of adaptation of the nematode to these two colinear resistant QTLs. Starting with two natural populations, an experimental evolution approach allowed, after 10 generations on resistant potato genotypes, selecting independent nematode lineages adapted to each QTL. These virulent lineages were analysed through a combination of phenotyping and genome scans approaches. Phenotyping enabled the quantification of virulence levels and confirmed resistance breakdowns. Pool-Seq whole genome sequencing followed by genome scan analyses identified genomic regions under selection, potentially involved in the adaptive mechanisms to each resistance factor. Candidate genes within these regions provided insights into the genetic basis of adaptation, revealing effectors known to suppress plant immunity. As genome scans highlighted distinct genomic regions for the adaptation to both resistant factors, we were able to predict and phenotypically confirm the absence of cross-virulence between nematode lineages evolving on GpaVvrn and GpaVspl. These findings have significant implications for the design of effective and sustainable resistance management strategies.

Animals

Dual regulation of the receptor-like kinase BIR1 involves site-directed transcript cleavage and 5'-leader-mediated translational control.

In Arabidopsis, BRASSINOSTEROID INSENSITIVE1-ASSOCIATED RECEPTOR KINASE 1 (BAK1)-INTERACTING RECEPTOR-LIKE KINASE 1 (BIR1) is a negative regulator of plant immunity and cell death. BIR1 was earlier described as a target of epigenetic and post-transcriptional degradation. During virus infections, degradome analysis of BIR1 transcripts mapped predominant mRNA cleavage sites at the 5'-untranslated leader region (site A) and the protein-coding sequence (sites B and C). Here, we identified another virus-associated cleavage site (D) within the BIR1 coding region and investigated the contribution of site-directed mRNA cleavage to BIR1 regulation. Mutations at B, C, and D sites enhanced mRNA stability by impairing transcript cleavage, resulting in increased BIR1 mRNA and protein accumulation. This regulation is disrupted in RNA silencing mutants, supporting a model of cis-directed small interfering RNA (siRNA)-mediated degradation. We next demonstrate that virus infection reduces BIR1 translation in Arabidopsis. Furthermore, our data reveal a repressive role for the 5'-leader in regulating BIR1 translation, potentially mediated by upstream open reading frames (uORFs) and a virus-responsive long non-coding RNA (lncRNA) derived from the natural antisense At4g39838 locus. Together, these findings reveal a multilayered regulatory mechanism that integrates sRNA-mediated cleavage with translational control, with broader implications for the fine-tuning of stress-responsive gene expression during infection.

Arabidopsis

Rapidly evolving aphid gall effector proteins exhibit saposin-like folds.

Many insects manipulate plants by injecting effector proteins. In one extreme example of this molecular "hijacking," Hormaphis cornu aphids inject bicycle proteins into Hamamelis virginiana, contributing to the development of novel organs called galls. Bicycle proteins share no amino acid sequence similarity with proteins of known function. Here, we report the crystal structures of two divergent bicycle proteins. Both proteins contain saposin-like folds: one with multiple disulfide bonds exhibits a swapped domain topology; the other has no disulfide bonds and possesses two distinct, tandem domains. To explore the structural evolution of bicycle proteins, we attempted to predict bicycle protein structures with Alphafold2 (AF2) and other deep learning programs. While AF2 did not recover the two experimental structures using existing databases, it succeeded when provided with multiple sequence alignments (MSAs) of protein sequences from newly sequenced closely related species. Using this approach, we generated 2,400 high-confidence bicycle protein predictions from seven aphid species. While all aphid bicycle proteins contain predicted saposin-like folds, they display a vast diversity of structural and physicochemical properties. While this diversity thwarts prediction of conserved functions encoded in structure, it suggests that bicycle proteins have evolved to target diverse plant processes and/or to evade plant immune surveillance. Our extension of AF2 with custom MSAs of proteins from closely related species provides a generalizable, powerful approach for predicting structures of rapidly evolving protein families.

Animals