Biology:Limb development

From HandWiki
Development of the limbs
6 weeks pregnant.png
Illustration of a human embryo at six weeks gestational age
9-Week Human Embryo from Ectopic Pregnancy.jpg
9-week human fetus from ectopic pregnancy
Anatomical terminology

Limb development in vertebrates is an area of active research in both developmental and evolutionary biology, with much of the latter work focused on the transition from fin to limb.[1]

Limb formation begins in the morphogenetic limb field, as mesenchymal cells from the lateral plate mesoderm proliferate to the point that they cause the ectoderm above to bulge out, forming a limb bud. Fibroblast growth factor (FGF) induces the formation of an organizer at the end of the limb bud, called the apical ectodermal ridge (AER), which guides further development and controls cell death. Programmed cell death is necessary to eliminate webbing between digits.

The limb field is a region specified by expression of certain Hox genes, a subset of homeotic genes, and T-box transcription factors – Tbx5 for forelimb or wing development, and Tbx4 for leg or hindlimb development. Establishment of the forelimb field (but not hindlimb field) requires retinoic acid signaling in the developing trunk of the embryo from which the limb buds emerge.[2][3] Also, although excess retinoic acid can alter limb patterning by ectopically activating Shh or Meis1/Meis2 expression, genetic studies in mouse that eliminate retinoic acid synthesis have shown that RA is not required for limb patterning.[4]

The limb bud remains active throughout much of limb development as it stimulates the creation and positive feedback retention of two signaling regions: the AER and its subsequent creation of the zone of polarizing activity (ZPA) with the mesenchymal cells.[5] In addition to the dorsal-ventral axis created by the ectodermal expression of competitive Wnt7a and BMP signals respectively, these AER and ZPA signaling centers are crucial to the proper formation of a limb that is correctly oriented with its corresponding axial polarity in the developing organism.[6][7] Because these signaling systems reciprocally sustain each other's activity, limb development is essentially autonomous after these signaling regions have been established.[5]

Limb formation

Limb bud

Main page: Biology:Limb bud

Limb formation begins in the morphogenetic limb field. Limb formation results from a series of reciprocal tissue interactions between the mesenchyme of the lateral plate mesoderm and the overlying ectodermally derived epithelial cells. Cells from the lateral plate mesoderm and the myotome migrate to the limb field and proliferate to the point that they cause the ectoderm above to bulge out, forming the limb bud. The lateral plate cells produce the cartilaginous and skeletal portions of the limb while the myotome cells produce the muscle components.

The lateral plate mesodermal cells secrete fibroblast growth factors (FGF7 and FGF10) to induce the overlying ectoderm to form an organizer at the end of the limb bud, called the apical ectodermal ridge (AER), which guides further development and controls cell death.[8] The AER secretes further growth factors FGF8 and FGF4 which maintain the FGF10 signal and induce proliferation in the mesoderm.[citation needed] The position of FGF10 expression is regulated by two Wnt signaling pathways: Wnt8c in the hindlimb and Wnt2b in the forelimb. The forelimb and the hindlimb are specified by their position along the anterior/posterior axis and possibly by two transcription factors: Tbx5 and Tbx4, respectively.[9][10]

Precartilage condensations

The limb's skeletal elements are prefigured by tight aggregates known as cellular condensations of the pre-cartilage mesenchymal cells.[11] Mesenchymal condensation is mediated by extracellular matrix and cell adhesion molecules.[12] In the process of chondrogenesis, chondrocytes differentiate from the condensations to form cartilage, giving rise to the skeletal primordia. In the development of most vertebrate limbs (though not in some amphibians), the cartilage skeleton is replaced by bone later in development.

Periodicities of the limb pattern

Vertebrate limbs are organized into stylopod, zeugopod, and autopod.
Vertebrate limbs are organized into stylopod, zeugopod, and autopod.

