Neural crest cells as a source of microevolutionary variation

Variable Craniofacial Shape and Development among Multiple Cave-Adapted Populations of Astyanax mexicanus

Astyanax mexicanus is a freshwater fish species with blind cave morphs and sighted surface morphs. Like other troglodytic species, independently evolved cave-dwelling A. mexicanus populations share several stereotypic phenotypes, including the expansion of certain sensory systems, as well as the loss of eyes and pigmentation. Here, we assess the extent to which there is also parallelism in craniofacial development across cave populations. Since multiple forces may be acting upon variation in the A. mexicanus system, including phylogenetic history, selection, and developmental constraint, several outcomes are possible. For example, eye regression may have triggered a conserved series of compensatory developmental events, in which case we would expect to observe highly similar craniofacial phenotypes across cave populations. Selection for cave-specific foraging may also lead to the evolution of a conserved craniofacial phenotype, especially in regions of the head directly associated with feeding. Alternatively, in the absence of a common axis of selection or strong developmental constraints, craniofacial shape may evolve under neutral processes such as gene flow, drift, and bottlenecking, in which case patterns of variation should reflect the evolutionary history of A. mexicanus. Our results found that cave-adapted populations do share certain anatomical features; however, they generally did not support the hypothesis of a conserved craniofacial phenotype across caves, as nearly every pairwise comparison was statistically significant, with greater effect sizes noted between more distantly related cave populations with little gene flow. A similar pattern was observed for developmental trajectories. We also found that morphological disparity was lower among all three cave populations versus surface fish, suggesting eye loss is not associated with increased variation, which would be consistent with a release of developmental constraint. Instead, this pattern reflects the relatively low genetic diversity within cave populations. Finally, magnitudes of craniofacial integration were found to be similar among all groups, meaning that coordinated development among anatomical units is robust to eye loss in A. mexicanus. We conclude that, in contrast to many conserved phenotypes across cave populations, global craniofacial shape is more variable, and patterns of shape variation are more in line with population structure than developmental architecture or selection.

Distinct genetic origins of eumelanin levels and barring patterns in cichlid fishes

Pigment patterns are incredibly diverse across vertebrates and are shaped by multiple selective pressures from predator avoidance to mate choice. A common pattern across fishes, but for which we know little about the underlying mechanisms, is repeated melanic vertical bars. To understand the genetic factors that modify the level or pattern of vertical barring, we generated a genetic cross of 322 F2 hybrids between two cichlid species with distinct barring patterns, Aulonocara koningsi and Metriaclima mbenjii. We identify 48 significant quantitative trait loci that underlie a series of seven phenotypes related to the relative pigmentation intensity, and four traits related to patterning of the vertical bars. We find that genomic regions that generate variation in the level of eumelanin produced are largely independent of those that control the spacing of vertical bars. Candidate genes within these intervals include novel genes and those newly-associated with vertical bars, which could affect melanophore survival, fate decisions, pigment biosynthesis, and pigment distribution. Together, this work provides insights into the regulation of pigment diversity, with direct implications for an animal's fitness and the speciation process.

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.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Associated congenital heart disease with Hirschsprung's disease: a retrospective cohort study on 2,174 children

Objective

To examine the incidence and phenotypes of congenital heart disease (CHD) in a large cohort of patients with Hirschsprung's disease (HSCR).

Study design

Retrospective data review of children with HSCR between 2003 and 2020 was conducted at the Provincial Key Laboratory for Structural Birth Defects in Guangzhou, Guangdong, China. HSCR was confirmed by pathological diagnosis. CHD was defined as a gross structural abnormality of the heart or intrathoracic great vessels that is of functional significance.

Results

A total of 2,174 HSCR patients (84.7% males) were studied and 306 of them underwent echocardiography. Overall, 27 children (1.2%) had associated CHD. Among them, CHDs mostly presented as atrial and ventricular septal defects (n = 5 and 12 respectively) and patent ductus arteriosus (n = 4). Three patients (1.4‰) presented as a severe CHD including complete atrioventricular canal, congenitally corrected transposition of the great arteries and double-outlet of right ventricle. Among 14 patients carrying a chromosomal abnormality, CHD was detected in 4 infants (28.6%), all being mild forms of septal defects.

Conclusions

Some new and severe types of CHD were found in patients with HSCR. Patients with syndromic features had higher incidence of CHD.

Single-Cell Transcriptome and Pigment Biochemistry Analysis Reveals the Potential for the High Nutritional and Medicinal Value of Purple Sea Cucumbers

