A nonsynonymous mutation in the transcriptional regulator lbh is associated with cichlid craniofacial adaptation and neural crest cell development

Molecular Innovations Shaping Beak Morphology in Birds

The beak, a pivotal evolutionary trait characterized by high morphological diversity and plasticity, has enabled birds to survive mass extinction events and subsequently radiate into diverse ecological niches worldwide. This remarkable ecological adaptability underscores the importance of uncovering the molecular mechanisms shaping avian beak morphology, particularly benefiting from the rapidly advancing archives of genomics and epigenomics. We review the latest advancements in understanding how genetic and epigenetic innovations control or regulate beak development and drive beak morphological adaptation and diversification over the past two decades. We conclude with several recommendations for future endeavors, expanding to more bird lineages, with a focus on beak shape and the lower beak, and conducting functional experiments. By directing research efforts toward these aspects and integrating advanced omics techniques, the complex molecular mechanisms involved in avian beak evolution and morphogenesis will be deeply interpreted.

Genetic and Developmental Divergence in the Neural Crest Program between Cichlid Fish Species

Neural crest (NC) is a vertebrate-specific embryonic progenitor cell population at the basis of important vertebrate features such as the craniofacial skeleton and pigmentation patterns. Despite the wide-ranging variation of NC-derived traits across vertebrates, the contribution of NC to species diversification remains underexplored. Here, leveraging the adaptive diversity of African Great Lakes' cichlid species, we combined comparative transcriptomics and population genomics to investigate the evolution of the NC genetic program in the context of their morphological divergence. Our analysis revealed substantial differences in transcriptional landscapes across somitogenesis, an embryonic period coinciding with NC development and migration. This included dozens of genes with described functions in the vertebrate NC gene regulatory network, several of which showed signatures of positive selection. Among candidates showing between-species expression divergence, we focused on teleost-specific paralogs of the NC-specifier sox10 (sox10a and sox10b) as prime candidates to influence NC development. These genes, expressed in NC cells, displayed remarkable spatio-temporal variation in cichlids, suggesting their contribution to interspecific morphological differences, such as craniofacial structures and pigmentation. Finally, through CRISPR/Cas9 mutagenesis, we demonstrated the functional divergence between cichlid sox10 paralogs, with the acquisition of a novel skeletogenic function by sox10a. When compared with teleost models zebrafish and medaka, our findings reveal that sox10 duplication, although retained in most teleost lineages, had variable functional fates across their phylogeny. Altogether, our study suggests that NC-related processes-particularly those controlled by sox10s-are involved in generating morphological diversification between species and lays the groundwork for further investigations into the mechanisms underpinning vertebrate NC diversification.

Single-cell transcriptome analysis of cavernous tissues reveals the key roles of pericytes in diabetic erectile dysfunction

Erectile dysfunction (ED) affects a significant proportion of men aged 40-70 and is caused by cavernous tissue dysfunction. Presently, the most common treatment for ED is phosphodiesterase 5 inhibitors; however, this is less effective in patients with severe vascular disease such as diabetic ED. Therefore, there is a need for development of new treatment, which requires a better understanding of the cavernous microenvironment and cell-cell communications under diabetic condition. Pericytes are vital in penile erection; however, their dysfunction due to diabetes remains unclear. In this study, we performed single-cell RNA sequencing to understand the cellular landscape of cavernous tissues and cell type-specific transcriptional changes in diabetic ED. We found a decreased expression of genes associated with collagen or extracellular matrix organization and angiogenesis in diabetic fibroblasts, chondrocytes, myofibroblasts, valve-related lymphatic endothelial cells, and pericytes. Moreover, the newly identified pericyte-specific marker, Limb Bud-Heart (Lbh), in mouse and human cavernous tissues, clearly distinguishing pericytes from smooth muscle cells. Cell-cell interaction analysis revealed that pericytes are involved in angiogenesis, adhesion, and migration by communicating with other cell types in the corpus cavernosum; however, these interactions were highly reduced under diabetic conditions. Lbh expression is low in diabetic pericytes, and overexpression of LBH prevents erectile function by regulating neurovascular regeneration. Furthermore, the LBH-interacting proteins (Crystallin Alpha B and Vimentin) were identified in mouse cavernous pericytes through LC-MS/MS analysis, indicating that their interactions were critical for maintaining pericyte function. Thus, our study reveals novel targets and insights into the pathogenesis of ED in patients with diabetes.

Genetic and developmental divergence in the neural crest programme between cichlid fish species

Neural crest (NC) is a vertebrate-specific embryonic progenitor cell population at the basis of important vertebrate features such as the craniofacial skeleton and pigmentation patterns. Despite the wide-ranging variation of NC-derived traits across vertebrates, the contribution of NC to species diversification remains underexplored. Here, leveraging the adaptive diversity of African Great Lakes' cichlid species, we combined comparative transcriptomics and population genomics to investigate the evolution of the NC genetic programme in the context of their morphological divergence. Our analysis revealed substantial differences in transcriptional landscapes across somitogenesis, an embryonic period coinciding with NC development and migration. This included dozens of genes with described functions in the vertebrate NC gene regulatory network, several of which showed signatures of positive selection. Among candidates showing between-species expression divergence, we focused on teleost-specific paralogs of the NC-specifier sox10 (sox10a and sox10b) as prime candidates to influence NC development. These genes, expressed in NC cells, displayed remarkable spatio-temporal variation in cichlids, suggesting their contribution to inter-specific morphological differences. Finally, through CRISPR/Cas9 mutagenesis, we demonstrated the functional divergence between cichlid sox10 paralogs, with the acquisition of a novel skeletogenic function by sox10a. When compared to the teleost models zebrafish and medaka, our findings reveal that sox10 duplication, although retained in most teleost lineages, had variable functional fates across their phylogeny. Altogether, our study suggests that NC-related processes - particularly those controlled by sox10s - might be involved in generating morphological diversification between species and lays the groundwork for further investigations into mechanisms underpinning vertebrate NC diversification.

Blocking LBH expression causes replication stress and sensitizes triple-negative breast cancer cells to ATR inhibitor treatment

Triple-negative (ER-PR-HER2-) breast cancers (TNBC) are highly aggressive and difficult to treat. TNBC exhibit high genomic instability, which enables them to adapt and become resistant to chemo/radiation therapy, leading to rapid disease relapse and mortality. The pro-survival factors that safeguard genome integrity in TNBC cells are poorly understood. LBH is an essential mammary stem cell-specific transcription regulator in the WNT pathway that is aberrantly overexpressed in TNBC, correlating with poor prognosis. Herein, we demonstrate a novel role for LBH in promoting TNBC cell survival. Depletion of LBH in multiple TNBC cell models triggered apoptotic cell death both in vitro and in vivo and led to S-G2M cell cycle delays. Mechanistically, LBH loss causes replication stress due to DNA replication fork stalling, leading to ssDNA breaks, ɣH2AX and 53BP1 nuclear foci formation, and activation of the ATR/CHK1 DNA damage response. Notably, ATR inhibition in combination with LBH downmodulation had a synergistic effect, boosting TNBC cell killing and blocking in vivo tumor growth. Our findings demonstrate, for the first time, that LBH protects the genome integrity of cancer cells by preventing replicative stress. Importantly, they uncover new synthetic lethal vulnerabilities in TNBC that could be exploited for future multi-modal precision medicine.

