<title>Introduction</title> <p>The genetic and developmental basis of the formation of organismal shape and form is a long-standing question in biology. The analysis of mutations has been essential in identifying the genes and regulatory networks underlying development. However, while the genetic basis of embryonic development has been extensively studied by systematic mutagenesis screens, we know little of the genes involved in the development of adult morphology. Yet, it is the heritable variation in adult form that natural selection primarily acts on during evolution. In order to understand the basis of variation, we need to know more about the genetic control of the development of adult form: which genes are involved, what are their function, and when are they required in development ##REF##10472187##[1]##,##UREF##0##[2]##. To identify genes important for development of adult structures, we initiated a large-scale mutagenesis screen in zebrafish and scored for mutants affected in the shape and pattern of adult structures. We isolated only adult viable mutants, therefore we selected for genes that have an increased probability to be involved in morphological change during evolution. Identification of zebrafish genes homologous to human genes associated with disease that arise during postnatal development into adulthood is also likely in this screen.</p> <p>We focused on mutants that exhibit defects in the dermal skeleton of the adult zebrafish. The dermal skeleton encompasses the external form of the adult fish. The most prominent dermal skeletal elements are the dermocranium of the skull and lateral bones of the opercular series, the scales, and the fin rays (or lepidotrichia). Additionally, the teeth and gill rakers (bones that support the gills in teleosts) are elements of the dermal skeleton ##REF##12803422##[3]##,##REF##9541261##[4]##. Unlike the ossification process that occurs during endochondrial bone development in which organic matrix is deposited by osteoblasts over a chondrogenic scaffold, dermal skeletal elements originate as direct mineralization of a collagenous matrix deposited by dermal fibroblasts. This process occurs in close association with the epidermis. The initiation and patterning of dermal elements are thought to be similar to epidermal appendages (e.g. hairs and feathers) and is controlled by reciprocal signaling between an epithelium and mesenchyme (see ##REF##14023393##[5]##,##REF##15071597##[6]##). Importantly, in zebrafish, as in most teleosts, the majority of dermal skeletal elements are not formed during larval development, rather through juvenile metamorphosis and development of the adult pattern. Those that begin to form in late larval development such as the teeth and gill rakers, do not fully attain their shape and pattern until juvenile metamorphosis.</p> <p>Variations in the shape of dermal skeletal elements of the fins, scales, cranium, and teeth play a significant role in adaptations of fish populations to new environments (e.g. dermal plate development and stickleback radiation ##UREF##1##[7]##). Additionally, integumentary appendages, such as hair and feathers, have been essential and defining traits of vertebrate classes. Early vertebrates, the conodonts, ostracoderms and placoderms, possessed a pronounced dermal skeleton, often in the absence of an ossified axial skeleton ##REF##12430166##[8]##. Through vertebrate evolution from fish to tetrapods, dermal structures such as lateral bones of the opercular series, scales, dermal plates and fin rays were either reduced or lost. This evolutionary transition was paralleled with the elaboration of the cartilaginous endoskeleton of the limbs and the evolution of specialized keratinized appendages of the integument such as epidermal scales, feathers and hairs. In contrast, the diversity of form in extant bony fishes involves modification in size, shape and number of the scales/dermal plates, fin rays, cranial dermal bones and teeth.</p> <p>Here, we describe a collection of mutants that have shared defects in the formation of the dermal skeletal elements of the skull, fins, scales and teeth of the adult zebrafish. The mutations disrupt the genes <italic>ectodysplasin</italic> (<italic>eda)</italic> and <italic>edar</italic> encoding the <italic>eda</italic> receptor. In mammals the EDA signaling pathway is involved in hair and teeth formation ##REF##12787560##[9]## and mutations affecting this pathway cause the human hereditary disease hypohidrotic ectodermal dysplasia (HED). Loss of Eda signaling in the zebrafish causes a spectrum of phenotypes corresponding to those described for HED in humans, and therefore the zebrafish mutants may serve as a genetic model of this disease. We describe the requirement of Eda signaling in the zebrafish epidermis for the formation of a structure resembling an epidermal placode seen in the early development of other vertebrate integumentary appendages. The mutations also result in defects of skeletal elements unique to fish and suggest an ancestral role of Eda signaling in the formation and patterning of the dermal skeleton. Lastly, whereas loss of function of Eda signaling causes a severe phenotype, the expressivity of dominant alleles is sensitive to background modifiers that buffer the phenotypic consequences of loss of Eda signaling. Additionally, we find that the response to reduction of Eda signaling is dose sensitive and organ specific. We suggest that such alleles may provide a basis for morphological variation in evolution.</p>
<title>Statistical Methods</title> <p>Analysis of cranial measurements were performed using Hotelling's T squared test for two dependent variables. For scale counts and size dimensions, a <italic>t</italic>-statistic for differential means was used to assess significance. Calculations and probability assessment were calculated using Biosoft 200 software (<ext-link ext-link-type="uri" xlink:href="http://www.biosoft.com">www.biosoft.com</ext-link>) and Excel statistical package.</p>
<title>Results</title> <title>
<italic>finless</italic> and <italic>Nackt</italic> Mutants Exhibit Defects in the Development of the Dermal Skeleton</title> <p>In a mutagenesis screen for mutations affecting adult zebrafish structures, we identified three mutants that showed nearly identical defects in the formation of scales, lepidotrichia, and shape of the skull of homozygous fish. These mutants fell into two complementation groups. The first is allelic to the <italic>finless (fls</italic>
<sup>te370f</sup>
<italic>)</italic> mutant that was previously isolated in the background of the Tübingen wild type stock (Tü) on the basis of the loss of fins in adults ##REF##24173565##[10]##. We isolated two new alleles of <italic>fls</italic> in the screen and found another in the background of the <italic>TLF</italic> wild type stock. The majority of the <italic>fls</italic> alleles isolated are recessive and have a strong phenotype (see below). However, the <italic>fls</italic>
<sup>dt3Tpl</sup> allele is dominant with a partial scale loss phenotype in heterozygotes (##FIG##0##Figure 1G, J##). One further <italic>fls</italic> allele was isolated in a screen for mutations that failed to complement the <italic>fls</italic>
<sup>te370f</sup> mutant (##FIG##0##Figure 1M##). We named this allele <italic>fang</italic> (<italic>fls</italic>
<sup>tfng</sup>) after its unique dental phenotype in homozygotes of having only one tooth on the fifth ceratobranchial (##FIG##0##Figure 1N##). The <italic>fls</italic>
<sup>tfng</sup> allele shows no effect on fin development and has a slight increase in the number of scales than the other <italic>fls</italic> alleles isolated.</p> <p>The second complementation group was comprised of a single gene, which we called <italic>Nackt</italic> (<italic>Nkt</italic>). This allele is dominant causing a slight defect in the patterning and shape of scales as heterozygotes (##FIG##1##Figure 2D##). The homozygous phenotype is more severe than that of strong <italic>fls</italic> alleles (##FIG##1##Figure 2A##).</p> <p>Phenotypic defects of <italic>fls</italic> and <italic>Nkt</italic> mutants become apparent in juvenile fish; as larvae, homozygous <italic>fls</italic> and <italic>Nkt</italic> mutants are visibly unaffected. Homozygous adults are viable, and of normal size. With the exception of the <italic>fls</italic>
<sup>tfng</sup> allele, the lepidotrichia that form during juvenile metamorphosis are defective, leading to fin loss in the adult (##FIG##0##Figure 1D##, ##FIG##1##2A##). The dermal bones of the pectoral girdle are present and patterned appropriately in both the <italic>fl</italic>s and <italic>Nkt</italic> mutants. By close examination of the visceral skeleton we found that neither the pharyngeal teeth, nor the bony substrates of the gills, the gill rakers, are formed (##FIG##0##Figure 1E and F##; ##FIG##1##Figure 2B, C##). In addition, scales are largely absent with infrequent formation of inappropriately shaped scales near the dorsal, anal and pectoral fins (##FIG##0##Figure 1D##, ##FIG##1##2A##). <italic>Nkt</italic> homozygous fish exhibit more severe defects in the formation of the dermal skeleton than <italic>fls</italic> alleles in the extent of lepidotrichial growth and number of residual scales formed (compare ##FIG##0##Figure 1D## and ##FIG##1##Figure 2A##). The skull of mutants has a normal appearance with all the bones being present, although the size, shape and relative proportion of the various bones differ compared to wild type individuals (##SUPPL##0##Figure S1##); no change in cranial shape was apparent in larvae.</p> <p>
<italic>Nkt</italic> heterozygous fish exhibit a dominant phenotype as they lack several scales on the flank and those present at the flank are elongated dorso-ventrally. The number of teeth and gill rakers is reduced, however lepidotrichia formation and growth of the fins are not affected (##FIG##1##Figure 2D–F##). The skulls of <italic>Nkt</italic> heterozygotes do not show increased size, but retain altered shape and proportion as seen in homozygotes (##SUPPL##0##Figure S1##).</p> <title>The Topless Allele of <italic>fls</italic> Uncovers Background Specific Modulation of <italic>fls</italic> Expressivity</title> <p>We isolated a dominant allele of <italic>fls</italic> that exhibits a distinct phenotype in heterozygotes that we named <italic>Topless</italic> (<italic>fls</italic>
<sup>dt3Tpl</sup>). Heterozygous <italic>fls</italic>
<sup>dt3Tpl</sup> have a reduction in the number of scales, teeth and gill rakers, but show little to no effect on lepidotrichia development (##FIG##0##Figure 1G–I##). Mutant <italic>fls</italic>
<sup>dt3Tpl</sup> fish exhibit the strong <italic>fls</italic> phenotype when homozygous or heterozygous with other <italic>fls</italic> alleles. Similar to <italic>Nkt</italic>, <italic>fls</italic>
<sup>dt3Tpl</sup> exhibits a dominant effect on skull shape as well (##SUPPL##0##Figure S1##).</p> <p>The expressivity of the dominant <italic>fls</italic>
<sup>dt3Tpl</sup> phenotype depends on the genetic background. Fish heterozygous for <italic>fls</italic>
<sup>dt3Tpl</sup> exhibited either a “strong” or a “weak” phenotype in the Tü background (##FIG##0##Figure 1G and J##, respectively). The “strong” phenotype shows loss of scales regionally in the midflank, loss of medial pharyngeal teeth along the fifth ceratobranchial and loss of posterior gill rakers of the anterior arches (##FIG##0##Figure 1G–I##). In contrast, the “weak” phenotype displays only subtle variation in scale patterning and no effect on the teeth or gill rakers could be detected (##FIG##0##Figure 1J–L##). The segregation pattern of the two phenotypic classes of <italic>fls</italic>
<sup>dt3Tpl</sup> suggests the presence of separate, unlinked, modifier loci in the Tü background affecting the number of scales (##TAB##0##Table 1## and data not shown). Additionally, we found that the <italic>fls</italic>
<sup>dt3Tpl</sup> “strong” phenotype was partially suppressed when crossed with the polymorphic WIK mapping strain indicating the presence of dominant modifier(s) in the WIK line (##TAB##0##Table 1##). The resulting heterozygous progeny showed reduced scale-loss compared to <italic>fls</italic>
<sup>dt3Tpl</sup> heterozygotes in a Tü background, but had similar defects in scale shape (##TAB##0##Table 1##). Therefore, dominant modifier loci are present in the WIK strain that buffer the expressivity of the <italic>fls</italic>
<sup>dt3Tpl</sup> dominant phenotype. None of the other <italic>fls</italic> alleles showed any dominance in the Tü, TLF, or WIK strains.</p> <title>Mutations in <italic>fls</italic> Disrupt the Ectodysplasin Receptor in Zebrafish</title> <p>We identified the affected loci of the <italic>fls</italic> mutants by positional cloning. The <italic>fls</italic>
<sup>te370f</sup> mutation was linked to SSLP markers on linkage group 9 (LG9). Due to similarity of the <italic>fls</italic> phenotype to ectodermal dysplasia phenotypes in mammals, we mapped several genes of the ectodysplasin pathway to the zebrafish radiation hybrid map to see if any of these genes were linked to <italic>fls</italic>. The <italic>edar</italic> gene is located on LG9 within the determined linkage interval for <italic>fls</italic> (see Methods). We cloned the full-length wild type cDNA of <italic>edar</italic> and found several polymorphisms in the Tü <italic>edar</italic> cDNA when compared to the WIK mapping strain; these polymorphisms were tightly linked with the <italic>fls</italic> mutation and did not show recombination in 238 meioses (##SUPPL##1##Figure S2##).</p> <p>The <italic>edar</italic> gene encodes a transmembrane protein with similarity to tumor necrosis factor receptor (TNFR). The Edar protein contains a conserved TNFR extracellular ligand binding domain and a cytoplasmic terminal death domain essential for protein interactions with signaling adaptor complexes. The <italic>fls</italic>
<sup>te370f</sup> mutation is an A to T transversion at a splice acceptor site, resulting in missplicing of the mRNA leading to a frame shift in translation and the generation of a premature stop codon (##FIG##2##Figure 3B## and ##SUPPL##1##Figure S2##). This allele is a likely molecular null mutation as only a fragment of the ligand-binding domain is present while the transmembrane and cytoplasmic death domains, which are essential for function of this protein, are both absent. The spontaneous mutation <italic>fls</italic>
<sup>t0sp212</sup> was found to have a splicing defect leading to the inclusion of intronic sequence. This is predicted to form a protein with incorrect amino acid sequence after residue 212, at the end of the transmembrane domain leading to a premature termination codon (##FIG##2##Figure 3B##, ##SUPPL##1##Figure S2##). The two alleles generated by the ENU mutagen both have missense mutations resulting in amino acid changes in the death domain (<italic>fls</italic>
<sup>t3R367W</sup>, R367W<sup>(C-T)</sup>; <italic>fls</italic>
<sup>dt3Tpl</sup>, I428F <sup>(A-T)</sup>). These mutations were found at identical positions as seen in familial cases of HED in humans (##FIG##2##Figure 3B, E##; ##REF##10431241##[11]##,##REF##17125505##[12]##.</p> <title>The <italic>fang</italic> Allele Uncovers Dose and Organ Specific Sensitivity to Levels of Eda Signaling</title> <p>The <italic>fang</italic> allele of <italic>fls</italic> was isolated in an allele screen for mutants that failed to complement <italic>fls<sup>te370f</sup></italic> (##FIG##0##Figure 1P##). <italic>fls<sup>tfng</sup></italic> homozygotes do not show any observable effect on lepidotrichia development yet have a reduction of scales and teeth/rakers as seen in other <italic>fls</italic> alleles (##FIG##0##Figure 1M–O##). The fang allele in <italic>trans</italic> to the te370f putative null allele shows an intermediate phenotype affecting lepidotrichial growth and a further reduction of teeth and scales suggesting that the <italic>fang</italic> allele is a hypomorph (##FIG##0##Figure 1P–R##); <italic>fls<sup>tfng</sup></italic> heterozygotes do not show any differences compared to wild type. The shape and number of the scales in <italic>fang</italic> is similar to the other homozygous <italic>fls</italic> alleles (##TAB##0##Table 1##). Analysis of <italic>edar</italic> RNA from homozygous <italic>fls</italic>
<sup>tfng</sup> showed the presence of two distinct transcripts with an additional larger isoform than seen in wildtype. Analysis of the sequence of the novel isoform showed the addition of intronic sequence leading to a premature termination codon (##FIG##2##Figure 3C##). The predicted protein would be similar to the <italic>fls</italic>
<sup>t0sp212</sup> allele having truncation just after the transmembrane domain at amino acid 218 (##FIG##2##Figure 3B## and ##SUPPL##1##Figure S2##). Analysis of the genomic sequence in the mutant revealed that the altered splicing is due to an A to G transition leading to the creation of a new splice donor site in the intron (##FIG##2##Figure 3C##). Given the presence of both isoforms in homozygous individuals, this novel splice site is used in addition to the normal splice junction. Using quantitative real time PCR we found that the <italic>fang</italic>-specific <italic>edar</italic> transcript represents 74% of the total pool of <italic>edar</italic> transcripts in homozygous mutants (##FIG##2##Figure 3D##). The dilution of wild type transcripts can explain the observed hypomorphic effect of the allele. From this unique allele of <italic>fls</italic>, it is clear that the phenotypic effect of loss of Eda signaling is dose dependent and that scales and teeth are more sensitive to alterations in the level of Eda signaling than fins.</p> <title>Ectodysplasin Is Mutated in the <italic>Nackt</italic> Mutant</title> <p>EDAR and its orthologue XEDAR recognize specific EDA isoforms that vary by two amino acids ##REF##11039935##[13]##–##REF##11309369##[16]##. The receptor-ligand complex signals though NF-κB using several adaptor proteins that are generally specific to each receptor. Together, mutations in <italic>Edar</italic> and <italic>Eda</italic> lead to the majority of cases of human HED in which the development of integumentary appendages (hairs, glands and teeth) are affected (OMIM 300451, ##REF##8696334##[17]##; OMIM 604095 ##REF##16435307##[18]##).</p> <p>We reasoned that, because of the similarity in phenotype to <italic>fls</italic>, the <italic>Nkt</italic> gene could be <italic>eda</italic>, encoding the ligand for Edar. We isolated the entire coding region for zebrafish <italic>eda</italic> by RACE (##SUPPL##2##Figure S3##). The <italic>eda</italic> transcript from the <italic>Nkt</italic> mutant shows a precocious stop codon predicting a truncation of the protein at the beginning of the TNF domain, which is involved in ligand-receptor binding (S243X<sup>(C-A)</sup>); ##FIG##2##Figure 3G## and ##SUPPL##2##Figure S3##). An analysis of the location of <italic>eda</italic> in the zebrafish radiation hybrid map placed <italic>eda</italic> on LG5<bold>.</bold> Subsequent linkage analysis of the <italic>Nkt</italic> mutation and markers indicated by radiation hybrid analysis demonstrated tight linkage of the mutant to this region (##FIG##2##Figure 3F##); the S243X mutation was always found in fish with the <italic>Nkt</italic> phenotype and served as a consistent genotypic marker.</p> <title>Role of Ectodysplasin Signaling in Regulating Epithelial Signaling Centers: Scale Placode Formation</title> <p>In fish, scales are bony elements that develop in the dermis underlying the epidermis. In amniotes, most integumentary organs affected by loss of Eda signaling have structural derivatives stemming from the epidermis (<italic>e.g</italic>. specific keratins of hair, feather and nail, secretory cells of glands). These integumentary organs develop from reciprocal signaling interactions between the basal epidermis and subjacent mesenchyme often controlled by a regional epithelial thickening called the epidermal placode. Eda signaling is necessary for the development and patterning of epithelial placodes of many integumentary organs in both the mouse and chick ##REF##10431226##[19]##–##REF##17362907##[22]##. Expression of <italic>Eda</italic> and <italic>Edar</italic> is found predominantly in the basal epidermal cells, but in the case of feathers <italic>Eda</italic> is detected in the subjacent mesenchyme as well ##REF##15673574##[21]##,##REF##17362907##[22]##. Whereas expression of developmental signaling genes such as <italic>sonic hedgehog (shh)</italic> in the development of integumentary appendages are comparable between vertebrates ##REF##15272389##[23]##, evidence for an early developmental role of the epidermis in induction or patterning of the teleost scale is lacking. The formation of an epithelial placode and signaling center in the development of amniote integumentary appendages is associated with histological changes in the basal cells of the epidermis; a similar structure has not been described in fish epidermis ##REF##9183678##[24]##. As early teleost scale development is quite different to that of other vertebrate integumentary organs, such as hairs and feathers, we addressed the question whether Eda signaling had a similar function in the epidermis of zebrafish during scale formation.</p> <p>We detected the expression of both <italic>edar</italic> and <italic>eda</italic> in the skin of juvenile fish by whole mount <italic>in situ</italic> analysis (WMISH). The expression of both genes presaged the formation of the initial scale row along the flank just ventral from the midline mysoseptum in the caudal peduncle (arrowheads, ##FIG##3##Figure 4A and C##; ##REF##9183678##[24]##). During scale formation, the expression of <italic>edar</italic> becomes progressively restricted to the posterior margin (##FIG##3##Figure 4B##) while <italic>eda</italic> expression persists throughout the scale primordia (##FIG##3##Figure 4D##). Developmental genes <italic>shh</italic> and <italic>bone morphogenic protein 2b</italic> (<italic>bmp2b</italic>), whose orthologues are known to be essential for placode development in the mouse and chick, show similar placodal expression as seen with <italic>edar</italic> (##FIG##3##Figure 4E and G##). Analysis of <italic>shh</italic> and <italic>bmp2b</italic> expression in <italic>fls</italic>
<sup>te370f</sup> indicated the necessity of <italic>edar</italic> function for their expression (##FIG##3##Figure 4F and H##).</p> <p>We investigated the development of scale primordia in wild type and mutant <italic>fls</italic>
<sup>te370f</sup> fish by light and transmission electron microscopy. Previous detailed histological work found evidence for raised signaling activity in the epidermis as measured by increased endoplasmic reticulum (ER), and secretory activity of the basal epidermal cells prior to scale formation ##REF##15272389##[23]##. However these changes in the basal epidermal cells were not associated regionally with sites of scale formation nor was there any indication of histological changes in basal cell morphology that are associated with placode formation in other vertebrates. To our surprise, in wild type juvenile fish, we discovered the formation of histologically defined, localized assemblies of cells of the basal epidermis that resemble early stages of the formation of hair and feather placodes.</p> <p>Prior to the development of the scale, the dermis consists of compact collagen layers (stratum compactum) and scattered dermal fibroblasts ##REF##9183678##[24]##. The epidermal basal cells have a uniform elongate morphology (black arrows, ##FIG##4##Figure 5A, D##) with high levels of basally located intermediate fibrils (##FIG##4##Figure 5G##). At the initiation of scale development, there is an accumulation of fibroblasts subjacent to the basal epidermis, associated with a reworking of the collagen strata ##REF##9183678##[24]##. We find a specific alteration in the morphology of the basal epidermal cells in wild type juvenile fish that coincides with the initial accumulation of fibroblasts at the sites of scale development (##FIG##4##Figure 5B, E##). These basal cells become cuboidal and have decreased width, such that they form a unit of closely packed cells (black arrows, ##FIG##4##Figure 5B, E##). This is observed above the localized accumulation of fibroblasts in the dermis (white arrows, ##FIG##4##Figure 5B, E##). In addition, in these placodal-like cells, the ER appears less prominent (data not shown), and hemidesmosomes, the cellular junctions involved in the attachment of the basal epidermal cells to the basal lamina, are almost completely absent (brackets, ##FIG##4##Figure 5H and I##). In contrast, the adjacent lateral epidermal cells show high levels of both ER and hemidesmosomes (data not shown and ##FIG##4##Figure 5G## brackets, respectively).</p> <p>In <italic>fls</italic>
<sup>te370f</sup> juvenile fish, at a corresponding site on the flank as in wild type, we detected the formation of similar aggregations of basal epidermal cells (black arrows. ##FIG##4##Figure 5C, F##). However, unlike the structures found in the wild type zebrafish, the epidermal cells of the placode were disorganized and showed histological evidence of cell death (##FIG##4##Figure 5F##). As is the case in wild type, the epidermal basal cells in the placode of <italic>fls</italic> display a reduced ER, however hemidesmosomes are present in the same high numbers as in adjacent cells in wildtype (brackets, ##FIG##4##Figure 5H## compared to brackets ##FIG##4##Fig. 5I##). Lateral basal epidermal cells in <italic>fls</italic>
<sup>te370f</sup>
<italic>/edar</italic> showed elongate morphology similar to those of their wild type siblings (data not shown). The expression of <italic>edar</italic> is seen in the basal cells of forming scale placodes (##FIG##4##Figure 5J–L##; stages 1–3 ##FIG##4##Figure 5M##) arising during early specification of the scale placode (s1; arrowhead ##FIG##4##Figure 5K##). We were unable to detect <italic>eda</italic> expression in sections due to the weak hybridization signal.</p> <p>These data support the notion that an epidermal placode is involved in dermal scale formation. Further we find that the epithelial organization and function of the developing scale placode is dependent on <italic>edar</italic>.</p> <title>Role of Ectodysplasin Signaling in Regulating Epithelial Signaling Centers:Maintenance of the Fin Fold and Establishing Anterior-Posterior Polarity of the Developing Fin</title> <p>The phenotype of both <italic>Nkt</italic> and <italic>fls</italic> demonstrate that Eda signaling is necessary for fin development. The growth and patterning of lepidotrichia are affected in all fins. Lepidotrichia are specified, however fail to maintain growth and elaboration of the fin rays (##FIG##5##Figure 6A–C, H–P##). Unpaired fins showed no defects in patterning of the endochondrial bones of the proximal and distal radials (##FIG##5##Figure 6K–P##); the dorsal pitch of the caudal fin is an indirect effect of the mutation on swimming without fin rays (amputated fins that fail to regenerate show similar morphology). In <italic>fls</italic> adults, fusions of the distal radials of the pectoral fin are seen at a low penetrance (##FIG##5##Figure 6G##, data not shown). In <italic>Nkt</italic>, there is an increase in the frequency of patterning and growth defects of the endochondrial components of the fin (##FIG##5##Figure 6G##). These alterations include the loss of the fourth proximal radial, altered growth patterns of anterior proximal radials 1 and 2, as well as lack of articulation of the distal radials (##FIG##5##Figure 6## D–F). <italic>Nkt</italic> causes a strong effect on lepidotrichial growth of both the pectoral and pelvic fins (##FIG##5##Figure 6E–F, J##). In contrast, a specific effect on the growth of anterior lepidotrichia of the pelvic fin is seen in <italic>fls</italic> where the dermal rays of the anterior (<italic>e.g.</italic> 1,2) are significantly shorter than rays at equivalent positions in wild type (##FIG##5##Figure 6I##). The asymmetry of lepidotrichial development suggests that, like the proximal endochondrial fin skeleton, the fin rays have a specific regional identity to provide the shape and form of the fin.</p> <p>Early limb development is driven by a localized organization of epithelial cells at the distal tip of the forming limb, termed the apical ectodermal ridge (AER) ##REF##18882505##[25]##. In zebrafish, the AER is involved in larval patterning of the paired fins, while the later stages of fin development are organized by an analogous epidermal formation of the fin fold in both paired and unpaired fins ##REF##10349624##[26]##,##REF##4031750##[27]##. In tetrapod limb development anterior-posterior specification is controlled by posterior mesenchyme expressing <italic>Shh</italic>. The function of this zone of polarizing activity (ZPA) expressing <italic>Shh</italic> is maintained by reciprocal signaling interactions between the ZPA and the AER. This interaction is necessary for proper patterning and growth of the tetrapod limb. In the zebrafish, <italic>shh</italic> and signals from the AER also orchestrate patterning and outgrowth of the early fin buds ##REF##7601313##[28]##–##REF##10976049##[30]##. In addition, genes functioning early in fin development, such as <italic>shh</italic>, play important roles during late fin development regulating growth and branching of lepidotrichia growth ##REF##9753672##[31]##.</p> <p>We investigated the regulation of <italic>edar</italic> in mid to late fin development focusing on the development of the paired fins. In early fin fold stage of pectoral fin development (8 mm), we detected <italic>edar</italic> expression in both the distal margin of the endochondrial radials (black arrowhead, ##FIG##6##Figure 7A##) as well as more distally in the forming lepidotrichial rays (##FIG##6##Figure 7A##). The expression of <italic>edar</italic> in the fin fold had a posterior bias in wild type fins (##FIG##6##Figure 7A##, white arrow). The pelvic fin showed similar expression of <italic>edar</italic> in forming lepidotrichial rays (##FIG##6##Figure 7E##). <italic>shh</italic> and <italic>bmp2b</italic> expression was observed in the forming lepidotrichia of both the pectoral and pelvic fins of wild type juveniles, having a similar distal bias in the leading margin of all rays (##FIG##6##Figure 7C, G and I, K##, respectively). Expression of all three genes in <italic>fls</italic> was decreased in the anterior portion of the pectoral fins (##FIG##6##Figure. 7B, D, J##). However, residual expression of all three genes was found in the posterior margin of the fin (##FIG##6##Figure 7B, D, J## arrow). In the pelvic fins, similar loss of anterior expression of <italic>edar</italic> (##FIG##6##Figure 7F##) and <italic>shh</italic> (##FIG##6##Figure 7H##) was seen in the <italic>fls</italic> mutant. We did not detect any difference in <italic>bmp2b</italic> expression in the pelvic fins even though obvious morphological differences in the developing rays of the samples could be seen (##FIG##6##Figure 7L##).</p> <p>We asked if the alteration of polarity of gene expression in the <italic>fls</italic> mutant was associated with regional cell death. Using acridine orange uptake as an assay for cell death (<italic>e.g.</italic>
##REF##8223253##[32]##), we analyzed fins of <italic>fls</italic> and siblings at size matched stages (7–9 mm) for regional patterns of cell death. In 7 mm juveniles, we detect a differential retention of acridine orange between <italic>fls</italic> and siblings in the anterior, and anterior-distal margin of the developing pectoral fin (##FIG##6##Figure 7M–N##). Similarly, in the pelvic fin of 8–9 mm juveniles, retention of the label was seen in the anterior distal margin of the fin (##FIG##6##Figure 7O–P##). Consistent with these data suggesting asymmetrical loss of the fin fold epidermis, we find that <italic>msxa</italic>, a marker of the distal epithelium, is differentially expressed in the mutant (##FIG##6##Figure 7Q and R##).</p> <p>We next analyzed gene expression during late development of the lepidotrichia. The expression of <italic>edar</italic> during late fin development was observed in forming lepidotrichia of all fins with a distal bias in its expression (##FIG##7##Figure 8C##). The expression in the forming dermal ray was similar in both location and timing to that of <italic>bmp2b</italic> and <italic>shh</italic> (##FIG##7##Figure 8A–B##, asterisk). In addition, <italic>edar</italic> was found expressed proximally between forming rays and at the distal margin (##FIG##7##Figure 8C## arrow). We were unable to resolve a clear signal for <italic>eda</italic> in the forming fins using WMISH. The expression of <italic>edar</italic> in the distal lepidotrichial tips suggests a late developmental role of Eda signaling in regulating formation of the lepidotrichia in concert with <italic>shh</italic> and <italic>bmp2b</italic>.</p> <p>Histological analysis of <italic>fls</italic> mutant fins at an early stage of lepidotrichial formation reveal a general deficiency of the development of the entire mesodermal component of the fin such as cartilage and muscle in both the paired (pectoral fin, ##FIG##7##Figure 8D, E##) and unpaired fins (anal and caudal fins, ##FIG##7##Figure 8F–G##; H–I, respectively). In contrast, the epidermis of the fin is formed and is similar to that of size matched siblings. However, close inspection of the distal tip of the fins showed disorganization of the epidermis and degeneration of distal epidermal nuclei (insets ##FIG##7##Figure 8F, J##). From these analyses, we hypothesize that loss of <italic>edar</italic>-mediated signaling leads to a defect in mesenchymal cell proliferation, muscle cell migration and defective lepidotrichial growth in the fin that correlates with degenerative defects seen in the distal epidermal fin fold.</p>
<title>Discussion</title> <p>We used a forward mutagenesis approach in the zebrafish to investigate the developmental mechanisms that underlie changes in adult form. Here, we identified a role of Eda signaling in the development of the dermal skeleton in the adult zebrafish. Mutations in either <italic>Eda</italic> or <italic>Edar</italic> have been shown previously to cause defects in integumentary appendages in several mammalian species. Additionally, Eda signaling genes have been associated with variation in morphology that occurred during the evolution of teleost fishes <italic>(eda)</italic> and in variation of human populations <italic>(Edar)</italic>
##REF##15790847##[33]##,##REF##17943131##[34]##; see below). Thus, through a forward genetic approach, we were successful in identifying genes that are important for the development and variation in adult form. We further show that the ENU generated alleles of <italic>fls</italic>, <italic>t3R357W</italic> and <italic>dt3Tpl</italic>, affect similar residues as those mutated in familial cases of HED ##REF##10431241##[11]##,##REF##17125505##[12]##,##REF##16435307##[18]## supporting the utility of adult zebrafish mutants as models for the investigation of heritable human disease.</p> <title>Conserved and Ancestral Role of Eda Signaling in Vertebrate Development</title> <p>We show that Eda signaling is necessary for the development and patterning of the dermal bones of the skull, scales, fin rays as well as teeth of the adult zebrafish. The correlated effect in these zebrafish structures is due to a developmental role of Eda signaling in organizing epithelial cells into signaling centers. In the case of scale development, Eda signaling is necessary for the basal epidermal cells to form a functional placode. Epidermal placodes are involved in the formation of integumentary appendages of other vertebrates such as hair, glands, feathers and teeth. These structures have been shown to act as signaling centers to orchestrate appendage development. We speculate that a primary function of Eda signaling in scale development is to promote cell-cell adhesion within the placode and that the coordinated signaling of the placode induces fibroblast assembly in the underlying dermis, an early step in scale formation.</p> <p>Schmidt-Ullrich et al. documents the formation of the hair placode and outline a stage series of placode formation in the mouse ##REF##16481354##[35]##. They report that the <italic>downless</italic> mouse mutant, which has a mutation that disrupts the mouse <italic>Edar</italic> gene ##REF##10431242##[36]##, causes arrest of placode formation at a pre-placode stage of development (P0–P1). This stage closely resembles the stage of scale placode formation that is affected in <italic>fls</italic> shown here. In agreement with our findings, Schmidt-Ullrich et al. further note a reduction of cell-to-cell adhesion within the placode and find increased apoptosis in the absence of Edar function. This suggests that there is a conserved developmental role of Edar between dermal scales of fish and mammalian hairs. During normal hair development, the hair placode invaginates to form the hair bulb. By contrast, the post-placodal events of scale formation in fish do not involve morphogenetic changes of the epidermis, rather the accumulation of mesenchymal cells subjacent to the epidermal placodal cells to form the scale pocket. Thus, Eda signaling in mammals and teleosts is conserved in the early phases of placode formation in controlling the functional continuity and signaling of the epidermal placode to orchestrate appendage formation. However, the downstream interpretation of the epithelial-mesenchymal signaling differs beyond this point leading to altered morphogenetic responses and histological differentiation to form diverse appendages such as scales and hair.</p> <p>In the fin, Eda signaling directs late stages of fin development such as the formation and growth of the dermal rays. The effect of loss of <italic>edar</italic> function on fin development uncovers an intrinsic developmental polarity of the late developing fish fin. This is seen both in the development of the proximal endochondrial bones as well as in the formation of the fin rays. We find that the change in patterning in the mutants is correlated with asymmetrical cell death of the distal marginal fin fold as well as a reduction of <italic>shh</italic> expression. This finding is similar to the effect of loss of AER function resulting in anterior-distal cell death and reduction of Shh activity in tetrapod limbs ##REF##18265010##[37]##–##REF##11476582##[39]##. While there has not been any previous indication of a role of Eda signaling in tetrapod limb development, both the expression of <italic>Edar</italic> and related receptor, <italic>Troy</italic>, have been detected in the AER of mice ##REF##12972005##[40]##,##REF##11023871##[41]##.</p> <p>A second developmental role of Eda signaling in the developing fin is observed in the outgrowth and patterning of the individual lepidotrichial rays evidenced by expression of <italic>edar</italic> in the distal tip of the forming rays and distal epidermis. The expression of <italic>edar</italic> is again associated with that of <italic>shh</italic> and <italic>bmp2b</italic>. The expression of <italic>shh</italic> and <italic>bmp2b</italic> has been shown to be within the basal epidermis overlying the forming lepidotrichia ##REF##9753672##[31]##. Given the expression of <italic>edar</italic> during fin development and the defects observed in the distal epidermis in the mutant, it is likely that the function of Eda signaling is to maintain the growth permissive function of the fin fold through its regulation of a distal signaling center of individual rays. The concomitant expression of <italic>edar</italic>, <italic>shh</italic> and <italic>bmp2b</italic> in both distal lepidotrichia development and during placode specification suggest that they work in concert to mediate the inductive and/or permissive effects of the epidermis – thus organizing signaling centers for the development of the dermal skeleton.</p> <p>While the nature of the defect in tooth formation or dermal bone patterning of the skull in the <italic>fls</italic> and <italic>Nkt</italic> mutants has not been characterized in detail, there is evidence that inductive signaling from the pharyngeal epithelium or cranial epidermis is necessary for appropriate development of both tooth ##REF##14550785##[42]## and skull ##REF##15071597##[6]##, respectively; Eda signaling likely shares a common role in inductive signaling in each of these diverse organs.</p> <title>Genetic Commonalities of Eda Signaling: From Fish to Man</title> <p>Mutations in the <italic>EDAR</italic> and <italic>EDA</italic> genes underlie a large percentage of autosomal and X-linked HED in humans, respectively ##REF##16435307##[18]##,##REF##11378824##[43]##. In the case of <italic>EDAR</italic>, both recessive and dominant mutations are associated with the HED phenotype in humans, however dominant mutations are found only within the death domain of the protein. These mutations are believed to act in a dominant negative fashion, although by unknown mechanisms ##REF##11570810##[44]##. We see similar dominance of a <italic>fls</italic> allele that affects the death domain of <italic>edar</italic> while all <italic>fls</italic> mutations outside this region do not show a dominant phenotype. Autosomal dominant HED in humans caused by mutation of <italic>EDAR</italic> within the death domain displays a large degree of phenotypic variability ##REF##17125505##[12]##. For example, the I418T mutation in human, which affects the same amino acid as <italic>fls</italic>
<sup>dt3Tpl</sup> (I327F), shows distinct phenotypic variability depending on genetic background ##REF##16435307##[18]##. Interestingly, the <italic>fls</italic>
<sup>dt3Tpl</sup> zebrafish mutant displays similar dominance and variation as the human allele affecting the same residue. These findings suggest that the molecular mechanisms of Edar function are similar between fish and humans.</p> <p>X-linked HED caused by mutations in the human <italic>EDA</italic> gene represents the majority of cases of this disease ##REF##11378824##[43]##,##UREF##2##[45]##. The zebrafish <italic>Nkt</italic> mutation described here is affected in the TNF domain and shows a mild dominant phenotype (S243X). As the <italic>EDA</italic> gene is sex linked in humans the molecular nature of different alleles can not be analyzed since the allele will be hemizygous in males and mosaic in female carriers. The zebrafish <italic>eda</italic> gene is autosomal in the zebrafish. Thus, <italic>Nkt</italic> exposes previously unknown dominant function of mutations in this gene since a true heterozygous condition is formed. Hemizygous wildtype condition in humans indicates that the dominance we see in <italic>Eda</italic> is probably not due to haploinsufficiency. Since EDA functions as a homotrimeric protein ##REF##10484778##[46]##, a plausible mechanistic explanation for the observed dominance of <italic>Nkt</italic> is that the C-terminal truncation inhibits the function of the wild type protein in binding to Edar.</p> <title>Eda Signaling and the Development and Variation of Adult Form</title> <p>Mutations affecting Eda signaling lead to impaired development of integumentary appendages of fish, birds, and man. These changes lead to viable changes in adult morphology. Mutations disrupting Eda signaling have been described for another teleost species. The spontaneous <italic>rs-3</italic> mutant in medaka (<italic>Oryzias latipes</italic>), is shown to have a transposon insertion in the 5′ UTR of <italic>edar</italic> resulting in the reduction of scales but no effect on fin or teeth development ##REF##11516953##[47]##. The zebrafish mutations described here show a previously undescribed role of Eda signaling in the development of the fins, teeth, as well as dermal bones of the skull – phenotypic traits observed in human alleles but not reported in the medaka mutant. As the phenotype of the <italic>rs-3</italic> mutant is similar to the <italic>fls</italic>
<sup>dfang</sup> allele in the zebrafish, it is likely that the more subtle phenotypes observed in the medaka mutant is due to partial loss of function of <italic>edar</italic> caused by a hypomorphic <italic>rs-3</italic> allele ##REF##11516953##[47]##.</p> <p>The graded effects seen in the expressivity of mutations affecting <italic>eda</italic> and <italic>edar</italic> points to a dose sensitive readout of the Eda signaling pathway that affects different organ systems with varied expressivity. In the dominant <italic>fls</italic>
<sup>dt3Tpl</sup> or <italic>Nkt</italic> heterozygotes, the shape and number of scales and teeth as well as patterning of the skull are affected, however there is no change in fin ray development. Similarly, the <italic>fang</italic> allele of <italic>fls</italic> clearly demonstrates this dose sensitivity as functional copies of Edar are titrated by the concomitant use of a new splice site in the mutant leading to the reduction in the amount of wild type transcript made (##FIG##2##Figure 3D##). This reduction in the amount of <italic>edar</italic> transcripts cause defects in scale and tooth development, however fins are normal. <italic>fang/te370f</italic>, in which the <italic>fang</italic> allele is in <italic>trans</italic> to a presumed null, further reduces the relative levels of wild type <italic>edar</italic> transcripts leading to further reduction of both teeth as well as fin lepidotrichia. Similar dose sensitive responses to levels of EDA signaling are seen in tooth development of the mouse regulating the number and shape of teeth ##REF##15031115##[48]##,##REF##15538367##[49]##. There are several reports of hypodontia in humans resulting from altered EDA function that do not show other phenotypes such as hypothrichosis or nail defects ##REF##18545687##[50]##–##REF##17256800##[52]##. Given our findings, it is likely that these particular alleles are hypomorphic and this is sufficient to explain the differential organ sensitivity to levels of EDA signaling during development. These data indicate that control of the level of Eda signaling in post-embryonic development is an essential component for the determination of the number and form of many different organ systems of the adult.</p> <p>Supporting this finding, we observed significant modification of expressivity of <italic>fls</italic>
<sup>dt3Tpl</sup> in different genetic backgrounds indicating the existence of genetic modifiers of Eda signaling. This sensitivity of Eda signaling to genetic modifiers occurs in other teleost fish as well. In our analysis of the medaka <italic>rs-3/edar</italic> mutant, we find a high degree of variability in the extent of scale formation (##SUPPL##3##Figure S4##) suggesting the existence of background modifiers of Edar function in this species. Additionally, evidence from the stickleback, <italic>Gasterosteus aculeatus</italic>, suggests that genetic variance at the <italic>eda</italic> locus underlies differences in the extent of dermal plate formation in diverged populations of this species ##REF##15790847##[33]##. A quantitative trait analysis (QTL) of lateral plate formation in a low-plated form of the stickleback indicates a significant modification of the reduced plate phenotype (<italic>eda</italic> locus) with modifying effects within and between loci affecting plate number and size ##REF##11780061##[53]##,##REF##15069472##[54]##. Interestingly, recent evidence also shows a significant association between the <italic>edar</italic> locus and dermal plate number in sticklebacks in addition to the predominant <italic>eda</italic> locus ##REF##17371397##[55]##. Thus variation at these gene loci may act in concert to regulate number of dermal plates/scales.</p> <p>Thus, while loss of Eda signaling can lead to severe phenotypes, the phenotypic consequences of variation in Eda signaling are graded and canalization of Eda signaling is prevalent. Therefore, buffering of the phenotypic outcome that results from defective Eda signaling could be a common mechanism that permits viable and diverse phenotypes. These viable phenotypic variations then could serve as a basis for selection. The lack of a coding change at the <italic>eda</italic> locus in sticklebacks that is associated with the loss of dermal plates has lead to the argument that, in this case, evolution of this trait is due to changes at <italic>cis</italic>-regulatory elements controlling <italic>eda</italic> expression ##REF##15790847##[33]##. Our findings on the dose and organ specific sensitivity of Eda signaling in different structures of the zebrafish argues that evolution of this trait could result from a regulation of absolute levels of expression.</p> <p>Interestingly, recent analysis of single nucleotide polymorphism (SNP) frequency in human populations supports the role of Eda signaling in causing phenotypic variation. Analysis of SNP variation between diverse human populations shows evidence of selection of the <italic>EDAR</italic> locus in East Asian and American populations ##REF##17356696##[56]##,##REF##17542651##[57]##. A defined allelic variant of <italic>EDAR</italic> within these populations leads to a coding change in the death domain of EDAR and is a candidate allele for altered gene function that could have lead to the region being fixed in these populations ##REF##17943131##[34]##. There is evidence from association data that this allele is associated with thick hair in these populations ##REF##18065779##[58]##, however the full extent of phenotypes that are affected in these populations that are related to EDA signaling has not been analyzed. It is interesting to note that recent work has identified this allele of <italic>EDAR</italic> as having an enhanced effect on Eda signaling in mouse models containing the altered human residue ##UREF##3##[59]##. Given that variation in the number and shape of integumentary derivatives of the dermal skeleton are a common morphological change in teleost evolution <italic>e.g.</italic>
##UREF##4##[60]##, it will be important to further investigate the prevalence and type of genetic changes in Eda signaling genes in cases of natural variation of these adult characters.</p> <title>Mutagenesis and Allele Designation</title> <p>Zebrafish mutagenesis was performed following ##UREF##5##[61]## with 5 treatments of 3.3–3.5 mM ethylnitrosourea. Screen design was similar to that described ##UREF##6##[62]##. Allele designation was determined using standard nomenclature with the addition of the molecular lesion or phenotypic description (when appropriate) to the designation. The serial numbers of the mutants found in the ZF models screen are as follows: <italic>fls</italic>
<sup>t3R367W</sup> (#0621); <italic>fls</italic>
<sup>dt3Tpl</sup> (#1248); <italic>Nkt</italic>
<sup>dt3S243X</sup> (#1261). Information on the screen can be found at <ext-link ext-link-type="uri" xlink:href="http://www.zf-models.org/">http://www.zf-models.org/</ext-link>. The screen for additional <italic>fls</italic> alleles used mutagenized <italic>TLF</italic> founder males treated similarly as Tü males used in the screen.</p> <title>Mapping</title> <p>Rough mapping of F2 progeny against a reference panel of SSLP markers ##UREF##7##[63]## indicated that <italic>fls</italic> was located on linkage group (chromosome) 9 (LG9) with loose linkage to z20031 (61.3cR; ##FIG##2##Figure 3A##). We found <italic>fls</italic> to be closely linked to markers z7001 and z11672. Results from radiation hybrid screening indicated linkage of zebrafish <italic>edar</italic> to markers positioned on LG9 in the region predicted by initial mapping analysis. Analysis of flanking markers and internal polymorphisms in <italic>edar</italic> showed tight linkage of the <italic>fls</italic> mutation to the <italic>edar</italic> gene. Using the defined molecular differences between WIK and Tü strains, we did not find recombination in 238 meioses indicating that the mutation was located less than 0.4 cM away.</p> <title>Cloning and Sequence Analysis</title> <p>We isolated the full-length cDNA of zebrafish <italic>edar</italic> and <italic>eda</italic> by reverse transcription (RT) PCR using sequences provided from genomic alignments and subsequent amplification of the cDNA ends by rapid amplification of DNA ends (RACE). cDNA was generated from RNA from blastemas of amputated caudal fins that had been allowed to regenerate for two days. cDNA sequences of zebrafish <italic>edar</italic> and <italic>eda</italic> genbank accession numbers are EF137867 and EF137866, respectively. Protein alignment of Edar and Eda were generated by ClustalW alignment (<ext-link ext-link-type="uri" xlink:href="http://www.ebi.ac.uk/clustalw/">http://www.ebi.ac.uk/clustalw/</ext-link>) and Box Shade software (<ext-link ext-link-type="uri" xlink:href="http://www.ch.embnet.org/software/BOX_form.html">http://www.ch.embnet.org/software/BOX_form.html</ext-link>) using a 0.4 identity threshold. <italic>Edar</italic> and <italic>eda</italic> sequences of other species were obtained from genomic databases at NCBI (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/">http://www.ncbi.nlm.nih.gov/</ext-link>), Sanger (<ext-link ext-link-type="uri" xlink:href="http://www.ensembl.org/index.html">http://www.ensembl.org/index.html</ext-link>), and Tigr (<ext-link ext-link-type="uri" xlink:href="http://www.tigr.org/tdb/tgi/">http://www.tigr.org/tdb/tgi/</ext-link>).</p> <p>Real Time PCR was performed on cDNA obtained from blastemas from two day old caudal fin regenerates. Calculations were made from three biological replicates and three technical replicates according to ##REF##11328886##[64]##. Primers were designed for wild type specific transcripts by using sequence from the neighbouring exon borders that are not adjacent to each other in Fang-transcripts. Fang specific primers were designed against the fang specific transcript sequence, which is spliced out in wild type. Crossing points of the control reaction, which were higher than in the water control, were set to the value for water. Normalization was done against the efficiency of primers to β-actin.</p> <title>Bone Stain and Measurements</title> <p>Adult bones were stained with alizarin red. Embryos were fixed in formalin (3.7% formaldehyde), briefly dehydrated in 70% ethanol, and placed in 1 g/l alizarin red; 0.5% KOH until bones suitably stained. Fish were destained in 1% KOH until background stain was lost and subsequently cleared in glycerol for analysis. For analysis of forming cartilage, fish were prestained in alcian blue from 4–24 hours. The fish were then destained, lightly trypsinized (3 g/l; 37°C) and processed for alzarin red staining.</p> <p>Skeletal measurements were made using digitizing software from Zeiss using a dissecting microscope. Measurements were made from fixed landmarks on each axis of the skull that did not vary depending on position of the suspensorium: the premaxilla was used for the distal most point on the length (L) axis, while the quadrate-anguloarticular joint was used as a ventral landmark for the height (H) axis. Raw measurement values are represented as normalized ratios of the distance along each axis in relation to the position of the center of the eye; values are normalized for standard length of the fish.</p> <title>Whole Mount In Situ Hybridisation and Immunostaining</title> <p>Probes for whole mount <italic>in situ</italic> hybridisation were generated by reverse transcription from cDNA made from regenerating caudal fin tissue. Digoxigenin labeled RNA probes were purified using P-30 micro bio-spin columns before use (BioRad). WMISH protocol was performed as described ##UREF##8##[65]## , at 70°C and with the addition of 0.1% CHAPS to hybridization and post hybridization wash buffers. Reactions were stopped in PBS, post-fixed and placed in methanol overnight to reduce non-specific staining.</p> <p>Acridine orange (Sigma) was used as a marker of apoptosis in developing tissue ##REF##8223253##[32]##. Juveniles were immersed in fish water containing 5 μg/ml acridine orange for 5 minutes and then washed with fish water, anesthetized and post-fixed in formalin to assist visualization of staining.</p> <title>Statistical Methods</title> <p>Analysis of cranial measurements were performed using Hotelling's T squared test for two dependent variables. For scale counts and size dimensions, a <italic>t</italic>-statistic for differential means was used to assess significance. Calculations and probability assessment were calculated using Biosoft 200 software (<ext-link ext-link-type="uri" xlink:href="http://www.biosoft.com">www.biosoft.com</ext-link>) and Excel statistical package.</p> <title>Electron Microscopy</title> <p>Specimens of 8–9 mm juvenile fish were fixed with a mixture of 4% formaldehyde in PBS and 1–2.5% glutaraldehyde at room temperature and subsequently placed at 4°C. After post-fixation with 1% osmium tetroxide in 100 mM PBS for 1 h on ice, samples were washed with H<sub>2</sub>O, treated with 1% aqueous uranyl acetate for 1 h at 4°C, dehydrated through a graded series of ethanol and embedded in Epon. Ultrathin sections were stained with uranyl acetate and lead citrate and viewed in a Philips CM10 electron microscope. In addition, toluidine blue stained Epon sections of 0.5 or 3 μm thickness were prepared for light microscopy.</p>
<p>Conceived and designed the experiments: MPH. Performed the experiments: MPH NR HS SP. Analyzed the data: MPH NR HS SP PK. Contributed reagents/materials/analysis tools: MPH NR. Wrote the paper: MPH NR CNV.</p> <p>The genetic basis of the development and variation of adult form of vertebrates is not well understood. To address this problem, we performed a mutant screen to identify genes essential for the formation of adult skeletal structures of the zebrafish. Here, we describe the phenotypic and molecular characterization of a set of mutants showing loss of adult structures of the dermal skeleton, such as the rays of the fins and the scales, as well as the pharyngeal teeth. The mutations represent adult-viable, loss of function alleles in the <italic>ectodysplasin</italic> (<italic>eda</italic>) and <italic>ectodysplasin receptor</italic> (<italic>edar</italic>) genes. These genes are frequently mutated in the human hereditary disease hypohidrotic ectodermal dysplasia (HED; OMIM 224900, 305100) that affects the development of integumentary appendages such as hair and teeth. We find mutations in zebrafish <italic>edar</italic> that affect similar residues as mutated in human cases of HED and show similar phenotypic consequences. <italic>eda</italic> and <italic>edar</italic> are not required for early zebrafish development, but are rather specific for the development of adult skeletal and dental structures. We find that the defects of the fins and scales are due to the role of Eda signaling in organizing epidermal cells into discrete signaling centers of the scale epidermal placode and fin fold. Our genetic analysis demonstrates dose-sensitive and organ-specific response to alteration in levels of Eda signaling. In addition, we show substantial buffering of the effect of loss of <italic>edar</italic> function in different genetic backgrounds, suggesting canalization of this developmental system. We uncover a previously unknown role of Eda signaling in teleosts and show conservation of the developmental mechanisms involved in the formation and variation of both integumentary appendages and limbs. Lastly, our findings point to the utility of adult genetic screens in the zebrafish in identifying essential developmental processes involved in human disease and in morphological evolution.</p> <title>Author Summary</title> <p>A major goal of the study of developmental genetics is to understand the genes and developmental mechanisms underlying the formation of organismal complexity and diversity. Here, we focus on genes controlling postembryonic development and describe mutations in genes of the ectodysplasin (Eda) pathway in regulating the formation of the scales, skull, fins, and teeth. Mutations in genes of this signaling pathway are common in humans with defects in ectodermal structures such as hair, glands, and teeth. We show that the similar phenotypes of loss of Eda signaling in fish and human are due to a conserved early developmental stage in the development of mammalian hair and fish scales; subsequent development of these two structures diverge. Our findings show that the Eda signaling pathway has an ancestral role in regulating the developmental interactions involved in patterning and growth of the dermal skeleton of fish. Recent work has shown that these genes are associated with morphological variation between humans and evolution within fish populations, suggesting that alteration in the function of these genes permits viable morphological change. Our data support the value of forward genetic studies on postembryonic development to reveal the genetic and developmental basis of both human disease and morphological evolution.</p>
<title>Supporting Information</title>
<p>The authors would like to thank Drs. M. Sonawane, K. Siegfried, C. Koehler, M. Levesque, and JY Sire for helpful comments on the manuscript, Dr. M. Hargrave for contributing the <italic>fls</italic>
<sup>t0sp213</sup> background allele, M. Akimenko for <italic>msxa</italic> cDNA, and Dr. Vu Nguyen for discussions and sharing of results before publication. We are grateful to Dr. Hiroshi Mitani in his gracious offer to provide medaka/<italic>rs-3</italic> mutant for anatomical analysis. In addition, the authors would like to thank the assistance of Dr. Robert Geisler and Ines Gehring for radiation hybrid analysis, Jennifer Zenker for assistance with characterization of the <italic>fls</italic>
<sup>te370f</sup> mutant, and Brigitte Sailer and Iris Koch for their expertise in histology. Lastly, we would like to thank the assistance and advice of three anonymous reviewers.</p>
<fig id="pgen-1000206-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000206.g001</object-id><label>Figure 1</label><caption><title>The formation of the adult dermal skeleton and pharyngeal teeth is affected in <italic>fls</italic> mutant zebrafish.</title><p>A) Alizarin red-stained wild type adult zebrafish shows staining of the scales, fin rays, dermal bones of the skull as well as the pharyngeal teeth along ceratobranchial 5 (bracket, B) and (C) gill rakers along both the anterior (antGr) and posterior edge (postGr) of non-teeth bearing ceratobranchials. D) <italic>fls</italic>
<sup>te370f</sup> shows loss of dermal skeletal structures of the fin rays, scales and alteration in the shape of the skull. Additionally, <italic>fls</italic>
<sup>te370f</sup> shows loss of pharyngeal teeth (E) and gill rakers (F). G–I) The Topless allele (<italic>fls<sup>dt3Tpl</sup></italic>) shows a dominant effect on scalation and tooth/gill raker formation while not affecting lepidotrichial growth. J–L) Expressivity of <italic>fls</italic>
<sup>dt3Tpl</sup> is sensitive to a modifier in the Tü strain leading to a “weak” <italic>fls</italic>
<sup>dt3Tpl</sup> phenotype; <italic>fls</italic>
<sup>dt3Tpl</sup> homozygotes were phenotypically identical to <italic>fls</italic>
<sup>te370f</sup> (not shown). The fang allele of <italic>fls</italic>, <italic>fls</italic>
<sup>tfang</sup>, isolated in a non-complementation screen with <italic>fls</italic>
<sup>te370f</sup>, shows no effect on fin development while exhibiting partial loss of scales (M), teeth (N), and gill rakers (O). Transallelic tfang/te370f zebrafish exhibit an intermediate phenotype between homozygous <italic>fls</italic>(te370f) and <italic>fls</italic>(tfang) (P–R).</p></caption></fig> <fig id="pgen-1000206-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000206.g002</object-id><label>Figure 2</label><caption><title>The dominant gene <italic>Nkt</italic> is phenotypically similar, however complements <italic>fls</italic> mutants.</title><p>
<italic>Nkt</italic> homozygotes show complete loss of scales, teeth and gill rakers resembling the <italic>fls</italic> phenotype (A–C). Heterozygous <italic>Nkt</italic> zebrafish show an intermediate phenotype of scale loss and patterning defect (arrows) while no effect on fin development is seen (D). Heterozygous <italic>Nkt</italic> also show a dominant effect on the number of teeth (arrows, E) and gill rakers (F), showing deficiencies along the posterior branchial arches and formation of rudimentary rakers along ceratobranchial 1 and 2 (arrows, F). <italic>Cb1-5</italic>, ceratobranchial bones.</p></caption></fig> <fig id="pgen-1000206-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000206.g003</object-id><label>Figure 3</label><caption><title>
<italic>fls</italic> and <italic>Nkt</italic> are mutations in genes encoding ectodysplasin receptor (<italic>edar</italic>) and its ligand ectodysplasin (<italic>eda)</italic>.</title><p>A) Mapping of <italic>fls</italic> using SSLP markers and placement of the <italic>edar</italic> gene within the candidate region on LG 9 by radiation hybrid mapping. The insert shows genetic linkage of the <italic>fls</italic> gene to local markers on LG 9. The numbers on the right of the insert indicate the number of recombinants seen in identified mutants per the number of meioses tested. B) Schematic of wild type Edar protein and mutant alleles. Polymorphisms seen in the WIK strain are noted above the wild type gene. Mutations that lead to premature termination are represented as truncated proteins showing the predicted residual fragment and position of the identified mutation. C) Analysis of the mutation in <italic>fls</italic>
<sup>dfang</sup>. A unique splice donor site (red) is generated leading to inclusion of additional coding sequence encoding a premature termination codon (underlined). D) Quantitative analysis of different <italic>edar</italic> transcript levels in <italic>fls</italic>
<sup>dfang</sup> homozygotes compared with wildtype. E) Similarity of altered residues in <italic>fls</italic>
<sup>t3R367W</sup>and <italic>fls</italic>
<sup>dt3Tpl</sup> with human HED shown in the death domain. The position of the mutated residues in <italic>fls</italic>
<sup>t3R367W</sup> (blue box) and in <italic>fls</italic>
<sup>dt3Tpl</sup> (red box) is identical to ones changed in cases of human autosomal dominant HED although the substitution is different. F) Linkage between the <italic>Nkt</italic> allele and <italic>eda</italic> on LG5. G) Schematic of wild type Eda protein and position of <italic>Nkt</italic> mutation. Numbers on gene diagrams represent amino acid length. <italic>TNF</italic>, tumor necrosis factor domain; <italic>TNFR</italic>, tumor necrosis factor receptor domain; <italic>TM</italic>, transmembrane domain; <italic>DD</italic>, death domain; the blue box in <italic>Eda</italic> is the furin binding site.</p></caption></fig> <fig id="pgen-1000206-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000206.g004</object-id><label>Figure 4</label><caption><title>The role of <italic>edar</italic> in expression of developmental patterning genes during early scale development.</title><p>Expression of <italic>edar</italic> (A, B) and <italic>eda</italic> (C, D) in early forming scales; arrowheads indicate site of expression of initial forming scales. A) <italic>edar</italic> expression above site of scale formation in 8 mm long (approximately 30 dpf juvenile fish) and in larger juveniles (9 mm; 30 dpf). C) <italic>eda</italic> expression during early scale development on the flank (8 mm) and in forming scales of older juvenile fish (10 mm, D). Expression of developmental genes <italic>bmp2b</italic> and <italic>shh</italic> in early scale development in wildtype (E, G) and <italic>fls</italic>
<sup>te370f</sup> (F, H) juveniles (9 mm).</p></caption></fig> <fig id="pgen-1000206-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000206.g005</object-id><label>Figure 5</label><caption><title>Eda signaling regulates the formation of an epidermal placode during scale development.</title><p>Histological analysis of wild type (A, D, B, E) and <italic>fls</italic>
<sup>te370f</sup> (C, F, I) integument of 8 mm standard length. In wild type juveniles (B, E), basal epidermal cells (black arrow heads) show a heightened, and cuboidal morphology at sites of scale development as indicated by an accumulation of migrating fibroblast-like cells (white arrowheads). (H) This morphology of the epidermis is associated with a reworking of the collagen layer of the stratum compactum (<italic>cpt</italic>; ##REF##9183678##[24]##). This is in contrast to the flattened morphology of basal epidermal cells lateral to those of the scale placode (A, D) and underlying dense stratum compactum (G). In <italic>fls</italic>
<sup>te370f</sup> this basal epidermal structure is disorganized and cell morphology is disrupted (C, F) including evidence of cell death (asterisk). The lack of reworking of the collagen of the stratum compactum in the <italic>fls</italic>
<sup>te370f</sup> mutant is associated with retention of hemidesmosomes (horizontal bracket G–I). <italic>edar</italic> is expressed in cells of the wildtype epidermis (J, K, L). Counterstaining of the same sections confirms the expression in basal cells overlying initial accumulating fibroblasts (white arrowheads; J′, K′, L′). Expression of <italic>edar</italic> is observed prior to organization of the placode and fibroblast aggregation and maintained in cells of the epidermal placode through early scale development (J–L). M) Schematic depicting scale development and <italic>edar</italic> expression. The stages of scale development are modeled using analogous stages as described for hair development ##REF##16481354##[35]##; stage 0, nascent epidermis; stage 1, placode specification; stage 2, scale pocket; stage 3, matrix deposition and ossification. Blue, <italic>edar</italic> expression; red, scale formation. <italic>ep</italic>, epidermis; <italic>cpt</italic> stratum compactum. The vertical bracket demarcates the extent of the epidermis in the sections. Measurement bar equals 10 μm.</p></caption></fig> <fig id="pgen-1000206-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000206.g006</object-id><label>Figure 6</label><caption><title>Fin development is defective in <italic>fls</italic> and <italic>Nkt</italic> mutant zebrafish.</title><p>Alizarin red stained adult zebrafish fins show a drastic effect of <italic>fls</italic>
<sup>te370f</sup> and <italic>Nkt</italic> on development of the lepidotrichial dermal rays of both the paired and unpaired fins. A–C,F) Pectoral fins, anterior-dorsal view; D, E) double staining developing pectoral fins with alcian blue and alizarin red show early patterning of the endochondrial bones of the pectoral fin of size matched wild type and <italic>Nkt</italic> homozygotes (asterisk indicates loss of fourth proximal radial). G) effect of <italic>fls</italic> and <italic>Nkt</italic> mutants on the patterning of the pectoral fin skeleton scored as the number of specimens showing alteration in pattern or form over total analyzed. The identity of the proximal radials is noted (I–IV). H–J) Analysis of pelvic fin development in <italic>fls</italic> and <italic>Nkt</italic> mutants. Numbers denote anterior-posterior identity of the lepidotrichia. K–M) Defects in the formation of the lepidotrichia in adult anal and (N–P) caudal fin of <italic>fls</italic> (L and O) and <italic>Nkt</italic> (M and P). <italic>ap</italic>, ascending process; <italic>cl</italic>, cleithrum; <italic>co</italic>, corticoid; <italic>dr</italic>, distal radial; <italic>sc</italic>, scalpula, <italic>le</italic>, lepidotrichia; <italic>pcl</italic>, postcleithrum; <italic>pr</italic>, proximal radial.</p></caption></fig> <fig id="pgen-1000206-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000206.g007</object-id><label>Figure 7</label><caption><title>Eda signaling and the maintenance of anterior-posterior pattern in late paired fin development.</title><p>Analysis of <italic>edar</italic> (A–B, E–F), <italic>shh</italic> (C–D, G–H), and <italic>bmp2b</italic> (I–L) expression in developing pectoral (A–D, –J) and pelvic fins (E–H, K–L) from 8 mm juvenile fish of wild type (A, E; C, G; I, K) and <italic>fls</italic>
<sup>te370f</sup> mutant fish (B, F; D, H; J, L). A–B) Arrowheads indicate two distinct patterns of <italic>edar</italic> expression in the pectoral fin: an expression that marks the posterior edge and distal region of the development of the proximal radials (black); and a posterior bias of <italic>edar</italic> expression in the forming lepidotrichia (white). Arrows point out the remaining posterior expression in mutant fins. Brackets in all panels outline anterior deficiencies in gene expression in <italic>fls</italic> mutant fins. M–P, analysis of patterns of cell death in the developing paired fins by retention of acridine orange stain. N, P) Arrows point out anterior distal regions of cell death in both pectoral and pelvic fins from the mutant; (M, N) pectoral fin and (O, P) pelvic fin respectively. Q, R) Expression of <italic>msxa</italic> in wild type and mutant pectoral fins. Region of expression outlined with brackets; asterisk marks an ectopic site of expression.</p></caption></fig> <fig id="pgen-1000206-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000206.g008</object-id><label>Figure 8</label><caption><title>Eda signaling is required for the function of the fin fold during late fin development.</title><p>Expression of <italic>bmp2b</italic> (A), <italic>shh</italic> (B) and <italic>edar</italic> (C) transcripts in developing juvenile (8 mm, 30 dpf) fin rays of the caudal fin; asterisks indicate regional expression within distal tip of developing ray; arrows in (C), expression in distal epidermis of the fin fold. D–I) Histological analysis of both paired (pectoral, D, E) and unpaired fins (anal F, G; and caudal, H, I) from <italic>fls</italic>
<sup>te370f</sup> and wild type siblings. <italic>fls</italic>
<sup>te370f</sup> fins showed a general deficiency in the maturation of the muscle and dermis of the fin (arrow G; acellular debris in anal fin of <italic>fls</italic>). Insets (D, E), tip of fin at higher magnification showing degeneration of the nuclei of the epidermis in the mutant fin. <italic>ffd</italic>, fin fold; <italic>le</italic>, lepidotrichia of the fin rays.</p></caption></fig>
<table-wrap id="pgen-1000206-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000206.t001</object-id><label>Table 1</label><caption><title>Quantitative effect of <italic>fls</italic> on scale number and shape and the effect of background modifiers in <italic>Danio rerio</italic> strains on <italic>fls</italic>
<sup>dt3Tpl</sup>.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Phenotype/Genotype</td><td align="left" rowspan="1" colspan="1">Scale #/ stl</td><td align="left" rowspan="1" colspan="1">n</td><td align="left" rowspan="1" colspan="1">Scale DV/AP</td><td align="left" rowspan="1" colspan="1">n</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">fish</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">scales</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">+/+</td><td align="left" rowspan="1" colspan="1">6.8±0.18</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">1.14±0.15</td><td align="left" rowspan="1" colspan="1">13</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>fls</italic>
<sup>dt3Tpl</sup> / Tü</td><td align="left" rowspan="1" colspan="1">3.0±0.20 ##</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1.52±0.29 #</td><td align="left" rowspan="1" colspan="1">8</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>fls</italic>
<sup>dt3Tpl</sup> / Tü; <italic>mod</italic>
</td><td align="left" rowspan="1" colspan="1">5.6±0.44 #</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">1.4±0.3 #</td><td align="left" rowspan="1" colspan="1">12</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>fls</italic>
<sup>dt3Tpl</sup> / WIK</td><td align="left" rowspan="1" colspan="1">5.84±0.66 #</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">1.43±0.35 #</td><td align="left" rowspan="1" colspan="1">32</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>fls</italic>
<sup>tfang</sup> / <italic>fls<sup>t</sup></italic>
<sup>fang</sup>
</td><td align="left" rowspan="1" colspan="1">0.97±0.50 ###</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1.57±0.18 #</td><td align="left" rowspan="1" colspan="1">7</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>fls</italic>
<sup>te370f</sup> / <italic>fls</italic>
<sup>te370f</sup>
</td><td align="left" rowspan="1" colspan="1">0.41±0.39 ###</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">1.8±0.64 #</td><td align="left" rowspan="1" colspan="1">16</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pgen.1000206.s001"><label>Figure S1</label><caption><p>
<italic>fls</italic> and <italic>Nkt</italic> alter size and proportion of the adult zebrafish skull. In addition to variations in integumentary structures, <italic>fls</italic> and <italic>Nkt</italic> exhibited a distinct change of the shape and size of the adult skull. Measurements of the absolute proportions of the adult skull, normalized for overall growth of the fish as determined by standard length, demonstrate that both <italic>fls</italic> and <italic>Nkt</italic> homozygous mutations result in overall larger skulls of the fish (<italic>fls</italic>
<sup>te370f/te370f</sup>, n = 16, T2 = 23.7, p&lt;0.001; <italic>Nkt</italic>, n = 10, T2 = 64.1, p&lt;0.001; <italic>fls</italic>
<sup>te370f/dt3Tpl</sup>, n = 6, T2 = 22.8, p&lt;0.005). The dominant effect of <italic>Nkt</italic> and <italic>fls</italic>
<sup>dt3Tpl</sup> seen in development of the scale pattern was not observed in the formation of skull size. However, an analysis of changes in the proportional development of the skull by measurements of the relative positioning of the eye within the skull (L1/L2, H1/H2; Panel A) showed a significant and dominant effect of <italic>Nkt</italic>, <italic>fls</italic>
<sup>dt3Tpl</sup> on the patterning of the skull (Panel C). This effect was seen in <italic>fls</italic>
<sup>te370f</sup> homozygotes as well and was not specific to particular alleles of <italic>fls</italic>. The alteration in skull size and shape in the mutants does not involve loss of a particular organ structure or specific bone, rather a change in proportions of the developing skull.</p><p>(3.39 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000206.s002"><label>Figure S2</label><caption><p>A comparison of <italic>edar</italic> sequence in representative vertebrates. Edar alleles <italic>fls</italic>
<sup>dt3Tpl</sup> and <italic>fls</italic>
<sup>t3R367W</sup> positioned above sequence. Sites of splicing defects of <italic>fls</italic>
<sup>te370f</sup> and <italic>fls</italic>
<sup>t0sp213</sup> alleles demarcated with ˆ marker. Yellow, TNFR domain; Grey, transmembrane domain; Green, death domain; Red, polymorphic sites in WIK mapping strain.</p><p>(0.04 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000206.s003"><label>Figure S3</label><caption><p>A comparison of <italic>eda</italic> sequence of representative vertebrates. Blue, transmembrane domain; Green, furin cleavage site; Yellow, TNF domain; asterisk <italic>Nkt</italic>
<sup>dtS238X</sup> allele; |, deleted residues in alternate spliced form of Eda-2.</p><p>(0.04 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000206.s004"><label>Figure S4</label><caption><p>Scale formation and variation in the rs3/<italic>edar</italic> medaka mutant on the cs-2 background. (A) Alizarin-red stained rs3 medaka showed substantial scale formation and variation of the extent of scalation. (B) Wild type cs-2 strain scalation pattern.</p><p>(5.72 MB TIF)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label/><p>The total number of scales on one side of alizarin red stained adults of different genotypes were counted and measured. Counts were normalized for standard length (stl) of individual fish as shape and number of scales in the mutants may vary as a measure of size. Shape characteristics of scales were quantified by measuring three to four scales from set positions across the flank of each fish and comparing the height (dorsal-ventral; DV) to length (anterior-posterior; AP) ratios. Results are presented as sample average and standard deviation around the mean. <italic>mod</italic>, inferred genotype of a modifier in Tü background leading to “weak” phenotype. The numerical symbol (#) indicates significant difference compared to wild type values (students <italic>t</italic>, p&lt;0.05). The different number of symbols signifies a significantly different phenotypic classes of scale development (#, ##, ###).</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>MPH was supported by a Max Planck Postdoctoral Fellowship.</p></fn></fn-group>
<graphic xlink:href="pgen.1000206.g001"/> <graphic xlink:href="pgen.1000206.g002"/> <graphic id="pgen-1000206-t001-1" xlink:href="pgen.1000206.t001"/> <graphic xlink:href="pgen.1000206.g003"/> <graphic xlink:href="pgen.1000206.g004"/> <graphic xlink:href="pgen.1000206.g005"/> <graphic xlink:href="pgen.1000206.g006"/> <graphic xlink:href="pgen.1000206.g007"/> <graphic xlink:href="pgen.1000206.g008"/>
<media xlink:href="pgen.1000206.s001.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000206.s002.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000206.s003.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000206.s004.tif"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["2"], "element-citation": ["\n"], "surname": ["Stern"], "given-names": ["DL"], "year": ["2000"], "article-title": ["Evolutionary developmental biology and the problem of variation."], "source": ["Evolution Int J Org Evolution"], "volume": ["54"], "fpage": ["1079"], "lpage": ["1091"]}, {"label": ["7"], "element-citation": ["\n"], "surname": ["Bell", "Foster"], "given-names": ["MA", "SA"], "year": ["1994"], "source": ["The Evolutionary Biology of the Threespine Stickleback"], "publisher-loc": ["Oxford"], "publisher-name": ["Oxford University Press"], "fpage": ["571"]}, {"label": ["45"], "element-citation": ["\n"], "surname": ["Paakkonen", "Cambiaghi", "Novelli", "Ouzts", "Penttinen"], "given-names": ["K", "S", "G", "LV", "M"], "year": ["2001"], "article-title": ["The mutation spectrum of the EDA gene in X-linked anhidrotic ectodermal dysplasia."], "source": ["Hum Mutat"], "volume": ["17"], "fpage": ["349"]}, {"label": ["59"], "element-citation": ["\n"], "surname": ["Mou", "Thomason", "Willan", "Clowes", "Harris"], "given-names": ["C", "HA", "PM", "C", "WE"], "year": ["2008"], "article-title": ["Enhanced ectodysplasin-A receptor (EDAR) signaling alters multiple fiber characteristics to produce the East Asian hair form."], "source": ["Hum Mutat. "], "italic": ["epub"]}, {"label": ["60"], "element-citation": ["\n"], "surname": ["Stock"], "given-names": ["DW"], "year": ["2007"], "article-title": ["Zebrafish dentition in comparative context."], "source": ["J Exp Zoolog B Mol Dev Evol"], "volume": ["308"], "fpage": ["523"], "lpage": ["549"]}, {"label": ["61"], "element-citation": ["\n"], "surname": ["Pelegri", "Nuesslein-Volhard", "Dahm"], "given-names": ["F", "C", "R"], "year": ["2002"], "article-title": ["Mutagenesis."], "source": ["Zebrafish"], "publisher-loc": ["Oxford"], "publisher-name": ["Oxford University Press"], "fpage": ["145"], "lpage": ["174"]}, {"label": ["62"], "element-citation": ["\n"], "surname": ["Haffter", "Granato", "Brand", "Mullins", "Hammerschmidt"], "given-names": ["P", "M", "M", "MC", "M"], "year": ["1996"], "article-title": ["The identification of genes with unique and essential functions in the development of the zebrafish, Danio rerio."], "source": ["Dev"], "volume": ["12"], "fpage": ["1"], "lpage": ["36"]}, {"label": ["63"], "element-citation": ["\n"], "surname": ["Geisler", "Nuesslein-Volhard", "Dahm"], "given-names": ["R", "C", "R"], "year": ["2002"], "article-title": ["Mapping and Cloning."], "source": ["Zebrafish"], "publisher-loc": ["Oxford"], "publisher-name": ["Oxford University Press"], "fpage": ["175"], "lpage": ["212"]}, {"label": ["65"], "element-citation": ["\n"], "surname": ["Shulte-Merker", "Nuesslein-Volhard", "Dahm"], "given-names": ["S", "C", "R"], "year": ["2002"], "article-title": ["Looking at Embryos."], "source": ["Zebrafish"], "publisher-loc": ["Oxford"], "publisher-name": ["Oxford University Press"], "fpage": ["39"], "lpage": ["58"]}]
PLoS Genet. 2008 Oct 3; 4(10):e1000206