Large-scale analysis of MYB genes in Cucurbitaceae identifies a novel gene regulating plant height

The MYB transcription factor (TF) family, which is involved in plant growth and development, is large and diverse. Previous studies on MYB family in Cucurbitaceae were mostly based on a single genome or focused on the R2R3 subfamily. Here, we analyzed 91 genomes of 11 Cucurbitaceae species and identified a total of 15 858 MYB genes. According to phylogenetic relationships, these genes were divided into 27 subgroups. The identified MYB genes were further classified into 121 MYB orthologous gene groups (OGGs), including 25 core, 57 softcore, 19 shell and 20 line-specific/cloud groups. Whole-genome duplication was the most common mechanism of MYB genes expansion. In core group, the higher proportions of MYB genes were found to be in the coexpression network constructed by the RNA-seq data. Through the comprehensive analysis including phylogeny and gene expression profile of cucumber MYB genes, as well as genetic variations in 103 cucumber germplasms, we identified a MYB gene CsRAX5, which may be related to cucumber plant height. We used gene editing technology to knockout and overexpress CsRAX5. In the knockout lines, Csrax5, the height was significantly increased compared with wild type (WT), whereas after overexpression the height of CsRAX5-OE plants was significantly decreased compared with WT. These results indicated that MYB gene CsRAX5 negatively regulated cucumber plant height. The large-scale analysis of MYB genes in Cucurbitaceae in this study provides insights for further investigating the evolution and function of MYB genes in Cucurbitaceae crops.