Epigenetics and the evolution of form: Experimental manipulation of a chromatin modification causes species-specific changes to the craniofacial skeleton

A central question in biology is the molecular origins of phenotypic diversity. While genetic changes are key to the genotype-phenotype relationship, alterations to chromatin structure and the physical packaging of histone proteins may also be important drivers of vertebrate divergence. We investigate the impact of such an epigenetic mechanism, histone acetylation, within a textbook example of an adaptive radiation. Cichlids of Lake Malawi have adapted diverse craniofacial structures, and here we investigate how histone acetylation influences morphological variation in these fishes. Specifically, we assessed the effect of inhibiting histone deacetylation using the drug trichostatin A (TSA) on developing facial structures. We examined this during three critical developmental windows in two cichlid species with alternate adult morphologies. Exposure to TSA during neural crest cell (NCC) migration and as postmigratory NCCs proliferate in the pharyngeal arches resulted in significant changes in lateral and ventral shape in Maylandia, but not in Tropheops. This included an overall shortening of the head, widening of the lower jaw, and steeper craniofacial profile, all of which are paedomorphic morphologies. In contrast, treatment with TSA during early chondrogenesis did not result in significant morphological changes in either species. Together, these data suggest a sensitivity to epigenetic alterations that are both time- and species-dependent. We find that morphologies are due to nonautonomous or potentially indirect effects on NCC development, including in part a global developmental delay. Our research bolsters the understanding that proper histone acetylation is essential for early craniofacial development and identifies a species-specific robustness to developmental change. Overall, this study demonstrates how epigenetic regulation may play an important role in both generating and buffering morphological variation.

Multi-cancer analysis reveals universal association of oncogenic LBH expression with DNA hypomethylation and WNT-Integrin signaling pathways

Limb-Bud and Heart (LBH) is a developmental transcription co-factor deregulated in cancer, with reported oncogenic and tumor suppressive effects. However, LBH expression in most cancer types remains unknown, impeding understanding of its mechanistic function Here, we performed systematic bioinformatic and TMA analysis for LBH in >20 different cancer types. LBH was overexpressed in most cancers compared to normal tissues (>1.5-fold; p < 0.05), including colon-rectal, pancreatic, esophageal, liver, stomach, bladder, kidney, prostate, testicular, brain, head & neck cancers, and sarcoma, correlating with poor prognosis. The cancer types showing LBH downregulation were lung, melanoma, ovarian, cervical, and uterine cancer, while both LBH over- and under-expression were observed in hematopoietic malignancies. In cancers with LBH overexpression, the LBH locus was frequently hypomethylated, identifying DNA hypomethylation as a potential mechanism for LBH dysregulation. Pathway analysis identified a universal, prognostically significant correlation between LBH overexpression and the WNT-Integrin signaling pathways. Validation of the clinical association of LBH with WNT activation in gastrointestinal cancer cell lines, and in colorectal patient samples by IHC uncovered that LBH is specifically expressed in tumor cells with nuclear beta-catenin at the invasive front. Collectively, these data reveal a high degree of LBH dysregulation in cancer and establish LBH as pan-cancer biomarker for detecting WNT hyperactivation in clinical specimens.

Neural crest cells as a source of microevolutionary variation

Vertebrates have some of the most complex and diverse features in animals, from varied craniofacial morphologies to colorful pigmentation patterns and elaborate social behaviors. All of these traits have their developmental origins in a multipotent embryonic lineage of neural crest cells. This "fourth germ layer" is a vertebrate innovation and the source of a wide range of adult cell types. While others have discussed the role of neural crest cells in human disease and animal domestication, less is known about their role in contributing to adaptive changes in wild populations. Here, we review how variation in the development of neural crest cells and their derivatives generates considerable phenotypic diversity in nature. We focus on the broad span of traits under natural and sexual selection whose variation may originate in the neural crest, with emphasis on behavioral factors such as intraspecies communication that are often overlooked. In all, we encourage the integration of evolutionary ecology with developmental biology and molecular genetics to gain a more complete understanding of the role of this single cell type in trait covariation, evolutionary trajectories, and vertebrate diversity.

Mutations in cdon and boc affect trunk neural crest cell migration and slow-twitch muscle development in zebrafish

The transmembrane proteins cdon and boc are implicated in regulating hedgehog signaling during vertebrate development. Recent work showing roles for these genes in axon guidance and neural crest cell migration suggest that cdon and boc may play additional functions in regulating directed cell movements. We use newly generated and existing mutants to investigate a role for cdon and boc in zebrafish neural crest cell migration. We find that single mutant embryos exhibit normal neural crest phenotypes, but that neural crest migration is strikingly disrupted in double cdon;boc mutant embryos. We further show that this migration phenotype is associated with defects in the differentiation of slow-twitch muscle cells, and the loss of a Col1a1a-containing extracellular matrix, suggesting that neural crest defects may be a secondary consequence to defects in mesoderm development. Combined, our data add to a growing literature showing that cdon and boc act synergistically to promote hedgehog signaling during vertebrate development, and suggest that the zebrafish can be used to study the function of hedgehog receptor paralogs.

Zebrafish crocc2 mutants exhibit divergent craniofacial shape, misregulated variability, and aberrant cartilage morphogenesis

BACKGROUND

Phenotypic variation is of paramount importance in development, evolution, and human health; however, the molecular mechanisms that influence organ shape and shape variability are not well understood. During craniofacial development, the behavior of skeletal precursors is regulated by both biochemical and environmental inputs, and the primary cilia play critical roles in transducing both types of signals. Here, we examine a gene that encodes a key constituent of the ciliary rootlets, crocc2, and its role in cartilage morphogenesis in larval zebrafish.

RESULTS

Geometric morphometric analysis of crocc2 mutants revealed altered craniofacial shapes and expanded variation. At the cellular level, we observed altered chondrocyte shapes and planar cell polarity across multiple stages in crocc2 mutants. Notably, cellular defects were specific to areas that experience direct mechanical input. Cartilage cell number, apoptosis, and bone patterning were not affected in crocc2 mutants.

CONCLUSIONS

Whereas "regulatory" genes are widely implicated in patterning the craniofacial skeleton, genes that encode "structural" aspects of the cell are increasingly implicated in shaping the face. Our results add crocc2 to this list, and demonstrate that it affects craniofacial geometry and canalizes phenotypic variation. We propose that it does so via mechanosensing, possibly through the ciliary rootlet. If true, this would implicate a new organelle in skeletal development and evolution.

The transcriptional state and chromatin landscape of cichlid jaw shape variation across species and environments

Adaptive phenotypes are shaped by a combination of genetic and environmental forces, but how they interact remains poorly understood. Here, we utilize the cichlid oral jaw apparatus to better understand these gene-by-environment effects. First, we employed RNA-seq in bony and ligamentous tissues important for jaw opening to identify differentially expressed genes between species and across foraging environments. We used two Lake Malawi species adapted to different foraging habitats along the pelagic-benthic ecomorphological axis. Our foraging treatments were designed to force animals to employ either suction or biting/scraping, which broadly mimic pelagic or benthic modes of feeding. We found a large number of differentially expressed genes between species, and while we identified relatively few differences between environments, species differences were far more pronounced when they were challenged with a pelagic versus benthic foraging mode. Expression data carried the signature of genetic assimilation, and implicated cell cycle regulation in shaping the jaw across species and environments. Next, we repeated the foraging experiment and performed ATAC-seq procedures on nuclei harvested from the same tissues. Cross-referencing results from both analyses revealed subsets of genes that were both differentially expressed and differentially accessible. This reduced dataset implicated notable candidate genes including the Hedgehog effector, KIAA0586 and the ETS transcription factor, etv4, which connects environmental stress and craniofacial morphogenesis. Taken together, these data provide novel insights into the epigenetic, genetic and cellular bases of species- and environment-specific bone shapes.

