# Patterns of Sequence and Expression Diversification Associate Members of the PADRE Gene Family With Response to Fungal Pathogens

Marie Didelon, Mehdi Khafif, Laurence Godiard, Adelin Barbacci, Sylvain Raffaele

► To cite this version:

Marie Didelon, Mehdi Khafif, Laurence Godiard, Adelin Barbacci, Sylvain Raffaele. Patterns of Sequence and Expression Diversification Associate Members of the PADRE Gene Family With Response to Fungal Pathogens. Frontiers in Genetics, 2020, 11, pp.1-16. 10.3389/fgene.2020.00491. hal-03051478

HAL Id: hal-03051478

https://hal.inrae.fr/hal-03051478

Submitted on 10 Dec 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Patterns of Sequence and Expression Diversification Associate Members of the PADRE Gene Family With Response to Fungal Pathogens

Marie Didelon, Mehdi Khafif, Laurence Godiard, Adelin Barbacci and Sylvain Raffaele*

Université de Toulouse, Laboratoire des Interactions Plantes Micro-organismes (LIPM), Institut National de Recherche pour l'Agriculture, l'Alimentation et l'Environnement (INRAE) - Centre National de la Recherche Scientifique (CNRS), Castanet-Tolosan, France

Pathogen infection triggers extensive reprogramming of the plant transcriptome, including numerous genes the function of which is unknown. Due to their wide taxonomic distribution, genes encoding proteins with Domains of Unknown Function (DUFs) activated upon pathogen challenge likely play important roles in disease. In *Arabidopsis thaliana*, we identified thirteen genes harboring a DUF4228 domain in the top 10% most induced genes after infection by the fungal pathogen *Sclerotinia sclerotiorum*. Based on functional information collected through homology and contextual searches, we propose to refer to this domain as the pathogen and abiotic stress response, cadmium tolerance, disordered region-containing (PADRE) domain. Genome-wide and phylogenetic analyses indicated that PADRE is specific to plants and diversified into 10 subfamilies early in the evolution of Angiosperms. PADRE typically occurs in small single-domain proteins with a bipartite architecture. PADRE N-terminus harbors conserved sequence motifs, while its C-terminus includes an intrinsically disordered region with multiple phosphorylation sites. A pangenomic survey of PADRE genes expression upon S. sclerotiorum inoculation in *Arabidopsis*, castor bean, and tomato indicated consistent expression across species within phylogenetic groups. Multi-stress expression profiling and co-expression network analyses associated AtPADRE genes with the induction of anthocyanin biosynthesis and responses to chitin and to hypoxia. Our analyses reveal patterns of sequence and expression diversification consistent with the evolution of a role in disease resistance for an uncharacterized family of plant genes. These findings highlight PADRE genes as prime candidates for the functional dissection of mechanisms underlying plant disease resistance to fungi.

**Keywords:** plant disease resistance, diversification, DUF4228, intrinsic disorder, pathogenesis-related, gene expression profiling

## INTRODUCTION

Wild plants and crops suffer from recurrent attacks by pathogenic microbes, threatening biodiversity and food production. Molecular and genetic studies revealed that plants possess an elaborate immune system able to detect pathogens and activate genetic pathways to mount effective defense responses (Dodds and Rathjen, 2010). Specific defense responses allow plants

**OPEN ACCESS**

**Edited by:**

Takaki Maekawa,
Max Planck Institute for Plant
Breeding Research, Germany

**Reviewed by:**

Thomas Griebel,
Freie Universität Berlin, Germany
Bilal Okmen,
University of Cologne, Germany
Andrea Ghelfi,
Kazusa DNA Research Institute,
Japan

**\*Correspondence:**

Sylvain Raffaele
sylvain.raffaele@inrae.fr

**Specialty section:**

This article was submitted to
Evolutionary and Population Genetics,
a section of the journal
Frontiers in Genetics

Received: 15 January 2020

Accepted: 20 April 2020

Published: 29 May 2020

**Citation:**

Didelon M, Khafif M, Godiard L, Barbacci A and Raffaele S (2020) Patterns of Sequence and Expression Diversification Associate Members of the PADRE Gene Family With Response to Fungal Pathogens. *Front. Genet.* 11:491. doi: 10.3389/fgene.2020.00491
to cope with microbial pathogens of diverse lifestyles and genotypes that target diverse plant organs (Glazebrook, 2005). In most cases, the activation of plant responses requires extensive transcriptional reprogramming, covering for instance up to 25% of the whole genome in *Arabidopsis thaliana* (Eulgem, 2005). In nature, one of the most frequent forms of plant immunity is designated as quantitative disease resistance (QDR) (Poland et al., 2009; Roux et al., 2014). QDR leads to a full continuum of disease resistance phenotypes in natural plant populations, from very susceptible to largely resistant, and generally involves a large number of genetic loci. Every gene adds a small contribution to form the overall resistance (Roux et al., 2014). Current knowledge of the molecular bases of QDR in plants remains very incomplete, but a few general properties have emerged. First, the molecular functions of QDR genes are very diverse, including for instance transporters (Krattinger et al., 2009), kinases (Derbyshire et al., 2019), proteases (Badet et al., 2017), and genes of unknown function (Fukuoka et al., 2009). Second, the function of QDR genes may not be limited to disease resistance and can include activity in cell morphology (Rajarammohan et al., 2018; Badet et al., 2019), metabolism (Rajarammohan et al., 2018), or embryogenesis (Derbyshire et al., 2019) in certain contexts. Third, QDR responses to a given pathogen species may involve hundreds or even thousands of genes (Corwin et al., 2016; Fordyce et al., 2018). Therefore, pathogen infection triggers extensive reprogramming of the plant transcriptome, including numerous genes the molecular function of which is currently unknown.

Recent progress in high-throughput omics techniques enabled the determination of the sequence of genes and proteins at an unprecedented pace. Homology relationships allow to rapidly transfer functional information from one sequence to another but suffer limitations (Pearson and Sierk, 2005), and our capacity to generate new sequences far exceeds our ability to interpret them. Sequence conservation across large evolutionary distances can identify previously unknown functional domains in proteins, such as in the case of the VAST domain (PF16016) (Khafif et al., 2014, 2017; Gatta et al., 2015). The Protein Family Database (Pfam) gathers protein families by their homology of sequence (El-Gebali et al., 2019). In 2019, the latest Pfam release (32.0) counted 17,929 entries, 3,961 (22%) of them being Domains of Unknown Function (DUFs). DUFs are protein families for which no member has an experimentally characterized function. Systematic structural analyses of DUF proteins revealed that a significant part of DUF proteins likely originate from extreme diversification and neofunctionalization of known protein domains (Jaroszewski et al., 2009). Due to their wide taxonomic distribution and their evolutionary sequence conservation, many DUFs are expected to compose essential proteins (Goodacre et al., 2013). Widely distributed genes encoding proteins with DUFs activated upon pathogen challenge are promising sources of new insight into the evolution and molecular mechanisms of plant disease resistance.

*Sclerotinia sclerotiorum* is a devastating fungal plant pathogen from the Ascomycota division with a necrotrophic lifestyle. It is responsible for the white and stem mold diseases on more than 400 plant species, including crops of high agricultural value like

sunflower, soybean, rapeseed, and tomato, among others (Boland and Hall, 1994; Hegedus and Rimmer, 2005). The host range of *S. sclerotiorum* also includes plants from the *Brassicaceae* family, such as the plant model *A. thaliana*. Resistance to *S. sclerotiorum* is typically quantitative with no complete resistance (Perchepied et al., 2010; Mbengue et al., 2016). The molecular bases of QDR to *S. sclerotiorum* are beginning to be elucidated, notably thanks to studies on *A. thaliana*, but remain very patchy (Mbengue et al., 2016). Global gene expression profiling by RNA sequencing revealed 4,703 *A. thaliana* genes significantly induced upon leaf inoculation with *S. sclerotiorum* (Badet et al., 2017), including several genes harboring a DUF4228 domain. Here, we took a survey of DUF4228 homologs across the plant kingdom and identified a few experimental insights into the function of these genes. We propose to refer to this domain as the pathogen and abiotic stress response, cadmium tolerance, disordered region-containing (PADRE) domain to facilitate future reference. We used phylogenetic analyses to document the extent diversity of PADRE sequences and infer scenarios for their evolution. PADRE proteins lack sequence similarity to characterized proteins but harbor a bipartite architecture with conserved motifs in the N-terminal region and a C-terminal region rich in phosphorylated residues and predicted to be intrinsically disordered. Pangenomic expression profiling in thale cress (*A. thaliana*), tomato (*Solanum lycopersicum*), and castor bean (*Ricinus communis*) plants inoculated by *S. sclerotiorum* identified groups of PADRE genes that respond to this fungal pathogen in a consistent manner across species. Finally, *AtPADRE* gene expression upon diverse stress treatments and co-expression network reconstruction suggests that several PADRE genes could function synergistically in plant defense. Our study reveals that responsiveness to fungal pathogen attack is conserved at the interspecific level in groups of PADRE genes and provides insights into the evolutionary history and functional diversification in this poorly characterized plant gene family.

## RESULTS

### Genes From the DUF4228 Family Are Over-Represented Among Genes Induced Upon *S. sclerotiorum* Inoculation

