The unicellular ancestor of animals had a complex repertoire of genes linked to multicellular processes. This suggests that changes in the regulatory genome, rather than in gene innovation, were key to the origin of animals. Here, we carry out multiple functional genomic assays in Capsaspora owczarzaki, the unicellular relative of animals with the largest known gene repertoire for transcriptional regulation. We show that changing chromatin states, differential lincRNA expression, and dynamic cis-regulatory sites are associated with life cycle transitions in Capsaspora. Moreover, we demonstrate conservation of animal developmental transcription-factor networks and extensive network interconnection in this premetazoan organism. In contrast, however, Capsaspora lacks animal promoter types, and its regulatory sites are small, proximal, and lack signatures of animal enhancers. Overall, our results indicate that the emergence of animal multicellularity was linked to a major shift in genome cis-regulatory complexity, most notably the appearance of distal enhancer regulation.

This work was supported by an Institució Catalana de Recerca i Estudis Avançats contract, a European Research Council Consolidator Grant (ERC-2012-Co-616960), and a grant from Ministerio de Economía y Competitividad (MINECO) (BFU-2011-23434) (to I.R.-T.). We also acknowledge financial support from Secretaria d’Universitats i Recerca del Departament d’Economia i Coneixement de la Generalitat de Catalunya (project 2014 SGR 619). The work in L.D.’s laboratory was supported by grants from the Spanish “Ministerio de Educación y Ciencia” (SAF2013-48926-P), AGAUR, and the European Commission’s 7th Framework Program 4DCellFate (277899). A.S.-P. is supported by an EMBO Long-Term Fellowship (ALTF 841-2014). J.L.G.-S. was funded by grants from Ministerio de Economía y Competitividad (BFU2013-41322-P) and the Andalusian Government (BIO-396). J.J.T. has a postdoctoral grant from the University Pablo de Olavide. The CRG/UPF Proteomics Unit is part of the “Plataforma de Recursos Biomoleculares y Bioinformáticos (ProteoRed)” supported by a grant from Instituto de Salud Carlos III (ISCIII) (PT13/0001). We thank Guadalupe Espadas for her support with the histone derivatization protocol and Núria Ros and Meritxell Antó for technical support. Finally, we thank the CRG Genomics Unit for helping with ChIP-seq and RNA-seq sequencing.