Morphological and temporal variation in early embryogenesis contributes to species divergence in Malawi cichlid fishes

The cichlid fishes comprise the largest extant vertebrate family and are the quintessential example of rapid "explosive" adaptive radiations and phenotypic diversification. Despite low genetic divergence, East African cichlids harbor a spectacular intra- and interspecific morphological diversity, including the hyper-variable, neural crest (NC)-derived traits such as coloration and craniofacial skeleton. Although the genetic and developmental basis of these phenotypes has been investigated, understanding of when, and specifically how early, in ontogeny species-specific differences emerge, remains limited. Since adult traits often originate during embryonic development, the processes of embryogenesis could serve as a potential source of species-specific variation. Consequently, we designed a staging system by which we compare the features of embryogenesis between three Malawi cichlid species-Astatotilapia calliptera, Tropheops sp. 'mauve' and Rhamphochromis sp. "chilingali"-representing a wide spectrum of variation in pigmentation and craniofacial morphologies. Our results showed fundamental differences in multiple aspects of embryogenesis that could underlie interspecific divergence in adult adaptive traits. First, we identified variation in the somite number and signatures of temporal variation, or heterochrony, in the rates of somite formation. The heterochrony was also evident within and between species throughout ontogeny, up to the juvenile stages. Finally, the identified interspecific differences in the development of pigmentation and craniofacial cartilages, present at the earliest stages of their overt formation, provide compelling evidence that the species-specific trajectories begin divergence during early embryogenesis, potentially during somitogenesis and NC development. Altogether, our results expand our understanding of fundamental cichlid biology and provide new insights into the developmental origins of vertebrate morphological diversity.

Non-synonymous variation and protein structure of candidate genes associated with selection in farm and wild populations of turbot (Scophthalmus maximus)

Non-synonymous variation (NSV) of protein coding genes represents raw material for selection to improve adaptation to the diverse environmental scenarios in wild and livestock populations. Many aquatic species face variations in temperature, salinity and biological factors throughout their distribution range that is reflected by the presence of allelic clines or local adaptation. The turbot (Scophthalmus maximus) is a flatfish of great commercial value with a flourishing aquaculture which has promoted the development of genomic resources. In this study, we developed the first atlas of NSVs in the turbot genome by resequencing 10 individuals from Northeast Atlantic Ocean. More than 50,000 NSVs where detected in the ~ 21,500 coding genes of the turbot genome, and we selected 18 NSVs to be genotyped using a single Mass ARRAY multiplex on 13 wild populations and three turbot farms. We detected signals of divergent selection on several genes related to growth, circadian rhythms, osmoregulation and oxygen binding in the different scenarios evaluated. Furthermore, we explored the impact of NSVs identified on the 3D structure and functional relationship of the correspondent proteins. In summary, our study provides a strategy to identify NSVs in species with consistently annotated and assembled genomes to ascertain their role in adaptation.

East African cichlid fishes

Cichlid fishes are a very diverse and species-rich family of teleost fishes that inhabit lakes and rivers of India, Africa, and South and Central America. Research has largely focused on East African cichlids of the Rift Lakes Tanganyika, Malawi, and Victoria that constitute the biodiversity hotspots of cichlid fishes. Here, we give an overview of the study system, research questions, and methodologies. Research on cichlid fishes spans many disciplines including ecology, evolution, physiology, genetics, development, and behavioral biology. In this review, we focus on a range of organismal traits, including coloration phenotypes, trophic adaptations, appendages like fins and scales, sensory systems, sex, brains, and behaviors. Moreover, we discuss studies on cichlid phylogenies, plasticity, and general evolutionary patterns, ranging from convergence to speciation rates and the proximate and ultimate mechanisms underlying these processes. From a methodological viewpoint, the last decade has brought great advances in cichlid fish research, particularly through the advent of affordable deep sequencing and advances in genetic manipulations. The ability to integrate across traits and research disciplines, ranging from developmental biology to ecology and evolution, makes cichlid fishes a fascinating research system.

Expression atlas of avian neural crest proteins: Neurulation to migration

Neural crest (NC) cells are a dynamic population of embryonic stem cells that create various adult tissues in vertebrate species including craniofacial bone and cartilage and the peripheral and enteric nervous systems. NC development is thought to be a conserved and complex process that is controlled by a tightly-regulated gene regulatory network (GRN) of morphogens, transcription factors, and cell adhesion proteins. While multiple studies have characterized the expression of several GRN factors in single species, a comprehensive protein analysis that directly compares expression across development is lacking. To address this lack in information, we used three closely related avian models, Gallus gallus (chicken), Coturnix japonica (Japanese quail), and Pavo cristatus (Indian peafowl), to compare the localization and timing of four GRN transcription factors, PAX7, SNAI2, SOX9, and SOX10, from the onset of neurulation to migration. While the spatial expression of these factors is largely conserved, we find that quail NC cells express SNAI2, SOX9, and SOX10 proteins at the equivalent of earlier developmental stages than chick and peafowl. In addition, quail NC cells migrate farther and more rapidly than the larger organisms. These data suggest that despite a conservation of NC GRN players, differences in the timing of NC development between species remain a significant frontier to be explored with functional studies.

Limb-bud and heart development (LBH) contributes to glioma progression in vitro and in vivo

Glioma is the predominant brain malignancy and is correlated with high mortality and severe morbidity. The transcription factor limb‐bud and heart (LBH) has been reported to be involved in the development of several cancers, although its role in glioma development remains elusive. Here, we examined the effect of LBH on glioma progression. The expression of LBH was increased in glioma samples from The Cancer Genome Atlas database, and upregulation of LBH was observed to be correlated with the poor survival of glioma patients. We also report that expression of LBH was elevated in clinical glioma tissues compared to adjacent normal tissues, and was also enhanced in glioma cell lines. LBH promotes proliferation and inhibits cell cycle arrest and apoptosis in glioma cells. In addition, LBH increased the migration and invasion of glioma cells in vitro. Moreover, tumorigenicity analysis revealed that LBH could promote the tumor growth of glioma cells in vivo. In conclusion, our findings suggest that LBH contributes to glioma progression in vitro and in vivo. Our findings provide new insights into the mechanism by which LBH promotes the development of glioma, improving our understanding of the correlation between LBH with cancer. LBH may have potential as a target for glioma therapy.

Complex genetic architecture of three-dimensional craniofacial shape variation in domestic pigeons