To get insights into plant processes activated during colonization by the fungal pathogen *S. sclerotiorum*, we analyzed RNA-Seq data for *A. thaliana* plants inoculated by *S. sclerotiorum*. Specifically, we focused on protein domains overrepresented among plant genes differentially expressed upon inoculation. To this end, we exploited the RNA sequencing data generated in Badet et al. (2017) (GSE106811). Differential expression analysis identified 4,703 genes significantly induced (log₂ fold change (LFC) > 1.5, adjusted *p*-value (padj) < 0.01) and 5,812 genes significantly down-regulated (LFC < 1.5, padj < 0.01) in *A. thaliana* during infection by *S. sclerotiorum*. We annotated genes by their protein domains using the Pfam database. Using a proportion Z-test (*p*-value < 0.01), we counted 53 protein
domains significantly overrepresented among induced genes with at least 10 occurrences in *A. thaliana* genome (Supplementary Table S1 and Figure 1). The ubiquitin-like domain PF14560 showed the strongest enrichment in induced genes (induced/total ratio = 0.83, *p*-val = 4.41e−05), and the protein kinase domain PF00069 had the most significant enrichment in induced genes (ratio 0.25, *p*-val 5.02e−10). Gene Ontology terms associated with the 74 overrepresented protein domains included defense mechanism and immune response in 64.5% of cases. For instance, 27 out of 72 genes harboring a WRKY domain (PF03106) were induced upon infection by *S. sclerotiorum* (ratio = 0.375, *p*-val 1.43 e−04). Other protein domains enriched in induced genes included ubiquitin-like domains (PF10302, PF14560, PF11976, and PF00240), transport-related domains (PF01105, PF08449, and PF03105), calcium binding (PF14658), and heat shock response (PF00011) (Supplementary Table S1).

One domain enriched in induced genes had no known molecular function and was identified as Domain of Unknown Function DUF4228 (ratio 0.5, *p*-val 3.86e−04). We identified 28 genes with a DUF4228 domain in the genome of *A. thaliana* (hmmscan e-value < 1E-10, Supplementary Table S2), 19 of them being differentially expressed upon infection by *S. sclerotiorum* (14 induced and 5 down-regulated). The DUF4228 gene AT5G37840 was induced over 1,000 times (LFC 10.22, *p*-val 1.45e−47), and 13 genes harboring a DUF4228 domain were in the top 10% most induced genes after infection by *S. sclerotiorum* in *A. thaliana* (LFC > 4.04, Supplementary Table S2). Because of their dramatic induction pattern and

although uncharacterized to date, some DUF4228 genes could function in plant defense responses.

## Taxonomic Distribution of the DUF4228/PADRE Domain

To document the taxonomic distribution of the DUF4228 domain across the tree of life, we performed a HMM search against the Refprot database of UniProtKB with an alignment of *A. thaliana* DUF4228 proteins as input (Supplementary Datasheet S1). We retrieved 3647 hits distributed in 98 species. As recently reported (Yang et al., 2020), DUF4228 appeared restricted to plants, including mosses, liverworts, and monocot and dicot species. The average size of DUF4228 domains detected in these proteins was 149.7 ± 40.3 amino acids, for proteins of 159.7 ± 37.6 amino acids long (Figure 2A). In good agreement, only 8.4% of proteins harboring a DUF4228 domain were multi-domain proteins. To identify the complete repertoire of DUF4228 in plant proteomes, we performed a HMM search against the Phytozome 12.1 database. Out of the 64 plant proteomes available at the time of our analysis, only the seven Chlorophyte proteomes did not show a single DUF4228 domain, indicating that the emergence of the DUF4228 domain occurred at least 450 million years ago (Figure 2B). Next, we used Timetree to relate the number of DUF4228-containing proteins with time of speciation in 49 plant species. In embryophytes, the number of DUF4228-containing proteins ranged from three (Selaginella moellendorffii) to 81 (Glycine max) (Figure 2B). A majority of embryophytes (28/45) had between 20 and 40 DUF4228-containing proteins, and there was no striking expansion of DUF4228 in a specific plant lineage. Recent whole genome duplication events were often associated with expanded DUF4228 repertoires, such as in Brassica rapa, Malus domestica, G. max, Zea mays, and Musa acuminata. Overall, the size of the DUF4228 family was well correlated (R² = 0.5066) with the total number of genes per genome across embryophyte species (Figure 2C).

Through homology and keyword searches, we found experimental insights into function for DUF4228 proteins. The *A. thaliana* AT4G37240 protein was identified as interacting with calmodulin proteins CAM4, 6, 7, 8, and 9 (Popescu et al., 2007). In addition, *A. thaliana* AT1G66480 was identified as interacting with Arabidopsis Response Regulator 14 (ARR14) in a yeast two-hybrid screen (Dortay et al., 2008). However, these protein-protein interactions have not been validated by independent approaches. In *Nicotiana tabacum*, the homolog of *A. thaliana* AT1G76600 was found responsive to tobacco mosaic virus and wounding and the corresponding protein designated as Pathogenesis-related protein of 23kDa (NtPRp23) (Akiyama et al., 2005). Its ortholog in *N. sylvestris* (LOC104235934) conferred tolerance to cadmium when expressed in yeast (Zhang et al., 2016). Recent work by Yang et al. (2020) revealed that several DUF4228 genes are responsive to drought, cold, or salt abiotic stress. Our analyses reported in this study indicated that several DUF4228 genes are responsive to infection by the fungal pathogen *S. sclerotiorum* and that *A. thaliana* DUF4228 proteins harbor intrinsically disordered regions. Based on this partial functional information and to facilitate further

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Protein domains enriched among *A. thaliana* genes upregulated upon *S. sclerotiorum* inoculation. Each bubble shows one of 54 PFAM domains significantly enriched (proportion test *p*-value < 0.01) in induced genes (LFC > 1.5, *p*-value < 0.01). Bubbles are sized according to the total number of genes containing the domain in the *A. thaliana* genome. Enrichment is shown as the *p*-value of a proportion Z-test for enrichment (X-axis), the ratio between the number of induced/total genes (Y-axis) and a composite enrichment score (color scale, see the section "Materials and Methods"). The DUF4228 domain is labeled in bold red. Associated raw data corresponds to *Arabidopsis thaliana* samples from GEO accession number GSE106811.</figcaption></figure>
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Taxonomic distribution of DUF4228 domains. (A) Relationship between total protein length and DUF4228 domain length in 3647 proteins from the Refprot database. (B) Timetree and DUF4228 family size for plant species from the Phytozome 12.1 database. The phylogenetic tree and divergence time estimates were obtained using Kumar et al. (2017). Polyploidization events described in the literature (Wang et al., 2019; Xu et al., 2019) are shown as stars with 4 spikes (tetraploidy), 6 spikes (hexaploidy), or 10 spikes (decaploidy). Greek letters correspond to common names attributed to polyploidy events in the literature, starting with events identified in the *Arabidopsis* genome. (C) Relationship between DUF4228 family size and total genome size (number of protein-coding loci) in 45 embryophyte species. Bubbles are sized according to genome size in Mbp and colored according to ploidy level.

reference, we propose to refer to this family as the pathogen and abiotic stress response, cadmium tolerance, disordered region-containing (PADRE) family.

## Sequence Diversification of the PADRE Domain

To analyze patterns of sequence diversification among PADRE proteins, we selected 13 plant genomes representative of the major Embryophyta lineages and constructed a phylogenetic tree of PADRE proteins from these species (see the section

"Materials and Methods"). For this, we generated a multiple protein alignment including 344 sequences and 116 informative sites located within the PADRE domain (Supplementary Datasheet S2, S3). We used maximum likelihood methods to represent phylogenetic relationships between these 354 PADRE domains as a tree (Figure 3A and Supplementary Datasheet S4). PADRE proteins classified into 10 monophyletic groups (a to j) supported by posterior probabilities ≥ 0.90 and encompassing 10 (clade e) to 55 (clade g) proteins. PADRE sequences diversified strongly since the divergence between Lycophytes and Angiosperms: groups a and i were restricted to Bryophytes
<figure><img src="image_5.png" /></figure>

**FIGURE 3** | Phylogenetic relationships between DUF4228 proteins in the complete proteome of 13 Embryophyta species. (A) Tree obtained by a maximum likelihood analysis, with the number of substitutions per site used as branch length, and branch support determined by an approximate likelihood ratio test (black labels if ≥0.90, blue otherwise). Terminal nodes are color-coded according to plant species (key shown in the upper panel). *A. thaliana* identifiers are labeled in red on the tree. Phylogenetic groups are labeled a to j on the outer circle. The upper panel shows the number of genes per species and per phylogenetic group as bubbles of increasing size. (B) Species tree showing rates of PADRE domain gain (green) and loss (red) in the evolution of Embryophyta as calculated with BadiRate (Librado et al., 2012). Neutral branches are shown in gray.
and Lycophytes, while groups b, c, e, f, g, h, and j were restricted to Angiosperms. Group d was represented in all species analyzed except *Sphagnum fallax*. Groups b, c, e, f, g, h, and j were represented in all Angiosperm species analyzed, with the exception of groups c and e that were absent from *Arabidopsis thaliana* and *Vitis vinifera*. This suggests that seven PADRE groups existed in the common Angiosperm ancestor and that groups c and e were lost in *A. thaliana* and *V. vinifera*. The number of PADRE groups expanded more rapidly in Angiosperms (reaching 6 and 8 distinct clades per species) than in Bryophytes and Lycophytes (2 or 3 clades per species), indicative of strong diversification of PADRE genes early in the evolution of Angiosperms. To estimate rates of domain birth and death in the PADRE family, we analyzed the species distribution of PADRE phylogenetic group with BadiRate (Librado et al., 2012) (Figure 3B). This revealed two major domain gain events during the emergence of Angiosperms and of core Eudicots and several lineage-specific gain events. Loss events mostly corresponded to the emergence of Tracheophytes and to terminal branches in the Fabids and Malvids clades.

## PADRE Is a Bipartite Domain Including Disordered and Phosphorylated C-Termini

