Evolution of developmental regulation in the vertebrate FgfD subfamily

Insights into Craniofacial Development and Anomalies: Exploring Fgf Signaling in Zebrafish Models

PURPOSE OF REVIEW

To illustrate the value of using zebrafish to understand the role of the Fgf signaling pathway during craniofacial skeletal development under normal and pathological conditions.

RECENT FINDINGS

Recent data obtained from studies on zebrafish have demonstrated the genetic redundancy of Fgf signaling pathway and have identified new molecular partners of this signaling during the early stages of craniofacial skeletal development. Studies on zebrafish models demonstrate the involvement of the Fgf signaling pathway at every stage of craniofacial development. They particularly emphasize the central role of Fgf signaling pathway during the early stages of the development, which significantly impacts the formation of the various structures making up the craniofacial skeleton. This partly explains the craniofacial abnormalities observed in disorders associated with FGF signaling. Future research efforts should focus on investigating zebrafish Fgf signaling during more advanced stages, notably by establishing zebrafish models expressing mutations responsible for diseases such as craniosynostoses.

Enigmatic Nodal and Lefty gene repertoire discrepancy: Latent evolutionary history revealed by vertebrate-wide phylogeny

Homology in vertebrate body plans is traditionally ascribed to the high-level conservation of regulatory components within the genetic programs governing them, particularly during the "phylotypic stage." However, advancements in embryology and molecular phylogeny have unveiled the dynamic nature of gene repertoires responsible for early development. Notably, the Nodal and Lefty genes, members of the transforming growth factor-beta superfamily producing intercellular signaling molecules and crucial for left-right (L-R) symmetry breaking, exhibit distinctive features within their gene repertoires. These features encompass among-species gene repertoire variations resulting from gene gain and loss, as well as gene conversion. Despite their significance, these features have been largely unexplored in a phylogenetic context, but accumulating genome-wide sequence information is allowing the scrutiny of these features. It has exposed hidden paralogy between Nodal1 and Nodal2 genes resulting from differential gene loss in amniotes. In parallel, the tandem cluster of Lefty1 and Lefty2 genes, which was thought to be confined to mammals, is observed in sharks and rays, with an unexpected phylogenetic pattern. This article provides a comprehensive review of the current understanding of the origins of these vertebrate gene repertoires and proposes a revised nomenclature based on the elucidated history of vertebrate genome evolution.

Elucidating the toxicity mechanisms of organophosphate esters by adverse outcome pathway network

With the widespread use of organophosphate esters (OPEs), the accumulation and toxicity effect of OPEs in biota are attracting more and more concern. In order to clarify the mechanism of toxicity of OPEs to organisms, this study reviewed the OPEs toxicity and systematically identified the mechanism of OPEs toxicity under the framework of adverse outcome pathway (AOP). OPEs were divided into three groups (alkyl-OPEs, aryl-OPEs, and halogenated-OPEs) and biota was divided into aquatic organism and mammals. The results showed that tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) and triphenyl phosphate (TPHP) mainly caused neurotoxicity, reproductive, and hepatotoxicity in different mechanisms. According to the constructed AOP network, the toxicity mechanism of OPEs on aquatic organisms and mammals is different, which is mainly attributed to the different biological metabolic systems of aquatic organisms and mammals. Interestingly, our results indicate that the toxicity effect of the three kinds of OPEs on aquatic organisms is different, while there was no obvious difference in the mechanism of toxicity of OPEs on mammals. This study provides a theoretical basis for OPEs risk assessment in the future.

Becoming female: Ovarian differentiation from an evolutionary perspective

