Custom methods to identify conserved genetic modules applied to novel transcriptomic data from Amborella trichopoda

We have devised a procedure for the inter-species comparison of transcriptomic data and used this procedure to reconstruct the expression dynamics of major genetic modules that were present at least 149 million years ago in the most recent common ancestor of living angiosperms. We began by using laser-assisted microdissection to generate novel transcriptomic data from female flower tissues of Amborella trichopoda, the likely sister to all other living angiosperms. We then employed a gene-expression clustering method, followed by a custom procedure to compare genetic modules on the basis of gene orthology between Amborella and the molecular-genetic model angiosperm Arabidopsis thaliana. Using this protocol, we succeeded in identifying nine major genetic modules that appear to have conserved their expression dynamics from an early stage in angiosperm evolution. The genes of these modules, representing over 5000 orthogroups, include around one third of those known to control female reproductive development in Arabidopsis. Our study constitutes a proof of concept for the comparison of transcriptomic data between widely diverged plant species and represents a first step in the large-scale analysis of gene expression dynamics in a macro-evolutionary context.