PADRE proteins do not display clear homology to functionally characterized proteins. In our alignment of PADRE proteins from 13 Tracheophyta species, sequence conservation appeared limited to four short motifs of 10 amino acids or less (Figure 4A). These conserved motifs correspond to motifs 1, 3, and 6 identified by Yang et al. (2020). As noted by Yang et al. (2020), additional short sequence motifs were restricted to specific PADRE groups. To get insights into PADRE protein sequence signatures and their potential functional implications, we scanned *A. thaliana* PADRE proteins with the ELM, PhosPhat, PrDOS, and Grantham Polarity calculation tools. First, we used the eukaryotic linear motif (ELM) resource to identify motifs similar to known functional sites in proteins (Gouw et al., 2018) (Figure 4B). Among motifs identified robustly in multiple AtPADRE proteins was an N-myristoylation motif, corresponding to the well-conserved GNXXX motif found at the very N-terminus of PADRE proteins. *In vitro* myristoylation provided experimental for N-myristoylation of AT4G37240 (group G) and AT1G10530 (group D) (Boisson et al., 2003, unpublished result available<fnref n="1" />). Furthermore, AT1G21010 (group G) and ATGG17350 (group J) were identified in plasma-membrane fractions as expected if N-myristoylated (Majeran et al., 2018). The conserved LXXG motif of PADRE proteins overlapped with a WH2 motif for interaction with actin (LIG_Actin_WH2). The conserved YFLLP motif overlapped with a Tyrosine-based signal for interaction with the adaptor protein complex (TRG_ENDOCYTIC_2), a LIR motif for binding to the autophagy protein Atg8 (LIG_LIR_Gen_1), and a protein phosphatase interacting motif (DOC_PP1_RVXF_1). Basic nuclear localization signals were detected at the C-terminus

of several PADRE proteins. The C-terminal WRPXLXXIXE motif overlapped with an APCC-binding Destruction motif required for targeting to ubiquitin-mediated proteasome-dependent degradation (DEG_APCC_DBOX_1). ELM also detected numerous putative phosphorylation sites at the C-terminus of PADRE proteins. We took advantage of the PhosPhAt 4.0 database to search for experimentally determined phospho-peptides in PADRE proteins (Durek et al., 2009). We retrieved phospho-peptides from seven AtPADRE proteins from group B (AT1G06980), D (AT1G60010), F (AT1G64700, AT1G66480, AT2G01340, AT5G37840), and G (AT1G76600) (Figure 4C). A large majority of the phosphorylated residues resided in the C-terminal half of the PADRE domain.

We used the PrDOS server (Ishida and Kinoshita, 2007) to predict natively disordered regions in the 28 AtPADRE proteins (Figure 4D). All AtPADRE proteins showed a relatively consistent pattern of disorder probability, indicating a short (~5 amino acids) N-terminal disordered region followed by an ordered region of ~60 amino acids and a C-terminal half with high probability of intrinsic disorder. To test whether the structural state of PADRE regions was associated with contrasted amino acid usage at the N and C terminus, we calculated the Grantham Polarity index (Grantham, 1974) along 28 AtPADRE proteins (Figure 4E). In average, the PADRE C-terminal region harbors more polar residues (average index 9.02) than the N-terminal region (average index 8.39).

## Responsiveness to S. sclerotiorum Varies Across PADRE Phylogenetic Groups

The clear delineation of phylogenetic groups in the PADRE family and recent investigations of *AtPADRE* gene expression upon abiotic stress (Yang et al., 2020) suggested that *PADRE* genes could have acquired several distinct functions over evolution. Here, we set to investigate whether responsiveness to the fungal pathogen *S. sclerotiorum* contrasts across PADRE phylogenetic groups and whether responsiveness to fungal infection is consistent across plant species. To this end, we analyzed the expression of the *PADRE* gene repertoire of *A. thaliana*, *Solanum lycopersicum*, and *Ricinus communis* by RNA-sequencing in leaves of healthy plants and upon inoculation by *S. sclerotiorum*. We detected the expression of 74 *PADRE* genes, including 28 *AtPADRE*, 26 *SlPADRE*, and 23 *RcPADRE* (Figure 5 and Supplementary Table S3). Among them, 31 were significantly induced, 11 were significantly down-regulated, and 32 were not differentially expressed. *PADRE* genes highly expressed in healthy leaves were frequent in groups d and h and a subgroup of group f. Group h was very homogenous with all four genes significantly induced upon *S. sclerotiorum* infection, and groups j and f included mostly induced genes, with group d that included a majority of down-regulated genes. To determine whether *PADRE* gene induction differs significantly between phylogenetic groups or between species, we performed ANOVA on expression LFC for the 74 *PADRE* genes. In a one-way ANOVA, the phylogenetic group effect was found highly significant (*p*-value 0.0052) while the species was not significant (*p*-value 0.68). In a two-way

<footnote n="1">¹https://www.i2bc.paris-saclay.fr/maturation/Myristoylome.html</footnote>
<figure><img src="image_6.png" /></figure>

ANOVA, the phylogenetic group effect was significant (*p*-value 0.018); the species effect and group × species interaction effect was not significant (*p*-value 0.72 and 0.99, respectively). We conclude that *PADRE* gene expression upon *S. sclerotiorum* inoculation differs between phylogenetic groups in a consistent manner across plant species.

Through a synteny analysis, we identified five pairs (AT1G06980/AT2G30230, AT1G10530/AT1G60010, AT1G210

10/AT1G76600, AT3G03280/AT5G17350, and AT3G10120/AT5G03890) and two quartets AT1G71015/AT2G01340/AT1G66480/AT5G37840, AT2G23690/AT4G37240/AT5G66580/AT3G50800) of *PADRE* genes associated as paralogs in the *A. thaliana* genome. In all instances, groups of paralogs belonged to the same phylogenetic group. The divergence in expression of the *PADRE* genes was limited within most of the paralog groups (Supplementary Figure S1). Only the pair of paralogs
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** | Genome-wide expression of *PADRE* genes in *Arabidopsis thaliana*, *Ricinus communis*, and *Solanum lycopersicum* in response to *S. sclerotiorum* determined by RNA sequencing. Genes are ordered according to phylogenetic relationships shown in Figure 3, with phylogenetic groups delimited by dark and light gray boxes. Bar plot shows LFC of gene expression upon *S. sclerotiorum* inoculation, with error bars showing standard error for LFC over three independent biological experiments. Bars are filled when *p*-value for differential expression determined by DESeq2 analysis is ≤0.01, empty otherwise. Bubbles are sized and colored according to normalized gene expression in healthy plants. Exp., gene expression; LFC, log₂ fold change of gene expression. The associated raw data is available from GEO accessions GSE106811 and GSE138039.</figcaption></figure>

AT3G10120-AT5G03890 (group f) showed significant divergence in expression (LFC-5.23 and 3.53, respectively), possibly due to their low basal level of expression. Altogether, our gene expression analysis suggests that responsiveness to pathogens was acquired early in the evolution of Angiosperms by specific groups or subgroups of *PADRE* genes.

**A. thaliana PADRE Genes Respond to Multiple Stress and Associate With Plant Defense Ontologies**

The finding that some *AtPADRE* genes are responsive to abiotic stresses (Yang et al., 2020) prompted us to investigate their expression under a range of biotic stresses. For this, we analyzed RNA sequencing data available in the Gene Expression Omnibus database (Figure 6A), collecting expression data for *A. thaliana* inoculated by the fungal pathogens *S. sclerotiorum* (Badet et al., 2017), *Botrytis cinerea* (Liu et al., 2015), *Alternaria brassicicola* (Rausch, 2016), and *Verticillium dahliae* (Scholz et al., 2018), the bacterial pathogen *Pseudomonas syringae* pv. *tomato* (*Pst*) DC3000 (Mine et al., 2018) and DC3000 expressing the effector *AvrRps4* (Bhandari et al., 2019), the Cabbage Leaf Curl Virus (CaLCuV) (Zorzatto et al., 2015), and the nematode *Heterodera schachtii* (Shanks et al., 2016). To serve as a reference, we also analyzed RNA sequencing data for plants inoculated with the endophytic fungus *Colletotrichum*