Differentiation of the bipotential gonadal primordium into ovaries and testes is a common process among vertebrate species. While vertebrate ovaries eventually share the same functions of producing oocytes and estrogens, ovarian differentiation relies on different morphogenetic, cellular, and molecular cues depending on species. The aim of this review is to highlight the conserved and divergent features of ovarian differentiation through an evolutionary perspective. From teleosts to mammals, each clade or species has a different story to tell. For this purpose, this review focuses on three specific aspects of ovarian differentiation: ovarian morphogenesis, the evolution of the role of estrogens on ovarian differentiation and the molecular pathways involved in granulosa cell determination and maintenance.

Genomic prediction of growth in a commercially, recreationally, and culturally important marine resource, the Australian snapper (Chrysophrys auratus)

Growth is one of the most important traits of an organism. For exploited species, this trait has ecological and evolutionary consequences as well as economical and conservation significance. Rapid changes in growth rate associated with anthropogenic stressors have been reported for several marine fishes, but little is known about the genetic basis of growth traits in teleosts. We used reduced genome representation data and genome-wide association approaches to identify growth-related genetic variation in the commercially, recreationally, and culturally important Australian snapper (Chrysophrys auratus, Sparidae). Based on 17,490 high-quality single-nucleotide polymorphisms and 363 individuals representing extreme growth phenotypes from 15,000 fish of the same age and reared under identical conditions in a sea pen, we identified 100 unique candidates that were annotated to 51 proteins. We documented a complex polygenic nature of growth in the species that included several loci with small effects and a few loci with larger effects. Overall heritability was high (75.7%), reflected in the high accuracy of the genomic prediction for the phenotype (small vs large). Although the single-nucleotide polymorphisms were distributed across the genome, most candidates (60%) clustered on chromosome 16, which also explains the largest proportion of heritability (16.4%). This study demonstrates that reduced genome representation single-nucleotide polymorphisms and the right bioinformatic tools provide a cost-efficient approach to identify growth-related loci and to describe genomic architectures of complex quantitative traits. Our results help to inform captive aquaculture breeding programs and are of relevance to monitor growth-related evolutionary shifts in wild populations in response to anthropogenic pressures.

Developmental constraints on fin diversity

The fish fin is a breathtaking repository full of evolutionary diversity, novelty, and convergence. Over 500 million years, the adaptation to novel habitats has provided landscapes of fin diversity. Although comparative anatomy of evolutionarily divergent patterns over centuries has highlighted the fundamental architectures and evolutionary trends of fins, including convergent evolution, the developmental constraints on fin evolution, which bias the evolutionary trajectories of fin morphology, largely remain elusive. Here, we review the evolutionary history, developmental mechanisms, and evolutionary underpinnings of paired fins, illuminating possible developmental constraints on fin evolution. Our compilation of anatomical and genetic knowledge of fin development sheds light on the canalized and the unpredictable aspects of fin shape in evolution. Leveraged by an arsenal of genomic and genetic tools within the working arena of spectacular fin diversity, evolutionary developmental biology embarks on the establishment of conceptual framework for developmental constraints, previously enigmatic properties of evolution.

Overactivity or blockade of transforming growth factor-β each generate a specific ureter malformation

