Constraint and diversification of developmental trajectories in cichlid facial morphologies

Background A major goal of evolutionary biology is to understand the origins of phenotypic diversity. Changes in development, for instance heterochrony, can be a potent source of phenotypic variation. On the other hand, development can also constrain the spectrum of phenotypes that can be produced. In order to understand these dual roles of development in evolution, we examined the developmental trajectory of a trait central to the extensive adaptive radiation of East African cichlid fishes: craniofacial adaptations that allow optimal exploitation of ecological niches. Specifically, we use geometric morphometric analysis to compare morphological ontogenies among six species of Lake Malawi cichlids (n \u3e 500 individuals) that span a major ecomorphological axis. We further evaluate how modulation of Wnt signaling impacts the long-term developmental trajectory of facial development. Results We find that, despite drastic differences in adult craniofacial morphologies, there are general similarities in the path of craniofacial ontogeny among species, suggesting that natural selection is working within a conserved developmental program. However, we also detect species-specific differences in the timing, direction, and/or duration of particular developmental trajectories, including evidence of heterochrony. Previous work in cichlids and other systems suggests that species-specific differences in adult morphology are due to changes in molecular signaling pathways that regulate early craniofacial development. In support of this, we demonstrate that modulation of Wnt signaling at early stages can shift a developmental trajectory into morphospace normally occupied by another species. However, without sustained modulation, craniofacial shape can recover by juvenile stages. This underscores the idea that craniofacial development is robust and that adult head shapes are the product of many molecular changes acting over extended periods of development. Conclusions Our results are consistent with the hypothesis that development acts to both constrain and promote morphological diversity. They also illustrate the modular nature of the craniofacial skeleton and hence the ability of selection to act upon distinct anatomical features in an independent manner. We propose that trophic diversity among cichlids has been achieved via shifts in both specific (e.g., stage-specific changes in gene expression) and global (e.g., heterochrony) ontogenetic processes acting within a conserved developmental program

Genetic analyses in Lake Malawi cichlids identify new roles for Fgf signaling in scale shape variation

Elasmoid scales are the most common epithelial appendage among vertebrates, however an understanding of the genetic mechanisms that underlie variation in scale shape is lacking. Using an F2 mapping cross between morphologically distinct cichlid species, we identified \u3e40 QTL for scale shape at different body positions. We show that while certain regions of the genome regulate variation in multiple scales, most are specific to scales at distinct positions. This suggests a degree of regional modularity in scale development. We also identified a single QTL for variation in scale shape disparity across the body. Finally, we screened a QTL hotspot for candidate loci, and identified the Fgf receptor fgfr1b as a prime target. Quantitative rtPCR and small molecule manipulation support a role for Fgf signaling in shaping cichlid scales. While Fgfs have previously been implicated in scale loss, these data reveal new roles for the pathway in scale shape variation

An RNA interference-based screen of transcription factor genes identifies pathways necessary for sensory regeneration in the avian inner ear

Sensory hair cells of the inner ear are the mechano-electric transducers of sound and head motion. In mammals, damage to sensory hair cells leads to hearing or balance deficits. Non-mammalian vertebrates such as birds can regenerate hair cells after injury. In a previous study, we characterized transcription factor gene expression during chicken hair cell regeneration. In those studies, a laser micro-beam or ototoxic antibiotics were used to damage the sensory epithelia (SE). The current study focused on 27 genes that were up-regulated in regenerating SE compared to untreated SE in the previous study. Those genes were knocked down by siRNA, to determine their requirement for supporting cell proliferation and to measure resulting changes in the larger network of gene expression. We identified 11 genes necessary for proliferation and also identified novel interactive relationships between many of them. Defined components of the WNT, PAX and AP1 pathways were shown to be required for supporting cell proliferation. These pathways intersect on WNT4, which is also necessary for proliferation. Among the required genes, the CCAAT enhancer binding protein, CEBPG, acts downstream of Jun Kinase and JUND in the AP1 pathway. The WNT co-receptor LRP5 acts downstream of CEBPG as does the transcription factor BTAF1. Both of these genes are also necessary for supporting cell proliferation. This is the first large scale screen of its type and suggests an important intersection between the AP1 pathway, the PAX pathway and WNT signaling in the regulation of supporting cell proliferation during inner ear hair cell regeneration

