Genome-Wide Survey of the RWP-RK Gene Family in Cassava (Manihot esculenta Crantz) and Functional Analysis

The plant-specific RWP-RK transcription factor family plays a central role in the regulation of nitrogen response and gametophyte development. However, little information is available regarding the evolutionary relationships and characteristics of the RWP-RK family genes in cassava, an important tropical crop. Herein, 13 RWP-RK proteins identified in cassava were unevenly distributed across 9 of the 18 chromosomes (Chr), and these proteins were divided into two clusters based on their phylogenetic distance. The NLP subfamily contained seven cassava proteins including GAF, RWP-RK, and PB1 domains; the RKD subfamily contained six cassava proteins including the RWP-RK domain. Genes of the NLP subfamily had a longer sequence and more introns than the RKD subfamily. A large number of hormone- and stress-related cis-acting elements were found in the analysis of RWP-RK promoters. Real-time quantitative PCR revealed that all MeNLP1-7 and MeRKD1/3/5 genes responded to different abiotic stressors (water deficit, cold temperature, mannitol, polyethylene glycol, NaCl, and H₂O₂), hormonal treatments (abscisic acid and methyl jasmonate), and nitrogen starvation. MeNLP3/4/5/6/7 and MeRKD3/5, which can quickly and efficiently respond to different stresses, were found to be important candidate genes for further functional assays in cassava. The MeRKD5 and MeNLP6 proteins were localized to the cell nucleus in tobacco leaf. Five and one candidate proteins interacting with MeRKD5 and MeNLP6, respectively, were screened from the cassava nitrogen starvation library, including agamous-like mads-box protein AGL14, metallothionein 2, Zine finger FYVE domain containing protein, glyceraldehyde-3-phosphate dehydrogenase, E3 Ubiquitin-protein ligase HUWE1, and PPR repeat family protein. These results provided a solid basis to understand abiotic stress responses and signal transduction mediated by RWP-RK genes in cassava.