Transforming growth factor-β (TGFβ) has been reported to be dysregulated in malformed ureters. There exists, however, little information on whether altered TGFβ levels actually perturb ureter development. We therefore hypothesised that TGFβ has functional effects on ureter morphogenesis. Tgfb1, Tgfb2 and Tgfb3 transcripts coding for TGFβ ligands, as well as Tgfbr1 and Tgfbr2 coding for TGFβ receptors, were detected by quantitative polymerase chain reaction in embryonic mouse ureters collected over a wide range of stages. As assessed by in situ hybridisation and immunohistochemistry, the two receptors were detected in embryonic urothelia. Next, TGFβ1 was added to serum-free cultures of embryonic day 15 mouse ureters. These organs contain immature smooth muscle and urothelial layers and their in vivo potential to grow and acquire peristaltic function can be replicated in serum-free organ culture. Such organs therefore constitute a suitable developmental stage with which to define roles of factors that affect ureter growth and functional differentiation. Exogenous TGFβ1 inhibited growth of the ureter tube and generated cocoon-like dysmorphogenesis. RNA sequencing suggested that altered levels of transcripts encoding certain fibroblast growth factors (FGFs) followed exposure to TGFβ. In serum-free organ culture exogenous FGF10 but not FGF18 abrogated certain dysmorphic effects mediated by exogenous TGFβ1. To assess whether an endogenous TGFβ axis functions in developing ureters, embryonic day 15 explants were exposed to TGFβ receptor chemical blockade; growth of the ureter was enhanced, and aberrant bud-like structures arose from the urothelial tube. The muscle layer was attenuated around these buds, and peristalsis was compromised. To determine whether TGFβ effects were limited to one stage, explants of mouse embryonic day 13 ureters, more primitive organs, were exposed to exogenous TGFβ1, again generating cocoon-like structures, and to TGFβ receptor blockade, again generating ectopic buds. As for the mouse studies, immunostaining of normal embryonic human ureters detected TGFβRI and TGFβRII in urothelia. Collectively, these observations reveal unsuspected regulatory roles for endogenous TGFβ in embryonic ureters, fine-tuning morphogenesis and functional differentiation. Our results also support the hypothesis that the TGFβ up-regulation reported in ureter malformations impacts on pathobiology. Further experiments are needed to unravel the intracellular signalling mechanisms involved in these dysmorphic responses. © 2019 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Double maternal-effect: duplicated nucleoplasmin 2 genes, npm2a and npm2b, with essential but distinct functions are shared by fish and tetrapods

BACKGROUND

Nucleoplasmin 2 (npm2) is an essential maternal-effect gene that mediates early embryonic events through its function as a histone chaperone that remodels chromatin. Recently, two npm2 (npm2a and npm2b) genes have been annotated in zebrafish. Thus, we examined the evolution of npm2a and npm2b in a variety of vertebrates, their potential phylogenetic relationships, and their biological functions using knockout models via the CRISPR/cas9 system.

RESULTS

We demonstrated that the two npm2 duplicates exist in a wide range of vertebrates, including sharks, ray-finned fish, amphibians, and sauropsids, while npm2a was lost in coelacanth and mammals, as well as some specific teleost lineages. Using phylogeny and synteny analyses, we traced their origins to the early stages of vertebrate evolution. Our findings suggested that npm2a and npm2b resulted from an ancient local gene duplication, and their functions diverged although key protein domains were conserved. We then investigated their functions by examining their tissue distribution in a wide variety of species and found that they shared ovarian-specific expression, a key feature of maternal-effect genes. We also demonstrated that both npm2a and npm2b are maternally-inherited transcripts in vertebrates, and that they play essential, but distinct, roles in early embryogenesis using zebrafish knockout models. Both npm2a and npm2b function early during oogenesis and may play a role in cortical granule function that impact egg activation and fertilization, while npm2b is also involved in early embryogenesis.

CONCLUSION

These novel findings will broaden our knowledge on the evolutionary history of maternal-effect genes and underlying mechanisms that contribute to vertebrate reproductive success. In addition, our results demonstrate the existence of a newly described maternal-effect gene, npm2a, that contributes to egg competence, an area that still requires further comprehension.

A Comparative View on Sex Differentiation and Gametogenesis Genes in Lungfish and Coelacanths

Gonadal sex differentiation and reproduction are the keys to the perpetuation of favorable gene combinations and positively selected traits. In vertebrates, several gonad development features that differentiate tetrapods and fishes are likely to be, at least in part, related to the water-to-land transition. The collection of information from basal sarcopterygians, coelacanths, and lungfishes, is crucial to improve our understanding of the molecular evolution of pathways involved in reproductive functions, since these organisms are generally regarded as “living fossils” and as the direct ancestors of tetrapods. Here, we report for the first time the characterization of >50 genes related to sex differentiation and gametogenesis in Latimeria menadoensis and Protopterus annectens. Although the expression profiles of most genes is consistent with the intermediate position of basal sarcopterygians between actinopterygian fish and tetrapods, their phylogenetic placement and presence/absence patterns often reveal a closer affinity to the tetrapod orthologs. On the other hand, particular genes, for example, the male gonad factor gsdf (Gonadal Soma-Derived Factor), provide examples of ancestral traits shared with actinopterygians, which disappeared in the tetrapod lineage.