Deciphering the genetic basis of vertebrate craniofacial variation is a longstanding biological problem with broad implications in evolution, development, and human pathology. One of the most stunning examples of craniofacial diversification is the adaptive radiation of birds, in which the beak serves essential roles in virtually every aspect of their life histories. The domestic pigeon (Columba livia) provides an exceptional opportunity to study the genetic underpinnings of craniofacial variation because of its unique balance of experimental accessibility and extraordinary phenotypic diversity within a single species. We used traditional and geometric morphometrics to quantify craniofacial variation in an F2 laboratory cross derived from the straight-beaked Pomeranian Pouter and curved-beak Scandaroon pigeon breeds. Using a combination of genome-wide quantitative trait locus scans and multi-locus modeling, we identified a set of genetic loci associated with complex shape variation in the craniofacial skeleton, including beak shape, braincase shape, and mandible shape. Some of these loci control coordinated changes between different structures, while others explain variation in the size and shape of specific skull and jaw regions. We find that in domestic pigeons, a complex blend of both independent and coupled genetic effects underlie three-dimensional craniofacial morphology.

Differences in Cell Proliferation and Craniofacial Phenotype of Closely Related Species in the Pupfish Genus Cyprinodon

Understanding the genetic basis for phenotypic differences is fundamental to the study of macroevolutionary patterns of biological diversity. While technological advances in DNA sequencing have made researching genetic variation in wild taxa routine, fully understanding how these variants affect phenotype requires taking the next step to investigate how genetic changes alter cell and tissue interactions that ultimately produce phenotypes. In this article, we investigate a role for cell proliferation as a developmental source of craniofacial diversity in a radiation of 3 species of Cyprinodon from San Salvador Island, Bahamas. Patterns of cell proliferation in the heads of hatching-age fish differ among species of Cyprinodon, and correlate with differences in allometric growth rate among the jaws of 3 distinct species. Regional patterns of cell proliferation in the head are complex, resulting in an unintuitive result in which lower levels of cell proliferation in the posterior head region are associated with the development of relatively larger jaws in one species. We combine these data with previously published morphological and genomic data to show how studying the mechanisms generating phenotype at the cellular and tissue levels of biological organization can help mechanistically link genomic studies with classic morphological studies.

Microcomputed tomography used to link head morphology and observed feeding behavior in cichlids of Lake Malaŵi

Cichlids inhabiting the Great Lakes of Africa have radiated extremely rapidly, with Lake Malaŵi harboring some 850 species. This rapid radiation may be linked to the diversity in behaviors, sexual selection, and phenotypic plasticity. To determine the relationships between morphology and behaviors, microcomputed tomography (microCT) was used to observe internal morphological structures. Observed morphological adaptations were linked with observed behavior of cichlids in Lake Malaŵi with respect to the various available food resources. Many of these adaptations have parallels, sometimes into the finest details, in other drainage systems and can thus be considered as variations of how cichlids in general respond to environmental opportunities and challenges. Variations in the structure and teeth of the pharyngeal jaws and the oral jaws allowed for fine tuning of specializations, so that various species can utilize the same source without direct competition. We suggested that high-resolution X-ray computed tomography will permit scientists to infer life history and behavior characters of rare or extinct taxa from a detailed examination of morphology and linkages between morphology and behavior found in extant species.

Ciliary Rootlet Coiled-Coil 2 (crocc2) Is Associated with Evolutionary Divergence and Plasticity of Cichlid Jaw Shape

Cichlid fishes exhibit rapid, extensive, and replicative adaptive radiation in feeding morphology. Plasticity of the cichlid jaw has also been well documented, and this combination of iterative evolution and developmental plasticity has led to the proposition that the cichlid feeding apparatus represents a morphological "flexible stem." Under this scenario, the fixation of environmentally sensitive genetic variation drives evolutionary divergence along a phenotypic axis established by the initial plastic response. Thus, if plasticity is predictable then so too should be the evolutionary response. We set out to explore these ideas at the molecular level by identifying genes that underlie both the evolution and plasticity of the cichlid jaw. As a first step, we fine-mapped an environment-specific quantitative trait loci for lower jaw shape in cichlids, and identified a nonsynonymous mutation in the ciliary rootlet coiled-coil 2 (crocc2), which encodes a major structural component of the primary cilium. Given that primary cilia play key roles in skeletal mechanosensing, we reasoned that this gene may confer its effects by regulating the sensitivity of bone to respond to mechanical input. Using both cichlids and zebrafish, we confirmed this prediction through a series of experiments targeting multiple levels of biological organization. Taken together, our results implicate crocc2 as a novel mediator of bone formation, plasticity, and evolution.

Microevolutionary change in viscerocranial bones under congeneric sympatry in the Lake Tanganyikan cichlid genus Tropheus

The endemic Lake Tanganyika cichlid genus Tropheus lives at rocky shores all around the lake and comprises six species which are subdivided into about 120 morphologically similar but color-wise distinct populations. Typically, they live without a second Tropheus species, but there are some regions where two or even three sister species live in sympatry. We previously showed that there are morphological differences concerning head shape, eye size and insertion of fins among populations living alone compared to those living in sympatry with a second Tropheus. This study goes one step further to test if sympatry affects the shape of viscerocranial bones. By means of geometric morphometrics, we compare the shape of four bones among thirteen Tropheus populations, some of which in sympatry and some living alone. We quantify patterns of shape variation and estimate morphological disparity among the four bony elements in the study species and populations. We found consistent differences in the shape of one bony element among non-sympatric and sympatric populations, besides an extensive variation in the shape of viscerocranial bones within and among species. Furthermore, sexual dimorphism in Tropheus is clearly evident in the viscerocranial bones analyzed. We suggest that the relatively subtle morphological signal in sympatric vs. non-sympatric Tropheus populations is owed to the fact that the depth segregation does not yet represent a full shift in the trophic niche, albeit our data confirm that differences in ecologically relevant traits, such as bones of the preorbital region, play an important role in the process of niche separation and in the context of explosive diversification of cichlid fishes.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10750-021-04536-7.

Trim-Away mediated knock down uncovers a new function for Lbh during gastrulation of Xenopus laevis

We previously identified the protein Lbh as necessary for cranial neural crest (CNC) cell migration in Xenopus through the use of morpholinos. However, Lbh is a maternally deposited protein and morpholinos achieve knockdowns through prevention of translation. In order to investigate the role of Lbh in earlier embryonic events, we employed the new technique "Trim-Away" to degrade this maternally deposited protein. Trim-Away utilizes the E3 ubiquitin ligase trim21 to degrade proteins targeted with an antibody and was developed in mammalian systems. Our results show that Xenopus is amenable to the Trim-Away technique. We also show that early knockdown of Lbh in Xenopus results in defects in gastrulation that present with a decrease in fibronectin matrix assembly, an increased in mesodermal cell migration and decrease in endodermal cell cohesion. We further show that the technique is also effective on a second abundant maternal protein PACSIN2. We discuss potential advantages and limit of the technique in Xenopus embryos as well as the mechanism of gastrulation inhibition.

Transcriptome sequencing reveals signatures of positive selection in the Spot-Tailed Earless Lizard