*tofieldiae* (Hacquard et al., 2016) and submitted to heat stress (Albihlal et al., 2018), cold stress (Zuther et al., 2019), and UV-B treatment (Tavridou et al., 2020). The analysis of PADRE gene differential expression revealed a cluster of six *PADRE* genes induced by multiple pathogens: *AT1G28190*, *AT5G12340*, *AT1G76600*, *AT1G21010*, *AT5G37840*, and *AT2G01340* are significantly induced in response to *S. sclerotiorum*, *B. cinerea*, *Pst* DC3000 *AvrRPS4*, and *V. dahliae*. Out of these six genes, five are also induced upon infection by *A. brassicicola*, three under heat stress. Three of them are down-regulated in root response to the non-pathogenic fungus *C. tofieldiae*. By contrast, a cluster of five PADRE genes (*AT4G37240*, *AT1G60010*, *AT1G06980*, *AT2G23690*, and *AT5G66580*) was down-regulated in response to pathogens and heat stress. The response of *PADRE* genes to heat stress shared more similarities with their response to pathogens than to other abiotic stimuli. *PADRE* genes were not responsive to all signals. Indeed, only *AT2G01340* was differentially expressed in response to the nematode *H. schachtii*, *AT3G50800* and *AT5G62900* upon infection by the virus *CaLCuV*, *AT1G76600*, *AT1G21010*, *AT2G01340*, and *AT3G10120* in response to *C. tofieldiae*, and *AT4G37240*, *AT3G61920*, *AT2G23690*, and *AT1G76600* to UV-B. To test for the relationship between phylogenetic clades and the response of *PADRE* genes to diverse stresses, we performed a two-way ANOVA. We found a significant effect of the phylogenetic group (*p*-value 0.042) and type of stress (*p*-value 4.15 10⁻⁵) on *PADRE*
<figure><img src="image_8.png" /><figcaption>**FIGURE 6** | Response to multiple pathogens and co-expression network for *AtPADRE* genes. (A) Expression profiles of *A. thaliana* *PADRE* genes under multiple biotic and abiotic stresses deduced from published RNA sequencing data. The expression levels of genes were normalized using min-max feature scaling to fit within the [-1; 1] range for all experiments. Non-significant LFCs are displayed as 0. The phylogenetic group of *AtPADRE* genes is given between square brackets. The associated raw data is available from GEO accessions GSE132169, GSE70094, GSE72548, GSE56922, GSE88798, GSE112225, GSE85653, GSE83478, GSE104590, GSE116269, GSE66290, and GSE106811. *Pst*, *Pseudomonas syringae* pv. *tomato*. (B) Co-expression network for *AtPADRE* genes deduced from an experiment of 14,668 microarrays. Nodes are color-coded according to subcellular localization predicted by WOLF-Psort, shown as hexagons for transcription factors and as circles otherwise. *AtPADRE* genes are outlined and labeled in red. Edge widths are scaled according to a mutual rank score index for co-expression. Gray-shaded areas show a subnetwork identified by network modularity analysis, with associated specific gene ontologies labeled in bold italics. (C) The same co-expression network as in (B) with nodes color-coded according to LFC of gene expression upon infection by *S. sclerotiorum* determined by RNA sequencing (GEO accession GSE106811). LFC, log₂ fold change of gene expression.</figcaption></figure>
gene expression. Using a Tukey HSD test, we found that the phylogenetic group effect is due to contrasted expression patterns of genes from groups h and d (p-value 0.022). A Tukey HSD test on the stress variable indicated that the stress effect is due to *S. sclerotiorum* infection triggering *PADRE* gene expression significantly different from every other stresses (p-value < 0.05), except for *Pst* DC3000 *AvrRPS4* infection and heat stress (p-values 0.361 and 0.097, respectively). Therefore, we detected an association (i) between two PADRE phylogenetic groups and responsiveness to multiple stresses and (ii) between the expression of *PADRE* genes and infection by *S. sclerotiorum*. In addition, visual inspection of **Figure 6A** suggested that *S. sclerotiorum*, *B. cinerea*, *Pst* DC3000 *AvrRPS4*, *V. dahliae*, *A. brassicicola*, and heat stress induced similar transcriptional responses in *PADRE* genes.