Identification of transcripts potentially involved in neural tube closure using RNA sequencing

Anencephaly is a fatal human neural tube defect (NTD) in which the anterior neural tube remains open. Zebrafish embryos with reduced Nodal signaling display an open anterior neural tube phenotype that is analogous to anencephaly. Previous work from our laboratory suggests that Nodal signaling acts through induction of the head mesendoderm and mesoderm. Head mesendoderm/mesoderm then, through an unknown mechanism, promotes formation of the polarized neuroepithelium that is capable of undergoing the movements required for closure. We compared the transcriptome of embryos treated with a Nodal signaling inhibitor at sphere stage, which causes NTDs, to embryos treated at 30% epiboly, which does not cause NTDs. This screen identified over 3,000 transcripts with potential roles in anterior neurulation. Expression of several genes encoding components of tight and adherens junctions was significantly reduced, supporting the model that Nodal signaling regulates formation of the neuroepithelium. mRNAs involved in Wnt, FGF, and BMP signaling were also differentially expressed, suggesting these pathways might regulate anterior neurulation. In support of this, we found that pharmacological inhibition of FGF-receptor function causes an open anterior NTD as well as loss of mesodermal derivatives. This suggests that Nodal and FGF signaling both promote anterior neurulation through induction of head mesoderm.

The fibroblast growth factor 8 family in the female reproductive tract

Several growth factor families have been shown to be involved in the function of the female reproductive tract. One subfamily of the fibroblast growth factor (FGF) superfamily, namely the FGF8 subfamily (including FGF17 and FGF18), has become important as Fgf8 has been described as an oocyte-derived factor essential for glycolysis in mouse cumulus cells and aberrant expression ofFGF18has been described in ovarian and endometrial cancers. In this review, we describe the pattern of expression of these factors in normal ovaries and uteri in rodents, ruminants and humans, as well as the expression of their receptors and intracellular negative feedback regulators. Expression of these molecules in gynaecological cancers is also reviewed. The role of FGF8 and FGF18 in ovarian and uterine function is described, and potential differences between rodents and ruminants have been highlighted especially with respect to FGF18 signalling within the ovarian follicle. Finally, we identify major questions about the reproductive biology of FGFs that remain to be answered, including (1) the physiological concentrations within the ovary and uterus, (2) which cell types within the endometrial stroma and theca layer express FGFs and (3) which receptors are activated by FGF8 subfamily members in reproductive tissues.

Tris(1,3-dichloro-2-propyl) phosphate disrupts dorsoventral patterning in zebrafish embryos

Tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) is a high-production volume organophosphate flame retardant widely used within the United States. Within zebrafish, initiation of TDCIPP exposure at 0.75 h post-fertilization (hpf) results in genome-wide alterations in methylation during cleavage (2 hpf) as well as epiboly delay or arrest (at higher concentrations) during late-blastula and early-gastrula (4–6 hpf). To determine whether these TDCIPP-induced effects were associated with impacts on the transcriptome, embryos were exposed to vehicle (0.1% DMSO) or 2 µM TDCIPP from 0.75 hpf to 6 hpf, and total RNA was extracted from triplicate embryo pools per treatment and hybridized onto duplicate Affymetrix Zebrafish Gene 1.0 ST Arrays per RNA sample. Based on transcriptome-wide profiling, TDCIPP resulted in a significant impact on biological processes involved in dorsoventral patterning and bone morphogenetic protein (BMP) signaling. Consistent with these responses, TDCIPP exposure also resulted in strongly dorsalized embryos by 24 hpf—a phenotype that mimicked the effects of dorsomorphin, a potent and selective BMP inhibitor. Moreover, the majority of dorsalized embryos were preceded by epiboly arrest at 6 hpf. Our microarray data also revealed that the expression of sizzled (szl)—a gene encoding a secreted Frizzled-related protein that limits BMP signaling—was significantly decreased by nearly 4-fold at 6 hpf. Therefore, we used a splice-blocking morpholino to test the hypothesis that knockdown of szl phenocopies TDCIPP-induced delays in epiboly progression. Interestingly, contrary to our hypothesis, injection of szl MOs did not affect epiboly progression but, similar to chordin (chd) morphants, resulted in mildly ventralized embryos by 24 hpf. Overall, our findings suggest that TDCIPP-induced epiboly delay may not be driven by decreased szl expression, and that TDCIPP-induced dorsalization may—similar to dorsomorphin—be due to interference with BMP signaling during early zebrafish development.

Fibroblast growth factor signaling is required for early somatic gonad development in zebrafish

The vertebrate ovary and testis develop from a sexually indifferent gonad. During early development of the organism, primordial germ cells (the gamete lineage) and somatic gonad cells coalesce and begin to undergo growth and morphogenesis to form this bipotential gonad. Although this aspect of development is requisite for a fertile adult, little is known about the genetic regulation of early gonadogenesis in any vertebrate. Here, we provide evidence that fibroblast growth factor (Fgf) signaling is required for the early growth phase of a vertebrate bipotential gonad. Based on mutational analysis in zebrafish, we show that the Fgf ligand 24 (Fgf24) is required for proliferation, differentiation, and morphogenesis of the early somatic gonad, and as a result, most fgf24 mutants are sterile as adults. Additionally, we describe the ultrastructural elements of the early zebrafish gonad and show that distinct somatic cell populations can be identified soon after the gonad forms. Specifically, we show that fgf24 is expressed in an epithelial population of early somatic gonad cells that surrounds an inner population of mesenchymal somatic gonad cells that are in direct contact with the germ cells, and that fgf24 is required for stratification of the somatic tissue. Furthermore, based on gene expression analysis, we find that differentiation of the inner mesenchymal somatic gonad cells into functional cell types in the larval and early juvenile-stage gonad is dependent on Fgf24 signaling. Finally, we argue that the role of Fgf24 in zebrafish is functionally analogous to the role of tetrapod FGF9 in early gonad development.

Incorporating tree-thinking and evolutionary time scale into developmental biology

Phylogenetic approaches are indispensable in any comparative molecular study involving multiple species. These approaches are in increasing demand as the amount and availability of DNA sequence information continues to increase exponentially, even for organisms that were previously not extensively studied. Without the sound application of phylogenetic concepts and knowledge, one can be misled when attempting to infer ancestral character states as well as the timing and order of evolutionary events, both of which are frequently exerted in evolutionary developmental biology. The ignorance of phylogenetic approaches can also impact non-evolutionary studies and cause misidentification of the target gene or protein to be examined in functional characterization. This review aims to promote tree-thinking in evolutionary conjecture and stress the importance of a sense of time scale in cross-species comparisons, in order to enhance the understanding of phylogenetics in all biological fields including developmental biology. To this end, molecular phylogenies of several developmental regulatory genes, including those denoted as "cryptic pan-vertebrate genes", are introduced as examples.

Lineage-specific loss of FGF17 within the avian orders Galliformes and Passeriformes