The limb is organized into three regions: stylopod, zeugopod, and autopod (in order from proximal to distal). The zeugopod and the autopod contain a number of periodic and quasi-periodic pattern motifs. The zeugopod consists of two parallel elements along the anteroposterior axis and the autopod contains 3-5 (in most cases) elements along the same axis. The digits also have a quasi-periodic arrangement along the proximodistal axis, consisting of tandem chains of skeletal elements. The generation of the basic limb plan during development results from the patterning of the mesenchyme by an interplay of factors that promote precartilage condensation and factors that inhibit it.[13]

The development of the basic limb plan is accompanied by the generation of local differences between the elements. For example, the radius and ulna of the forelimb, and the tibia and fibula of the hindlimb of the zeugopod are distinct from one another, as are the different fingers or toes in the autopod. These differences can be treated schematically by considering how they are reflected in each of the limb's three main axes.

A general consensus is that the patterning of the limb skeleton involves one or more Turing-type reaction–diffusion mechanisms.[1]

Evolution and development

The evolution of limbs from paired fins has been an area of much focus.[1] There have been many studies aimed at elucidating specific genes and transcription factors that are found responsible for limb development (See Table 1 below). In particular, studies have shown that SHH, DACH1, and the BMP, HOX, T-box, FGF, and WNT gene families all play a critical role in successful limb development and formation.[14][15] In order to study the genes involved in limb development (and thus evolution) limb reduction and limb loss in snakes is a complementary approach.[16] Conserved sequences involved in limb development are retained in the genomes of snakes. Certain limb-enhancer sequences are also conserved between different types of appendage, such as limbs the phallus.[16][17] For instance, limb-development signalling plays a role both in the development of the limbs and of the genital tubercle in mice.[16][17] The study of limb reduction and limb loss is unravelling the genetic pathways that control limb development.[16] The Turing system has enabled a number of possible outcomes in the evolutionary steps of patterning networks.[1]

Table 1: Various genes known to be responsible for limb development (separated by gene family)
Bmp2, Bmp4, Bmp7
Dach1
En1
Fgf4, Fgf8, Fgf9, Fgf10, Fgf17, Fgfr1
Gli3
Gremlin1
Hand2
Hoxa13, Hoxd13
Msx1, Msx2
Pitx1
Shh
Tbx4, Tbx5
Wnt3, Wnt5a

Many of the genes listed in Table 1 play an important role in embryonic development, specifically during skeletal patterning and limb bud formation.[18] The Shh gene, and genes belonging to the BMP, Hox, T-box, FGF, and Wnt families, all play a pivotal role in cell signaling and differentiation to regulate and promote successful limb formation. Various other genes listed above, one example being Dach1, are DNA-binding proteins that regulate gene expression. The intricate combination of gene expression, regulation, activation, and de-activation allows these genes to produce limbs during embryonic development. Interestingly, many of these genes remain present even in animals that do not have limbs, such as snakes.

Snake evolution and limb loss

Depiction of the spatiotemporal regulation of the Shh gene in coordinating vertebrate digits (fingers) [19]

An interesting aspect in understanding limb development is addressing the question of how snakes lost their legs. Snakes are a particularly good example for studying limb loss, as they underwent limb loss and regeneration multiple times throughout their evolution before they finally lost their legs for good. Much of the gene expression during embryonic development is regulated via spatiotemporal and chemotactic signaling,[20] as depicted by the image to the right. Recent evidence suggests that the highly conserved genes responsible for limb development (Table 1) still remain present in limbless vertebrates,[21] indicating that during embryonic development, the production of limbs, or lack thereof, may best be explained by gene regulation.

Possible role of Shh enhancer in snake limb loss

One theory suggests that the degradation of enhancer sequences may have contributed to the progressive limb loss in snake evolution. In particular, many studies have focused on the ZPA Regulatory Sequence (ZRS) - the enhancer of the Sonic Hedgehog gene (Shh). This long-range enhancer is required for proper limb formation in several vertebrate species, with mutations in this sequence known to cause limb deformities.[17] As such, this sequence is highly conserved across a variety of vertebrate species.

