Aerial plant architecture is predominantly determined by shoot branching and leaf morphology, which are governed by apparently unrelated developmental processes, axillary meristem formation, and leaf dissection. Here, we show that in tomato (Solanum lycopersicum), these processes share essential functions in boundary establishment. Potato leaf (C), a key regulator of leaf dissection, was identified to be the closest paralog of the shoot branching regulator Blind (Bl). Comparative genomics revealed that these two R2R3 MYB genes are orthologs of the Arabidopsis thaliana branching regulator REGULATOR OF AXILLARY MERISTEMS1 (RAX1). Expression studies and complementation analyses indicate that these genes have undergone sub- or neofunctionalization due to promoter differentiation. C acts in a pathway independent of other identified leaf dissection regulators. Furthermore, the known leaf complexity regulator Goblet (Gob) is crucial for axillary meristem initiation and acts in parallel to C and Bl. Finally, RNA in situ hybridization revealed that the branching regulator Lateral suppressor (Ls) is also expressed in leaves. All four boundary genes, C, Bl, Gob, and Ls, may act by suppressing growth, as indicated by gain-of-function plants. Thus, leaf architecture and shoot architecture rely on a conserved mechanism of boundary formation preceding the initiation of leaflets and axillary meristems.

The R2R3MYB proteins comprise one of the largest families of transcription factors and play regulatory roles in developmental processes and defense responses in plants. However, there has been relatively little effort to systematically carry out comprehensive genomic and functional analyses of these genes in tomato (Solanum lycopersicum L.), a reference species for Solanaceae plants, and the model plant for fruit development. In this study, a total of 121 R2R3MYB genes were identified in the tomato genome released recently and further classified into 29 subgroups based on the phylogenetic analysis of the complete protein sequences. Phylogenetic comparison of the members of this superfamily among tomato, Arabidopsis, grape, rice, poplar, soybean, cucumber and apple revealed that the putative functions of some tomato R2R3MYB proteins were clustered into the Arabidopsis functional clades. The chromosome distribution pattern revealed that tomato R2R3MYB genes were enriched on several chromosomes and 52 % of the family members were tandemly duplicated genes. Tissue specificity or different expression levels of SlR2R3MYBs in different tissues suggested differential regulation of tissue development as well as metabolic regulation. The transcript abundance level analysis during abiotic conditions identified a group of R2R3MYB genes that responded to one or more treatments suggesting that the SlR2R3MYBs played major roles in the plant response to abiotic conditions and involved in signal transduction pathways. This study not only provides a solid foundation for further functional dissection of tomato R2R3MYB family genes, but may also be profitable for, in the future, the improvement of tomato stress tolerance and fruit quality.