The genomic and developmental complexity of vertebrates is commonly attributed to two rounds of whole genome duplications which occurred at the base of the vertebrate radiation. These duplications led to the rise of several, multi-gene families of developmental proteins like the fibroblast growth factors (FGFs); a signaling protein family which functions at various stages of embryonic development. One of the major FGF assemblages arising from these duplications is the FGF8 subfamily, which includes FGF8, FGF17, and FGF18 in tetrapods. While FGF8 and FGF18 are found in all tetrapods and are critical for embryonic survival, genomic analyses suggest putative loss of FGF17 in various lineages ranging from frogs and fish, to the chicken. This study utilizes 27 avian genomes in conjunction with molecular analyses of chicken embryos to confirm the loss of FGF17 in chicken as a true, biological occurrence. FGF17 is also missing in the turkey, black grouse, Japanese quail and northern bobwhite genomes. These species, along with chicken, form a monophyletic clade in the order Galliformes. Four additional species, members of the clade Passeroidea, within the order Passeriformes, are also missing FGF17. Additionally, analysis of intact FGF17 in other avian lineages reveals that it is still under strong purifying selection, despite being seemingly dispensable. Thus, FGF17 likely represents a molecular spandrel arising from a genome duplication event and due to its high connectivity with FGF8/FGF18, and potential for interference with their function, is retained under strong purifying selection, despite itself not having a strong selective advantage.

Dynamic expression of a Hydra FGF at boundaries and termini

Guidance of cells and tissue sheets is an essential function in developing and differentiating animal tissues. In Hydra, where cells and tissue move dynamically due to constant cell proliferation towards the termini or into lateral, vegetative buds, factors essential for guidance are still unknown. Good candidates to take over this function are fibroblast growth factors (FGFs). We present the phylogeny of several Hydra FGFs and analysis of their expression patterns. One of the FGFs is expressed in all terminal regions targeted by tissue movement and at boundaries crossed by moving tissue and cells with an expression pattern slightly differing in two Hydra strains. A model addressing an involvement of this FGF in cell movement and morphogenesis is proposed: Hydra FGFf-expressing cells might serve as sources to attract tissue and cells towards the termini of the body column and across morphological boundaries. Moreover, a function in morphogenesis and/or differentiation of cells and tissue is suggested.

Pleiotropic constraints, expression level, and the evolution of miRNA sequences

Post-transcriptional gene regulation mediated by microRNAs (miRNAs) plays critical roles during development by modulating gene expression and conferring robustness to stochastic errors. Phylogenetic analyses suggest that miRNA acquisition could play a role in phenotypic innovation. Moreover, miRNA-induced regulation strongly impacts genome evolution, increasing selective constraints on 3'UTRs, protein sequences, and expression level divergence. Thus, it is essential to understand the factors governing sequence evolution for this important class of regulatory molecules. Investigation of the patterns of molecular evolution at miRNA loci have been limited in Caenorhabditis elegans because of the lack of a close outgroup. Instead, I used Caenorhabditis briggsae as the focus point of this study because of its close relationship to Caenorhabditis sp. 9. I also corroborated the patterns of sequence evolution in Caenorhabditis using published orthologous relationships among miRNAs in Drosophila. In nematodes and in flies, miRNA sequence divergence is not influenced by the genomic neighborhood (i.e., intronic or intergenic) but is nevertheless affected by the genomic context because X-linked miRNAs evolve faster than autosomal miRNAs. However, this effect of chromosomal linkage can be explained by differential expression levels rather than a fast-X effect. The results presented here support a universal negative relationship between rates of molecular evolution and expression level, and suggest that mutations in highly expressed miRNAs are more likely to be deleterious because they potentially affect a larger number of target genes. Finally, I show that many single family member miRNAs evolve faster than miRNAs from multigene families and have limited functional scope, suggesting that they are not strongly integrated in gene regulatory networks.

Toxicogenomic responses of zebrafish embryos/larvae to tris(1,3-dichloro-2-propyl) phosphate (TDCPP) reveal possible molecular mechanisms of developmental toxicity