Changes in the Shh enhancer ZRS (~800bp) throughout snake evolution contributed to progressive loss of enhancer function and thus limb development.[22]

Comparisons of the core ZRS in several snake species to the mouse and lizard sequences shows the presence of snake-specific alterations.[23] The core ZRS proved to be mainly conserved in basal snakes like the boa constrictor and python, which still contain pelvic girdle bones.[22] In contrast, advanced snakes such as the viper and cobra, in which no skeletal limb structures remain, have a much higher rate of nucleotide changes when compared to the mouse and lizard ZRS.[22] It is thought that these cumulative changes in the snake ZRS are indicative of a progressive loss of function in this enhancer throughout snake evolution.[22]

Alignment of partial ZPA Regulatory Sequence (ZRS) in vertebrates show increased substitutions in advanced snakes compared to limbed vertebrates and earlier basal snakes. Genomes from the UCSC Genome Browser and GigaDB, and orthologous ZRS enhancer sequences were compared by BLAST. After Kvon et al.[22]

Further investigation into these changes showed an increased rate of substitution in binding sites for transcription factors such as ETS1, whose binding to ZRS has been shown to activate Shh transcription.[24] This degradation in ZRS suggests that this enhancer may be important in further exploring the molecular mechanisms that propelled the morphological evolution of snakes.

Current conclusions and limitations

Snakes are not a common model organism, i.e. they are not easily genetically tractable. In addition, their genome sequence data is incomplete and suffers from poor annotation and quality. These factors make it difficult to understand the mechanism of snake limb loss using a genetic approach, targeting and observing the presence and activity of these genes and their regulatory enhancers. Many of the genes necessary for limb formation are still retained in snakes, hence limb loss can probably not be explained by gene loss.[citation needed]

Axial patterning and related issues

The developing limb has to align itself in relation to three axes of symmetry.[25] These are the craniocaudal (head to tail), dorsoventral (back to front), and proximodistal (near to far) axes.[25]

Many investigations into the development of the limb skeletal pattern have been influenced by the positional information concept proposed by Lewis Wolpert in 1971.[26] In tune with this idea, efforts have been made to identify diffusive signaling molecules (morphogens) that traverse orthogonal axes of developing limbs and determine locations and identities of skeletal elements in a concentration-dependent fashion.

Proximodistal patterning

Hox genes contribute to the specification of the stylopod, zeugopod and autopod. Mutations in Hox genes lead to proximodistal losses or abnormalities.[27] Three different models have been advanced for explaining the patterning of these regions.

Progress zone model

The apical ectodermal ridge (AER) creates and maintains a zone of cell proliferation known as the progress zone.[28] It is thought that cells here gain the positional information they need to travel to their destined position.[28] It was proposed that their positional value was determined by the length of time that the cells were in the progress zone but this has yet to be proved (as of 2001).[28] Proximal structures were proposed to be formed by the first cells to leave the zone and distal ones, by cells that left later.[28]

The Progress Zone model was proposed 30 years ago but recent evidence has conflicted with this model.[29]

Experimental evidence:

  • Removing the AER at a later period of development results in less disruption of distal structures than if the AER was removed early in development.
  • Grafting an early limb bud tip onto a late wing results in duplication of structures, while grafting a late wing bud tip onto an early limb results in a deletion of structures.

Early allocation and progenitor expansion model (or prespecification model)

Cells are specified for each segment in the early limb bud and this population of cells expand out as the limb bud grows. This model is consistent with the following observations. Cell division is seen throughout the limb bud. Cell death occurs within a 200 μm zone subjacent to the AER when it is removed; cell death forecloses some patterning. FGF-releasing beads are able to rescue limb development when the AER is removed by preventing this cell death.