Large Scale Gene Expression Profiles of Regenerating Inner Ear Sensory Epithelia

Loss of inner ear sensory hair cells (HC) is a leading cause of human hearing loss and balance disorders. Unlike mammals, many lower vertebrates can regenerate these cells. We used cross-species microarrays to examine this process in the avian inner ear. Specifically, changes in expression of over 1700 transcription factor (TF) genes were investigated in hair cells of auditory and vestibular organs following treatment with two different damaging agents and regeneration in vitro. Multiple components of seven distinct known signaling pathways were clearly identifiable: TGFβ, PAX, NOTCH, WNT, NFKappaB, INSULIN/IGF1 and AP1. Numerous components of apoptotic and cell cycle control pathways were differentially expressed, including p27KIP and TFs that regulate its expression. A comparison of expression trends across tissues and treatments revealed identical patterns of expression that occurred at identical times during regenerative proliferation. Network analysis of the patterns of gene expression in this large dataset also revealed the additional presence of many components (and possible network interactions) of estrogen receptor signaling, circadian rhythm genes and parts of the polycomb complex (among others). Equal numbers of differentially expressed genes were identified that have not yet been placed into any known pathway. Specific time points and tissues also exhibited interesting differences: For example, 45 zinc finger genes were specifically up-regulated at later stages of cochlear regeneration. These results are the first of their kind and should provide the starting point for more detailed investigations of the role of these many pathways in HC recovery, and for a description of their possible interactions

A Cross-Species Analysis of MicroRNAs in the Developing Avian Face

Higher vertebrates use similar genetic tools to derive very different facial features. This diversity is believed to occur through temporal, spatial and species-specific changes in gene expression within cranial neural crest (NC) cells. These contribute to the facial skeleton and contain species-specific information that drives morphological variation. A few signaling molecules and transcription factors are known to play important roles in these processes, but little is known regarding the role of micro-RNAs (miRNAs). We have identified and compared all miRNAs expressed in cranial NC cells from three avian species (chicken, duck, and quail) before and after species-specific facial distinctions occur. We identified 170 differentially expressed miRNAs. These include thirty-five novel chicken orthologs of previously described miRNAs, and six avian-specific miRNAs. Five of these avian-specific miRNAs are conserved over 120 million years of avian evolution, from ratites to galliforms, and their predicted target mRNAs include many components of Wnt signaling. Previous work indicates that mRNA gene expression in NC cells is relatively static during stages when the beak acquires species-specific morphologies. However, miRNA expression is remarkably dynamic within this timeframe, suggesting that the timing of specific developmental transitions is altered in birds with different beak shapes. We evaluated one miRNA:mRNA target pair and found that the cell cycle regulator p27KIP1 is a likely target of miR-222 in frontonasal NC cells, and that the timing of this interaction correlates with the onset of phenotypic variation. Our comparative genomic approach is the first comprehensive analysis of miRNAs in the developing facial primordial, and in species-specific facial development

Data from: Genetic analyses in Lake Malawi cichlids identify new roles for Fgf signaling in scale shape variation

Elasmoid scales are the most common epithelial appendage among vertebrates, however an understanding of the genetic mechanisms that underlie variation in scale shape is lacking. Using an F2 mapping cross between morphologically distinct cichlid species, we identified \u3e40 QTL for scale shape at different body positions. We show that while certain regions of the genome regulate variation in multiple scales, most are specific to scales at distinct positions. This suggests a degree of regional modularity in scale development. We also identified a single QTL for variation in scale shape disparity across the body. Finally, we screened a QTL hotspot for candidate loci, and identified the Fgf receptor fgfr1b as a prime target. Quantitative rtPCR and small molecule manipulation support a role for Fgf signaling in shaping cichlid scales. While Fgfs have previously been implicated in scale loss, these data reveal new roles for the pathway in scale shape variation.,Parental_F2_Residuals_dataraw data for scale shape measures for parental species and a subset of F2 hybrids.LFxTRC_F2_scale_pheno_dataF2 phenotypes for QTL mapping. Genotypic data are under a different DOI.F2_scale1-6raw landmark data for all F2 hybrids in TPS format.LF_TRC_parental_scale1-6raw landmark data for parental scales in TPS format.SU5402_treatmentraw phenotypic data for scales post small molecule treatment.Cichlid_Scale_qPCR_DataqPCR data used to compare scale types between and within species.Parental_F2_sex_SLSex and standard length of fish used in this studyparental_F2scale_numberScale number in the fish used in this study