To get insights into genes functioning in the same processes as PADRE genes, we retrieved a co-expression network for *AtPADRE* genes from the ATTED-II database covering 14,668 microarray samples (Obayashi et al., 2018) (**Supplementary Datasheet S5** and **Figure 6B**). The network was composed of 225 nodes and 523 undirected edges, including 19 *AtPADRE* genes. We mapped LFC of gene expression upon *S. sclerotiorum* inoculation obtained from our RNA sequencing analysis onto this network, revealing one major sector including predominantly highly induced genes and another sector including mostly down-regulated genes (**Figure 6C**). To emphasize biological processes involving *AtPADRE* genes, we performed a modularity analysis based on the network topology (Blondel et al., 2008) to compute subnetworks and test next if every subnetwork corresponded to gene ontology. The modularity analysis identified 12 subnetworks, four of which were significantly associated with a specific biological function (**Figures 6B,C**). A subnetwork, strongly overexpressed during *S. sclerotiorum* infection (**Figure 6C**), was involved in the perception of the fungus cell wall and in the response to chitin (FDR = 5.3E-11). Response to hypoxia (FDR = 9.45E-3, FDR = 9.29E-27) was overexpressed during infection by *S. sclerotiorum* whereas genes involved in the cell-cell junction assembly were downregulated (**Figure 6C**). The subnetwork grouping genes associated with the biosynthesis of anthocyanin (FDR = 2.7E-4), secondary metabolites with antifungal activity (Kumar Sudheeran et al., 2019), appeared overexpressed during infection. To test further the role played by *AtPADRE* genes in the topology of the network, we computed the local centrality (or degree) of every gene (**Supplementary Figure S2**). Despite the high centrality of the *At5g17350* gene, centralities of *AtPADRE* genes did not differ significantly from other genes of the network (mean degree *AtPadre* 4.63, others 4.65, Wilcoxon's test *p*-value = 0.53).

DISCUSSION

Plant genomes harbor a remarkably large number of gene families that are not found in other life kingdoms, several of which function in cell signaling (Yamasaki et al., 2013) and defense (Tenhaken et al., 2005; Raffaele et al., 2007; Weidenbach et al., 2016). Through contextual searches, we identified experimental

evidence that DUF4228 genes are involved in pathogen and abiotic stress response, cadmium tolerance, disordered region-containing family (Zhang et al., 2016; Yang et al., 2020), and we propose to refer to this family as the PADRE family to reflect these functional and architectural information. Naming domains based on the first functional clues is unlikely to reflect all or the most prominent function of gene families but can foster further research on the function of these genes (Doerks et al., 2000; Habermann, 2004; Tenhaken et al., 2005). Some PADRE genes are responsive to environmental stimuli such as wounding and viruses (Akiyama et al., 2005), drought, cold, and salt (Yang et al., 2020), pointing toward yet uncharacterized molecular functions.

Insights Into the Evolutionary History of
the PADRE Domain

We identified 344 high-quality PADRE protein sequences across 13 plant genomes and used this information in a phylogenetic analysis to explore the dynamics of the PADRE domain evolution. Analysis of the extent diversity of PADRE proteins suggests that they originated before the divergence between Bryophyta and Tracheophyta, like an estimated ~50% of plant-specific domains (Kersting et al., 2012). We classified PADRE proteins into 10 phylogenetic groups, corresponding approximately to subdivisions of the three groups proposed by Yang et al. (2020). In our analysis, the phylogenetic signal was too weak to infer a common ancestor to several groups and combine them with confidence. The BadiRate analysis highlighted a strong radiation of the PADRE domain at the base of the Angiosperms, around ~350 to 175 million years ago. It should be noted that our dataset does not include sequences from the Pinophyta and Pteridophyta lineages, so that the burst of PADRE diversification may date back to the divergence of these groups or to the Angiosperm most recent common ancestor. The recent duplication of PADRE genes from groups b, d, f, g, and j in *A. thaliana* is well supported by the phylogeny and synteny analysis and consistent with Yang et al. (2020). Recent duplications in these groups are also likely in *M. truncatula*, *S. lycopersicum*, and *A. coerulea*. However, PADRE domain births remained limited or null within the core Eudicot clade, where domain loss seemed predominant. This could indicate selection toward some degree of functional specialization in the PADRE family, favoring the expansion of a few clades to the detriment of the overall domain diversity. Our pangenomic expression analysis supported somewhat consistent patterns of PADRE gene expression upon *S. sclerotiorum* inoculation within phylogenetic groups and across species. This suggests that responsiveness to fungal infection was acquired by PADRE groups f, g, h, and j early in the evolution of core Eudicots. Nevertheless, there was striking contrast in expression within groups f and g, which may indicate some degree of neo- or subfunctionalization.

A Probable Bipartite Architecture With
Structured and Disordered Regions

Sequence analysis pointed toward a bipartite architecture
for the PADRE domain, with a combination of structured
and intrinsically disordered regions. Intrinsically disordered
<table><caption>**TABLE 1** | List of experiments used to analyze the response of *A. thaliana* PADRE to multiple biotic and abiotic stresses, with the corresponding GEO accession numbers.</caption><thead><tr><th>Plants and treatment</th><th>Plant age</th><th>Time post treatment</th><th>GEO accession</th></tr></thead><tbody><tr><td>A. thaliana and S. lycopersicum inoculated by S. sclerotiorum</td><td>4 weeks</td><td>48 h</td><td>GSE106811</td></tr><tr><td>R. communis inoculated by S. sclerotiorum</td><td>4 weeks</td><td>50 h</td><td>GSE138039</td></tr><tr><td>A. thaliana inoculated by Botrytis cinerea</td><td>4 weeks</td><td>14 h</td><td>GSE66290</td></tr><tr><td>A. thaliana inoculated by Alternaria brassicicola</td><td>NA</td><td>NA</td><td>GSE83478</td></tr><tr><td>A. thaliana inoculated by Verticillium dahlia</td><td>14 days</td><td>24 h</td><td>GSE104590</td></tr><tr><td>A. thaliana inoculated by Colletotrichum tofieldiae</td><td>10 days</td><td>10 days</td><td>GSE70094</td></tr><tr><td>A. thaliana inoculated by Pseudomonas syringae pv. tomato DC3000</td><td>30 days</td><td>24 h</td><td>GSE88798</td></tr><tr><td>A. thaliana inoculated by Pseudomonas syringae pv. tomato DC3000 AvrRPS4</td><td>4 weeks</td><td>24 h</td><td>GSE116269</td></tr><tr><td>A. thaliana inoculated by Heterodera schachtii</td><td>10 days</td><td>10 days</td><td>GSE72548</td></tr><tr><td>A. thaliana inoculated by CaLCuV</td><td>7 leaf stage</td><td>20 days</td><td>GSE56922</td></tr><tr><td>Heat-treated (37°C) A. thaliana Col-0</td><td>5 weeks</td><td>30 min</td><td>GSE85653</td></tr><tr><td>Cold-treated (4°C) A. thaliana Col-0</td><td>28 days</td><td>3 days</td><td>GSE112225</td></tr><tr><td>UV-B-treated A. thaliana Col-0</td><td>4 days</td><td>6 h</td><td>GSE132169</td></tr></tbody></table>

NA, not available.

regions (IDRs) are flexible protein regions lacking a stable 3D fold in solution, which may transition to an ordered state upon binding to natural ligands (Uversky, 2013). Proteins with IDRs are abundant in eukaryotic genomes and are depleted in hydrophobic residues and enriched in polar and charged residues. We found higher amino acid polarity at the C-terminus of PADRE proteins, in agreement with high disorder probability in this region. The peculiar composition and folding properties of IDRs confer specific functional properties (Sun et al., 2013; Uversky, 2013). First, IDRs are generally able to establish protein-protein interactions with multiple partners and are commonly found in hub proteins in eukaryotic networks. One paradigmatic example in plant immunity is RPM1-interacting protein 4 (RIN4) which interacts with multiple plant resistance proteins and bacterial effectors (Sun et al., 2014). In line with this property, PADRE proteins were shown experimentally to interact with calmodulins (Popescu et al., 2007) and response regulators (Dortay et al., 2008). Our co-expression network also suggests a high degree of connectivity for PADRE genes. Screening for protein-protein interactions involving PADRE proteins should prove an insightful avenue for future research. Second, IDRs are highly accessible regions and can therefore undergo complex regulations by post-translational modifications. For instance, Remorins are plant-specific proteins with a role in plant immunity (Raffaele et al., 2009; Bozkurt et al., 2014) containing structured and disordered regions, with their IDRs harboring multiple phosphorylation sites (Marín and Ott, 2012; Marín et al., 2012; Perraki et al., 2018). RIN4 also undergoes multiple post-translational modifications and regulation by proteolysis (Toruño et al., 2019). Similarly, we identified multiple phosphorylated residues in the C-terminal region of PADRE proteins, as well as degradation signals. Third, the ability to undergo a disorder-to-order transition can confer transient functionality to IDRs, such as membrane binding in Remorins (Perraki et al., 2012), cytotoxic activity of Bordetella CyaA toxin (O'Brien et al., 2018), and protein complex formation by cAMP response element-binding (CREB)

protein (Arai et al., 2015). We could then speculate that every PADRE protein could adopt several functions according to their cellular environment.

## Toward a Functional Understanding of PADRE Family

We report the significant induction of 31 PADRE genes upon inoculation by S. sclerotiorum, including 14 AtPADRE, 7 RcPADRE, and 10 SIPADRE genes. Radiation of the PADRE family into 10 phylogenetic groups could provide the basis for some degree of functional diversification. In line with hypothesis, Yang et al. (2020) identified 3 AtPADRE genes induced upon osmotic stress, 4 upon salt stress, and 5 upon cold stress. Our work revealed an intrinsically disordered region in PADRE proteins, suggesting that PADRE gene function could be context-dependent. This could explain why Yang et al. (2020) found several AtPADRE genes mis-regulated by salt while none were significantly responding to NaCl in the RNA sequencing dataset we analyzed (Suzuki et al., 2016). The identification of multiple subcellular localization signals in PADRE proteins (N-myristoylation, NLS, endocytic vesicles) prevents predictions regarding the site of PADRE action. The use of fluorescent protein reporter fusions in a structure-function analysis will be required to this end. This approach shall be challenging given the presence of targeting signals at both ends of the PADRE domain. We found six AtPADRE genes induced upon inoculation by several fungal pathogens with a necrotrophic lifestyle (S. sclerotiorum, B. cinerea, and A. brassicicola), a bacterial pathogen (P. syringae pv. tomato), and a hemibiotrophic root-infecting fungus (V. dahliae), indicating that pathogens are very potent inducers of AtPADRE genes. The PADRE co-expression network included several important players in plant immunity such as the syntaxin SYP122 (Zhang et al., 2007), the C2-domain protein BAP1 (Yang et al., 2006), the patatin-like protein 2 PLP2 (La Camera et al., 2009), members of the RPM1-interacting protein 4 RIN4 (At3g48450), wall-associated kinase-like WAKL10 (At1g79680), and the nematode resistance protein-like HSPRO2 (At2g40000).
These findings are consistent with a role for members of the PADRE family in disease resistance.

# MATERIALS AND METHODS

## Pfam Domain Annotation and Enrichment Analyses

Pfam domains were annotated using hmmscan 3.1b1 with e-value threshold 1E-10 against the Pfam-A 32.0 database. Enrichment of Pfam domains among genes induced after S. sclerotiorum infection was analyzed using a two-proportion Z-test in R. Arabidopsis thaliana gene expression from GEO accession GSE106811 (Badet et al., 2017) was used in this analysis. Briefly, total RNA was extracted from the edge of developed necrotic lesions of leaves from 4-week-old plants inoculated by S. sclerotiorum strain 1980, as described in Peyraud et al. (2019). Samples were collected in triplicates from three plants in independent inoculation experiments. RNA sequencing was performed on an Illumina HiSeq 2500 instrument as described in Badet et al. (2017). A composite enrichment score taking into account the significance of the Z-test and the enrichment ratio was calculated with the formula $R_Z(i) * R_r(i)$, where $R_Z(i)$ is the normalized rank of domain i for the Z-test p-value and $R_r(i)$ is the normalized rank of domain i for the enrichment ratio.

## RNA Sequencing Data Analysis

Raw data for RNA sequencing experiments used in this work is available in the NCBI Gene Expression Omnibus (GEO) database with accession numbers provided in Table 1. All raw datasets were processed separately with DESeq2 to calculate normalized read counts (Basemean) and log2 fold change (LFC) of expression and identify genes differentially expressed between control and treated samples. Genes were considered differentially expressed for LFC ≥ 1.5 and adjusted p-value ≤ 0.01. In the multiple stress analysis, raw data were used for statistical analysis and LFC values were normalized for the heatmap, as follows: (±)log2(1 + |LFC|).

## Taxonomic Distribution of PADRE Proteins

We used MAFFT Version 7.407 (Katoh et al., 2002) to align the *Arabidopsis thaliana* DUF4228 protein sequences using default parameters. After manual curation, 24 *A. thaliana* DUF4228 proteins expressed in our RNA sequencing data (GSE106811) were kept for further analysis. This alignment (**Supplementary Datasheet S1**) was used in a phmm器 search on the HMMER webserver<fnref n="2" /> (Potter et al., 2018) against the UniProt References Proteomes in UniProtKB (The Uniprot Consortium, 2019). The search was carried out using parameters -E 1e-10 -domE 1 -incE 1e-10 -incdomE 0.03 -seqdb uniprotrefprot identifying 3467 significant sequence hits. The ‘target length’ and length of the target alignment from the output of the HMM search were used to compare total protein length and DUF4228 domain length (**Figure 2A**). Genome sizes and total number of genes per genomes were obtained from the Phytozome 12.1 database

(Goodstein et al., 2011). Ploidy levels were obtained from the Plant DNA C-values Database on the Kew Royal Botanic Gardens website (Plant DNA C-values Database | Royal Botanic Gardens, Kew) and from the original genome papers. The timetree was generated on http://www.timetree.org/ (Kumar et al., 2017) using species names as input. Polyploidization events described in the literature were collected from Wang et al. (2019) and Xu et al. (2019). DUF4228 proteins in complete plant genomes were identified using hmmsearch against a local instance of the Phytozome 12.1 proteome database, using the same parameters as previously.

## Phylogenetic Analysis of PADRE Proteins

We extracted DUF4228 proteins from *Marchantia polymorpha*, *Physcomitrella patens*, *Sphagnum fallax*, *Selaginella moellendorffii*, *Amborella trichopoda*, *Brachypodium distachyon*, *Setaria italica*, *Aquilegia coerulea*, *Solanum lycopersicum*, *Medicago truncatula*, *Arabidopsis thaliana*, *Ricinus communis*, *Vitis vinifera*, and *Theobroma cacao* from our hmmsearch against Phytozome 12.1. Prior to alignment, we removed sequence Pp3c24_13210 for having <40 amino acids and truncated the 650, 650, and 1650 N-terminal amino acids from Solyc05g013500, Thecc1EG010515, and Medtr8g069400, respectively. A first sequence alignment was performed in ClustalO (Madeira et al., 2019), and 10 sequences were removed for being too divergent, leaving 344 sequences (**Supplementary Datasheet S2**). These sequences were aligned with ClustalO, and the alignment was manually edited in Jalview to keep positions with no gap in at least 172/344 sequences, yielding a final alignment of 116 amino acids long (**Supplementary Datasheet S3**). Phylogenetic relationships were determined by a maximum likelihood approach using PhyML (Guindon et al., 2010) with aLRT branch support in phylogeny.fr (Dereeper et al., 2008), with no alignment and no alignment curation steps (**Supplementary Datasheet S4**), using the LG substitution model (Le and Gascuel, 2008) and a gamma distribution with four categories. The resulting tree had a log-likelihood of -48028.9 and gamma shape parameter 1.720. The tree was rooted on *M. polymorpha* Mapoly0024s0004 and rendered with FigTree<fnref n="3" /> v1.4.3. Phylogenetic groups were defined based on the most ancestral branch with support ≥ 0.9. The rates of birth and death of PADRE domains were calculated using BadiRate 1.35 (Librado et al., 2012) with parameters -bmodel FR -ep CML -family.

## Bioinformatics Analyses of PADRE Sequence Features

Conserved motifs were identified in the 28 *A. thaliana* PADRE proteins using the alignment provided in **Supplementary Datasheet S3** and rendered using WebLogo 3 (Crooks et al., 2004). ELMs were identified using the ELM webserver (Gouw et al., 2018) with *A. thaliana* as species and subcellular localization not specified. Phosphorylated peptides were identified with a 'Basic search' in 'Experiment data' in the PhosPhAt 4.0 database (Durek et al., 2009). Intrinsic disorder probability was calculated using the PrDOS webserver

<footnote n="2">²https://www.ebi.ac.uk/Tools/hmmer/</footnote>

<footnote n="3">³http://tree.bio.ed.ac.uk/software/figtree/</footnote>
(Ishida and Kinoshita, 2007) with a false-positive rate of 5%. Grantham residue polarity was determined using the ProtScale tool in ExPASy (Gasteiger et al., 2005) with a window size of 9.

## Reconstruction of PADRE Co-expression Network

The co-expression network was built using the NetworkDrawer tool in ATTED-II version 9.2 (Obayashi et al., 2018) using the Ath-m version C7.1 platform including 14,668 microarray samples, with the Coex option "Add many genes" and PPI option "Add a few genes." The resulting network was rendered in Cytoscape 3.6.1 (Shannon et al., 2003). Gene expression LFC upon S. sclerotiorum corresponds to the A. thaliana RNA sequencing data from Badet et al. (2017) (GSE106811), with LFC values provided as node attribute table in Cytoscape (Shannon et al., 2003). The modularity of the network was computed by the algorithm proposed by Blondel et al. (2008). Gene ontologies associated with subnetworks were determined using the GO enrichment analysis online tools<fnref n="4" />. Cutoff on FDR was set at 1E-2.

## DATA AVAILABILITY STATEMENT

RNA sequencing read datasets are available from the NCBI Gene Expression Omnibus (GEO) database with accession numbers GSE106811, GSE138039, GSE66290, GSE83478, GSE104590, GSE70094, GSE116269, GSE72548, GSE56922, and GSE72806. The datasets generated by the analyses presented in this study are included in the article **Supplementary Material**.

## AUTHOR CONTRIBUTIONS

MD, MK, and SR performed phylogenetic analyses. MD, LG, and SR performed gene expression analyses. AB and SR performed co-expression network analysis. SR conceived and designed the study. All authors contributed to writing the manuscript draft, reviewed the manuscript, and approved the final article.

## FUNDING

This work was supported by a Ph.D. grant from INRAE SPE Division to MD, and the French Laboratory of Excellence Project TULIP (Grants ANR-10-LABX-41 and ANR-11-IDEX-0002-02).

## REFERENCES

* Akiyama, K., Ueda, H., Yoda, H., and Sano, H. (2005). A novel basic pathogenesis-related protein from tobacco plants. *Plant Biotechnol.* 22, 157–159. doi: 10.5511/plantbiotechnology.22.157

* Albihlal, W. S., Obomighie, I., Blein, T., Persad, R., Chernukhin, I., Crespi, M., et al. (2018). *Arabidopsis* HEAT SHOCK TRANSCRIPTION FACTOR1b regulates multiple developmental genes under benign and stress conditions. *J. Exp. Bot.* 69, 2847–2862. doi: 10.1093/jxb/ery142

* Arai, M., Sugase, K., Dyson, H. J., and Wright, P. E. (2015). Conformational propensities of intrinsically disordered proteins influence the mechanism of

The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## ACKNOWLEDGMENTS

We thank the INRAE SPE division and Labex TULIP community for stimulating discussions and support.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00491/full#supplementary-material

**FIGURE S1** | Genomic distribution and synteny analysis of DUF4228 genes in A. thaliana. For each gene, the phylogenetic clade and the level of expression (log2 fold change) upon infection by S. sclerotiorum is provided. Genes significantly differentially expressed upon infection are indicated with a star (\*). Lines in the center of the graph show PADRE synteny blocks. The quartet AT2G23690/AT4G37240/AT5G66580/AT3G50800 is reported as two gene pairs in Yang et al. (2020), the quartet AT1G71015/AT2G01340/AT1G66480/AT5G37840 was not associated with gene duplication by Yang et al. (2020).

**FIGURE S2** | Modularity of genes included in the AtPADRE co-expression network. Genes with values higher than average are labeled on the graph. Boxplots show first and third quartiles (box), median (thick line) and the most dispersed values within 1.5 times the interquartile range (whiskers).

**TABLE S1** | PFAM domains enriched in genes induced upon S. sclerotiorum inoculation.

**TABLE S2** | Expression of 24 PADRE genes from A. thaliana under a range of biotic and abiotic stress, with mention of their inclusion in the list of genes upregulated upon S. sclerotiorum inoculation, down-regulated upon S. sclerotiorum inoculation and top 10% most induced genes after infection by S. sclerotiorum.

**TABLE S3** | Expression of 74 PADRE genes from A. thaliana, S. lycopersicum and R. communis upon S. sclerotiorum inoculation determined by RNA sequencing.

**DATASHEET S1** | Sequence alignment of 24 A. thaliana DUF4228 proteins produced by MAFFT and used for phmmer search, in fasta format.

**DATASHEET S2** | Full-length sequence of 344 PADRE proteins from 13 Embryophyte species in fasta format.

**DATASHEET S3** | Sequence alignment of 344 PADRE proteins from 13 Embryophyte species produced by ClustalO and used for phylogenetic analyses, in fasta format.

**DATASHEET S4** | Maximum likelihood phylogenetic tree of 344 PADRE proteins from 13 Embryophyte species, in Newick format.

**DATASHEET S5** | Co-expression gene network shown in **Figures 6B,C**, in .xgmm1 format.

* binding and folding. *Proc. Natl. Acad. Sci. U.S.A.* 112, 9614–9619. doi: 10.1073/pnas.1512799112

* Badet, T., Léger, O., Barascud, M., Voisin, D., Sadon, P., Vincent, R., et al. (2019). Expression polymorphism at the ARPC 4 locus links the actin cytoskeleton with quantitative disease resistance to Sclerotinia sclerotiorum in *Arabidopsis thaliana*. *New Phytol.* 222, 480–496. doi: 10.1111/nph.15580

* Badet, T., Voisin, D., Mbengue, M., Barascud, M., Sucher, J., Sadon, P., et al. (2017). Parallel evolution of the POQR prolyl oligo peptidase gene conferring plant quantitative disease resistance. *PLoS Genet.* 13:e1007143. doi: 10.1371/journal.pgen.1007143

<footnote n="4">⁴geneontology.org</footnote>
Bhandari, D. D., Lapin, D., Kracher, B., von Born, P., Bautor, J., Niefind, K., et al. (2019). An ED1 heterodimer signalling surface enforces timely reprogramming of immunity genes in *Arabidopsis*. *Nat. Commun.* 10:772. doi: 10.1038/s41467-019-08783-0

Blondel, V. D., Guillaume, J.-L., Lambiotte, R., and Lefebvre, E. (2008). Fast unfolding of communities in large networks. *J. Stat. Mech. Theory Exp.* 2008:10008. doi: 10.1088/1742-5468/2008/10/P10008

Boisson, B., Giglione, C., and Meinnel, T. (2003). Unexpected protein families including cell defense components feature in the N-myristoylome of a higher eukaryote. *J. Biol. Chem.* 278, 43418–43429. doi: 10.1074/jbc.M307321200

Boland, G. J., and Hall, R. (1994). Index of plant hosts of *Sclerotinia sclerotiorum*. *Can. J. Plant Pathol.* 16, 93–108.

Bozkurt, T. O., Richardson, A., Dagdas, Y. F., Mongrand, S., Kamoun, S., and Raffaele, S. (2014). The plant membrane-associated REMORIN1.3 accumulates in discrete perihaustorial domains and enhances susceptibility to phytophthora infestans. *Plant Physiol.* 165, 1005–1018. doi: 10.1104/pp.114.235804

Corwin, J. A., Copeland, D., Feusier, J., Subedy, A., Eshbaugh, R., Palmer, C., et al. (2016). The quantitative basis of the *Arabidopsis* innate immune system to endemic pathogens depends on pathogen genetics. *PLoS Genet.* 12:e1005789. doi: 10.1371/journal.pgen.1005789

Crooks, G. E., Hon, G., Chandonia, J. M., and Brenner, S. E. (2004). WebLogo: a sequence logo generator. *Genome Res.* 14, 1188–1190. doi: 10.1101/gr.849004

Derbyshire, M., Mbengue, M., Barascud, M., Navaud, O., and Raffaele, S. (2019). Small RNAs from the plant pathogenic fungus *Sclerotinia sclerotiorum* highlight host candidate genes associated with quantitative disease resistance. *Mol. Plant Pathol.* 20, 1279–1297. doi: 10.1111/mpp.12841

Dereeper, A., Guignon, V., Blanc, G., Audic, S., Buffet, S., Chevenet, F., et al. (2008). Phylogeny. *fr*: robust phylogenetic analysis for the non-specialist. *Nucleic Acids Res.* 36, W465–W469. doi: 10.1093/nar/gkn180

Dodds, P. N., and Rathjen, J. P. (2010). Plant immunity: towards an integrated view of plant-pathogen interactions. *Nat. Rev. Genet.* 11, 539–548. doi: 10.1038/nrg2812

Doerks, T., Strauss, M., Brendel, M., and Bork, P. (2000). GRAM, a novel domain in glucosyltransferases, myotubularins and other putative membrane-associated proteins. *Trends Biochem. Sci.* 25, 483–485. doi: 10.1016/s0968-0004(00)01 664-9

Dortay, H., Gruhn, N., Pfeifer, A., Schwerdtner, M., Schmülling, T., and Heyl, A. (2008). Toward an interaction map of the two-component signaling pathway of *Arabidopsis thaliana*. *J. Proteome Res.* 7, 3649–3660. doi: 10.1021/pr0703831

Durek, P., Schmidt, R., Heazlewood, J. L., Jones, A., MacLean, D., Nagel, A., et al. (2009). PhosPhAt: The *Arabidopsis thaliana* phosphorylation site database. An update. *Nucleic Acids Res.* 38, D828–D834. doi: 10.1093/nar/gkp810

El-Gebali, S., Mistry, J., Bateman, A., Eddy, S. R., Luciani, A., Potter, S. C., et al. (2019). The Pfam protein families database in 2019. *Nucleic Acids Res.* 47, D427–D432. doi: 10.1093/nar/gky995

Eulgem, T. (2005). Regulation of the *Arabidopsis* defense transcriptome. *Trends Plant Sci.* 10, 71–78. doi: 10.1016/j.tplants.2004.12.006

Fordyce, R. F., Soltis, N. E., Caseys, C., Gwinner, R., and Corwin, J. A. (2018). Digital imaging combined with genome-wide association mapping links loci to plant-pathogen interaction traits. *Plant Physiol.* 178, 1406–1422. doi: 10.1104/ pp.18.00851

Fukuoka, S., Saka, N., Koga, H., Ono, K., Shimizu, T., Ebana, K., et al. (2009). Loss of function of a proline-containing protein confers durable disease resistance in rice. *Science* 325, 998–1001. doi: 10.1126/science.1175550

Gasteiger, E., Hoogland, C., Gattiker, A., Duvaud, S., Wilkins, M. R., Appel, R. D., et al. (2005). "Protein identification and analysis tools on the ExPASy server," in *The Proteomics Protocols Handbook*, ed. J. M. Walker (Totowa, NJ: Humana Press), 571–607. doi: 10.1385/1-59259-584-7:531

Gatta, A. T., Wong, L. H., Sere, Y. Y., Calderón-Noreña, D. M., Cockcroft, S., Menon, A. K., et al. (2015). A new family of StART domain proteins at membrane contact sites has a role in ER-PM sterol transport. *eLife* 4, 1–46. doi: 10.7554/eLife.07253

Glazebrook, J. (2005). Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. *Annu. Rev. Phytopathol.* 43, 205–227. doi: 10.1146/ annurev.phyto.43.040204.135923

Goodacre, N. F., Gerloff, D. L., and Uetz, P. (2013). Protein domains of unknown function are essential in bacteria. *mBio* 5:e00744-13. doi: 10.1128/mBio. 00744-13

Goodstein, D. M., Shu, S., Howson, R., Neupane, R., Hayes, R. D., Fazo, J., et al. (2011). Phytozome: a comparative platform for green plant genomics. *Nucleic Acids Res.* 40, D1178–D1186. doi: 10.1093/nar/gkr944

Gouw, M., Michael, S., Sámano-Sánchez, H., Kumar, M., Zeke, A., Lang, B., et al. (2018). The eukaryotic linear motif resource - 2018 update. *Nucleic Acids Res.* 46, D428–D434. doi: 10.1093/nar/gkx1077

Grantham, R. (1974). Amino acid difference formula to help explain protein evolution. *Science* 185, 862–864. doi: 10.1126/science.185.415 4.862

Guindon, S., Dufayard, J.-F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. *Syst. Biol.* 59, 307–321. doi: 10.1093/sysbio/syq010

Habermann, B. (2004). The BAR-domain family of proteins: a case of bending and binding? *EMBO Rep.* 5, 250–255. doi: 10.1038/sj.embor.7400105

Hacquard, S., Kracher, B., Hiruma, K., Münch, P. C., Garrido-Oter, R., Thon, M. R., et al. (2016). Survival trade-offs in plant roots during colonization by closely related beneficial and pathogenic fungi. *Nat. Commun.* 7:11362. doi: 10.1038/ncomms11362

Hegedus, D. D., and Rimmer, S. R. (2005). Sclerotinia sclerotiorum: when “to be or not to be” a pathogen? *FEMS Microbiol. Lett.* 251, 177–184. doi: 10.1016/j. femsle.2005.07.040

Ishida, T., and Kinoshita, K. (2007). PrDOS: prediction of disordered protein regions from amino acid sequence. *Nucleic Acids Res.* 35, W460–W464. doi: 10.1093/nar/gkm363

Jaroszewski, L., Li, Z., Krishna, S. S., Bakolitsa, C., Wooley, J., Deacon, A. M., et al. (2009). Exploration of uncharted regions of the protein universe. *PLoS Biol.* 7:e1000205. doi: 10.1371/journal.pbio.1000205

Katoh, K., Misawa, K., Kuma, K., and Miyata, T. (2002). MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. *Nucleic Acids Res.* 30, 3059–3066. doi: 10.1093/nar/gkf436

Kersting, A. R., Bornberg-Bauer, E., Moore, A. D., and Grath, S. (2012). Dynamics and adaptive benefits of protein domain emergence and arrangements during plant genome evolution. *Genome Biol. Evol.* 4, 316–329. doi: 10.1093/gbe/ evs004

Khapff, M., Balaguè, C., Huard-Chauveau, C., and Roby, D. (2017). An essential role for the VASt domain of the *Arabidopsis* VAD1 protein in the regulation of defense and cell death in response to pathogens. *PLoS One* 12:e0179782. doi: 10.1371/journal.pone.0179782

Khapff, M., Cottret, L., Balagué, C., and Raffaele, S. (2014). Identification and phylogenetic analyses of VASt, an uncharacterized protein domain associated with lipid-binding domains in Eukaryotes. *BMC Bioinformatics* 15:222. doi: 10.1186/1471-2105-15-222

Krattinger, S. G., Lagudah, E. S., Spielmeyer, W., Singh, R. P., Huerta-Espino, J., McFadden, H., et al. (2009). A Putative ABC Transporter Confers Durable Resistance to Multiple Fungal Pathogens in Wheat. *Science* 323, 1360–1363. doi: 10.1126/science.1166453

Kumar, S., Stecher, G., Suleski, M., and Hedges, S. B. (2017). Timetree: a resource for timelines, timetrees, and divergence times. *Mol. Biol. Evol.* 34, 1812–1819. doi: 10.1093/molbev/msx116

Kumar Sudheeran, P., Ovadia, R., Galsarker, O., Maoz, I., Sela, N., Maurer, D., et al. (2019). Glycosylated flavonoids: fruit's concealed antifungal arsenal. *New Phytol.* 225, 1788–1798. doi: 10.1111/nph.16251

La Camera, S., Balagué, C., Göbel, C., Geoffroy, P., Legrand, M., Feussner, I., et al. (2009). The *Arabidopsis* patatin-like protein 2 (PLP2) plays an essential role in cell death execution and differentially affects biosynthesis of oxylipins and resistance to pathogens. *Mol. Plant Microbe Interact.* 22, 469–481. doi: 10.1094/MPMI-22-4-0469

Le, S. Q., and Gascuel, O. (2008). An improved general amino acid replacement matrix. *Mol. Biol. Evol.* 25, 1307–1320. doi: 10.1093/molbev/msn067

Leitch, I. J., Johnston, E., Pellicer, J., Hidalgo, O., and Bennett, M. D. (2019). Plant DNA C-values database (release 7.1, Apr 2019) Available online at: https://cvalues.science.kew.org/

Librado, P., Vieira, F. G., and Rozas, J. (2012). BadiRate: estimating family turnover rates by likelihood-based methods. *Bioinformatics* 28, 279–281. doi: 10.1093/ bioinformatics/btr623

Liu, S., Kracher, B., Ziegler, J., Birkenbihl, R. P., and Somssich, I. E. (2015). Negative regulation of ABA signaling by WRKY33 is critical for *Arabidopsis*
immunity towards Botrytis cinerea 2100. *eLife* 4:e07295. doi: 10.7554/eLife.07295 and functional perspectives. *Plant Physiol.* 145, 593–600. doi: 10.1104/pp.107.108639

Madeira, F., Park, Y. M., Lee, J., Buso, N., Gur, T., Madhusoodanan, N., et al. (2019). The EMBL-EBI search and sequence analysis tools APIs in 2019. *Nucleic Acids Res.* 47, W636–W641. doi: 10.1093/nar/gkz268

Majeran, W., Le Caer, J. P., Ponnala, L., Meinnel, T., and Giglione, C. (2018). Targeted profiling of *Arabidopsis thaliana* subproteomes illuminates co- and posttranslationally N-terminal myristoylated proteins. *Plant Cell* 30, 543–562. doi: 10.1105/tpc.17.00523

Marín, M., and Ott, T. (2012). Phosphorylation of intrinsically disordered regions in remorin proteins. *Front. Plant Sci.* 3:86. doi: 10.3389/fpls.2012.00086

Marín, M., Thallmair, V., Ott, T., Marin, M., Thallmair, V., and Ott, T. (2012). The intrinsically disordered N-terminal region of AtREM1.3 remorin protein mediates protein-protein interactions. *J. Biol. Chem.* 287, 39982–39991. doi: 10.1074/jbc.M112.414292

Mbengue, M., Navaud, O., Peyraud, R., Barascud, M., Badet, T., Vincent, R., et al. (2016). Emerging trends in molecular interactions between plants and the broad host range fungal pathogens botrytis cinerea and *Sclerotinia sclerotiorum*. *Front. Plant Sci.* 7:422. doi: 10.3389/fpls.2016.00422

Mine, A., Seyfferth, C., Kracher, B., Berens, M. L., Becker, D., and Tsuda, K. (2018). The defense phytohormone signaling network enables rapid, high-amplitude transcriptional reprogramming during effector-triggered immunity. *Plant Cell* 30, 1199–1219. doi: 10.1105/tpc.17.00970

Obayashi, T., Aoki, Y., Tadaka, S., Kagaya, Y., and Kinoshita, K. (2018). ATTED-II in 2018: a plant coexpression database based on investigation of the statistical property of the mutual rank index. *Plant Cell Physiol.* 59:e3. doi: 10.1093/pcp/pcx191

O'Brien, D. P., Perez, A. C. S., Karst, J., Cannella, S. E., Enguéné, V. Y. N., Hessel, A., et al. (2018). Calcium-dependent disorder-to-order transitions are central to the secretion and folding of the CyaA toxin of Bordetella pertussis, the causative agent of whooping cough. *Toxicon* 149, 37–44. doi: 10.1016/j.toxicon.2018.01.007

Pearson, W. R., and Sierk, M. L. (2005). The limits of protein sequence comparison? *Curr. Opin. Struct. Biol.* 15, 254–260. doi: 10.1016/j.sbi.2005.05.005

Perchepied, L., Balagué, C., Riou, C., Claudel-Renard, C., Rivière, N., Grezes-Besset, B., et al. (2010). Nitric oxide participates in the complex interplay of defense-related signaling pathways controlling disease resistance to *Sclerotinia sclerotiorum* in *Arabidopsis thaliana*. *Mol. Plant Microbe Interact.* 23, 846–860. doi: 10.1094/MPMI-23-7-0846

Perraki, A., Cacas, J.-L., Crowet, J.-M., Lins, L., Castroviejo, M., German-Retana, S., et al. (2012). Plasma membrane localization of *Solanum tuberosum* Remorin from group 1, homolog 3 is mediated by conformational changes in a novel C-terminal anchor and required for the restriction of potato virus X movement. *Plant Physiol.* 160, 624–637. doi: 10.1104/pp.112.200519

Perraki, A., Gronnier, J., Gouguet, P., Boudsocq, M., Deroubaix, A. F., Simon, V., et al. (2018). REM1.3's phospho-status defines its plasma membrane nanodomain organization and activity in restricting PVX cell-to-cell movement. *PLoS Pathog.* 14:e1007378. doi: 10.1371/journal.ppat.1007378

Peyraud, R., Mbengue, M., Barbacci, A., and Raffaele, S. (2019). Intercellular cooperation in a fungal plant pathogen facilitates host colonization. *Proc. Natl. Acad. Sci. U.S.A.* 116, 3193–3201. doi: 10.1073/pnas.1811267116

Poland, J. A., Balint-Kurti, P. J., Wisser, R. J., Pratt, R. C., and Nelson, R. J. (2009). Shades of gray: the world of quantitative disease resistance. *Trends Plant Sci.* 14, 21–29. doi: 10.1016/j.tplants.2008.10.006

Popescu, S. C., Popescu, G. V., Bachan, S., Zhang, Z., Seay, M., Gerstein, M., et al. (2007). Differential binding of calmodulin-related proteins to their targets revealed through high-density *Arabidopsis* protein microarrays. *Proc. Natl. Acad. Sci. U.S.A.* 104, 4730–4735. doi: 10.1073/pnas.0611615104

Potter, S. C., Luciani, A., Eddy, S. R., Park, Y., Lopez, R., and Finn, R. D. (2018). HMMER web server: 2018 update. *Nucleic Acids Res.* 46, W200–W204. doi: 10.1093/nar/gky448

Raffaele, S., Bayer, E., Lafarge, D., Cluzet, S., Retana, S. G., Boubekeur, T., et al. (2009). Remorin, a solanaceae protein resident in membrane shafts and plasmodesmata, impairs potato virus X movement. *Plant Cell Online* 21, 1541–1555. doi: 10.1105/tpc.108.064279

Raffaele, S., Mongrand, S., Gamas, P., Niebel, A., and Ott, T. (2007). Genome-wide annotation of remorins, a plant-specific protein family: evolutionary

and functional perspectives. *Plant Physiol.* 145, 593–600. doi: 10.1104/pp.107.108639

Rajarammohan, S., Pradhan, A. K., Pental, D., and Kaur, J. (2018). Genome-wide association mapping in *Arabidopsis* identifies novel genes underlying quantitative disease resistance to Alternaria brassicica. *Mol. Plant Pathol.* 19, 1719–1732. doi: 10.1111/mpp.12654

Rausch, S. (2016). The function of *Arabidopsis* microRNAs in defense against the necrotrophic fungal pathogen *Alternaria brassicicola*. *Biologie* 7, 1–126. doi: 10.15496/publikation-16913

Roux, F., Voisin, D., Badet, T., Balagué, C., Barlet, X., Huard-Chauveau, C., et al. (2014). Resistance to phytopathogens e tutti quanti: placing plant quantitative disease resistance on the map. *Mol. Plant Pathol.* 15, 427–432. doi: 10.1111/mpp.12138

Scholz, S. S., Schmidt-Heck, W., Guthke, R., Furch, A. C. U., Reichelt, M., Gershenzon, J., et al. (2018). Verticillium dahliae-*Arabidopsis* interaction causes changes in gene expression profiles and jasmonate levels on different time scales. *Front. Microbiol.* 9:217. doi: 10.3389/fmicb.2018.00217

Shanks, C. M., Rice, J. H., Zubo, Y., Schaller, G. E., Hewezi, T., and Kieber, J. J. (2016). The Role of cytokinin during infection of *Arabidopsis thaliana* by the cyst nematode heterodera schachtii. *Mol. Plant Microbe Interact.* 29, 57–68. doi: 10.1094/MPMI-07-15-0156-R

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software Environment for integrated models of biomolecular interaction networks. *Genome Res.* 13, 2498–2504. doi: 10.1101/gr.1239303

Sun, X., Greenwood, D. R., Templeton, M. D., Libich, D. S., McGhie, T. K., Xue, B., et al. (2014). The intrinsically disordered structural platform of the plant defence hub protein RPM1-interacting protein 4 provides insights into its mode of action in the host-pathogen interface and evolution of the nitrate-induced domain protein family. *FEBS J.* 281, 3955–3979. doi: 10.1111/febs.12937

Sun, X., Rikkerink, E. H. A., Jones, W. T., and Uversky, V. N. (2013). Multifarious roles of intrinsic disorder in proteins illustrate its broad impact on plant biology. *Plant Cell* 25, 38–55. doi: 10.1105/tpc.112.106062

Suzuki, N., Bassil, E., Hamilton, J. S., Inupakutika, M. A., Zandalinas, S. I., Tripathy, D., et al. (2016). ABA is required for plant acclimation to a combination of salt and heat stress. *PLoS One* 11:e0147625. doi: 10.1371/journal.pone.0147625

Tavridou, E., Pireyre, M., and Ulm, R. (2020). Degradation of the transcription factors PIF4 and PIF5 under UV-B promotes UVR8-mediated inhibition of hypocotyl growth in *Arabidopsis*. *Plant J.* 101, 507–517. doi: 10.1111/tpj.14556

Tenhaken, R., Doerks, T., and Bork, P. (2005). DCD - A novel plant specific domain in proteins involved in development and programmed cell death. *BMC Bioinformatics* 6:169. doi: 10.1186/1471-2105-6-169

The Uniprot Consortium. (2019). UniProt: a worldwide hub of protein knowledge. *Nucleic Acids Res.* 47, D506–D515. doi: 10.1093/nar/gky1049

Toruño, T. Y., Shen, M., Coaker, G., and Mackey, D. (2019). Regulated disorder: posttranslational modifications control the RIN4 plant immune signaling hub. *Mol. Plant Microbe Interact.* 32, 56–64. doi: 10.1094/MPMI-07-18-0212-FI

Uversky, V. N. (2013). A decade and a half of protein intrinsic disorder: biology still waits for physics. *Protein Sci.* 22, 693–724. doi: 10.1002/pro.2261

Wang, J., Qin, J., Sun, P., Ma, X., Yu, J., Li, Y., et al. (2019). Polyploidy index and its implications for the evolution of polyploids. *Front. Genet.* 10:807. doi: 10.3389/fgene.2019.00807

Weidenbach, D., Esch, L., Möller, C., Hensel, G., Kumlehn, J., Höfle, C., et al. (2016). Polarized defense against fungal pathogens is mediated by the jacalin-related lectin domain of modular poaceae-specific proteins. *Mol. Plant* 9, 514–527. doi: 10.1016/j.molp.2015.12.009

Xu, X., Yang, Y., Liu, C., Sun, Y., Zhang, T., Hou, M., et al. (2019). The evolutionary history of the sucrose synthase gene family in higher plants. *BMC Plant Biol.* 19:566. doi: 10.1186/s12870-019-2181-4

Yamasaki, K., Kigawa, T., Seki, M., Shinozaki, K., and Yokoyama, S. (2013). DNA-binding domains of plant-specific transcription factors: structure, function, and evolution. *Trends Plant Sci.* 18, 267–276. doi: 10.1016/j.tplants.2012.09.001

Yang, H., Li, Y., and Hua, J. (2006). The C2 domain protein BAP1 negatively regulates defense responses in *Arabidopsis*. *Plant J.* 48, 238–248. doi: 10.1111/j.1365-313X.2006.02869.x
Yang, Q., Niu, X., Tian, X., Zhang, X., Cong, J., Wang, R., et al. (2020). Comprehensive genomic analysis of the DUF4228 gene family in land plants and expression profiling of ATDUF4228 under abiotic stresses. *BMC Genomics* 21:12. doi: 10.1186/s12864-019-6389-3

Zhang, M., Mo, H., Sun, W., Guo, Y., and Li, J. (2016). Systematic isolation and characterization of cadmium tolerant genes in tobacco: a cDNA library construction and screening approach. *PLoS One* 11:e0161147. doi: 10.1371/journal.pone.0161147

Zhang, Z., Feechan, A., Pedersen, C., Newman, M. A., Qiu, J. L., Olesen, K. L., et al. (2007). A SNARE-protein has opposing functions in penetration resistance and defence signalling pathways. *Plant J.* 49, 302–312. doi: 10.1111/j.1365-313X.2006.02961.x

Zorzatto, C., MacHado, J. P. B., Lopes, K. V. G., Nascimento, K. J. T., Pereira, W. A., Brustolini, O. J. B., et al. (2015). NIK1-mediated translation suppression functions as a plant antiviral immunity mechanism. *Nature* 520, 679–682. doi: 10.1038/nature14171

Zuther, E., Schaarschmidt, S., Fischer, A., Erban, A., Pagter, M., Mubeen, U., et al. (2019). Molecular signatures associated with increased freezing tolerance due to low temperature memory in *Arabidopsis*. *Plant Cell Environ.* 42, 854–873. doi: 10.1111/pce.13502

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Didelon, Khafif, Godiard, Barbacci and Raffaele. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.