The continual loss of threatened biodiversity is occurring at an accelerated pace. High-throughput sequencing technologies are now providing opportunities to address this issue by aiding in the generation of molecular data for many understudied species of high conservation interest. Our overall goal of this study was to begin building the genomic resources to continue investigations and conservation of the Spot-Tailed Earless lizard. Here we leverage the power of high-throughput sequencing to generate the liver transcriptome for the Northern Spot-Tailed Earless Lizard (Holbrookia lacerata) and Southern Spot-Tailed Earless Lizard (Holbrookia subcaudalis), which have declined in abundance in the past decades, and their sister species, the Common Lesser Earless Lizard (Holbrookia maculata). Our efforts produced high quality and robust transcriptome assemblies validated by 1) quantifying the number of processed reads represented in the transcriptome assembly and 2) quantifying the number of highly conserved single-copy orthologs that are present in our transcript set using the BUSCO pipeline. We found 1,361 1-to-1 orthologs among the three Holbrookia species, Anolis carolinensis, and Sceloporus undulatus. We carried out dN/dS selection tests using a branch-sites model and identified a dozen genes that experienced positive selection in the Holbrookia lineage with functions in development, immunity, and metabolism. Our single-copy orthologous sequences additionally revealed significant pairwise sequence divergence (~.73%) between the Northern H. lacerata and Southern H. subcaudalis that further supports the recent elevation of the Southern Spot-Tailed Earless Lizard to full species.

CHD9 upregulates RUNX2 and has a potential role in skeletal evolution

BACKGROUND

Changes in gene regulation are widely recognized as an important driver of adaptive phenotypic evolution. However, the specific molecular mechanisms that underpin such changes are still poorly understood. Chromatin state plays an essential role in gene regulation, by influencing the accessibility of coding loci to the transcriptional machinery. Changes in the function of chromatin remodellers are therefore strong candidates to drive changes in gene expression associated with phenotypic adaptation. Here, we identify amino acid homoplasies in the chromatin remodeller CHD9, shared between the extinct marsupial thylacine and eutherian wolf which show remarkable skull convergence. CHD9 is involved in osteogenesis, though its role in the process is still poorly understood. We examine whether CHD9 is able to regulate the expression of osteogenic target genes and examine the function of a key substitution in the CHD9 DNA binding domain.

RESULTS

We examined whether CHD9 was able to upregulate its osteogenic target genes, RUNX2, Osteocalcin (OC) and ALP in HEK293T cells. We found that overexpression of CHD9 upregulated RUNX2, the master regulator of osteoblast cell fate, but not the downstream genes OC or ALP, supporting the idea that CHD9 regulates osteogenic progenitors rather than terminal osteoblasts. We also found that the evolutionary substitution in the CHD9 DNA binding domain does not alter protein secondary structure, but was able to drive a small but insignificant increase in RUNX2 activation. Finally, CHD9 was unable to activate an episomal RUNX2 promoter-reporter construct, suggesting that CHD9 requires the full chromatin complement for its function.

CONCLUSIONS

We provide new evidence to the role of CHD9 in osteogenic differentiation through its newly observed ability to upregulate the expression of RUNX2. Though we were unable to identify significant functional consequences of the evolutionary substitution in HEK293T cells, our study provides important steps forward in the functional investigation of protein homoplasy and its role in developmental processes. Mutations in coding genes may be a mechanism for driving adaptive changes in gene expression, and their validation is essential towards determining the functional consequences of evolutionary homoplasy.

Overexpression of Limb-Bud and Heart (LBH) promotes angiogenesis in human glioma via VEGFA-mediated ERK signalling under hypoxia

Background

Glioma is the most common primary malignant tumor in the central nervous system with frequent hypoxia and angiogenesis. Limb-Bud and Heart (LBH) is a highly conserved transcription cofactor that participates in embryonic development and tumorigenesis.

Methods

The conditioned media from LBH regulated human glioma cell lines and patient-derived glioma stem cells (GSCs) were used to treat the human brain microvessel endothelial cells (hBMECs). The function of LBH on angiogenesis were examined through methods of MTS assay, Edu assay, TUNEL assay, western blotting analysis, qPCR analysis, luciferase reporter assay and xenograft experiment.

Findings

Our study found for the first time that LBH was overexpressed in gliomas and was associated with a poor prognosis. LBH overexpression participated in the angiogenesis of gliomas via the vascular endothelial growth factor A (VEGFA)-mediated extracellular signal-regulated kinase (ERK) signalling pathway in human brain microvessel endothelial cells (hBMECs). Rapid proliferation of gliomas can lead to tissue hypoxia and hypoxia inducible factor-1 (HIF-1) activation, while HIF-1 can directly transcriptionally regulate the expression of LBH and result in a self-reinforcing cycle.

Interpretation

LBH may be a possible treatment target to break the vicious cycle in glioma treatment.

 

 

Divergence in larval jaw gene expression reflects differential trophic adaptation in haplochromine cichlids prior to foraging

Background

Understanding how variation in gene expression contributes to morphological diversity is a major goal in evolutionary biology. Cichlid fishes from the East African Great lakes exhibit striking diversity in trophic adaptations predicated on the functional modularity of their two sets of jaws (oral and pharyngeal). However, the transcriptional basis of this modularity is not so well understood, as no studies thus far have directly compared the expression of genes in the oral and pharyngeal jaws. Nor is it well understood how gene expression may have contributed to the parallel evolution of trophic morphologies across the replicate cichlid adaptive radiations in Lake Tanganyika, Malawi and Victoria.

Results

We set out to investigate the role of gene expression divergence in cichlid fishes from these three lakes adapted to herbivorous and carnivorous trophic niches. We focused on the development stage prior to the onset of exogenous feeding that is critical for understanding patterns of gene expression after oral and pharyngeal jaw skeletogenesis, anticipating environmental cues. This framework permitted us for the first time to test for signatures of gene expression underlying jaw modularity in convergent eco-morphologies across three independent adaptive radiations. We validated a set of reference genes, with stable expression between the two jaw types and across species, which can be important for future studies of gene expression in cichlid jaws. Next we found evidence of modular and non-modular gene expression between the two jaws, across different trophic niches and lakes. For instance, prdm1a, a skeletogenic gene with modular anterior-posterior expression, displayed higher pharyngeal jaw expression and modular expression pattern only in carnivorous species. Furthermore, we found the expression of genes in cichlids jaws from the youngest Lake Victoria to exhibit low modularity compared to the older lakes.

Conclusion

Overall, our results provide cross-species transcriptional comparisons of modularly-regulated skeletogenic genes in the two jaw types, implicating expression differences which might contribute to the formation of divergent trophic morphologies at the stage of larval independence prior to foraging.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1483-3) contains supplementary material, which is available to authorized users.

Elevated limb-bud and heart development (LBH) expression indicates poor prognosis and promotes gastric cancer cell proliferation and invasion via upregulating Integrin/FAK/Akt pathway

