Deep phylogenomics of a tandem-repeat galectin regulating appendicular skeletal pattern formation
© The Author(s). 2016
Received: 8 February 2016
Accepted: 26 July 2016
Published: 18 August 2016
A multiscale network of two galectins Galectin-1 (Gal-1) and Galectin-8 (Gal-8) patterns the avian limb skeleton. Among vertebrates with paired appendages, chondrichthyan fins typically have one or more cartilage plates and many repeating parallel endoskeletal elements, actinopterygian fins have more varied patterns of nodules, bars and plates, while tetrapod limbs exhibit tandem arrays of few, proximodistally increasing numbers of elements. We applied a comparative genomic and protein evolution approach to understand the origin of the galectin patterning network. Having previously observed a phylogenetic constraint on Gal-1 structure across vertebrates, we asked whether evolutionary changes of Gal-8 could have critically contributed to the origin of the tetrapod pattern.
Translocations, duplications, and losses of Gal-8 genes in Actinopterygii established them in different genomic locations from those that the Sarcopterygii (including the tetrapods) share with chondrichthyans. The sarcopterygian Gal-8 genes acquired a potentially regulatory non-coding motif and underwent purifying selection. The actinopterygian Gal-8 genes, in contrast, did not acquire the non-coding motif and underwent positive selection.
These observations interpreted through the lens of a reaction-diffusion-adhesion model based on avian experimental findings can account for the distinct endoskeletal patterns of cartilaginous, ray-finned, and lobe-finned fishes, and the stereotypical limb skeletons of tetrapods.
KeywordsGalectin-8 Limb skeleton Pattern formation Mathematical modeling Homology Phylogeny
Galectin-8 (Gal-8), encoded by the gene lgals8, belongs to the family of β-galactoside-binding proteins [1–3]. Alternative splicing of lgals8 results in two protein isoforms [4, 5]: prototype galectins, which contain a single carbohydrate recognition domain, and tandem-repeat galectins, which contain two carbohydrate recognition domains (CRDs) with distinct binding affinities [6–8] and different evolutionary origins . Gal-8 is an important regulator of cell adhesion in adult tissues [10, 11] and is differentially expressed in normal and cancer tissues . During avian embryogenesis, Gal-8 is expressed in the limb bud and mediates the patterning of the precartilage mesenchymal condensations that constitute the primordia of the appendicular skeleton [13, 14]. Specifically, Gal-8 upregulates expression of Gal-1A, the cell adhesive homolog of Gal-1, through a mutually reinforcing feedback loop while also inhibiting cell adhesion by competing with the binding of Gal-1A to its cognate ligand/receptor .
Represented in the form of a mathematical model , these findings suggest that the two galectins participate in a reaction-diffusion-type mechanism [16, 17] of the kind that best integrates the patterning and morphogenesis of skeletal elements during limb skeletal pattern formation [18–20]. Such empirically based models allow for testable hypotheses about the mechanisms that underlie the evolution of endoskeletal diversity in tetrapod appendages. Specifically, they can be used to explore how the modulation of parameters of these patterning networks may have been responsible for differences observed in limb skeletal anatomy between major gnathostome clades [18, 20].
The paired fins and limbs of gnathostomes are characterized by endoskeletal elements (cartilages and the endochondral bones that arise from them) . Gal-1A is the Gal-1 homolog that mediates precartilage condensation formation in the chicken. The Gal-1 s of actinopterygians (ray-finned fishes) resemble Gal-1A more closely in sequence and fold than they do the non-skeletogenic homolog of Gal-1 (Gal-1B), which evolved in the sauropsid lineage (which comprises birds and reptiles) . Furthermore, genes encoding some of the Gal-1 homologs in amphibians and the single gene encoding Gal-1 in mouse specify proteins with the Gal-1A-type fold structure seen in the ray-finned fish and sauropsids. Therefore, a potentially cartilage-inducing Gal-1 homolog is likely to have evolved before the origin of digits and thus was not the key factor responsible for innovating the tetrapod limb skeletal patterning network.
To trace the origin of the tetrapod skeletal patterning network we therefore turned our attention to Gal-8, which in the chicken limb regulates the number and spacing of condensations, not their initiation and morphogenesis [13, 14]. Here we used a combination of phylogenetic methods to compare the evolution of sarcopterygian and actinopterygian Gal-8s relative to their chondrichthyan ortholog. With respect to synteny, selected sequence signatures at the gene level, and residue conservation at the protein level, actinopterygian Gal-8s differ more extensively from their chondrichthyan orthologs than do sarcopterygian Gal-8s. Employing a previously described mathematical model of the galectin-based patterning network for avian limb skeletogenesis , we show how changes in both the regulation of gene expression and in the coding sequence of Gal-8 could have enabled the transformation of a precartilage condensation pattern like that of chondrichthyan fins to one characteristic of tetrapod limbs.
Phylogenetic analysis of Gal-8 protein sequence shows a deep split between Actinopterygii and Sarcopterygii
Actinopterygian and sarcopterygian lgals8s show distinct synteny
To determine whether one of the lgals8 syntenies was ancestral to the other, we mapped the synteny of the lgals8 homolog in elephant shark. Orthologs of two genes actn2b and heatr1 were observed to flank the lgals8 ortholog of elephant shark. Both these genes are also proximal to lgals8 of every sarcopterygian species examined (along with edaradd) but not part of actinopterygian lgals8 syntenies. We then searched for heatr1 in the spotted gar genome and found its ortholog proximal to those for actin2b and edaradd. The proximity in positions of these three genes was absent in teleost genomes. In spotted gar, we observed lgals8 in a new location flanked by kif11, rrm2, tmem242 and arid1b. In teleosts, all of which were observed to have two lgals8 paralogs, the flanking genes for each overlapped partially with their spotted gar paralogs but not with each other. In sarcopterygian genomes, the orthologs of actin2b and heatr1 were found to flank lgals8. Our findings suggest that a single genome transposition of lgals8 took place very early in the actinopterygian lineage, before the divergence of gar and teleosts. This transposition was followed by the teleost genome duplication. Post-duplication genome evolutionary dynamics is probably responsible for the largely non-overlapping syntenies of the actinopterygian lgals8 paralogs. The synteny of lgals8 in the elephant shark was therefore retained in the sarcopterygian lineage.
Sarcopterygian Gal-8 evolved more slowly than actinopterygian Gal-8
dN/dS for different site classes in Gal-8 using branch site test (Zhang et al. 2005)
In the Gal-8 C-CRD, we identified one potential residue under possible positive selection in Actinopterygii: Glu251, which is conserved in sarcopterygians but substituted with Pro, His, Ser or Gln in actinopterygians. We next sought to determine whether sarcopterygian Gal-8s were more similar to shark Gal-8 or to their actinopterygian counterparts. We aligned actinopterygian and sarcopterygian sequences, taking into account the crystallographically elucidated secondary structure of the human Gal-8 CRD  and quantified shared invariant and variant residues specific to each class. We then overlaid both sets with the shark Gal-8 sequence and performed a similar analysis (Additional file 1: Figure S1). We found that sarcopterygian Gal-8s shared a higher percentage of strongly and weakly conserved residues with shark Gal-8 than the latter did with actinopterygian Gal-8s (Additional file 2: Figure S2). Both CRDs of sarcopterygian Gal-8 show higher residue conservation than their actinopterygian counterparts, with N-CRD showing relatively greater conservation. We also found greater conservation of residues in the short sequence that precedes the sarcopterygian Gal-8 N-CRD sequence relative to those of actinopterygians (Additional file 2: Figure S2).
A broad range of Gal-8 structural and expression parameters is consistent with endoskeletal patterning by the Gal-1-Gal-8 network
We had previously used a mathematical representation of the Gal-1-Gal-8 network to identify Gal-8-related parameters whose variations modulate condensation patterning . Patterns capable of being produced by the described mechanism include repeated cartilage elements with regular spacing comparable to the element widths . For the purposes of the present paper, we define a system as “patterned” when it exhibits two or more such elements. One such modulatory network parameter is β: a function of the binding affinity of this galectin to its receptors. This parameter has an obvious relationship to the 3D fold structure, and hence the sequence, of Gal-8. The evidence described above for purifying selection and sequence conservation of sarcopterygian Gal-8s, including regions known to be involved in binding to its carbohydrate ligands, suggests a phylogenetic constraint on β values during tetrapod evolution. This assumes, as suggested by previous evidence , that the relevant folds of Gal-1 have been subject to purifying selection since their origination.
A comparison of the structure of human Gal-8 with those of chicken Gal-1A and Gal-1B (the -1B paralog being a non-skeletogenic protein that arose from the ancestral Gal-1 after its duplication in the sauropsids and divergence from its skeletogenic paralog Gal-1A ), shows a greater similarity between the folds of each of the CRDs of human Gal-8 (the only Gal-8 CRD structures to be experimentally elucidated) and the fold of chicken Gal-1A, than with the fold of chicken Gal-1B (Additional file 4: Figure S3). Gal-8 in Sarcopterygii therefore appears to have evolved under pressure to remain similar in fold to the basal (skeletogenic) isoform of Gal-1. Such purifying selection could have ensured that Gal-8 binding to its shared ligand with Gal-1 was in a range that made it neither negligible nor too avid. The crystallographic elucidation of the tertiary folds of sarcopterygian Gal-8 CRDs in addition to those of human Gal-8 would provide a better understanding of how the structure of Gal-8 evolved in the context of the split between Actinopterygii and Sarcopterygii. As a corollary, our model predicts that a progressive decrease in element number in the face of a phylogenetic constraint on β could take place through a monotonic increase in values of μ (i.e., expression levels of Gal-8). Consistent with this, knocking down the expression of chicken Gal-8 using RNAi in developing chicken limbs led to ectopic digit formation (data not shown).
A conserved non-coding motif upstream of lgals8 is present exclusively in sarcopterygians
In addition to requiring a Gal-8 with the capacity to interfere with cell-cell adhesion mediated by Gal-1, our model predicts that a change in the regulatory regime of the lgals8 that would permit its regulated elevated expression may have been instrumental in the emergence of a tetrapod-type patterning network from an ancestral gnathostome one. We therefore analyzed the non-coding regions upstream of a broad selection of lgals8 orthologs to identify potential evolutionary changes in determinants of Gal-8 expression.
We have provided evidence for significant differences between major gnathostome clades in the evolution of galectin Gal-8 and its specifying gene lgals8. These differences include conservation within Sarcopterygii of the synteny of lgals8 in its shared ancestor with cartilaginous fish, and retention of key residues in sarcopterygian Gal-8s, presumably by purifying selection, after their divergence from the actinopterygians. Furthermore, we also detected the presence of a conserved non-coding motif (CNM) containing binding sites for transcription factors preferentially expressed in embryonic tetrapod limb buds upstream of the sarcopterygian Gal-8 encoding genes.
In light of the skeletogenic two-galectin network we have inferred from investigations of an avian system, we suggest that the new findings presented here on Gal-8, in conjunction with earlier work on the phylogenetics and fold-structure of Gal-1  can provide insight into the evolution and developmental regulation of a sarcopterygian appendicular patterning network in which there is regulated elaboration of small numbers of elements (generally increasing proximodistally). The elaboration is ultimately refined to the stylopod, zeugopod and autopod of tetrapod vertebrates. In contrast, the endoskeletons of the paired fins of cartilaginous and ray-finned fishes exhibit a wide variety of plates, nodules, and multiple parallel bars of cartilage or endochondral bone (Fig. 4). If, as suggested by our analysis, the fins of actinopterygians (for example) evolved unconstrained by a tetrapod-type patterning network (incorporating Gal-8 with certain structural refinements and a CNM upstream of its gene), their adaptive radiation could have followed less stereotypical pathways, reflected in both their divergent fin endoskeletal patterns and positive selection on Gal-8. In contrast, the stereotypical tetrapod pattern persisted in the fins of vertebrates that became secondarily aquatic.
The constraint on the range of protein conformations of Gal-8 that appears to have accompanied the galectin-mediated transformation of the appendicular skeletal pattern was based on purifying selection which, according to our residue conservation analysis, was already underway in crown gnathostomes. In terms of our mathematical model, this evolutionary process corresponds to centering the value of the parameter β in the skeletogenic galectin network around 0.30 in the relative units of the parameter space (Fig. 3). According to the experiments that motivated the model , these values of β should correspond to a Gal-8 conformation that allows it to bind to its common receptor with Gal-1, but not so strongly that it completely displaces the latter. By the geometry of the parameter space these β values are well within the range permissive for pattern formation, but represent the locus within this space in which the number of elements can be precisely controlled developmentally by changes in the expression level of Gal-8.
While the model predicts that a wide range of expression levels of lgals8 will be consistent with pattern formation, only elevated levels of the protein lead to the consistently small numbers of elements characteristic of sarcopterygian fin endoskeletons and tetrapod limbs. This is particularly true for values of β around 0.30 units (Fig. 3), although there is nothing in our analysis that precludes individual chondrichthyan and actinopterygian species from having acquired β and μ values that are compatible with small numbers of elements, as is sometimes seen in the skeletal anatomies of these groups [29, 30]. If our interpretation is correct that the purifying selection on Gal-8 in the sarcopterygian lineage (Table 1) preserved a value of β in the patterning network for which the expression levels of the protein can directly calibrate the number of condensations, then the appearance of a conserved non-coding motif (CNM) with multiple potential transcription factor binding sites immediately upstream the promoter of the sarcopterygian lgals8 gene becomes of great interest. While it is not presently known which factors (apart from Gal-1A ) in the developing tetrapod limb regulate the production of Gal-8, evidence that this protein came under a novel regulatory regime in sarcopterygians supports the general outline of our model.
In addition to having a small number of discrete skeletal elements, the paired appendages of the vast majority of tetrapods (and some lobe-finned fish) exhibit a stereotypical proximodistal increase in the number of parallel skeletal elements. This arrangement, as well as the proximodistal order of their generation, are predictable consequences (based on reaction-diffusion schemes like the two-galectin one discussed here) of the distal suppression of precartilage condensation by the fibroblast growth factor-8 (FGF8) produced by the ectoderm at the limb bud tip [18, 31, 32]. Indeed, we have observed the in vitro patterning function of this network to be markedly inhibited by FGF8 .
In a recent publication, Hoxa and Hoxd enhancers that specify digit and wrist identity in murine limbs were found to be utilized in the distal radial regions of pectoral fins of an actinopterygian fish . Rather than taking this, as suggested, as evidence of homology of these elements across these distant groups , we are led by our results to consider the associated Hox proteins as modulators of an ancestral skeletogenic system that only took on a specific digit-related role after the sarcopterygian two-galectin network was in place [18–20].
Evolutionary changes in the gnathostome genes specifying the animal lectin Gal-8 at the protein and noncoding DNA levels, when analyzed in terms of an experimentally based mathematical model, suggest a phylogenetic route of emergence of a reaction-diffusion network capable of generating sarcopterygian-type limb skeletal patterns. In the context of the structural conservation of Gal-1, the skeletogenic component of the network, across the vertebrates, the structural and implied functional differences in Gal-8, the pattern regulatory component, among the cartilaginous, ray-finned and lobe-finned fishes/tetrapods, support a galectin-based evolutionary-developmental hypothesis for the fin-limb transition.
Protein and nucleic acid sequence search
Peptide sequences of Gal-8 and its homologs were retrieved from Ensembl (http://www.ensembl.org; Release 78), NCBI (http://www.ncbi.nlm.nih.gov/protein/) and Xenbase (http://www.xenbase.org; Version 3.0; Xenopus tropicalis v7.1 and Xenopus laevis v7.1). We used Basic Local Alignment Search Tool (BLAST)/BLAT (BLAST-like Alignment Tool) algorithm to identify vertebrate Gal-8 sequences using chicken Gal-8 (ENSGALT00000006738) peptide sequence as input. Nucleotide sequences of lgals8 genes were accessed from Ensembl (http://www.ensembl.org; Release 78) and verified by translating them using ExPASY Translate tool . Gene and protein sequences used in the study have been deposited in the Dryad repository .
Sequence alignment and phylogenetic tree construction
A rapid inference of peptide sequence phylogeny was carried out by aligning them using MUSCLE (MUltiple Sequence Comparison by Log-Expectation)  followed by tree construction using neighborhood joining method (BIONJ; Poisson distribution)  (with bootstrap analysis: 100000 replicates) using SeaView (V4.5.2) phylogenetic analysis software . This was followed by a tree construction using maximum likelihood method with PhyML . We used LG, a model of amino acid replacement matrix with improved performance over other models such as JTT and Whelan and Goldman, and optimized for both invariant sites and across-the-tree variation in rate of evolution. Posterior branch support was computed using both approximate Likelihood Ratio test (aLRT)  and bootstrap analysis (with 100 replicates). The tree searching operation was set to Nearest-Neighbor Interchange.
The location of lgals8 genes were ascertained using Ensembl and the Genomicus Browser  (http://www.genomicus.biologie.ens.fr/genomicus-78.01/cgi-bin/search.pl; version 78.01) was used to obtain a simple visual representation of gene syntenies. For the chromosome-level analysis, we used the dotplots option from the Synteny Database  to compare the spatial maps of chromosome 17 and 20 of Danio rerio with all chromosomes of Mus musculus (http://teleost.cs.uoregon.edu/dotplots/; Ensembl version 70).
Analysis of conservation of residues
Primary structures (amino acid sequences) of sarcopterygian Gal-8s, and actinopterygian Gal-8s were aligned in separate subsets using MUSCLE and overlaid with secondary structure of Gal-8 (locations of β-strands [S1–S6b, F1–F5]). The pre N terminal CRD (pre-N-CRD) N-terminal CRD (N-CRD) and C-terminal CRD (C-CRD) domains were demarcated. The percentage of conserved (identical amino acids and amino acids with strong or weak similar properties) and the percentage of strongly conserved (identical amino acids) were ascertained for the whole sequence as well as for individual domains.
Fold prediction and comparative analysis
The PDB files for H. sapiens Gal-8 C-CRD (3OJB, 2YRO), H. sapiens Gal-8 N-CRD (2YV8, 3BMB, 3VKN), G. gallus Gal-1A (1QMJ) and Gal-1B (3DUI) were retrieved from the RCSB Protein Data Bank  (http://www.rcsb.org/pdb/home/home.do, last accessed February 5, 2015). Each Gal-8 CRD PDB was compared with the tertiary folds of chicken Gal-1A and chicken Gal-1B, using PDBeFold (http://www.ebi.ac.uk/msd-srv/ssm/, last accessed on February 5, 2015) , which uses the Secondary Structure Matching algorithm to achieve the best Cα alignment of amino acids. The metric used for comparing topological similarity was Q score, which takes into account Nalign (the maximum number of aligned residues) as well as a measure of the distance between the Cα atoms of the matched residues (RMSD) when the target and query sequences are superposed in 3D. Q scores from alignment comparisons between two crystal structures were computed using “A” chain identifiers of both PDB files.
Test for rate of protein evolution
The PAML4 package was used to assess the of clade-specific Gal-8 evolution by quantifying the rate of non-synonymous substitutions of lgals8 . Guided by the amino acid alignment, the codons of the genes were aligned in TranslatorX . The free ratios model was used to calculate the maximum likelihood estimate of non-synonymous substitution (dN) for each branch of the given tree. The null model assumes that all sites are evolving under stochastic forces or under purifying selection. If there is an increase in the substitution rate for reasons other than selection, the likelihood ratio test will not reject the null hypothesis.
Non-coding DNA motif search
The near-promoter regions (2000 bp upstream of start codon) of sarcopterygian lgals8 were searched for conserved non-coding motifs (CNMs) (sequences deposited in Dryad repository ). The list of near-promoter regions were used as input for MEME Suite  (http://meme.nbcr.net/meme/; Version 4.9.0), which represents motifs as position-dependent probability matrices and uses an Expectation Maximization algorithm to fit a two-component finite mixture model to the set, searching for motifs of given length and number of occurrences (minimum width = 2, maximum width = 300, any number of occurrences). The P-value (statistical significance for the presence of the motif) was measured for each sarcopterygian. We then used MAST  (http://meme.nbcr.net/meme/cgi-bin/mast.cgi) to search for candidate motifs within sarcopterygian lgals8 near-promoter regions. TOMTOM  (http://meme.nbcr.net/meme/cgi-bin/tomtom.cgi), was used to search the CNM against the JASPAR Core Vertebrate database  (http://jaspar.genereg.net/) which contains motifs that are known binding sites of transcription factors (TFBSs). Possible limb-tissue expression of selected transcription factors was searched using the now dysfunct EMBRYS database (http://embrys.jp/embrys/html/MainMenu.html and have been since been confirmed to show similar results using the Mouse Genome informatics databasehttp://www.informatics.jax.org/
aLRT, approximate likelihood ratio test; CNM, conserved non-coding motif; CRD, carbohydrate recognition domain; Gal-1, Galectin-1; Gal-8, Galectin-8; MUSCLE, Multiple sequence comparison by log-expectation; PAML, phylogenetic analysis by maximum likelihood; RMSD, root mean square deviation; TFBS, transcription factor binding sites.
We thank Hans Larsson and an anonymous reviewer for helping ensure the rigor and clarity of our presentation.
Availability of data and materials
The datasets supporting the study are available in the Dryad repository (doi:10.5061/dryad.n5095) (http://datadryad.org/review?doi=doi:10.5061/dryad.n5095).
RB designed and performed the phylogenetic analyses. TG designed the in silico experiments and performed the simulations. MC performed selection analysis. RB, TG, MC, TS and SAN interpreted the results and wrote the manuscript. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interest.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Varki A. Essentials of glycobiology. 2nd ed. Cold Spring Harbor: Cold Spring Harbor Laboratory Press; 2009.Google Scholar
- Cooper DN. Galectinomics: finding themes in complexity. Biochim Biophys Acta. 2002;1572:209–31.View ArticlePubMedGoogle Scholar
- Gabius HJ. The sugar code : fundamentals of glycosciences. Weinheim Chichester: Wiley-VCH; John Wiley distributor; 2009.Google Scholar
- Bidon N, Brichory F, Bourguet P, Le Pennec JP, Dazord L. Galectin-8: a complex sub-family of galectins (Review). Int J Mol Med. 2001;8:245–50.PubMedGoogle Scholar
- Lahm H, Andre S, Hoeflich A, Kaltner H, Siebert HC, Sordat B, von der Lieth CW, Wolf E, Gabius HJ. Tumor galectinology: insights into the complex network of a family of endogenous lectins. Glycoconj J. 2004;20:227–38.View ArticlePubMedGoogle Scholar
- Stowell SR, Arthur CM, Slanina KA, Horton JR, Smith DF, Cummings RD. Dimeric Galectin-8 induces phosphatidylserine exposure in leukocytes through polylactosamine recognition by the C-terminal domain. J Biol Chem. 2008;283:20547–59.View ArticlePubMedPubMed CentralGoogle Scholar
- Carlsson S, Oberg CT, Carlsson MC, Sundin A, Nilsson UJ, Smith D, Cummings RD, Almkvist J, Karlsson A, Leffler H. Affinity of galectin-8 and its carbohydrate recognition domains for ligands in solution and at the cell surface. Glycobiology. 2007;17:663–76.View ArticlePubMedGoogle Scholar
- Kumar S, Frank M, Schwartz-Albiez R. Understanding the specificity of human Galectin-8C domain interactions with its glycan ligands based on molecular dynamics simulations. PLoS One. 2013;8:e59761.View ArticlePubMedPubMed CentralGoogle Scholar
- Houzelstein D, Goncalves IR, Fadden AJ, Sidhu SS, Cooper DN, Drickamer K, Leffler H, Poirier F. Phylogenetic analysis of the vertebrate galectin family. Mol Biol Evol. 2004;21:1177–87.View ArticlePubMedGoogle Scholar
- Levy Y, Arbel-Goren R, Hadari YR, Eshhar S, Ronen D, Elhanany E, Geiger B, Zick Y. Galectin-8 functions as a matricellular modulator of cell adhesion. J Biol Chem. 2001;276:31285–95.View ArticlePubMedGoogle Scholar
- Hadari YR, Arbel-Goren R, Levy Y, Amsterdam A, Alon R, Zakut R, Zick Y. Galectin-8 binding to integrins inhibits cell adhesion and induces apoptosis. J Cell Sci. 2000;113(Pt 13):2385–97.PubMedGoogle Scholar
- Zick Y, Eisenstein M, Goren RA, Hadari YR, Levy Y, Ronen D. Role of galectin-8 as a modulator of cell adhesion and cell growth. Glycoconj J. 2004;19:517–26.View ArticlePubMedGoogle Scholar
- Bhat R, Lerea KM, Peng H, Kaltner H, Gabius HJ, Newman SA. A regulatory network of two galectins mediates the earliest steps of avian limb skeletal morphogenesis. BMC Dev Biol. 2011;11:6.View ArticlePubMedPubMed CentralGoogle Scholar
- Lorda-Diez CI, Montero JA, Diaz-Mendoza MJ, Garcia-Porrero JA, Hurle JM. Defining the earliest transcriptional steps of chondrogenic progenitor specification during the formation of the digits in the embryonic limb. PLoS One. 2011;6:e24546.View ArticlePubMedPubMed CentralGoogle Scholar
- Glimm T, Bhat R, Newman SA. Modeling the morphodynamic galectin patterning network of the developing avian limb skeleton. J Theor Biol. 2014;346:86–108.View ArticlePubMedGoogle Scholar
- Turing AM. The chemical basis of morphogenesis. Philos Trans R Soc Lond B Biol Sci. 1952;237:37–72.View ArticleGoogle Scholar
- Kondo S, Miura T. Reaction–diffusion model as a framework for understanding biological pattern formation. Science. 2010;329:1616–20.View ArticlePubMedGoogle Scholar
- Zhu J, Zhang YT, Alber MS, Newman SA. Bare bones pattern formation: a core regulatory network in varying geometries reproduces major features of vertebrate limb development and evolution. PLoS One. 2010;5:e10892.View ArticlePubMedPubMed CentralGoogle Scholar
- Sheth R, Marcon L, Bastida MF, Junco M, Quintana L, Dahn R, Kmita M, Sharpe J, Ros MA. Hox genes regulate digit patterning by controlling the wavelength of a Turing-type mechanism. Science. 2012;338:1476–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Raspopovic J, Marcon L, Russo L, Sharpe J. Modeling digits. Digit patterning is controlled by a Bmp-Sox9-Wnt Turing network modulated by morphogen gradients. Science. 2014;345:566–70.View ArticlePubMedGoogle Scholar
- Janvier P. Early vertebrates. Oxford; New York: Clarendon Press; Oxford University Press; 1996.Google Scholar
- Bhat R, Chakraborty M, Mian IS, Newman SA. Structural divergence in vertebrate phylogeny of a duplicated prototype galectin. Genome Biol Evol. 2014;6:2721–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Hahn MW. Distinguishing among evolutionary models for the maintenance of gene duplicates. J Hered. 2009;100:605–17.View ArticlePubMedGoogle Scholar
- Ideo H, Matsuzaka T, Nonaka T, Seko A, Yamashita K. Galectin-8-N-domain recognition mechanism for sialylated and sulfated glycans. J Biol Chem. 2011;286:11346–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Yoshida H, Yamashita S, Teraoka M, Itoh A, Nakakita S, Nishi N, Kamitori S. X-ray structure of a protease-resistant mutant form of human galectin-8 with two carbohydrate recognition domains. FEBS J. 2012;279:3937–51.View ArticlePubMedGoogle Scholar
- Hentschel HG, Glimm T, Glazier JA, Newman SA. Dynamical mechanisms for skeletal pattern formation in the vertebrate limb. Proc Biol Sci. 2004;271:1713–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Mercader N, Leonardo E, Azpiazu N, Serrano A, Morata G, Martinez C, Torres M. Conserved regulation of proximodistal limb axis development by Meis1/Hth. Nature. 1999;402:425–9.View ArticlePubMedGoogle Scholar
- Yu L, Liu H, Yan M, Yang J, Long F, Muneoka K, Chen Y. Shox2 is required for chondrocyte proliferation and maturation in proximal limb skeleton. Dev Biol. 2007;306:549–59.View ArticlePubMedPubMed CentralGoogle Scholar
- Zangerl R, Case GR. Iniopterygia: a new order of Chondrichthyan fishes from the Pennsylvanian of North America. In: Fieldiana Geology Memoirs, vol. 6. Chicago: Field Museum of Natural History; 1973.Google Scholar
- Shubin NH, Daeschler EB, Jenkins Jr FA. The pectoral fin of Tiktaalik roseae and the origin of the tetrapod limb. Nature. 2006;440:764–71.View ArticlePubMedGoogle Scholar
- Moftah MZ, Downie SA, Bronstein NB, Mezentseva N, Pu J, Maher PA, Newman SA. Ectodermal FGFs induce perinodular inhibition of limb chondrogenesis in vitro and in vivo via FGF receptor 2. Dev Biol. 2002;249:270-282.View ArticlePubMedGoogle Scholar
- Mercader N, Leonardo E, Piedra ME, Martinez AC, Ros MA, Torres M. Opposing RA and FGF signals control proximodistal vertebrate limb development through regulation of Meis genes. Development. 2000;127:3961–70.PubMedGoogle Scholar
- Bhat R. Molecular and dynamical mediators of avian limb pattern formation: unpublished doctoral dissertation. New York Medical College, Valhalla, NY; 2010.Google Scholar
- Gehrke AR, Schneider I, de la Calle-Mustienes E, Tena JJ, Gomez-Marin C, Chandran M, Nakamura T, Braasch I, Postlethwait JH, Gomez-Skarmeta JL, Shubin NH. Deep conservation of wrist and digit enhancers in fish. Proc Natl Acad Sci U S A. 2015;112:803–8.View ArticlePubMedGoogle Scholar
- Ahn D, Ho RK. Tri-phasic expression of posterior Hox genes during development of pectoral fins in zebrafish: implications for the evolution of vertebrate paired appendages. Dev Biol. 2008;322:220–33.View ArticlePubMedGoogle Scholar
- Stockinger H, Altenhoff AM, Arnold K, Bairoch A, Bastian F, Bergmann S, Bougueleret L, Bucher P, Delorenzi M, Lane L, et al. Fifteen years SIB Swiss Institute of Bioinformatics: life science databases, tools and support. Nucleic Acids Res. 2014;42:W436–441.View ArticlePubMedPubMed CentralGoogle Scholar
- Bhat R, Chakraborty M, Glimm T, Stewart T, Newman SA. Data from: Deep phylogenomics of a tandem-repeat galectin regulating appendicular skeletal pattern formation. In: Dryad digital repository. 2016.Google Scholar
- Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113.View ArticlePubMedPubMed CentralGoogle Scholar
- Gascuel O. BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol Biol Evol. 1997;14:685–95.View ArticlePubMedGoogle Scholar
- Gouy M, Guindon S, Gascuel O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27:221–4.View ArticlePubMedGoogle Scholar
- Le SQ, Gascuel O. An improved general amino acid replacement matrix. Mol Biol Evol. 2008;25:1307–20.View ArticlePubMedGoogle Scholar
- Anisimova M, Gascuel O. Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol. 2006;55:539–52.View ArticlePubMedGoogle Scholar
- Muffato M, Louis A, Poisnel CE, Roest Crollius H. Genomicus: a database and a browser to study gene synteny in modern and ancestral genomes. Bioinformatics. 2010;26:1119–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Catchen JM, Conery JS, Postlethwait JH. Automated identification of conserved synteny after whole-genome duplication. Genome Res. 2009;19:1497–505.View ArticlePubMedPubMed CentralGoogle Scholar
- Rose PW, Bi C, Bluhm WF, Christie CH, Dimitropoulos D, Dutta S, Green RK, Goodsell DS, Prlic A, Quesada M, et al. The RCSB Protein Data Bank: new resources for research and education. Nucleic Acids Res. 2013;41:D475–482.View ArticlePubMedGoogle Scholar
- Krissinel E, Henrick K. Secondary-structure matching (SSM), a new tool for fast protein structure alignment in three dimensions. Acta Crystallogr D Biol Crystallogr. 2004;60:2256–68.View ArticlePubMedGoogle Scholar
- Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.View ArticlePubMedGoogle Scholar
- Abascal F, Zardoya R, Telford MJ. TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Res. 2010;38:W7–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Bailey TL, Elkan C. Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994;2:28–36.PubMedGoogle Scholar
- Bailey TL, Gribskov M. Combining evidence using p-values: application to sequence homology searches. Bioinformatics. 1998;14:48–54.View ArticlePubMedGoogle Scholar
- Gupta S, Stamatoyannopoulos JA, Bailey TL, Noble WS. Quantifying similarity between motifs. Genome Biol. 2007;8:R24.View ArticlePubMedPubMed CentralGoogle Scholar
- Wasserman WW, Sandelin A. Applied bioinformatics for the identification of regulatory elements. Nat Rev Genet. 2004;5:276–87.View ArticlePubMedGoogle Scholar
- Tulenko FJ, Augustus GJ, Massey JL, Sims SE, Mazan S, Davis MC. HoxD expression in the fin-fold compartment of basal gnathostomes and implications for paired appendage evolution. Sci Rep. 2016;6:22720.View ArticlePubMedPubMed CentralGoogle Scholar
- Silva JP, Carvalho MR. Systematics and morphology of Potamotrygon orbignyi (Castelnau, 1855) and allied forms (Chondrichthyes: Myliobatiformes: Potamotrygonidae). Zootaxa. 2015;3982:1–82.View ArticlePubMedGoogle Scholar
- Johanson Z, Joss J, Boisvert CA, Ericsson R, Sutija M, Ahlberg PE. Fish fingers: digit homologues in sarcopterygian fish fins. J Exp Zool B Mol Dev Evol. 2007;308:757–68.View ArticlePubMedGoogle Scholar
- Mellor LCL, Torre J, Brownell RL. Paedomorphic ossification in porpoises with an emphasis on the vaquita (Phocoena sinus). Aquat Mamm. 2009;35:193–202.View ArticleGoogle Scholar
- Cuervo R, Hernandez-Martinez R, Chimal-Monroy J, Merchant-Larios H, Covarrubias L. Full regeneration of the tribasal Polypterus fin. Proc Natl Acad Sci U S A. 2012;109:3838–43.View ArticlePubMedPubMed CentralGoogle Scholar