Genome-Wide Identification of GYF-Domain Encoding Genes in Three Brassica Species and Their Expression Responding to Sclerotinia sclerotiorum in Brassica napus

GYF (glycine-tyrosine-phenylalanine)-domain-containing proteins, which were reported to participate in many aspects of biological processes in yeast and animals, are highly conserved adaptor proteins existing in almost all eukaryotes. Our previous study revealed that GYF protein MUSE11/EXA1 is involved in nucleotide-binding leucine-rich repeat (NLR) receptor-mediated defense in Arabidopsis thaliana. However, the GYF-domain encoding homologous genes are still not clear in other plants. Here, we performed genome-wide identification of GYF-domain encoding genes (GYFs) from Brassica napus and its parental species, Brassica rapa and Brassica oleracea. As a result, 26 GYFs of B. napus (BnaGYFs), 11 GYFs of B. rapa (BraGYFs), and 14 GYFs of B. oleracea (BolGYFs) together with 10 A. thaliana (AtGYFs) were identified, respectively. We, then, conducted gene structure, motif, cis-acting elements, duplication, chromosome localization, and phylogenetic analysis of these genes. Gene structure analysis indicated the diversity of the exon numbers of these genes. We found that the defense and stress responsiveness element existed in 23 genes and also identified 10 motifs in these GYF proteins. Chromosome localization exhibited a similar distribution of BnaGYFs with BraGYFs or BolGYFs in their respective genomes. The phylogenetic and gene collinearity analysis showed the evolutionary conservation of GYFs among B. napus and its parental species as well as Arabidopsis. These 61 identified GYF domain proteins can be classified into seven groups according to their sequence similarity. Expression of BnaGYFs induced by Sclerotinia sclerotiorum provided five highly upregulated genes and five highly downregulated genes, which might be candidates for further research of plant–fungal interaction in B. napus.