[Identification and expression analysis of R2R3-MYB gene family in Andrographis paniculata]

The plant growth, development, and secondary metabolism are regulated by R2 R3-MYB transcription factors. This study identified the R2 R3-MYB genes in the genome of Andrographis paniculata and analyzed the chromosomal localization, gene structure, and conserved domains, phylogenetic relationship, and promoter cis-acting elements of these R2 R3-MYB genes. Moreover, the gene expression profiles of R2 R3-MYB genes under abiotic stress and hormone treatments were generated by RNA-seq and validated by qRT-PCR. The results showed that A. paniculata contained 73 R2 R3-MYB genes on 21 chromosomes. These members belonged to 34 subfamilies, 19 of which could be classified into the known subfamilies in Arabidopsis thaliana. The 73 R2 R3-MYB members included 36 acidic proteins and 37 basic proteins, with the lengths of 148-887 aa. The domains, motifs, and gene structures of R2 R3-MYBs in A. paniculata were conserved. The promoter regions of these genes contains a variety of cis-acting elements related to the responses to environmental factors and plant hormones including light, ABA, MeJA, and drought. Based on the similarity of functions of R2 R3-MYBs in the same subfamily and the transcription profiles, ApMYB13/21/35/67/73(S22) may regulate drought stress through ABA pathway; ApMYB20(S11) and ApMYB55(S2) may play a role in the response of A. paniculata to high temperature and UV-C stress; ApMYB5(S7) and ApMYB33(S20) may affect the accumulation of andrographolide by regulating the expression of key enzymes in the MEP pathway. This study provides theoretical reference for further research on the functions of R2 R3-MYB genes in A. paniculata and breeding of A. paniculata varieties with high andrographolide content.