Tris(1,3-dichloro-2-propyl) phosphate (TDCPP) is frequently present in indoor dust and can be detected in human milk. In order to evaluate the effects of TDCPP on vertebrate development, zebrafish embryos/larvae were used as an animal model to examine developmental phenotypes and explore possible mechanisms of toxicity by employing microarrays and iTRAQ labeling quantitative proteomics. The results demonstrated that treatment with TDCPP (3 μM) from 0.75 h postfertilization (hpf) inhibited cell rearrangement at 4 hpf, caused delay in epiboly at 5.7 and 8.5 hpf, and led to abnormal development (e.g., short tail, reduced body size) and lethality between 14 and 45 hpf, which might be related with altered expression of genes regulating embryogenesis. Furthermore, trunk curvature was observed as the main phenotype in 96 hpf zebrafish larvae exposed to 1 or 3 μM TDCPP, possibly by changing somite formation and expression of proteins related to fast muscle and cartilage development. Collectively, our results suggest that exposure to TDCPP causes developmental toxicity in vertebrates and warrant the need for studies to evaluate the potential health risks of TDCPP to developing human embryos/infants/children, due to its frequent presence in indoor dust and potential for human exposure.

The African coelacanth genome provides insights into tetrapod evolution

It was a zoological sensation when a living specimen of the coelacanth was first discovered in 1938, as this lineage of lobe-finned fish was thought to have gone extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features . Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain, and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues demonstrate the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

Evolution of the tbx6/16 subfamily genes in vertebrates: insights from zebrafish

In any comparative studies striving to understand the similarities and differences of the living organisms at the molecular genetic level, the crucial first step is to establish the homology (orthology and paralogy) of genes between different organisms. Determination of the homology of genes becomes complicated when the genes have undergone a rapid divergence in sequence or when the involved genes are members of a gene family that has experienced a differential gain or loss of its constituents in different taxonomic groups. Organisms with duplicated genomes such as teleost fishes might have been especially prone to these problems because the functional redundancies provided by the duplicate copies of genes would have allowed a rapid divergence or loss of genes during evolution. In this study, we will demonstrate that much of the ambiguities in the determination of the homology between fish and tetrapod genes resulting from the problems like these can be eliminated by complementing the sequence-based phylogenies with nonsequence information, such as the exon-intron structure of a gene or the composition of a gene's genomic neighbors. We will use the Tbx6/16 subfamily genes of zebrafish (tbx6, tbx16, tbx24, and mga genes), which have been well known for the ambiguity of their evolutionary relationships to the Tbx6/16 subfamily genes of tetrapods, as an illustrative example. We will show that, despite the similarity of sequence and expression to the tetrapod Tbx6 genes, zebrafish tbx6 gene is actually a novel T-box gene more closely related to the tetrapod Tbx16 genes, whereas the zebrafish tbx24 gene, hitherto considered to be a novel gene due to the high level of sequence divergence, is actually an ortholog of tetrapod Tbx6 genes. We will also show that, after their initial appearance by the multiplication of a common ancestral gene at the beginning of vertebrate evolution, the Tbx6/16 subfamily of vertebrate T-box genes might have experienced differential losses of member genes in different vertebrate groups and gradual pooling of member gene's functions in surviving members, which might have prevented the revelation of the true identity of member genes by way of the comparison of sequence and function.

Vgll2a is required for neural crest cell survival during zebrafish craniofacial development

Invertebrate and vertebrate vestigial (vg) and vestigial-like (VGLL) genes are involved in embryonic patterning and cell fate determination. These genes encode cofactors that interact with members of the Scalloped/TEAD family of transcription factors and modulate their activity. We have previously shown that, in mice, Vgll2 is differentially expressed in the developing facial prominences. In this study, we show that the zebrafish ortholog vgll2a is expressed in the pharyngeal endoderm and ectoderm surrounding the neural crest derived mesenchyme of the pharyngeal arches. Moreover, both the FGF and retinoic acid (RA) signaling pathways, which are critical components of the hierarchy controlling craniofacial patterning, regulate this domain of vgll2a expression. Consistent with these observations, vgll2a is required within the pharyngeal endoderm for NCC survival and pharyngeal cartilage development. Specifically, knockdown of Vgll2a in zebrafish embryos using Morpholino injection results in increased cell death within the pharyngeal arches, aberrant endodermal pouch morphogenesis, and hypoplastic cranial cartilages. Overall, our data reveal a novel non-cell autonomous role for Vgll2a in development of the NCC-derived vertebrate craniofacial skeleton.

