Transcriptome dynamics in developing leaves from C and C Flaveria species.

C species have evolved more than 60 times independently from C ancestors. This multiple and parallel evolution of the complex C trait suggests common underlying evolutionary mechanisms, which could be identified by comparative analysis of closely related C and C species. Efficient C function depends on a distinctive leaf anatomy that is characterised by enlarged, chloroplast-rich bundle sheath cells and narrow vein spacing. To elucidate the molecular mechanisms that generate the Kranz anatomy, we analysed a developmental series of leaves from the C plant Flaveria bidentis and the closely related C species Flaveria robusta by comparing anatomies and transcriptomes. Vascular density measurements of all nine leaf developmental stages identified three leaf anatomical zones whose proportions vary with respect to the developmental stage. We then deconvoluted the transcriptome datasets using non-negative matrix factorisation, which identified four distinct transcriptome patterns in the growing leaves of both species. By integrating the leaf anatomy and transcriptome data, we were able to correlate the different transcriptional profiles with different developmental zones in the leaves. These comparisons revealed an important role for auxin metabolism, in particular auxin homeostasis (conjugation and deconjugation), in establishing the high vein density typical of C species.