The origin of bmp16, a novel Bmp2/4relative, retained in teleost fish genomes
© Feiner et al; licensee BioMed Central Ltd. 2009
Received: 22 July 2009
Accepted: 1 December 2009
Published: 1 December 2009
Whole genome sequences have allowed us to have an overview of the evolution of gene repertoires. The target of the present study, the TGFβ superfamily, contains many genes involved in vertebrate development, and provides an ideal system to explore the relationships between evolution of gene repertoires and that of developmental programs.
As a result of a bioinformatic survey of sequenced vertebrate genomes, we identified an uncharacterized member of the TGFβ superfamily, designated bmp16, which is confined to teleost fish species. Our molecular phylogenetic study revealed a high affinity of bmp16 to the Bmp2/4 subfamily. Importantly, further analyses based on the maximum-likelihood method unambiguously ruled out the possibility that this teleost-specific gene is a product of teleost-specific genome duplication. This suggests that the absence of a bmp16 ortholog in tetrapods is due to a secondary loss. In situ hybridization showed embryonic expression of the zebrafish bmp16 in the developing swim bladder, heart, tail bud, and ectoderm of pectoral and median fin folds in pharyngula stages, as well as gut-associated expression in 5-day embryos.
Comparisons of expression patterns revealed (1) the redundancy of bmp16 expression with its homologs in presumably plesiomorphic expression domains, such as the fin fold, heart, and tail bud, which might have permitted its loss in the tetrapod lineage, and (2) the loss of craniofacial expression and gain of swim bladder expression of bmp16 after the gene duplication between Bmp2, -4 and -16. Our findings highlight the importance of documenting secondary changes of gene repertoires and expression patterns in other gene families.
In the pre-genomic era, genes were identified and sequenced in a heuristic manner to address specific biological questions. After a long period of such heuristic exploration, many genes are now comprehensively identified as parts of genomic sequences determined in large-scale projects. With several relatively well-assembled genome sequences both for teleost fishes and tetrapods, one might expect that we have discovered almost all genes commonly possessed by vertebrates (see ). This should be proven with care, however, especially for gene families whose members were identified and functionally characterized intensively in the very early phase of development of molecular biology. This caution is emphasized when we take into account a generalized view that diverse species possess similar gene repertoires, for well-characterized gene families, such as those containing many regulatory genes involved in animal development, namely 'toolkit genes' (; also see ), such as Hox genes .
One of those is a group of genes containing vertebrate bone morphogenetic proteins (bmp) 2, bmp4, and fly decapentaplegic (dpp) genes. These genes are members of the transforming growth factor (TGF) β superfamily, which encode proteins that are cleaved to be active extracellular signaling ligands [5, 6]. The TGFβ superfamily contains other genes encoding signal proteins, such as Bmp2-15, growth differentiation factors (GDFs), TGFβ, inhibins and activins (note that Bmp1 [tolloid] is a protocollagen C-end peptidase that belongs to the Zinc-dependent metalloprotease family ). In vertebrates, Bmp2 and -4 are involved in various developmental processes, such as axial patterning, tissue specification and organogenesis (e.g. [8–11]; reviewed in ). In particular, there are several pioneering studies that highlighted the importance of this group of genes in understanding the molecular mechanisms that underlie morphological changes during evolution (e.g. [13–17]).
In silico identification of a novel bmp2/4homolog in teleost genomes
Ensembl entries for bmp16 genes
Ensembl Gene ID
Danio rerio (zebrafish)
Gasterosteus aculeatus (stickleback)
Oryzias latipes (medaka)
Structural properties of bmp16
By means of RT-PCR, we amplified 5' and 3' cDNA fragments of the zebrafish bmp16 (see Methods for details). By assembling the two cDNA fragments, we obtained a contiguous cDNA sequence containing a small fraction of 5' UTR, full protein-coding region and a full 3' UTR, which was submitted to EMBL Nucleotide Sequence Database under the accession number [EMBL:FN400947]. In an alignment between two predicted transcript sequences currently available in Ensembl (ENSDART00000098425 and ENSDART00000098426) and our assembled bmp16 cDNA, we found distinct differences in nucleotide sequences, probably due to putative mis-annotation of exons in the Ensembl entries (data not shown).
We then aligned our cDNA with genomic DNA sequences retrieved from Ensembl, to find out correct exon-intron boundaries (Additional file 1: Figure S1). It turned out that the differences between our and predicted Ensembl transcripts originates partly from an incorrect assembly of genomic sequences, which is verified by our PCR amplification of the corresponding genomic region (Additional file 1: Figure S1). Our genomic sequence was also submitted to EMBL Nucleotide Sequence Database under the accession number [EMBL:FN400946].
In the deduced amino acid sequences of bmp16, it was found that six amino acid residues, out of nine that for human Bmp2 interact with its receptor , are rigidly conserved between all teleost fishes analyzed here (Figure 2).
In silico scan for locations of the bmp16 gene in the genomes of teleost fishes, except zebrafish, highlighted its linkage with leucine-rich repeat containing (lrrc) 68, gem-associated protein (gemin) 7, splicing factor arginine/serine-rich (sfrs) 16, reticulon (rtn) 2 and relB genes in a approximately 500-kb long genomic stretch (Additional file 2: Figure S2). In the human genome, their orthologs were found tightly clustered in a genomic stretch with a similar size on chromosome 19 (Additional file 2: Figure S2).
Molecular phylogenetic analyses: in or out?
Statistical supports for alternative tree topologies for relationships among Bmp2/4/16 genes
(((((PmB, PmC), (PmA, bmp16)), Bmp2), Bmp4), Out)
((((((PmB, PmC), bmp16), PmA), Bmp2), Bmp4), Out)
1.49 ± 3.77
((((((PmB, PmC), bmp16), PmA), Bmp4), Bmp2), Out)
3.72 ± 4.68
(((((PmB, PmC), bmp16), PmA),(Bmp4, Bmp2)), Out)
3.77 ± 4.59
(((((PmB, PmC), (PmA, bmp16)), Bmp4), Bmp2), Out)
2.22 ± 2.55
((((((bmp16, PmC), PmB), PmA), Bmp2), Bmp4), Out)
5.33 ± 5.81
((((((bmp16, PmB), PmC), PmA), Bmp2), Bmp4), Out)
5.33 ± 5.81
((((PmB, PmC), (PmA, bmp16)),(Bmp2, Bmp4)), Out)
2.27 ± 2.42
(((Bmp4, (Bmp2, bmp16)), (PmA, (PmB, PmC))), Out)
13.81 ± 5.99
((((Bmp2, PmC), (bmp16, PmA)), (Bmp4, PmB)), Out)
30.84 ± 10.87
Gene expression analysis of bmp16
Between 36 and 48 hpf, bmp16 is expressed in the apical ectodermal fold, a specialized epithelial structure homologous to the tetrapod apical ectodermal ridge (AER), of the pectoral fins (Figure 6F and 6H). Around 36 hpf, bmp16 also starts to be expressed in a region of the midline endoderm at the level of the anterior extent of the pectoral fins (Figure 6H). To better characterize these cells we compared bmp16 expression with that of the marker for the dorsal pancreatic bud, islet-1 (isl-1; ), and found that bmp16 localizes to endodermal cells at the boundary between somites 2 and 3, while isl-1 is expressed slightly more caudally at the somite 3-4 boundary (Figure 6H and 6I). At 52 hpf, bmp16 expression continues to align with the midline (Figure 6G), and by 72 hpf is found in an additional, more anterior patch of cells that is associated with the foregut endoderm, whereas the posterior expression domain marks the developing swim bladder  (Figure 6J). Accordingly, at 5 dpf bmp16 is strongly expressed in the swim bladder as well as in the gut epithelium (Figure 6K and 6L).
Timings of gene duplications between Bmp2, -4 and -16
Since the teleost-specific genome duplication was first proposed, many genes unique to teleost genomes have been attributed to this event. This event per se is one of the best systems to characterize the evolutionary impact of whole genome duplications [4, 23]. However, our finding of the novel bmp16 genes unique to teleost genomes does not align with this framework. All teleost genomes examined so far have a single ortholog (bmp4) to tetrapod Bmp4, while two Bmp2 orthologs (bmp2a and -2b) are found in the zebrafish. Notably, our molecular phylogenetic analyses did not place the bmp16 into either of Bmp2 or Bmp4 subfamilies (Figure 1B and 4). Instead, the closer relationship of bmp16 to three lamprey Bmp2/4 homologs was supported (Figure 5A). Based on statistical evaluation of likelihood values, we conclude that the origin of bmp16 dates back to the earliest phase of vertebrate evolution. Despite our additional attempt to characterize the timing of the gene duplication by comparing gene order in the flanking genomic regions (syntenies) between teleost bmp16, -2, and -4, no convincing evidence for chromosomal duplications and their timings was obtained, probably because of frequent chromosomal rearrangements in their adjacent genomic regions (data not shown). Overall, our analysis favors the scenario that all the gene duplications giving rise to Bmp2, -4, and -16 occurred before the cyclostome-gnathostome split. If these three subtypes are products of the so-called two-round whole genome duplications, their timing is consistent with an emerging hypothesis that the common ancestor of all extant vertebrates already experienced these events [24, 25].
According to a previous study, the coverage of the available sea lamprey genome sequences should be sufficient for detecting at least one exon for any given gene --consistently, fragments of the previously reported three sea lamprey Bmp2/4 homologs were detected. Additionally, our effort to survey publicly available sea lamprey genomic and transcriptomic sequences failed to detect additional cyclostome Bmp2/4/16 homologs. Therefore, after our identification of bmp16, both the lamprey and gnathostomes possess three Bmp2/4/16-related homologs. However, the most parsimonious phylogenetic scenario that postulates the smallest number of secondary gene duplications and losses (Figure 5C) was significantly rejected (Table 2). If we assume the topology of the ML tree in Figure 4 to be correct, it suggests that only P. marinus Bmp2/4-A is the genuine ortholog of bmp16, and that the gnathostome ortholog for P. marinus Bmp2/4-B and -C was secondarily lost in its early evolution. Accordingly, we should also assume that cyclostome orthologs for gnathostome Bmp2 and/or Bmp4 are unidentified or secondarily lost. In the light of low statistical support for the ML tree, we cannot rule out other similar hypotheses (Table 2). However, with the currently available dataset and traditional framework of molecular phylogenetics, the scenario for gene duplications and retentions of duplicates during early vertebrate evolution seems to be much more complicated than previously reported . Gene sampling for other early vertebrates will be a future task. In particular, the group of hagfishes, which diverged from the future lamprey lineage just after the cyclostome-gnathostome split , might provide information for resolving deep divergence of bmp16 from others.
Bmp2/4/16gene repertoires: teleost fishes and tetrapods
As well as between gnathostomes and the lamprey, our analysis has highlighted the contrast of gene repertoires between teleosts and tetrapods. In any tetrapod species for which whole genome sequences are available, bmp16 orthologs have not been found. The absence of this gene in the entire taxon Tetrapoda was also confirmed by our intensive survey of ESTs in NCBI. Judging from the maintenance of syntenic organization of genes surrounding the bmp16 gene even in species without its ortholog (Additional file 2: Figure S2), the possible loss of this gene doesn't seem to have involved deletion of a large genomic portion. The most parsimonious explanation is that the bmp16 ortholog was lost from the genome in a small-scale event once at the basal lineage of tetrapods. A careful scan of sarcopterygian fish genomes (namely, of lungfish and coelacanth) could narrow the time range of the gene loss.
Another interesting observation is that even though the zebrafish genome contains two bmp2 orthologs (namely bmp2a and -2b), other model teleost fishes categorized in Acanthomorpha seem to possess only one, namely bmp2b (Figure 4). Our additional survey in all available EST sequences for Acanthomorpha did not detect any more unidentified bmp2 relative. This suggests that bmp2a is confined to at least some of species in Otocephala (e.g. zebrafish, Mexican cavefish), and the bmp2a ortholog was secondarily lost in the lineage leading to Acanthomorpha. This line of evidence for gene losses underlines an importance of elaborate comparisons of gene repertoires even between relatively closely related animals. It is likely that more 'toolkit' genes involved in vertebrate development, which are normally considered to be present universally in all relevant genomes, are revealed to have experienced this type of evolutionary change, as shown in some previous studies (e.g. ).
Roles of bmp16during zebrafish embryogenesis: did redundancy permit the loss?
In contrast to the divergent nature of bmp16 from other Bmp2/4 members, the deduced amino acid sequences of bmp16 genes seem to be relatively conserved among teleosts, especially within Acanthomorpha (Figure 2). We confirmed transcription of this uncharacterized gene in zebrafish with RT-PCR and in situ hybridization, and in medaka with a survey of publicly available ESTs. These data, together with the conservation of the majority of amino acid residues that are crucial for recognition by the BMP type I receptor , suggest that bmp16 functions as an active ligand in the TGFβ signaling pathway in particular biological processes.
Our analysis revealed bmp16 expressions in several different tissues and cell types during zebrafish embryogenesis, some of which are known to express other members of the TGFβ superfamily in tetrapods. In light of our molecular phylogenetic evidence of the loss of the bmp16 ortholog in the tetrapod lineage, one can ask whether its original functions might have been taken over by other members of the TGFβ superfamily in the tetrapods.
This assumption holds true for pectoral fins, because in addition to bmp2b and bmp4 [29–31], several other bmp genes (e.g. bmp6, bmp7b [32, 33]) are co-expressed in the zebrafish apical ectodermal fold. In addition, expression of bmp2a, bmp2b, bmp4 and bmp7 is up-regulated in the epithelium of the developing caudal fin fold [29–32]. The AER of tetrapod limbs also expresses a number of Bmp genes (reviewed in ). Therefore, individual Bmp subtypes are likely to act in a redundant fashion. The loss of the bmp16 ortholog in the tetrapod ancestor might have been compensated in the developing specialized epithelia of median fins and paired appendages by a redundancy of Bmp gene function in these tissues.
Similarly, bmp2b and bmp4, as well as bmp16, are commonly expressed in the developing zebrafish heart [35–37]. Possible deleterious effects associated with the loss of the bmp16 ortholog in the tetrapod lineage could eventually have been compensated by the expression of Bmp2 and Bmp4, and the recruitment of additional homologs along with the evolution of an advanced chambering system. Accordingly, tetrapods not only express Bmp2 and Bmp4 during heart development, but also Bmp5 and Bmp7 [9, 10, 38, 39].
Because the swim bladder and mammalian lungs share a common evolutionary origin , it is likely that their development is controlled by shared regulatory pathways. In the mouse, Bmp4 is expressed during lung development where it is initially found at high levels in the lung buds, followed by expression in the distal epithelial tips and with lower levels in the adjacent mesenchyme . Bmp4 promotes the development of distal epithelial cell types at the expense of proximal cell types [42, 43]. Accordingly, the proximal lung epithelium expresses the Bmp4 antagonist, Gremlin1, which determines proximal lung epithelial cell fate through inhibition of Bmp signaling. By regulating the proliferation, survival and morphogenetic activity of distal lung epithelial cells, Bmp signaling plays an essential role in lung morphogenesis [41, 42, 44–49].
In zebrafish, it has been shown that signaling through Bmp2a/2b regulates the decision between hepatic versus pancreatic fate in the endodermal midline [50, 51]. Bmp signaling in adjacent lateral plate mesodermal tissues patterns the endoderm through direct contact, as disruption of coordinated cell movements between these tissues causes patterning defects in the endoderm and ceased proliferation and apoptosis of cells in the liver primordium that results from a failure of bmp2a expressing cells to properly contact the endoderm .
Comparably little is known about the roles of Bmps during swim bladder development. However, in the zebrafish, the swim bladder expresses the protein related to Dan and Cerberus (prdc) gene (alias: gremlin2; ), which encodes an antagonist of signaling through bmp2a/2b and bmp4 . It will therefore be interesting to see whether bmp16 and prdc play comparable antagonistic roles during swim bladder development like Bmp4/Gremlin1 in the mouse lung.
Interestingly, we observed bmp16 expression in the tail bud between 24 and 36 hpf. Although bmp2a and bmp4 are transiently expressed in a broad domain of the caudal trunk-forming mesoderm after the end of gastrulation , neither gene is expressed at late stages of somitogenesis. Due to its dynamic expression in the developing neural keel, it is possible that bmp16 plays a role during the ongoing neurulation process in the tail, a process that probably also involves bmp7a . Its expression in the zebrafish tail may correspond to that in the apical surface of the neural tube in tetrapods. In the mouse embryo, Bmp2 together with Bmp7, are expressed in the surface ectoderm overlying the neural folds . Faure et al. showed that Smad-1, which transduces Bmp signaling, is activated in cells of the opened neural tube as well as in the ectodermal cells covering the neural folds in 10-day mouse embryos, i.e. during neurulation of the trunk region . These observations support the idea that bmp16 in teleosts and Bmp2 and -7 in the mouse may serve similar functions during tail neurulation. Interestingly, the zebrafish bmp16 is expressed in partially overlapping domains with notch1b and notch5 in the tail bud mesenchyme that forms the neural rod during segmentation stages . It has been shown that ectopic expression of Bmp4 in the Xenopus tail bud leads to the formation of an extra tail-like outgrowth . The formation of a neural tube in this structure is Notch signaling-dependent . It will therefore be of interest to define the possible genetic interactions that might exist between these signaling pathways during caudal neural keel development in the zebrafish.
Shuffling of expression domains among Bmp2/4/16genes
The bmp16 expression patterns, scrutinized above from the developmental point of view, can also be dissected from the phylogenetic point of view to shed light on how embryonic expressions for this group of genes have been shuffled during vertebrate evolution. As for the swim bladder-associated expression of bmp16 (Figure 6G, H and 6J-L), our observation is the first evidence of involvement of a member of the TGFβ superfamily in its development. While neither zebrafish bmp2a/2b nor bmp4 are currently known to be expressed in the developing swim bladder , a focused analysis on their expression patterns in the developing endoderm at the time of swim bladder formation might be required to determine whether or not one of these genes is involved in this process. This should reveal whether other members of the bmp2/4/16 subfamily that duplicated prior to the evolution of the swim bladder were recruited for its development. Nonetheless, as described above, mouse Bmp4 is expressed in the developing lung, which is a homologous structure to the swim bladder . So far, there is no report of lung-associated expression for Bmp2 or -4 genes in non-mammalian vertebrates. Although this might be a result of recent lineage-specific recruitment of paralogous genes in the morphogenesis of homologous structures, this endodermal expression might have already been possessed by the ancestral gene before the split of Bmp4 and bmp16 genes.
In contrast, Bmp2 and -4 expression in fin/limb buds, heart and tail bud, where zebrafish bmp16 is expressed (Figure 6), are found in a more taxonomically broad set of vertebrates (see  for limb bud expressions for various vertebrates). It is highly likely that bmp16 expression in those tissues was inherited from the ancestral Bmp2/4/16 gene, with some lineage-specific exceptions in which some expression domains are absent (e.g. absence of tail bud expression for zebrafish bmp4).
On the other hand, zebrafish bmp16 expression is absent in many of the tissues where Bmp2 and/or -4 are frequently expressed, such as in the central nervous system, otic vesicle, olfactory placode, neural crest-derived ectomesenchyme (Figure 6). Moreover, the zebrafish bmp16 was not expressed in pre-gastrulation stages in which Bmp2 and -4 are expressed in the ventral mesoderm to regulate dorsoventral axis specification . This implies that after the gene duplication between Bmp2, -4 and -16, bmp16 expression was secondarily lost in those tissues, and was instead recruited for swim bladder development.
Even though in the molecular phylogenetic trees bmp16 is placed more closely to sea lamprey Bmp2/4 homologs than to jawed-vertebrate Bmp2/4 (Figure 4), a direct comparison of expression patterns between them shows that the only expression domain they have in common is the heart (between bmp16 and PmBmp2/4A) . In addition to the heart, PmBmp2/4A is expressed in a wide variety of tissues where Bmp2/4 expression is usually observed in bony vertebrates, while PmBmp2/4B and -C expression is restricted to branchial arches . In exploring the evolutionary history of Bmp2/4/16 expression patterns, this deviation in expression patterns of the sea lamprey homologs prevents us from drawing a clear-cut scenario--it would not be so surprising if early-branching lineages have accumulated a considerable number of secondary changes (see  for an example of cyclostome lineage-specific changes in gene order; see also  for a review). Information for missing early-branching lineages (e.g. chondrichthyes) might provide more clues about whether the variety of expression patterns originates from (1) frequent lineage-specific losses of expression domains from an ancestral template with a full set of domains, or (2) independent gains of similar expression domains in different lineages.
Our scan in a public genome database identified a novel Bmp2/4 relative, bmp16, confined to teleost fish genomes. We suggest that this gene duplicated from Bmp2 and -4 early in vertebrate evolution, and that the absence of the tetrapod ortholog is due to the secondary loss in its basal lineage. This loss might have been permitted by the redundant expression patterns shared with other Bmp subtypes in fin folds, heart, and tail bud, seen in the zebrafish. At the same time, this gene was revealed to be expressed in the swim bladder, which so far has not been shown to express any other Bmp homolog. Our study demonstrates the power of bioinformatic, phylogenetic approaches to newly identify unexplored genes responsible for crucial biological processes. In particular, an importance of making a similar type of efforts for other gene families to document unidentified genes and secondary changes of their expression patterns is underlined.
Retrieval of sequences
Sequences for members of the TGFβ superfamily were retrieved from the Ensembl genome database (version 53; ) and NCBI Protein database, by running Blastp  using mammalian Bmp2 peptide sequences as queries. In parallel, peptide sequences of Branchiostoma floridae were also downloaded from Joint Genome institute (JGI; http://genome.jgi-psf.org/Brafl1/Brafl1.download.ftp.html), and used as a target of a Blastp search. An optimal multiple alignment of the retrieved amino acid sequences and those sequenced in this study was constructed using the alignment editor XCED, in which the MAFFT program is implemented .
Molecular phylogenetic analysis
Molecular phylogenetic trees were inferred using the regions that were unambiguously aligned with no gaps. The neighbor-joining tree reconstruction  was processed using XCED. To calculate the support values for the molecular phylogenetic tree shown in Figures 3 and 4, we used PhyML  and MrBayes 3.1 . Bootstrap probabilities were obtained with 100 resamplings. To investigate the phylogenetic relationship between Bmp2/4/16 and P. marinus genes (Table 2), the maximum-likelihood (ML) tree was inferred using Tree-Puzzle , assuming JTT+I+Γ4 model (shape parameter of the gamma distribution α = 0.47). In this ML analysis, we performed an exhaustive search of the ML tree in the 'user defined trees' mode of Tree-Puzzle with all 945 possible topologies consisting of seven OTUs, namely, gnathostome Bmp2, gnathostome Bmp4, teleost bmp16, P. marinus PmBmp2/4-A, P. marinus PmBmp2/4-B, P. marinus PmBmp2/4-C and outgroup. Phylogenetic relationships within individual OTUs were constrained according to generally accepted phylogenetic relationships of relevant species. Statistical tests for evaluation of alternative tree topologies were performed using CONSEL .
Total RNA was extracted from a whole 52 hpf zebrafish embryo. This RNA was reverse transcribed into cDNA using SuperScript III (Invitrogen) following the instruction of 3'RACE System (Invitrogen). This cDNA was used as a template for PCR with sequence specific primers designed based on the zebrafish bmp16 transcript sequence found in Ensembl (ID: ENSDART00000098426). The first 3'RACE PCR was performed using the forward primer 5'-CAAGAACTGTGTTGTTTCAGCTC-3', and the product of this PCR was used as a template of the nested PCR reaction with forward primer 5'-GCTAGTCAGCTCTGTAAACGGA-3'. In amplifying the 5'-end of the bmp16 cDNA, we designed the reverse primer 5'-CCTCCTGGTCCAGAAAGAGT-3' based on the afore-mentioned 3'-end cDNA sequence, while the forward primer 5'-CTCCCGTTATTCCAGAGATCCA-3' was designed based on expressed sequence tags (ESTs) in NCBI dbEST (accession numbers [NCBI:CK709830], [NCBI:AW421680] and [NCBI:AW305992]) containing the 5' UTR of this gene. The assembled zebrafish bmp16 cDNA sequence is deposited in EMBL Nucleotide Sequence Database under the accession number [EMBL:FN400947].
Genomic DNA was extracted from 24 hpf embryos. Sequence specific primers designed based on the cDNA sequence of zebrafish bmp16 we obtained as above were used to amplify a 527-bp fragment containing a region where we found an incorrect assembly of genomic sequences in Ensembl (zebrafish genome assemblies, version Zv7 and Zv8) (see Additional file 1: Figure S1). This partial genomic sequence for the zebrafish bmp16 is deposited in EMBL Nucleotide Sequence Database under the accession number [EMBL:FN400947].
The bmp16 riboprobes were prepared using the two cDNA fragments (5'- and 3'- fragments) amplified as above, between which we confirmed positive expression signals observed in common. In Figure 6, only results of in situ hybridization using the 3' probe are shown. isl-1 probe synthesis employed the cDNA template previously reported . Whole mount in situ hybridization was performed as previously described , using the riboprobes labeled with digoxigenin-UTP (Roche Applied Science). Hybridization was detected with alkaline phosphate-conjugated anti-digoxigenin antibody followed by incubation with NBT (nitroblue tetrazolium) and BCIP (5-bromo- 4-chloro-3-indolyl phosphate). Stained embryos were examined with a Zeiss Axiophot microscope. Images were processed using Zeiss Axiovision and Adobe Photoshop software.
This study was supported by the Young Scholar Fund, University of Konstanz to S. Kuraku. We thank Nicola Blum, Silke Pittlik, Ursula Topel, and Elke Hespeler for technical support in cDNA cloning, handling of zebrafish embryos and in situ hybridization, and Helen Gunter, Tereza Manousaki and Huan Qiu for insightful discussion and comments to the manuscript. Our acknowledgement extends to the Zebrafish Nomenclature Committee which provided us with suggestions in our gene naming.
- Muffato M, Crollius HR: Paleogenomics in vertebrates, or the recovery of lost genomes from the mist of time. Bioessays. 2008, 30: 122-134. 10.1002/bies.20707.View ArticlePubMedGoogle Scholar
- Carroll SB, Grenier JK, Weatherbee SD: From DNA to diversity: molecular genetics and the evolution of animal design. 2005, Malden, MA: Blackwell Pub, 2Google Scholar
- Hoekstra HE, Coyne JA: The locus of evolution: evo devo and the genetics of adaptation. Evolution. 2007, 61: 995-1016. 10.1111/j.1558-5646.2007.00105.x.View ArticlePubMedGoogle Scholar
- Kuraku S, Meyer A: The evolution and maintenance of Hox gene clusters in vertebrates and the teleost-specific genome duplication. Int J Dev Biol. 2009, 53: 765-773. 10.1387/ijdb.072533km.View ArticlePubMedGoogle Scholar
- Chen D, Zhao M, Mundy GR: Bone morphogenetic proteins. Growth Factors. 2004, 22: 233-241. 10.1080/08977190412331279890.View ArticlePubMedGoogle Scholar
- Wu MY, Hill CS: Tgf-Beta superfamily signaling in embryonic development and homeostasis. Dev Cell. 2009, 16: 329-343. 10.1016/j.devcel.2009.02.012.View ArticlePubMedGoogle Scholar
- Mac Sweeney A, Gil-Parrado S, Vinzenz D, Bernardi A, Hein A, Bodendorf U, Erbel P, Logel C, Gerhartz B: Structural basis for the substrate specificity of bone morphogenetic protein 1/tolloid-like metalloproteases. J Mol Biol. 2008, 384: 228-239. 10.1016/j.jmb.2008.09.029.View ArticlePubMedGoogle Scholar
- Francis PH, Richardson MK, Brickell PM, Tickle C: Bone morphogenetic proteins and a signalling pathway that controls patterning in the developing chick limb. Development. 1994, 120: 209-218.PubMedGoogle Scholar
- Jones CM, Lyons KM, Hogan BL: Involvement of Bone Morphogenetic Protein-4 (BMP-4) and Vgr-1 in morphogenesis and neurogenesis in the mouse. Development. 1991, 111: 531-542.PubMedGoogle Scholar
- Lyons KM, Pelton RW, Hogan BL: Organogenesis and pattern formation in the mouse: RNA distribution patterns suggest a role for bone morphogenetic protein-2A (BMP-2A). Development. 1990, 109: 833-844.PubMedGoogle Scholar
- Padgett RW, Wozney JM, Gelbart WM: Human BMP sequences can confer normal dorsal-ventral patterning in the Drosophila embryo. Proc Natl Acad Sci USA. 1993, 90: 2905-2909. 10.1073/pnas.90.7.2905.PubMed CentralView ArticlePubMedGoogle Scholar
- Mizutani CM, Bier E: EvoD/Vo: the origins of BMP signalling in the neuroectoderm. Nat Rev Genet. 2008, 9: 663-677. 10.1038/nrg2417.PubMed CentralView ArticlePubMedGoogle Scholar
- Abzhanov A, Kuo WP, Hartmann C, Grant BR, Grant PR, Tabin CJ: The calmodulin pathway and evolution of elongated beak morphology in Darwin's finches. Nature. 2006, 442: 563-567. 10.1038/nature04843.View ArticlePubMedGoogle Scholar
- Abzhanov A, Protas M, Grant BR, Grant PR, Tabin CJ: Bmp4 and morphological variation of beaks in Darwin's finches. Science. 2004, 305: 1462-1465. 10.1126/science.1098095.View ArticlePubMedGoogle Scholar
- Fraser GJ, Hulsey CD, Bloomquist RF, Uyesugi K, Manley NR, Streelman JT: An ancient gene network is co-opted for teeth on old and new jaws. PLoS Biol. 2009, 7: e31-10.1371/journal.pbio.1000031.View ArticlePubMedGoogle Scholar
- Shigetani Y, Sugahara F, Kawakami Y, Murakami Y, Hirano S, Kuratani S: Heterotopic shift of epithelial-mesenchymal interactions in vertebrate jaw evolution. Science. 2002, 296: 1316-1319. 10.1126/science.1068310.View ArticlePubMedGoogle Scholar
- Yu M, Wu P, Widelitz RB, Chuong CM: The morphogenesis of feathers. Nature. 2002, 420: 308-312. 10.1038/nature01196.PubMed CentralView ArticlePubMedGoogle Scholar
- Allendorph GP, Vale WW, Choe S: Structure of the ternary signaling complex of a TGF-beta superfamily member. Proc Natl Acad Sci USA. 2006, 103: 7643-7648. 10.1073/pnas.0602558103.PubMed CentralView ArticlePubMedGoogle Scholar
- McCauley DW, Bronner-Fraser M: Conservation and divergence of BMP2/4 genes in the lamprey: expression and phylogenetic analysis suggest a single ancestral vertebrate gene. Evol Dev. 2004, 6: 411-422. 10.1111/j.1525-142X.2004.04054.x.View ArticlePubMedGoogle Scholar
- Kanki JP, Ho RK: The development of the posterior body in zebrafish. Development. 1997, 124: 881-893.PubMedGoogle Scholar
- Biemar F, Argenton F, Schmidtke R, Epperlein S, Peers B, Driever W: Pancreas development in zebrafish: early dispersed appearance of endocrine hormone expressing cells and their convergence to form the definitive islet. Dev Biol. 2001, 230: 189-203. 10.1006/dbio.2000.0103.View ArticlePubMedGoogle Scholar
- Winata CL, Korzh S, Kondrychyn I, Zheng W, Korzh V, Gong Z: Development of zebrafish swimbladder: The requirement of Hedgehog signaling in specification and organization of the three tissue layers. Dev Biol. 2009, 331: 222-236. 10.1016/j.ydbio.2009.04.035.View ArticlePubMedGoogle Scholar
- Meyer A, Peer Van de Y: From 2R to 3R: evidence for a fish-specific genome duplication (FSGD). Bioessays. 2005, 27: 937-945. 10.1002/bies.20293.View ArticlePubMedGoogle Scholar
- Kuraku S: Insights into cyclostome phylogenomics: pre-2R or post-2R. Zool Sci. 2008, 25: 960-968. 10.2108/zsj.25.960.View ArticlePubMedGoogle Scholar
- Kuraku S, Meyer A, Kuratani S: Timing of genome duplications relative to the origin of the vertebrates: did cyclostomes diverge before or after?. Mol Biol Evol. 2009, 26: 47-59. 10.1093/molbev/msn222.View ArticlePubMedGoogle Scholar
- Grus WE, Zhang J: Origin of the genetic components of the vomeronasal system in the common ancestor of all extant vertebrates. Mol Biol Evol. 2009, 26: 407-419. 10.1093/molbev/msn262.PubMed CentralView ArticlePubMedGoogle Scholar
- Kuraku S, Ota KG, Kuratani S: Jawless fishes (Cyclostomata). The Timetree of Life. Edited by: Hedges SB, Kumar S. 2009, New York: Oxford Univ Press, 317-319.Google Scholar
- Kuraku S, Takio Y, Tamura K, Aono H, Meyer A, Kuratani S: Noncanonical role of Hox14 revealed by its expression patterns in lamprey and shark. Proc Natl Acad Sci USA. 2008, 105: 6679-6683. 10.1073/pnas.0710947105.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee KH, Marden JJ, Thompson MS, MacLennan H, Kishimoto Y, Pratt SJ, Schulte-Merker S, Hammerschmidt M, Johnson SL, Postlethwait JH, Beier DC, Zon LI: Cloning and genetic mapping of zebrafish BMP-2. Dev Genet. 1998, 23: 97-103. 10.1002/(SICI)1520-6408(1998)23:2<97::AID-DVG1>3.0.CO;2-0.View ArticlePubMedGoogle Scholar
- Martinez-Barbera JP, Toresson H, Da Rocha S, Krauss S: Cloning and expression of three members of the zebrafish Bmp family: Bmp2a, Bmp2b and Bmp4. Gene. 1997, 198: 53-59. 10.1016/S0378-1119(97)00292-8.View ArticlePubMedGoogle Scholar
- Kishimoto Y, Lee KH, Zon L, Hammerschmidt M, Schulte-Merker S: The molecular nature of zebrafish swirl: BMP2 function is essential during early dorsoventral patterning. Development. 1997, 124: 4457-4466.PubMedGoogle Scholar
- Shawi M, Serluca F: Identification of a BMP7 homolog in zebrafish expressed in developing organ systems. Gene Expr Patterns. 2008, 8: 369-375. 10.1016/j.gep.2008.05.004.View ArticlePubMedGoogle Scholar
- Thisse B, Thisse C: Fast Release Clones: A High Throughput Expression Analysis. ZFIN Direct Data Submission. 2004, [http://zfin.org]Google Scholar
- Niswander L: Pattern formation: old models out on a limb. Nat Rev Genet. 2003, 4: 133-143. 10.1038/nrg1001.View ArticlePubMedGoogle Scholar
- Chen JN, van Eeden FJ, Warren KS, Chin A, Nusslein-Volhard C, Haffter P, Fishman MC: Left-right pattern of cardiac BMP4 may drive asymmetry of the heart in zebrafish. Development. 1997, 124: 4373-4382.PubMedGoogle Scholar
- Chin AJ, Chen J-N, Weinberg ES: Bone morphogenetic protein-4 expression characterizes inductive boundaries in organs of developing zebrafish. Dev Genes Evol. 1997, 207: 107-114. 10.1007/s004270050097.View ArticleGoogle Scholar
- Wang YX, Qian LX, Liu D, Yao LL, Jiang Q, Yu Z, Gui YH, Zhong TP, Song HY: Bone morphogenetic protein-2 acts upstream of myocyte-specific enhancer factor 2a to control embryonic cardiac contractility. Cardiovasc Res. 2007, 74: 290-303. 10.1016/j.cardiores.2007.02.007.View ArticlePubMedGoogle Scholar
- Danesh SM, Villasenor A, D C, Soukup C, Cleaver O: BMP and BMP receptor expression during murine organogenesis. Gene Expr Patterns. 2009, 9: 255-265. 10.1016/j.gep.2009.04.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Dudley AT, Robertson EJ: Overlapping expression domains of bone morphogenetic protein family members potentially account for limited tissue defects in BMP7 deficient embryos. Dev Dyn. 1997, 208: 349-362. 10.1002/(SICI)1097-0177(199703)208:3<349::AID-AJA6>3.0.CO;2-I.View ArticlePubMedGoogle Scholar
- Perry SF, Wilson RJ, Straus C, Harris MB, Remmers JE: Which came first, the lung or the breath?. Comp Biochem Physiol A Mol Integr Physiol. 2001, 129: 37-47. 10.1016/S1095-6433(01)00304-X.View ArticlePubMedGoogle Scholar
- Bellusci S, Henderson R, Winnier G, Oikawa T, Hogan B: Evidence from normal expression and targeted misexpression that bone morphogenetic protein (Bmp-4) plays a role in mouse embryonic lung morphogenesis. Development. 1996, 122: 1693-1702.PubMedGoogle Scholar
- Lu M, Yang H, Zhang L, Shu W, Blair D, Morrisey E: The bone morphogenic protein antagonist gremlin regulates proximal-distal patterning of the lung. Dev Dyn. 2001, 222: 667-680. 10.1002/dvdy.1231.View ArticlePubMedGoogle Scholar
- Weaver M, Yingling J, Dunn N, Bellusci S, Hogan B: Bmp signaling regulates proximal-distal differentiation of endoderm in mouse lung development. Development. 1999, 126: 4005-4015.PubMedGoogle Scholar
- Eblaghie M, Reedy M, Oliver T, Mishina Y, Hogan B: Evidence that autocrine signaling through Bmpr1a regulates the proliferation, survival and morphogenetic behavior of distal lung epithelial cells. Dev Biol. 2006, 291: 67-82. 10.1016/j.ydbio.2005.12.006.View ArticlePubMedGoogle Scholar
- Hsu D, Economides A, Wang X, Eimon P, Harland R: The Xenopus dorsalizing factor Gremlin identifies a novel family of secreted proteins that antagonize BMP activities. Mol Cell. 1998, 1: 673-683. 10.1016/S1097-2765(00)80067-2.View ArticlePubMedGoogle Scholar
- Khokha M, Hsu D, Brunet L, Dionne M, Harland R: Gremlin is the BMP antagonist required for maintenance of Shh and Fgf signals during limb patterning. Nat Genet. 2003, 34: 303-307. 10.1038/ng1178.View ArticlePubMedGoogle Scholar
- Michos O, Panman L, Vintersten K, Beier K, Zeller R, Zuniga A: Gremlin-mediated BMP antagonism induces the epithelial-mesenchymal feedback signaling controlling metanephric kidney and limb organogenesis. Development. 2004, 131: 3401-3410. 10.1242/dev.01251.View ArticlePubMedGoogle Scholar
- Shi W, Zhao J, Anderson K, Warburton D: Gremlin negatively modulates BMP-4 induction of embryonic mouse lung branching morphogenesis. Am J Physiol Lung Cell Mol Physiol. 2001, 280: L1030-1039.PubMedGoogle Scholar
- Topol L, Bardot B, Zhang Q, Resau J, Huillard E, Marx M, Calothy G, Blair D: Biosynthesis, post-translation modification, and functional characterization of Drm/Gremlin. J Biol Chem. 2000, 275: 8785-8793. 10.1074/jbc.275.12.8785.View ArticlePubMedGoogle Scholar
- Chung W, Shin C, Stainier D: Bmp2 signaling regulates the hepatic versus pancreatic fate decision. Dev Cell. 2008, 15: 738-748. 10.1016/j.devcel.2008.08.019.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang H, Ruan H, Aw M, Hussain A, Guo L, Gao C, Qian F, Leung T, Song H, Kimelman D, Wen J, Peng Z: Mypt1-mediated spatial positioning of Bmp2-producing cells is essential for liver organogenesis. Development. 2008, 135: 3209-3218. 10.1242/dev.024406.View ArticlePubMedGoogle Scholar
- Muller I, Knapik E, Hatzopoulos A: Expression of the protein related to Dan and Cerberus gene-prdc-During eye, pharyngeal arch, somite, and swim bladder development in zebrafish. Dev Dyn. 2006, 235: 2881-2888. 10.1002/dvdy.20925.View ArticlePubMedGoogle Scholar
- Sudo S, Avsian-Kretchmer O, Wang L, Hsueh A: Protein related to DAN and cerberus is a bone morphogenetic protein antagonist that participates in ovarian paracrine regulation. J Biol Chem. 2004, 279: 23134-23141. 10.1074/jbc.M402376200.View ArticlePubMedGoogle Scholar
- Schmid B, Furthauer M, Connors SA, Trout J, Thisse B, Thisse C, Mullins MC: Equivalent genetic roles for bmp7/snailhouse and bmp2b/swirl in dorsoventral pattern formation. Development. 2000, 127: 957-967.PubMedGoogle Scholar
- Ybot-Gonzalez P, Gaston-Massuet C, Girdler G, Klingensmith J, Arkell R, Greene N, Copp A: Neural plate morphogenesis during mouse neurulation is regulated by antagonism of Bmp signalling. Development. 2007, 134: 3203-3211. 10.1242/dev.008177.View ArticlePubMedGoogle Scholar
- Faure S, de Santa Barbara P, Roberts D, Whitman M: Endogenous patterns of BMP signaling during early chick development. Dev Biol. 2002, 244: 44-65. 10.1006/dbio.2002.0579.View ArticlePubMedGoogle Scholar
- Westin J, Lardelli M: Three novel Notch genes in zebrafish: implications for vertebrate Notch gene evolution and function. Dev Genes Evol. 1997, 207: 51-63. 10.1007/s004270050091.View ArticlePubMedGoogle Scholar
- Beck CW, Whitman M, Slack JM: The role of BMP signaling in outgrowth and patterning of the Xenopus tail bud. Dev Biol. 2001, 238: 303-314. 10.1006/dbio.2001.0407.View ArticlePubMedGoogle Scholar
- Kondo M: Bone morphogenetic proteins in the early development of zebrafish. FEBS J. 2007, 274: 2960-2967. 10.1111/j.1742-4658.2007.05838.x.View ArticlePubMedGoogle Scholar
- Furlong RF, Younger R, Kasahara M, Reinhardt R, Thorndyke M, Holland PW: A degenerate ParaHox gene cluster in a degenerate vertebrate. Mol Biol Evol. 2007, 24: 2681-2686. 10.1093/molbev/msm194.View ArticlePubMedGoogle Scholar
- Hubbard TJ, Aken BL, Ayling S, Ballester B, Beal K, Bragin E, Brent S, Chen Y, Clapham P, Clarke L, Coates G, Fairley S, Fitzgerald S, Fernandez-Banet J, Gordon L, Graf S, Haider S, Hammond M, Holland R, Howe K, Jenkinson A, Johnson N, Kahari A, Keefe D, Keenan S, Kinsella R, Kokocinski F, Kulesha E, Lawson D, Longden I, et al: Ensembl 2009. Nucleic Acids Res. 2009, 37: D690-697. 10.1093/nar/gkn828.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33: 511-518. 10.1093/nar/gki198.PubMed CentralView ArticlePubMedGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.View ArticlePubMedGoogle Scholar
- Shimodaira H, Hasegawa M: CONSEL: for assessing the confidence of phylogenetic tree selection. Bioinformatics. 2001, 17: 1246-1247. 10.1093/bioinformatics/17.12.1246.View ArticlePubMedGoogle Scholar
- Korzh V, Edlund T, Thor S: Zebrafish primary neurons initiate expression of the LIM homeodomain protein Isl-1 at the end of gastrulation. Development. 1993, 118: 417-425.PubMedGoogle Scholar
- Begemann G, Schilling TF, Rauch GJ, Geisler R, Ingham PW: The zebrafish neckless mutation reveals a requirement for raldh2 in mesodermal signals that pattern the hindbrain. Development. 2001, 128: 3081-3094.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.