Seminal studies in vertebrate protein evolution concluded that gene regulatory changes likely drive anatomical innovations. Even in the unparalleled East African cichlid fish radiation, we demonstrate cis-regulatory divergence as a contributor to phenotypic diversity. To further investigate this mechanism, we extended the Arboretum algorithm, initially designed for yeast adaptation to study the evolution of regulatory expression divergence in complex vertebrate species. For the first time, we reconstruct tissuespecific gene regulatory networks underpinning evolutionary adaptations from multiple vertebrate species in a phylogeny. Our framework consists of identifying ancestral reconstructed and extant species gene co-expression modules and integrating their associated regulators to investigate gene regulatory network evolution contributing to traits of phenotypic diversity in the East African cichlids. Along the phylogeny, we identify tissue-specific co-expression patterns across six tissues of five cichlids species that were predicted to be regulated by divergent suites of regulators. We report striking cases of rapid network rewiring for adaptive trait genes, such as the visual system. In regulatory regions of visual opsin genes e.g. sws1, in vitro assays confirm that single nucleotide polymorphisms (SNPs) in transcription factor binding sites (TFBSs) have driven network rewiring between species sharing the same visual palette. SNPs overlapping TFBSs segregate according to phylogeny and ecology suggesting ecotype-associated network rewiring in East African cichlid radiations. Our unique integrative approach inferring multispecies regulatory networks allowed us to identify regulatory changes associated with traits of phenotypic diversity in radiating cichlids.These species occupy a large diversity of ecological niches and differ dramatically in phenotypic traits, including skeletal morphology, dentition, color patterning, and a range of behavioral traits. Such explosive phenotypic diversification of East African cichlids is unparalleled among vertebrates making cichlid fishes an ideal model to investigate the evolution of gene regulatory networks associated with phenotypes under selection.Sequencing and analyses of representative East African cichlid genomes and transcriptomes suggested that bursts of gene duplication events combined with the rapid evolution of regulatory elements and selected protein coding genes contributed to evolutionary innovations 17 . This, coupled with low levels of nucleotide diversity between cichlid species pairs (Lake Malawi cichlids, 0.1-0.25%) 18 , implies high impact of low level regulatory sequence divergence to gene expression diversity. This likely fuels GRN rewiring, serving as a substrate for the evolution of phenotypic diversity in cichlids.To investigate the functional implications of changes in gene regulation on phenotypic and ecotypic diversity in cichlids, we developed a framework to characterize changes in gene regulation at the sequence, regulatory edge, expression and...