Experimental evidence:

  • Labeled cells in different position of an early limb bud were restricted to single segments of the limb.[30]
  • Limbs lacking expression of required FGF4 & FGF8 showed all structures of the limb and not just the proximal parts.[31]

More recently, however, the investigators primarily responsible for both the Progress Zone and Prespecification models have acknowledged that neither of these models accounts adequately for the available experimental data.[29]

Turing-type reaction–diffusion model

This model, a reaction–diffusion model first proposed in 1979,[32] is based on the self-organizing properties of excitable media described by Alan Turing in 1952.[33] The excitable medium is the limb bud mesenchyme, in which cells interact by positively autoregulatory morphogens such as transforming growth factor beta (TGF-β) and inhibitory signaling pathways involving fibroblast growth factor (FGF) and Notch. Proximodistal and craniocaudal axes are not considered to be independently specified, but instead emerge by transitions in the number of parallel elements as the undifferentiated apical zone of the growing limb bud undergoes reshaping.[34] This model only specifies a "bare bones" pattern. Other factors like sonic hedgehog (Shh) and Hox proteins, primary informational molecules in the other models, are proposed instead to play a fine-tuning role.

Experimental evidence:

  • Limb mesenchymal cells, when dissociated and grown in culture or reintroduced within ectodermal "hulls" can recapitulate essential aspects of pattern formation, morphogenesis and differentiation.[35][36]
  • Peculiarities of the limb skeletal pattern in the mouse Doublefoot mutant are predicted outcomes of a Turing-type mechanism.[37]
  • Progressive reduction in distal Hox genes in a Gli3-null background results in progressively more severe polydactyly, displaying thinner and densely packed digits, suggesting (with the aid of computer modeling) that the dose of distal Hox genes modulates the period or wavelength of digits specified by a Turing-type mechanism.[38]

Craniocaudal patterning

Early signals that define the craniocaudal and proximodistal axis in vertebrate limb development.
Early signals that define the craniocaudal (anterior-posterior), and proximodistal axes in vertebrate limb development.

The discovery, in 1957, of the zone of polarizing activity (ZPA) in the limb bud provided a model for understanding the patterning activity by the action of a morphogenic gradient of sonic hedgehog (Shh).[39] Shh is recognised as a limb-specific enhancer.[40] Shh is both sufficient and necessary to create the ZPA and specify the craniocaudal pattern in the distal limb (Shh is not necessary for the polarity of the stylopod). Shh is turned on in the posterior through the early expression of Hoxd genes, the expression of Hoxb8, and the expression dHAND. Shh is maintained in the posterior through a feedback loop between the ZPA and the AER. Shh induces the AER to produce FGF4 and FGF8 which maintains the expression of Shh.

Digits 3,4 and 5 are specified by a temporal gradient of Shh. Digit 2 is specified by a long-range diffusible form of Shh and Digit 1 does not require Shh. Shh cleaves the Ci/Gli3 transcriptional repressor complex to convert the transcription factor Gli3 to an activator which activates the transcription of HoxD genes along the craniocaudal. Loss of the Gli3 repressor leads to the formation of generic (non-individualized) digits in extra quantities.[41]

Dorsoventral patterning

Dorsoventral patterning is mediated by Wnt7a signals in the overlying ectoderm not the mesoderm. Wnt7a is both necessary and sufficient to dorsalize the limb. Wnt7a also influences the craniocaudal and loss of Wnt7a causes the dorsal side of limbs to become ventral sides and causes missing posterior digits. Replacing Wnt7a signals rescues this defect. Wnt7a is also required to maintain expression of Shh.

Wnt7a also causes Lmx1b, a LIM Homebox gene (and thus a transcription factor), to be expressed. Lmx1b is involved in dorsalisation of the limb, which was shown by knocking out the Lmx1b gene in mice.[42] The mice lacking the Lmx1b produced ventral skin on both sides of their paws. There are other factors thought to control the DV patterning; Engrailed-1 represses the dorsalizing effect of Wnt7a on the ventral side of the limbs.[43]

