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

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.

The phylogeographical pattern of the Amur minnow Rhynchocypris lagowskii (Cypriniformes: Cyprinidae) in the Qinling Mountains

In this study, the phylogeographical pattern of the Amur minnow (Rhynchocypris lagowskii) widely distributed in the cold freshwaters of the Qinling Mountains was examined. A total of 464 specimens from 48 localities were sequenced at a 540-bp region of the mitochondrial cytochrome b (Cytb) gene, and 69 haplotypes were obtained. The mean ratio of the number of synonymous and nonsynonymous substitutions per site (dN/dS) was 0.028 and indicated purifying selection. Haplotype diversity (h) and nucleotide diversity (π) of natural populations of R. lagowskii varied widely between distinct localities. Phylogenetic trees based on Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods, and network analysis showed five well-differentiated lineages, but these did not completely correspond to localities and geographic distribution. Meanwhile, analysis of molecular variances (AMOVA) indicated the highest proportion of genetic variation was attributed to the differentiation between populations rather than by our defined lineages. In addition, there was no significant correlation between the pairwise Fst values and geographic distance (p > .05). Based on the molecular clock calibration, the time to the most recent common ancestor (TMRCA) was estimated to have emerged from the Late Miocene to the Early Pleistocene. Finally, the results of demographic history based on the neutrality test, mismatch distribution, and Bayesian skyline plot (BSP) analyses showed that collectively, the populations were stable during the Pleistocene while one lineage (lineage E) probably underwent a slight contraction during the Middle Pleistocene and a rapid expansion from the Middle to the Late Pleistocene. Therefore, the study suggests the current phylogeographical pattern of R. lagowskii was likely shaped by geological events that led to vicariance followed by dispersal and secondary contact, river capture, and climatic oscillation during the Late Miocene to the Early Pleistocene in the Qinling Mountains.

Contrasting patterns of molecular evolution in metazoan germ line genes

BACKGROUND

Germ lines are the cell lineages that give rise to the sperm and eggs in animals. The germ lines first arise from primordial germ cells (PGCs) during embryogenesis: these form from either a presumed derived mode of preformed germ plasm (inheritance) or from an ancestral mechanism of inductive cell-cell signalling (induction). Numerous genes involved in germ line specification and development have been identified and functionally studied. However, little is known about the molecular evolutionary dynamics of germ line genes in metazoan model systems.

RESULTS

Here, we studied the molecular evolution of germ line genes within three metazoan model systems. These include the genus Drosophila (N=34 genes, inheritance), the fellow insect Apis (N=30, induction), and their more distant relative Caenorhabditis (N=23, inheritance). Using multiple species and established phylogenies in each genus, we report that germ line genes exhibited marked variation in the constraint on protein sequence divergence (dN/dS) and codon usage bias (CUB) within each genus. Importantly, we found that de novo lineage-specific inheritance (LSI) genes in Drosophila (osk, pgc) and in Caenorhabditis (pie-1, pgl-1), which are essential to germ plasm functions under the derived inheritance mode, displayed rapid protein sequence divergence relative to the other germ line genes within each respective genus. We show this may reflect the evolution of specialized germ plasm functions and/or low pleiotropy of LSI genes, features not shared with other germ line genes. In addition, we observed that the relative ranking of dN/dS and of CUB between genera were each more strongly correlated between Drosophila and Caenorhabditis, from different phyla, than between Drosophila and its insect relative Apis, suggesting taxonomic differences in how germ line genes have evolved.

CONCLUSIONS

Taken together, the present results advance our understanding of the evolution of animal germ line genes within three well-known metazoan models. Further, the findings provide insights to the molecular evolution of germ line genes with respect to LSI status, pleiotropy, adaptive evolution as well as PGC-specification mode.

Towards a gene regulatory network shaping the fins of the Princess cichlid

Variation in fin shape and size contributes to the outstanding morphological diversity of teleost fishes, but the regulation of fin growth has not yet been studied extensively outside the zebrafish model. A previous gene expression study addressing the ornamental elongations of unpaired fins in the African cichlid fish Neolamprologus brichardi identified three genes (cx43, mmp9 and sema3d) with strong and consistent expression differences between short and elongated fin regions. Remarkably, the expression patterns of these genes were not consistent with inferences on their regulatory interactions in zebrafish. Here, we identify a gene expression network (GRN) comprising cx43, mmp9, and possibly also sema3d by a stepwise approach of identifying co-expression modules and predicting their upstream regulators. Among the transcription factors (TFs) predicted as potential upstream regulators of 11 co-expressed genes, six TFs (foxc1, foxp1, foxd3, myc, egr2, irf8) showed expression patterns consistent with their cooperative transcriptional regulation of the gene network. Some of these TFs have already been implicated in teleost fish fin regeneration and formation. We particularly discuss the potential function of foxd3 as driver of the network and its role in the unexpected gene expression correlations observed in N. brichardi.

Germ layers, the neural crest and emergent organization in development and evolution

Discovered in chick embryos by Wilhelm His in 1868 and named the neural crest by Arthur Milnes Marshall in 1879, the neural crest cells that arise from the neural folds have since been shown to differentiate into almost two dozen vertebrate cell types and to have played major roles in the evolution of such vertebrate features as bone, jaws, teeth, visceral (pharyngeal) arches, and sense organs. I discuss the discovery that ectodermal neural crest gave rise to mesenchyme and the controversy generated by that finding; the germ layer theory maintained that only mesoderm could give rise to mesenchyme. A second topic of discussion is germ layers (including the neural crest) as emergent levels of organization in animal development and evolution that facilitated major developmental and evolutionary change. The third topic is gene networks, gene co-option, and the evolution of gene-signaling pathways as key to developmental and evolutionary transitions associated with the origin and evolution of the neural crest and neural crest cells.

A gene expression study of ornamental fin shape in Neolamprologus brichardi, an African cichlid species

The diversity of fin morphology within and across fish taxa offers great, but still largely unexplored, opportunities to investigate the proximate mechanisms underlying fin shape variation. Relying on available genetic knowledge brought forth mainly by the comprehensive study of the zebrafish caudal fin, we explored candidate molecular mechanisms for the maintenance and formation of the conspicuously elongated filaments adorning the unpaired fins of the East African "princess cichlid" Neolamprologus brichardi. Via qPCR assays, we detected expression differences of candidate genes between elongated and short regions of intact and regenerating fins. The identified genes include skeletogenic and growth factors (igf2b, fgf3, bmp2 and bmp4), components of the WNT pathway (lef1, wnt5b and wnt10) and a regulatory network determining fin ray segment size and junction (cx43, esco2 and sema3d), as well as other genes with different roles (mmp9, msxb and pea3). Interestingly, some of these genes showed fin specific expression differences which are often neglected in studies of model fish that focus on the caudal fin. Moreover, while the observed expression patterns were generally consistent with zebrafish results, we also detected deviating expression correlations and gene functions.