The limb-bud and heart development (LBH) gene is a highly conserved, tissue-specific transcription cofactor in vertebrates that regulates multiple key genes in embryonic development. The role of LBH in various cancer types is still controversial, and its specific role and molecular mechanism in the oncogenesis of gastric cancer (GC) remains largely unexplored. In the present study, the prognostic significance and clinicopathological characteristics of LBH in GC was determined. The LBH mRNA expression was first investigated in four independent public datasets (TCGA-STAD, GSE15459, GSE29272, and GSE62254) and then validated with our samples at the protein level. LBH was overexpressed at both the mRNA and protein levels in cancer compared with normal tissues. High LBH expression was correlated with advanced T, N, and M stages. Kaplan–Meier analysis and log-rank test indicated that higher LBH expression was statistically correlated with shorter overall survival (OS) in the public datasets and our study samples. Univariate and multivariate Cox regression analysis showed that LBH was an independent prognostic biomarker for survival in TCGA-STAD, GSE15459, GSE62254 cohorts, and our GC patients. In vitro experiments showed that knockdown of LBH can significantly inhibit the proliferation and invasion of HGC-27 cells, while overexpression of LBH can significantly enhance the proliferation and invasion of BGC-823 cells. Gene Set Enrichment Analysis (GSEA), Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomics (KEGG) indicated that high LBH expression is associated with the PI3K-Akt pathway, focal adhesion, and extracellular matrix (ECM)-receptor interaction. Western blot analysis showed that knockdown of LBH significantly inhibited the expression of integrin α5, integrin β1, p-FAK, and p-Akt. Therefore, results from the present study indicate that LBH is a potential independent prognostic biomarker and promotes proliferation and invasion of GC cells by activating the integrin/FAK/Akt pathway.

Genetic architecture and sex-specific selection govern modular, male-biased evolution of doublesex

doublesex regulates early embryonic sex differentiation in holometabolous insects, along with the development of species-, sex-, and morph-specific adaptations during pupal stages. How does a highly conserved gene with a critical developmental role also remain functionally dynamic enough to gain ecologically important adaptations that are divergent in sister species? We analyzed patterns of exon-level molecular evolution and protein structural homology of doublesex from 145 species of four insect orders representing 350 million years of divergence. This analysis revealed that evolution of doublesex was governed by a modular architecture: Functional domains and female-specific regions were highly conserved, whereas male-specific sequences and protein structures evolved up to thousand-fold faster, with sites under pervasive and/or episodic positive selection. This pattern of sex bias was reversed in Hymenoptera. Thus, highly conserved yet dynamic master regulators such as doublesex may partition specific conserved and novel functions in different genic modules at deep evolutionary time scales.

A Natural Animal Model System of Craniofacial Anomalies: The Blind Mexican Cavefish

Natural model systems evolving under extreme environmental pressures provide the opportunity to advance our knowledge of how the craniofacial complex evolves in nature. Unlike traditional models, natural systems are less inbred, and, therefore, better model the complex variation of the human population. Owing to the nature of certain craniofacial aberrations in blind Mexican cavefish, we suggest that this organism can provide new insights to a variety of craniofacial changes. Diverse cranial features have evolved in natural cave-dwelling Astyanax fish, which have thrived in the unforgiving darkness and nutrient-poor environment of the cave for countless generations. While the genetic and environmental underpinnings of various cranial anomalies have been investigated for decades, a comprehensive characterization of their molecular and developmental origins remains incomplete. Cavefish provide numerous advantages given the availability of genomic resources, developmental and molecular tools, and the presence of a normative surface-dwelling "ancestral" surrogate for comparative studies. By leveraging the frequency of abnormal and asymmetric cranial features in cavefish, we anticipate advances in our knowledge of the etiologies of irregular cranial features. Extreme adaptations in cavefish are expected to offer new insights into the complex and multifactorial nature of craniofacial disorders and facial asymmetry. Anat Rec, 2018. © 2018 American Association for Anatomy.

Genetic Dissection of a Supergene Implicates Tfap2a in Craniofacial Evolution of Threespine Sticklebacks

In nature, multiple adaptive phenotypes often coevolve and can be controlled by tightly linked genetic loci known as supergenes. Dissecting the genetic basis of these linked phenotypes is a major challenge in evolutionary genetics. Multiple freshwater populations of threespine stickleback fish (Gasterosteus aculeatus) have convergently evolved two constructive craniofacial traits, longer branchial bones and increased pharyngeal tooth number, likely as adaptations to dietary differences between marine and freshwater environments. Prior QTL mapping showed that both traits are partially controlled by overlapping genomic regions on chromosome 21 and that a regulatory change in Bmp6 likely underlies the tooth number QTL. Here, we mapped the branchial bone length QTL to a 155 kb, eight-gene interval tightly linked to, but excluding the coding regions of Bmp6 and containing the candidate gene Tfap2a Further recombinant mapping revealed this bone length QTL is separable into at least two loci. During embryonic and larval development, Tfap2a was expressed in the branchial bone primordia, where allele specific expression assays revealed the freshwater allele of Tfap2a was expressed at lower levels relative to the marine allele in hybrid fish. Induced loss-of-function mutations in Tfap2a revealed an essential role in stickleback craniofacial development and show that bone length is sensitive to Tfap2a dosage in heterozygotes. Combined, these results suggest that closely linked but genetically separable changes in Bmp6 and Tfap2a contribute to a supergene underlying evolved skeletal gain in multiple freshwater stickleback populations.

The Role of Alternative Splicing and Differential Gene Expression in Cichlid Adaptive Radiation

Species diverge eco-morphologically through the continuous action of natural selection on functionally important structures, producing alternative adaptive morphologies. In cichlid fishes, the oral and pharyngeal jaws are such key structures. Adaptive variation in jaw morphology contributes to trophic specialization, which is hypothesized to fuel their rapid speciation in the East African Great Lakes. Much is known about the genes involved in cichlid jaw and craniofacial development. However, it is still unclear what salient sources of variation gave rise to trophic-niche specialization, facilitating adaptive radiation. Here, we explore two sources of transcriptional variation that may underlie species-specific disparities in jaw morphology. Using whole transcriptome RNA-sequencing, we analyze differences in gene expression and alternative splicing, at the end of postlarval development, in fully functional jaws of six species of cichlids from the Lake Tanganyika tribe Tropheini. Our data reveal a surprisingly high degree of alternative splicing events compared with gene expression differences among species and trophic types. This suggests that differential trophic adaptation of the jaw apparatus may have been shaped by transcriptional rewiring of splicing as well as gene expression variation during the rapid radiation of the Tropheini. Specifically, genes undergoing splicing across most species were found to be enriched for pharyngeal jaw gene ontology terms. Overall, jaw transcriptional patterns at postlarval developmental stage were highly dynamic and species-specific. In conclusion, this work indicates that shifts in alternative splicing could have played a more important role in cichlid adaptive radiation, and possibly adaptive radiation in general, than currently recognized.

Effect of craniofacial genotype on the relationship between morphology and feeding performance in cichlid fishes

The relationship between morphology and performance is complex, but important for understanding the adaptive nature of morphological variation. Recent studies have sought to better understand this system by illuminating the interconnectedness of different functional systems; however, the role of genetics is often overlooked. In this study, we attempt to gain insights into this relationship by examining the effect of genotypic variation at putative craniofacial loci on the relationship between morphology and feeding performance in cichlids. We studied two morphologically disparate species, as well as a morphologically intermediate hybrid population. We assessed feeding performance, jaw protrusion, and general facial morphology for each fish. We also genotyped hybrid animals at six previously identified craniofacial loci. Cichlid species were found to differ in facial geometry, kinematic morphology, and performance. Significant correlations were also noted between these variables; however, the explanatory power of facial geometry in predicting performance was relatively poor. Notably, when hybrids were grouped by genotype, the relationship between shape and performance improved. This relationship was especially robust in animals with the specialist allele at sox9b, a well-characterized regulator of craniofacial development. These data suggest a novel role for genotype in influencing complex relationships between form and function.

