Genome-Wide Identification and Expression Analysis of the R2R3-MYB Gene Family in Rubber Trees

The plant MYB transcription factor family featured as highly conserved DNA-binding domains consisting of 1 to 4 imperfect repeats (R). Increasing evidence indicates that MYB genes participates in growth, differentiation, metabolism, and biotic and abiotic stress responses. However, the functions of MYB genes in the rubber tree remain to be deeply elucidated, especially R2R3-MYB gene family. In this study, molecular biology, bioinformatics, and qRT-PCR were used to identify and analyze HbR2R3-MYB gene family members in the rubber tree. A total of 132 members of the R2R3-MYB gene family were identified in the rubber tree based on genome-wide level. Most of the HbR2R3-MYBs were mapped to 17 rubber tree chromosomes except four genes. A phylogenetic analysis divided all the HbR2R3-MYBs into 20 subgroups with Arabidopsis thaliana. MEME analysis showed that the protein of HbR2R3-MYBs was characterized by 9 conserved motifs. Twenty-six representative R2R3 HbMYBs from different subgroups were selected for expression profiles analysis and the results revealed that the HbR2R3-MYBs members showed various expression patterns in different tissues, powdery mildew-infected and ethylene treatment, implying the diversity of their functions in rubber trees. These results provide fundamental knowledge for further studying the response of the HbR2R3-MYB family to stress and regulation latex flow in rubber tree.