See also

References

  1. 1.0 1.1 1.2 1.3 Stewart, TA; Bhat, R; Newman, SA (2017). "The evolutionary origin of digit patterning.". EvoDevo 8: 21. doi:10.1186/s13227-017-0084-8. PMID 29201343. 
  2. "Retinoic acid is required for the initiation of outgrowth in the chick limb bud". Curr Biol 6 (9): 1124–33. 1996. doi:10.1016/S0960-9822(02)70679-9. PMID 8805369. 
  3. "Retinoic acid promotes limb induction through effects on body axis extension but is unnecessary for limb patterning". Curr. Biol. 19 (12): 1050–7. June 2009. doi:10.1016/j.cub.2009.04.059. PMID 19464179. 
  4. Cunningham, T.J.; Duester, G. (2015). "Mechanisms of retinoic acid signalling and its roles in organ and limb development". Nat. Rev. Mol. Cell Biol. 16 (2): 110–123. doi:10.1038/nrm3932. PMID 25560970. 
  5. 5.0 5.1 Tickle, C (October 2015). "How the embryo makes a limb: determination, polarity and identity.". Journal of Anatomy 227 (4): 418–30. doi:10.1111/joa.12361. PMID 26249743. 
  6. Parr, BA; McMahon, AP (23 March 1995). "Dorsalizing signal Wnt-7a required for normal polarity of D-V and A-P axes of mouse limb.". Nature 374 (6520): 350–3. doi:10.1038/374350a0. PMID 7885472. Bibcode1995Natur.374..350P. 
  7. Pizette, S; Abate-Shen, C; Niswander, L (November 2001). "BMP controls proximodistal outgrowth, via induction of the apical ectodermal ridge, and dorsoventral patterning in the vertebrate limb.". Development 128 (22): 4463–74. doi:10.1242/dev.128.22.4463. PMID 11714672. 
  8. "FGF7 and FGF10 directly induce the apical ectodermal ridge in chick embryos". Dev. Biol. 211 (1): 133–43. 1999. doi:10.1006/dbio.1999.9290. PMID 10373311. 
  9. "Correlation of wing-leg identity in ectopic FGF-induced chimeric limbs with the differential expression of chick Tbx5 and Tbx4". Development 125 (1): 51–60. 1998. doi:10.1242/dev.125.1.51. PMID 9389663. 
  10. "The T-box genes Tbx4 and Tbx5 regulate limb outgrowth and identity". Nature 398 (6730): 814–8. 1999. doi:10.1038/19769. PMID 10235264. Bibcode1999Natur.398..814R. 
  11. DeLise, AM; Fischer, L; Tuan, RS (September 2000). "Cellular interactions and signaling in cartilage development.". Osteoarthritis and Cartilage 8 (5): 309–34. doi:10.1053/joca.1999.0306. PMID 10966838. 
  12. "All for one and one for all: condensations and the initiation of skeletal development". BioEssays 22 (2): 138–47. 2000. doi:10.1002/(SICI)1521-1878(200002)22:2<138::AID-BIES5>3.0.CO;2-4. PMID 10655033. 
  13. "Activator-inhibitor dynamics of vertebrate limb pattern formation". Birth Defects Research Part C: Embryo Today: Reviews 81 (4): 305–19. 2007. doi:10.1002/bdrc.20112. PMID 18228262. 
  14. Leal, Francisca; Cohn, Martin J. (2016-11-07). "Loss and Re-emergence of Legs in Snakes by Modular Evolution of Sonic hedgehog and HOXD Enhancers". Current Biology 26 (21): 2966–2973. doi:10.1016/j.cub.2016.09.020. ISSN 0960-9822. PMID 27773569. 
  15. Cohn, Martin J. (2008-10-07), Cardew, Gail; Goode, Jamie A., eds., "Developmental Mechanisms of Vertebrate Limb Evolution", Novartis Foundation Symposium, Novartis Foundation Symposia (Chichester, UK: John Wiley & Sons, Ltd) 232: 47–62, doi:10.1002/0470846658.ch5, ISBN 978-0-471-49433-1, PMID 11277086 
  16. 16.0 16.1 16.2 16.3 Infante, CR; Rasys, AM; Menke, DB (January 2018). "Appendages and gene regulatory networks: Lessons from the limbless.". Genesis 56 (1): e23078. doi:10.1002/dvg.23078. PMID 29076617. 
  17. 17.0 17.1 17.2 "Loss and Re-emergence of Legs in Snakes by Modular Evolution of Sonic hedgehog and HOXD Enhancers.". Curr Biol 26 (21): 2966–2973. 2016. doi:10.1016/j.cub.2016.09.020. PMID 27773569. 
  18. Tickle, Cheryll; Towers, Matthew (2017). "Sonic Hedgehog Signaling in Limb Development". Frontiers in Cell and Developmental Biology 5: 14. doi:10.3389/fcell.2017.00014. ISSN 2296-634X. PMID 28293554. 
  19. Harfe, Brian D.; Scherz, Paul J.; Nissim, Sahar; Tian, Hua; McMahon, Andrew P.; Tabin, Clifford J. (2004-08-20). "Evidence for an Expansion-Based Temporal Shh Gradient in Specifying Vertebrate Digit Identities" (in en). Cell 118 (4): 517–528. doi:10.1016/j.cell.2004.07.024. ISSN 0092-8674. PMID 15315763. 
  20. Harfe, Brian D.; Scherz, Paul J.; Nissim, Sahar; Tian, Hua; McMahon, Andrew P.; Tabin, Clifford J. (2004-08-20). "Evidence for an Expansion-Based Temporal Shh Gradient in Specifying Vertebrate Digit Identities" (in en). Cell 118 (4): 517–528. doi:10.1016/j.cell.2004.07.024. ISSN 0092-8674. PMID 15315763. 
  21. Kvon, Evgeny Z.; Kamneva, Olga K.; Melo, Uirá S.; Barozzi, Iros; Osterwalder, Marco; Mannion, Brandon J.; Tissières, Virginie; Pickle, Catherine S. et al. (2016-10-20). "Progressive Loss of Function in a Limb Enhancer during Snake Evolution" (in English). Cell 167 (3): 633–642.e11. doi:10.1016/j.cell.2016.09.028. ISSN 0092-8674. PMID 27768887. 
  22. 22.0 22.1 22.2 22.3 22.4 "Progressive Loss of Function in a Limb Enhancer during Snake Evolution.". Cell 167 (3): 633–642.e11. 2016. doi:10.1016/j.cell.2016.09.028. PMID 27768887. 
  23. "Development: Sonic snakes and regulation of limb formation.". Nat Rev Genet 17 (12): 715. 2016. doi:10.1038/nrg.2016.149. PMID 27843141. 
  24. "A novel ZRS variant causes preaxial polydactyly type I by increased sonic hedgehog expression in the developing limb bud.". Genet Med 22 (1): 189–198. 2020. doi:10.1038/s41436-019-0626-7. PMID 31395945. 
  25. 25.0 25.1 Larsen, William (2001). Human embryology (3rd ed.). Philadelphia, Pa.: Churchill Livingstone. p. 335. ISBN 978-0-443-06583-5. 
  26. Positional information and pattern formation. Current Topics in Developmental Biology. 6. 1971. 183–224. doi:10.1016/S0070-2153(08)60641-9. ISBN 9780121531065. 
  27. "Hox10 and Hox11 genes are required to globally pattern the mammalian skeleton". Science 301 (5631): 363–7. 2003. doi:10.1126/science.1085672. PMID 12869760. Bibcode2003Sci...301..363W. 
  28. 28.0 28.1 28.2 28.3 Larsen, William (2001). Human embryology (3rd ed.). Philadelphia, Pa.: Churchill Livingstone. p. 337. ISBN 978-0-443-06583-5. 
  29. 29.0 29.1 "Rethinking the proximodistal axis of the vertebrate limb in the molecular era". Genes Dev. 21 (12): 1433–42. 2007. doi:10.1101/gad.1547407. PMID 17575045. 
  30. "A re-examination of proximodistal patterning during vertebrate limb development". Nature 418 (6897): 539–44. 2002. doi:10.1038/nature00945. PMID 12152081. Bibcode2002Natur.418..539D. 
  31. "Functions of FGF signalling from the apical ectodermal ridge in limb development". Nature 418 (6897): 501–8. 2002. doi:10.1038/nature00902. PMID 12152071. Bibcode2002Natur.418..501S. 
  32. "Dynamics of skeletal pattern formation in developing chick limb". Science 205 (4407): 662–668. 1979. doi:10.1126/science.462174. PMID 462174. Bibcode1979Sci...205..662N. 
  33. Turing, AM (1952). "The Chemical Basis of Morphogenesis". Philosophical Transactions of the Royal Society B 237 (641): 37–72. doi:10.1098/rstb.1952.0012. Bibcode1952RSPTB.237...37T. 
  34. "Bare bones pattern formation: a core regulatory network in varying geometries reproduces major features of vertebrate limb development and evolution". PLOS ONE 5 (5): e:10892. 2010. doi:10.1371/journal.pone.0010892. PMID 20531940. Bibcode2010PLoSO...510892Z. 
  35. "Ectodermal FGFs induce perinodular inhibition of limb chondrogenesis in vitro and in vivo via FGF receptor 2". Dev. Biol. 249 (2): 270–82. 2002. doi:10.1006/dbio.2002.0766. PMID 12221006. 
  36. "Recombinant limbs as a model to study homeobox gene regulation during limb development". Dev. Biol. 166 (1): 59–72. 1994. doi:10.1006/dbio.1994.1296. PMID 7958460. 
  37. "Mixed-mode pattern in Doublefoot mutant mouse limb−Turing reaction-diffusion model on a growing domain during limb development". J Theor Biol 240 (4): 562–73. 2006. doi:10.1016/j.jtbi.2005.10.016. PMID 16364368. Bibcode2006JThBi.240..562M. https://ora.ox.ac.uk/objects/uuid:bc2d0b21-69b3-4871-a1ed-c981fd73f6eb. 
  38. "Hox genes regulate digit patterning by controlling the wavelength of a Turing-type mechanism". Science 338 (6113): 1476–1480. 2012. doi:10.1126/science.1226804. PMID 23239739. Bibcode2012Sci...338.1476S. 
  39. Zhu, J; Mackem, S (15 September 2017). "John Saunders' ZPA, Sonic hedgehog and digit identity - How does it really all work?". Developmental Biology 429 (2): 391–400. doi:10.1016/j.ydbio.2017.02.001. PMID 28161524. 
  40. "Developmental, genetic, and genomic insights into the evolutionary loss of limbs in snakes.". Genesis 56 (1): e23077. January 2018. doi:10.1002/dvg.23077. PMID 29095557. 
  41. "Manifestation of the limb prepattern: limb development in the absence of sonic hedgehog function". Dev Biol 236 (2): 421–35. 2001. doi:10.1006/dbio.2001.0346. PMID 11476582. 
  42. "'Induction of the LIM homeobox gene Lmx1b by WNT7a establishes dorsoventral pattern in the vertebrate limb". Cell 83 (2): 631–40. 1995. doi:10.1016/0092-8674(95)90103-5. PMID 7585966. 
  43. "Molecular Models for Vertebrate Limb Development". Cell 90 (6): 979–990. September 19, 1997. doi:10.1016/S0092-8674(00)80364-5. PMID 9323126.