The Cyprinodon variegatus genome reveals gene expression changes underlying differences in skull morphology among closely related species

BACKGROUND

Understanding the genetic and developmental origins of phenotypic novelty is central to the study of biological diversity. In this study we identify modifications to the expression of genes at four developmental stages that may underlie jaw morphological differences among three closely related species of pupfish (genus Cyprinodon) from San Salvador Island, Bahamas. Pupfishes on San Salvador Island are trophically differentiated and include two endemic species that have evolved jaw morphologies unlike that of any other species in the genus Cyprinodon.

RESULTS

We find that gene expression differs significantly across recently diverged species of pupfish. Genes such as Bmp4 and calmodulin, previously implicated in jaw diversification in African cichlid fishes and Galapagos finches, were not found to be differentially expressed among species of pupfish. Instead we find multiple growth factors and cytokine/chemokine genes to be differentially expressed among these pupfish taxa. These include both genes and pathways known to affect craniofacial development, such as Wnt signaling, as well as novel genes and pathways not previously implicated in craniofacial development. These data highlight both shared and potentially unique sources of jaw diversity in pupfish and those identified in other evolutionary model systems such as Galapagos finches and African cichlids.

CONCLUSIONS

We identify modifications to the expression of genes involved in Wnt signaling, Igf signaling, and the inflammation response as promising avenues for future research. Our project provides insight into the magnitude of gene expression changes contributing to the evolution of morphological novelties, such as jaw structure, in recently diverged pupfish species.

Signalling pathways in trophic skeletal development and morphogenesis: Insights from studies on teleost fish

During the development of the vertebrate feeding apparatus, a variety of complicated cellular and molecular processes participate in the formation and integration of individual skeletal elements. The molecular mechanisms regulating the formation of skeletal primordia and their development into specific morphological structures are tightly controlled by a set of interconnected signalling pathways. Some of these pathways, such as Bmp, Hedgehog, Notch and Wnt, are long known for their pivotal roles in craniofacial skeletogenesis. Studies addressing the functional details of their components and downstream targets, the mechanisms of their interactions with other signals as well as their potential roles in adaptive morphological divergence, are currently attracting considerable attention. An increasing number of signalling pathways that had previously been described in different biological contexts have been shown to be important in the regulation of jaw skeletal development and morphogenesis. In this review, I provide an overview of signalling pathways involved in trophic skeletogenesis emphasizing studies of the most species-rich group of vertebrates, the teleost fish, which through their evolutionary history have undergone repeated episodes of spectacular trophic diversification.

Cichlid fishes as a model to understand normal and clinical craniofacial variation

We have made great strides towards understanding the etiology of craniofacial disorders, especially for 'simple' Mendelian traits. However, the facial skeleton is a complex trait, and the full spectrum of genetic, developmental, and environmental factors that contribute to its final geometry remain unresolved. Forward genetic screens are constrained with respect to complex traits due to the types of genes and alleles commonly identified, developmental pleiotropy, and limited information about the impact of environmental interactions. Here, we discuss how studies in an evolutionary model - African cichlid fishes - can complement traditional approaches to understand the genetic and developmental origins of complex shape. Cichlids exhibit an unparalleled range of natural craniofacial morphologies that model normal human variation, and in certain instances mimic human facial dysmorphologies. Moreover, the evolutionary history and genomic architecture of cichlids make them an ideal system to identify the genetic basis of these phenotypes via quantitative trait loci (QTL) mapping and population genomics. Given the molecular conservation of developmental genes and pathways, insights from cichlids are applicable to human facial variation and disease. We review recent work in this system, which has identified lbh as a novel regulator of neural crest cell migration, determined the Wnt and Hedgehog pathways mediate species-specific bone morphologies, and examined how plastic responses to diet modulate adult facial shapes. These studies have not only revealed new roles for existing pathways in craniofacial development, but have identified new genes and mechanisms involved in shaping the craniofacial skeleton. In all, we suggest that combining work in traditional laboratory and evolutionary models offers significant potential to provide a more complete and comprehensive picture of the myriad factors that are involved in the development of complex traits.

Quantitative trait locus analysis of body shape divergence in nine-spined sticklebacks based on high-density SNP-panel

Heritable phenotypic differences between populations, caused by the selective effects of distinct environmental conditions, are of commonplace occurrence in nature. However, the actual genomic targets of this kind of selection are still poorly understood. We conducted a quantitative trait locus (QTL) mapping study to identify genomic regions responsible for morphometric differentiation between genetically and phenotypically divergent marine and freshwater nine-spined stickleback (Pungitius pungitius) populations. Using a dense panel of SNP-markers obtained by restriction site associated DNA sequencing of an F2 recombinant cross, we found 22 QTL that explained 3.5-12.9% of phenotypic variance in the traits under investigation. We detected one fairly large-effect (PVE = 9.6%) QTL for caudal peduncle length-a trait with a well-established adaptive function showing clear differentiation among marine and freshwater populations. We also identified two large-effect QTL for lateral plate numbers, which are different from the lateral plate QTL reported in earlier studies of this and related species. Hence, apart from identifying several large-effect QTL in shape traits showing adaptive differentiation in response to different environmental conditions, the results suggest intra- and interspecific heterogeneity in the genomic basis of lateral plate number variation.

Developmental basis of phenotypic integration in two Lake Malawi cichlids

Background

Cichlid fishes from the Rift Lakes of East Africa have undergone the most spectacular adaptive radiations in vertebrate history. Eco-morphological adaptations in lakes Victoria, Malawi and Tanganyika have resulted in a vast array of skull shapes and sizes, yet primary axes of morphological variation are conserved in all three radiations, prominently including the size of the preorbital region of the skull. This conserved pattern suggests that development may constrain the trajectories of cichlid head morphological evolution.

Results

Here, we (1) present a comparative analysis of adult head morphology in two sand-dweller cichlids from Lake Malawi with preorbital size differences representative of the main axis of variation among the three lakes and (2) analyze the ontogeny of shape and size differences by focusing on known developmental modules throughout the head. We find that (1) developmental differences between the two species correlate with known developmental modules; (2) differences in embryonic cartilage development result in phenotypically integrated changes among all bones derived from a single cartilage, while differences in dermal bone development tend to influence isolated regions within a bone; and lastly (3) species-specific morphologies appear in the embryo as subtle differences, which become progressively amplified throughout ontogeny. We propose that this amplification takes place at skeletal growth zones, the locations and shapes of which are patterned during embryogenesis.

Conclusions

This study is the most anatomically comprehensive analysis of the developmental differences underlying cichlid skull evolution in the Rift Lakes of East Africa. The scale of our analysis reveals previously unnoticed correlations between developmental modules and patterns of phenotypic integration. We propose that the primary axes of morphological variation among East African cichlid adaptive radiations are constrained by the hierarchical modularity of the teleost head skeleton.

Molecular Evolution of the Neural Crest Regulatory Network in Ray-Finned Fish