Evolution of pleiotropy: epistatic interaction pattern supports a mechanistic model underlying variation in genotype-phenotype map

The genotype-phenotype (GP) map consists of developmental and physiological mechanisms mapping genetic onto phenotypic variation. It determines the distribution of heritable phenotypic variance on which selection can act. Comparative studies of morphology as well as of gene regulatory networks show that the GP map itself evolves, yet little is known about the actual evolutionary mechanisms involved. The study of such mechanisms requires exploring the variation in GP maps at the population level, which presently is easier to quantify by statistical genetic methods rather than by regulatory network structures. We focus on the evolution of pleiotropy, a major structural aspect of the GP map. Pleiotropic genes affect multiple traits and underlie genetic covariance between traits, often causing evolutionary constraints. Previous quantitative genetic studies have demonstrated population-level variation in pleiotropy in the form of loci, at which genotypes differ in the genetic covariation between traits. This variation can potentially fuel evolution of the GP map under selection and/or drift. Here, we propose a developmental mechanism underlying population genetic variation in covariance and test its predictions. Specifically, the mechanism predicts that the loci identified as responsible for genetic variation in pleiotropy are involved in trait-specific epistatic interactions. We test this prediction for loci affecting allometric relationships between traits in an advanced intercross between inbred mouse strains. The results consistently support the prediction. We further find a high degree of sign epistasis in these interactions, which we interpret as an indication of adaptive gene complexes within the diverged parental lines.

Conservation of gene linkage in dispersed vertebrate NK homeobox clusters

Nk homeobox genes are important regulators of many different developmental processes including muscle, heart, central nervous system and sensory organ development. They are thought to have arisen as part of the ANTP megacluster, which also gave rise to Hox and ParaHox genes, and at least some NK genes remain tightly linked in all animals examined so far. The protostome-deuterostome ancestor probably contained a cluster of nine Nk genes: (Msx)-(Nk4/tinman)-(Nk3/bagpipe)-(Lbx/ladybird)-(Tlx/c15)-(Nk7)-(Nk6/hgtx)-(Nk1/slouch)-(Nk5/Hmx). Of these genes, only NKX2.6-NKX3.1, LBX1-TLX1 and LBX2-TLX2 remain tightly linked in humans. However, it is currently unclear whether this is unique to the human genome as we do not know which of these Nk genes are clustered in other vertebrates. This makes it difficult to assess whether the remaining linkages are due to selective pressures or because chance rearrangements have "missed" certain genes. In this paper, we identify all of the paralogs of these ancestrally clustered NK genes in several distinct vertebrates. We demonstrate that tight linkages of Lbx1-Tlx1, Lbx2-Tlx2 and Nkx3.1-Nkx2.6 have been widely maintained in both the ray-finned and lobe-finned fish lineages. Moreover, the recently duplicated Hmx2-Hmx3 genes are also tightly linked. Finally, we show that Lbx1-Tlx1 and Hmx2-Hmx3 are flanked by highly conserved noncoding elements, suggesting that shared regulatory regions may have resulted in evolutionary pressure to maintain these linkages. Consistent with this, these pairs of genes have overlapping expression domains. In contrast, Lbx2-Tlx2 and Nkx3.1-Nkx2.6, which do not seem to be coexpressed, are also not associated with conserved noncoding sequences, suggesting that an alternative mechanism may be responsible for the continued clustering of these genes.