Genetic basis of continuous variation in the levels and modular inheritance of pigmentation in cichlid fishes

Variation in pigmentation type and levels is a hallmark of myriad evolutionary radiations, and biologists have long been fascinated by the factors that promote and maintain variation in coloration across populations. Here, we provide insights into the genetic basis of complex and continuous patterns of colour variation in cichlid fishes, which offer a vast diversity of pigmentation patterns that have evolved in response to both natural and sexual selection. Specifically, we crossed two divergent cichlid species to generate an F2 mapping population that exhibited extensive variation in pigmentation levels and patterns. Our experimental design is robust in that it combines traditional quantitative trait locus (QTL) analysis with population genomics, which has allowed us to move efficiently from QTL interval to candidate gene. In total, we detected 41 QTL and 13 epistatic interactions that underlie melanocyte- and xanthophore-based coloration across the fins and flanks of these fishes. We also identified 2 QTL and 1 interaction for variation in the magnitude of integration among these colour traits. This finding in particular is notable as there are marked differences both within and between species with respect to the complexity of pigmentation patterns. While certain individuals are characterized by more uniform ‘integrated’ colour patterns, others exhibit many more degrees of freedom with respect to the distribution of colour ‘modules’ across the fins and flank. Our data reveal, for the first time, a genetic basis for this difference. Finally, we implicate pax3a as a mediator of continuous variation in the levels of xanthophore-based colour along the cichlid flank

Genetic basis of continuous variation in the levels and modular inheritance of pigmentation in cichlid fishes

Variation in pigmentation type and levels is a hallmark of myriad evolutionary radiations, and biologists have long been fascinated by the factors that promote and maintain variation in coloration across populations. Here, we provide insights into the genetic basis of complex and continuous patterns of colour variation in cichlid fishes, which offer a vast diversity of pigmentation patterns that have evolved in response to both natural and sexual selection. Specifically, we crossed two divergent cichlid species to generate an F2 mapping population that exhibited extensive variation in pigmentation levels and patterns. Our experimental design is robust in that it combines traditional quantitative trait locus (QTL) analysis with population genomics, which has allowed us to move efficiently from QTL interval to candidate gene. In total, we detected 41 QTL and 13 epistatic interactions that underlie melanocyte- and xanthophore-based coloration across the fins and flanks of these fishes. We also identified 2 QTL and 1 interaction for variation in the magnitude of integration among these colour traits. This finding in particular is notable as there are marked differences both within and between species with respect to the complexity of pigmentation patterns. While certain individuals are characterized by more uniform ‘integrated’ colour patterns, others exhibit many more degrees of freedom with respect to the distribution of colour ‘modules’ across the fins and flank. Our data reveal, for the first time, a genetic basis for this difference. Finally, we implicate pax3a as a mediator of continuous variation in the levels of xanthophore-based colour along the cichlid flank

allele specific expression data

Data based on Sanger sequencing for estimating the relative expression of different alleles

Interactions among multiple selective pressures on the form–function relationship in insular stream fishes

Relationships between body shape and escape performance are well established for many species. However, organisms can face multiple selection pressures that might impose competing demands. Many fishes use fast starts for escaping predator attacks, whereas some species of gobiid fishes have evolved the ability to climb waterfalls out of predator-dense habitats. The ancestral ‘powerburst’ climbing mechanism uses lateral body undulations to move up waterfalls, whereas a derived ‘inching’ mechanism uses rectilinear locomotion. We examined whether fast-start performance is impacted by selection imposed from the new functional demands of climbing. We predicted that non-climbing species would show morphology and fast-start performance that facilitate predator evasion, because these fish live consistently with predators and are not constrained by the demands of climbing. We also predicted that, by using lateral undulations, powerburst climbers would show escape performance superior to that of inchers. We compared fast starts and body shape across six goby species. As predicted, non-climbing fish exhibited distinct morphology and responded more frequently to an attack stimulus than climbing species. Contrary to our predictions, we found no differences in escape performance among climbing styles. These results indicate that selection for a competing pressure need not limit the ability of prey to escape predator attacks