Gene regulatory networks (GRN) are central to developmental processes. They are composed of transcription factors and signaling molecules orchestrating gene expression modules that tightly regulate the development of organisms. The neural crest (NC) is a multipotent cell population that is considered a key innovation of vertebrates. Its derivatives contribute to shaping the astounding morphological diversity of jaws, teeth, head skeleton, or pigmentation. Here, we study the molecular evolution of the NC GRN by analyzing patterns of molecular divergence for a total of 36 genes in 16 species of bony fishes. Analyses of nonsynonymous to synonymous substitution rate ratios (dN/dS) support patterns of variable selective pressures among genes deployed at different stages of NC development, consistent with the developmental hourglass model. Model-based clustering techniques of sequence features support the notion of extreme conservation of NC-genes across the entire network. Our data show that most genes are under strong purifying selection that is maintained throughout ray-finned fish evolution. Late NC development genes reveal a pattern of increased constraints in more recent lineages. Additionally, seven of the NC-genes showed signs of relaxation of purifying selection in the famously species-rich lineage of cichlid fishes. This suggests that NC genes might have played a role in the adaptive radiation of cichlids by granting flexibility in the development of NC-derived traits-suggesting an important role for NC network architecture during the diversification in vertebrates.

Differential expression of the aryl hydrocarbon receptor pathway associates with craniofacial polymorphism in sympatric Arctic charr

Background

The developmental basis of craniofacial morphology hinges on interactions of numerous signalling systems. Extensive craniofacial variation in the polymorphic Arctic charr, a member of the salmonid family, from Lake Thingvallavatn (Iceland), offers opportunities to find and study such signalling pathways and their key regulators, thereby shedding light on the developmental pathways, and the genetics of trophic divergence.

Results

To identify genes involved in the craniofacial differences between benthic and limnetic Arctic charr, we used transcriptome data from different morphs, spanning early development, together with data on craniofacial expression patterns and skeletogenesis in model vertebrate species. Out of 20 genes identified, 7 showed lower gene expression in benthic than in limnetic charr morphs. We had previously identified a conserved gene network involved in extracellular matrix (ECM) organization and skeletogenesis, showing higher expression in developing craniofacial elements of benthic than in limnetic Arctic charr morphs. The present study adds a second set of genes constituting an expanded gene network with strong, benthic–limnetic differential expression. To identify putative upstream regulators, we performed knowledge-based motif enrichment analyses on the regulatory sequences of the identified genes which yielded potential binding sites for a set of known transcription factors (TFs). Of the 8 TFs that we examined using qPCR, two (Ahr2b and Ap2) were found to be differentially expressed between benthic and limnetic charr. Expression analysis of several known AhR targets indicated higher activity of the AhR pathway during craniofacial development in benthic charr morphotypes.

Conclusion

These results suggest a key role of the aryl hydrocarbon receptor (AhR) pathway in the observed craniofacial differences between distinct charr morphotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0022-6) contains supplementary material, which is available to authorized users.

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 > 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.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0020-8) contains supplementary material, which is available to authorized users.

Zebrafish Lbh-like Is Required for Otx2-mediated Photoreceptor Differentiation

The homeobox transcription factor orthodenticle homolog 2 (otx2) is supposed as an organizer that orchestrates a transcription factor network during photoreceptor development. However, its regulation in the process remains unclear. In this study, we have identified a zebrafish limb bud and heart-like gene (lbh-like), which is expressed initially at 30 hours post fertilization (hpf) in the developing brain and eyes. Lbh-like knockdown by morpholinos specifically inhibits expression of multiple photoreceptor-specific genes, such as opsins, gnat1, gnat2 and irbp. Interestingly, otx2 expression in the morphants is not significantly reduced until 32 hpf when lbh-like begins to express, but its expression level in 72 hpf morphants is higher than that in wild type embryos. Co-injection of otx2 and its downstream target neuroD mRNAs can rescue the faults in eyes of Lbh-like morphants. Combined with the results of promoter-reporter assay, we suggest that lbh-like is a new regulator of photoreceptor differentiation directly through affecting otx2 expression in zebrafish. Furthermore, knockdown of lbh-like increases the activity of Notch pathway and perturbs the balance among proliferation, differentiation and survival of photoreceptor precursors.

Long-term expandable SOX9+ chondrogenic ectomesenchymal cells from human pluripotent stem cells

Here we report the successful generation and long-term expansion of SOX9-expressing CD271⁺PDGFRα⁺CD73⁺ chondrogenic ectomesenchymal cells from the PAX3/SOX10/FOXD3-expressing MIXL1⁻CD271^hiPDGFRα^loCD73⁻ neural crest-like progeny of human pluripotent stem cells in a chemically defined medium supplemented with Nodal/Activin/transforming growth factorβ (TGFβ) inhibitor and fibroblast growth factor (FGF). When “primed” with TGFβ, such cells efficiently formed translucent cartilage particles, which were completely mineralized in 12 weeks in immunocompromized mice. The ectomesenchymal cells were expandable without loss of chondrogenic potential for at least 16 passages. They maintained normal karyotype for at least 10 passages and expressed genes representing embryonic progenitors (SOX4/12, LIN28A/B), cranial mesenchyme (ALX1/3/4), and chondroprogenitors (SOX9, COL2A1) of neural crest origin (SOX8/9, NGFR, NES). Ectomesenchyme is a source of many craniofacial bone and cartilage structures. The method we describe for obtaining a large quantity of human ectomesenchymal cells will help to model craniofacial disorders in vitro and potentially provide cells for the repair of craniofacial damage.

•

Generation of SOX9⁺ chondrogenic ectomesenchymal cells from hPSCs

•

Long-term expandability of the cells without loss of in vitro chondrogenicity

•

Need for TGFβ priming for cartilage particle formation

•

In vivo ossification capacity of cartilage developed in vitro

Embryonic and larval development in the Midas cichlid fish species flock (Amphilophus spp.): a new evo-devo model for the investigation of adaptive novelties and species differences

Background

Central American crater lake cichlid fish of the Midas species complex (Amphilophus spp.) are a model system for sympatric speciation and fast ecological diversification and specialization. Midas cichlids have been intensively analyzed from an ecological and morphological perspective. Genomic resources such as transcriptomic and genomic data sets, and a high-quality draft genome are available now. Many ecologically relevant species-specific traits and differences such as pigmentation and cranial morphology arise during development. Detailed descriptions of the early development of the Midas cichlid in particular, will help to investigate the ontogeny of species differences and adaptations.

Results

We describe the embryonic and larval development of the crater lake cichlid, Amphilophus xiloaensis, until seven days after fertilization. Similar to previous studies on teleost development, we describe six periods of embryogenesis - the zygote, cleavage, blastula, gastrula, segmentation, and post-hatching period. Furthermore, we define homologous stages to well-described teleost models such as medaka and zebrafish, as well as other cichlid species such as the Nile tilapia and the South American cichlid Cichlasoma dimerus. Key morphological differences between the embryos of Midas cichlids and other teleosts are highlighted and discussed, including the presence of adhesive glands and different early chromatophore patterns, as well as variation in developmental timing.

Conclusions

The developmental staging of the Midas cichlid will aid researchers in the comparative investigation of teleost ontogenies. It will facilitate comparative developmental biological studies of Neotropical and African cichlid fish in particular. In the past, the species flocks of the African Great Lakes have received the most attention from researchers, but some lineages of the 300–400 species of Central American lakes are fascinating model systems for adaptive radiation and rapid phenotypic evolution. The availability of genetic resources, their status as a model system for evolutionary research, and the possibility to perform functional experiments including transgenesis makes the Midas cichlid complex a very attractive model for evolutionary-developmental research.