The sea cucumber Apostichopus japonicus has important nutritional and medicinal value. Unfortunately, we know little of the source of active chemicals in this animal, but the plentiful pigments of these animals are thought to function in intriguing ways for translation into clinical and food chemistry usage. Here, we found key cell groups with the gene activity predicted for the color morphology of sea cucumber body using single-cell RNA-seq. We refer to these cell populations as melanocytes and quinocytes, which are responsible for the synthesis of melanin and quinone pigments, respectively. We integrated analysis of pigment biochemistry with the transcript profiles to illuminate the molecular mechanisms regulating distinct pigment formation in echinoderms. In concert with the correlated pigment analysis from each color morph, this study expands our understanding of medically important pigment production, as well as the genetic mechanisms for color morphs, and provides deep datasets for exploring advancements in the fields of bioactives and nutraceuticals.

Developmental origin underlies evolutionary rate variation across the placental skull

The placental skull has evolved into myriad forms, from longirostrine whales to globular primates, and with a diverse array of appendages from antlers to tusks. This disparity has recently been studied from the perspective of the whole skull, but the skull is composed of numerous elements that have distinct developmental origins and varied functions. Here, we assess the evolution of the skull's major skeletal elements, decomposed into 17 individual regions. Using a high-dimensional morphometric approach for a dataset of 322 living and extinct eutherians (placental mammals and their stem relatives), we quantify patterns of variation and estimate phylogenetic, allometric and ecological signal across the skull. We further compare rates of evolution across ecological categories and ordinal-level clades and reconstruct rates of evolution along lineages and through time to assess whether developmental origin or function discriminate the evolutionary trajectories of individual cranial elements. Our results demonstrate distinct macroevolutionary patterns across cranial elements that reflect the ecological adaptations of major clades. Elements derived from neural crest show the fastest rates of evolution, but ecological signal is equally pronounced in bones derived from neural crest and paraxial mesoderm, suggesting that developmental origin may influence evolutionary tempo, but not capacity for specialisation. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Distinct genetic origins of eumelanin intensity and barring patterns in cichlid fishes

Pigment patterns are incredibly diverse across vertebrates and are shaped by multiple selective pressures from predator avoidance to mate choice. A common pattern across fishes, but for which we know little about the underlying mechanisms, is repeated melanic vertical bars. In order to understand genetic factors that modify the level or pattern of vertical barring, we generated a genetic cross of 322 F2 hybrids between two cichlid species with distinct barring patterns, Aulonocara koningsi and Metriaclima mbenjii. We identify 48 significant quantitative trait loci that underlie a series of seven phenotypes related to the relative pigmentation intensity, and four traits related to patterning of the vertical bars. We find that genomic regions that generate variation in the level of eumelanin produced are largely independent of those that control the spacing of vertical bars. Candidate genes within these intervals include novel genes and those newly-associated with vertical bars, which could affect melanophore survival, fate decisions, pigment biosynthesis, and pigment distribution. Together, this work provides insights into the regulation of pigment diversity, with direct implications for an animal's fitness and the speciation process.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underly facial shape variation, yet how those in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest even during the late migratory phase results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both Tfap2 members dysregulated numerous midface GRN components involved in midface fusion, patterning, and differentiation. Notably, Alx1/3/4 (Alx) transcript levels are reduced, while ChIP-seq analyses suggest TFAP2 directly and positively regulates Alx gene expression. TFAP2 and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish further implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a mutant zebrafish present abnormal alx3 expression patterns, and the two genes display a genetic interaction in this species. Together, these data demonstrate a critical role for TFAP2 in regulating vertebrate midfacial development in part through ALX transcription factor gene expression.

Shared reproductive disruption, not neural crest or tameness, explains the domestication syndrome

Altered neural crest cell (NCC) behaviour is an increasingly cited explanation for the domestication syndrome in animals. However, recent authors have questioned this explanation, while others cast doubt on whether domestication syndrome even exists. Here, we review published literature concerning this syndrome and the NCC hypothesis, together with recent critiques of both. We synthesize these contributions and propose a novel interpretation, arguing shared trait changes under ancient domestication resulted primarily from shared disruption of wild reproductive regimes. We detail four primary selective pathways for 'reproductive disruption' under domestication and contrast these succinct and demonstrable mechanisms with cryptic genetic associations posited by the NCC hypothesis. In support of our perspective, we illustrate numerous important ways in which NCCs contribute to vertebrate reproductive phenotypes, and argue it is not surprising that features derived from these cells would be coincidentally altered under major selective regime changes, as occur in domestication. We then illustrate several pertinent examples of Darwin's 'unconscious selection' in action, and compare applied selection and phenotypic responses in each case. Lastly, we explore the ramifications of reproductive disruption for wider evolutionary discourse, including links to wild 'self-domestication' and 'island effect', and discuss outstanding questions.

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.