Abstract

Birds display remarkable diversity in the distribution and morphology of scales and feathers on their feet, yet the genetic and developmental mechanisms governing this diversity remain unknown. Domestic pigeons have striking variation in foot feathering within a single species, providing a tractable model to investigate the molecular basis of skin appendage differences. We found that feathered feet in pigeons result from a partial transformation from hindlimb to forelimb identity mediated by cis-regulatory changes in the genes encoding the hindlimb-specific transcription factor Pitx1 and forelimb-specific transcription factor Tbx5. We also found that ectopic expression of Tbx5 is associated with foot feathers in chickens, suggesting similar molecular pathways underlie phenotypic convergence between these two species. These results show how changes in expression of regional patterning genes can generate localized changes in organ fate and morphology, and provide viable molecular mechanisms for diversity in hindlimb scale and feather distribution.

https://doi.org/10.7554/eLife.12115.001

eLife digest

Animals ranging from fish to birds display dramatic diversity within and among species; yet remarkably little is known about the genetic and developmental mechanisms that underlie this variation. In birds and their extinct dinosaur relatives, the distribution of scales and feathers on the feet is a highly variable trait.

Different breeds of domestic pigeon all belong to the same species but have feet that can be feathery or scaly to different extents. Classical genetics experiments suggested that only a few genes are involved in the transition from scaled to feathered skin on the feet of pigeons. However, the molecular basis for this transition was unknown.

Domyan et al. set out to identify the genes involved in the transition from scaled to feathered feet by mating different breeds of pigeon in the laboratory and then sequencing the birds’ DNA. They also surveyed the entire DNA sequences of many additional pigeon breeds with and without feathered feet. This combined approach showed that two regions of the pigeon genome have a profound effect on the number and size of feathers on the feet of domestic pigeons. These regions contain genes that are known to play key roles in controlling the growth of a limb and whether it develops into a leg or a wing. In developing pigeon embryos, Domyan et al. found that a gene called Pitx1, which is typically considered a hindlimb gene, is expressed at lower levels in the developing legs of breeds with feathered feet than in a breed with scaled feet. The experiments also showed that Tbx5, a gene that is expressed in the forelimbs of many animals, is expressed abnormally in the embryonic hindlimbs of breeds of pigeon and chicken with feathery feet.

Together, these findings suggest that the hindlimbs of domestic birds with feathery feet are more like wings at the molecular level, which results in them being covered in feathers rather than scales. Future work will now aim to discover the specific DNA sequences that alter the expression of Pitx1 and Tbx5 in feather-footed breeds, and whether the same genes control the foot feathers of other species of birds.

https://doi.org/10.7554/eLife.12115.002

Introduction

In birds, the genetic and developmental mechanisms that control the decision between scale and feather development remain poorly understood. Most birds possess scales on the foot (tarsometatarsus and toes) and feathers elsewhere. Exceptions to this pattern can provide insights into the evolutionary and developmental basis of skin appendage diversity. Some raptors and boreal birds evolved foot feathers instead of scales ('ptilopody'; Danforth, 1919; Lucas and Stettenheim, 1972), but the lack of appendage variation within these species precludes their use as genetic models. Likewise, paravians (birds and their close non-avian theropod dinosaur relatives) and other dinosaurs evolved diverse feather coverings on their legs and feet that sometimes resemble flight-like feathers (Xu et al., 2003; Hu et al., 2009; Turner et al., 2012; Zheng et al., 2013; Foth et al., 2014; Godefroit et al., 2014), but the absence of living specimens preclude mechanistic molecular studies.

In contrast, domestic pigeons (Columba livia) exhibit stunning variation within a single extant species (Shapiro and Domyan, 2013). Most breeds have feet covered by scaled epidermis (wild-type), but scales are replaced by small feathers in grouse (gr) mutants, and by larger feather 'muffs' in birds that also carry mutant alleles at the Slipper (Sl) locus (Doncaster, 1912; Wexelsen, 1934; Hollander, 1937; Levi, 1986) (Figure 1A). In muffed breeds, scutellate scales are generally absent or poorly developed on the feathered epidermis covering the tarsometatarsus and toes, and feathers are surrounded by soft integument. The molecular identities of both gr and Sl are unknown, and additional loci probably control quantitative variation in the muff phenotype. Because both scale-footed and feather-footed pigeon breeds belong to the same species, we can use traditional genetic crosses and whole-genome resequencing to map the genes that control this striking variation (Shapiro et al., 2013; Domyan et al., 2014). Therefore, we can study diversity of the magnitude usually observed among different species without the roadblock of hybrid incompatibility that often eliminates the possibility of genetic mapping studies.

Figure 1 with 1 supplement see all
Two QTL differentiate scale- and feather-footed domestic pigeons.

(A) Common phenotypes of domestic rock pigeon, in order of increasing foot feathering (left to right): scaled, groused, small- and large-muffed. (B-F) QTL scans and effect plots: proportion of tarsus feathered (B,C), number of toe feathers (D,E), and length of toe feathers (F,G). Mean phenotypes ± S.E. are plotted in (C,E,G). For (E) and (G), genotypes at the QTL with the higher LOD score are on the x-axis, and genotypes at the other QTL are inset. See Tables 1 and 2 for detailed QTL statistics. S, allele from scale-footed grandparent; F, allele from feather-footed grandparent.

https://doi.org/10.7554/eLife.12115.003

During development in vertebrates, skin appendages form through interactions between the ectoderm-derived epidermis and the mesoderm-derived dermis, and signals from the dermis determine epidermal appendage fate (Dhouailly, 2009; Hughes et al., 2011). Previous analyses of mutants and gene misexpression in chickens suggest candidates for feathered feet in the Hedgehog (Harris et al., 2002), BMP (Zou and Niswander, 1996; Harris et al., 2002; 2004), Delta-Notch (Crowe et al., 1998), and Wnt (Chang et al., 2004) pathways. Our study of genetic variation and embryonic development in pigeons, however, reveals a surprisingly different mechanism with broad implications for limb identity and patterning.

Results

Two loci of major effect control foot feathering in pigeons

To identify chromosome regions that contribute to feathered feet, we generated an F2 intercross between muffed (Pomeranian pouter) and scaled (Scandaroon) breeds (Figure 1—figure supplement 1). F1 hybrids displayed intermediate foot feathering, demonstrating a semi-dominant inheritance pattern. Scaled, muffed, and intermediate phenotypes were recovered in the F2 population, confirming that a small number of major-effect loci contribute to this trait. Among F2 offspring, digit 3 bore the largest and greatest number of feathers (digit 1, 7.04 ± 5.08 feathers; digit 2, 7.64 ± 7.54; digit 3, 15.46 ± 10.43; digit 4, 6.13 ± 5.70). Using quantitative trait locus (QTL) mapping with 130 F2 offspring genotyped at 3803 polymorphic markers (Broman et al., 2003), we identified two linkage groups (LG11 and LG20) that had significant effects on three different aspects of foot feathering (log10 odds ratio (LOD) > 4.6; Figure 1B–G, Table 1). LG11 had the largest effects on the proportion of the tarsometatarsal epidermis that was transformed from scaled to feathered (LOD = 7.69, percent variance explained (PVE) = 28.4%) and toe feather number (LOD = 6.72, PVE = 21.3%), and a smaller effect on toe feather length (LOD = 8.51, PVE = 15.8%). LG20 had the largest effect on toe feather length (LOD = 20.9, PVE = 52.2%), and a smaller effect on toe feather number (LOD = 5.36, PVE = 16.5%). When toe feather number was analyzed for each digit individually, the same two QTL were identified and had the most pronounced effects on digits 3 and 4 (Table 2). In summary, two major QTL have marked and separable effects on qualitative and quantitative variation in epidermal appendages.

Table 1

Summary of foot feathering QTL.

https://doi.org/10.7554/eLife.12115.005
TraitLGLoc (cM)ScaffoldPositionLODPVEMean ± S.D.
SSSFFF
Proportion tarsus feathered11117799,205,2867.6928.40.46 ± 0.040.58 ± 0.03**0.80 ± 0.04***
Number of toe feathers111247912,325,9776.7321.343.6 ± 8.967.3 ± 6.5*105.3 ± 8.1***
Number of toe feathers2015951,451,1275.3616.545.3 ± 8.178.2 ± 6.5**100.3 ± 9.9***
Toe feather length (mm)20070136,74620.8952.25.2 ± 2.319.9 ± 2.0***37.3 ± 2.5***
Toe feather length (mm)111247912,325,9778.5115.811.4 ± 3.318.7 ± 2.4*28.5 ± 3.0**
  1. LG, linkage group; Loc, genetic location of peak LOD score in centimorgans; S, allele from scaled parent; PVE, percent variance explained; F, allele from feathered parent; LOD, log10 odds ratio. (Welch two sample t-test of means compared to SS homozygote; *p≤ 0.05, **p≤0.005, ***p≤0.0005.)

Table 2

Summary of QTL for numbers of feathers on individual toes.

https://doi.org/10.7554/eLife.12115.006
DigitLGLoc (cM)ScaffoldPositionLODPVEMean ± S.D.
SSSFFF
Digit 2, left foot111487911,624,7015.2020.243.75 ± 1.136.95 ± 1.12*12.19 ± 1.16***
Digit 3, right foot20070136,7467.7124.19.06 ± 1.7017.35 ± 1.48**19.86 ± 1.83***
Digit 3, right foot111487911,624,7016.4719.59.98 ± 1.6015.60 ± 1.62*21.02 ± 1.65***
Digit 3, left foot20070136,74610.7929.848.66 ± 1.5917.06 ± 1.41**21.09 ± 1.86***
Digit 3, left foot111487911,624,7019.4825.529.37 ± 1.5216.26 ± 1.52**21.30 ± 1.57***
Digit 4, right foot2032952,464,7887.3321.362.44 ± 0.907.24 ± 0.75**10.14 ± 1.10***
Digit 4, right foot11118795,475,4746.7819.53.17 ± 0.896.91 ± 0.92*9.95 ± 1.00***
Digit 4, left foot2030952,464,7889.3726.832.44 ± 0.796.63 ± 0.66**9.65 ± 0.96***
Digit 4, left foot111487911,624,7017.5820.833.75 ± 1.146.95 ± 1.13*12.19 ± 1.17***
  1. LG, linkage group; Loc, genetic location of peak LOD score in centimorgans; S, allele from scaled parent; PVE, percent variance explained; F, allele from feathered parent; LOD, log10 odds ratio. (Welch two sample t-test of means compared to SS homozygote; *p≤0.05, **p≤0.005, ***p≤0.0005.)

In parallel to our experimental cross, we performed probabilistic whole-genome scans of allele-frequency differentiation (pFst; see Kronenberg et al., 2014) across a genetically and phenotypically diverse panel of breeds by comparing 15 feather-footed birds (4 groused and 11 muffed) to 28 scale-footed birds (Shapiro et al., 2013). Using this independent approach across breeds, the two most-highly differentiated pFst signals implicate the same genomic regions as the QTL study: genomic scaffold 79 is located on LG11 (p=4.44 x 10–16, genome-wide significance threshold = 2.11 x 10–9), and scaffolds 70 and 95 are adjacent to one another on LG20 (p=9.81 x 10–13) (Figure 2A, Figure 2—figure supplement 1A,B)

Figure 2 with 3 supplements see all
Two regions of genomic differentiation and H3K27ac enrichment distinguish scale- and feather-footed pigeons.

(A) Whole-genome pFst comparisons between genomes of feather-footed and scale-footed pigeons. Scaffolds are ordered by genetic position in a linkage map from an F2 cross (see Figure 1). Dashed line, genome-wide significance threshold. (B) pFst and extended haplotype homozygosity (EHH) plots for region of high differentiation on scaffold 79. Feather-footed birds (n=10, red in EHH plot) homozygous for a 44-kb deletion are differentiated from scale-footed birds (n=28, black) and show extended haplotype homozygosity in this region. Smoothed lines follow a generalized additive model (Wickham, 2009). (C) H3K27ac ChIP-seq enrichment differed significantly between embryonic wing and leg buds of the scale-footed racing homer (RH) in several regions (blue shading), including within a 44-kb interval that is deleted in the muffed Indian fantail (IF; blue arrowheads). This deleted region is orthologous to a known human limb enhancer (hs1473). (D) Selection scans show similar patterns of differentiation on scaffold 70 between muffed (n=11, red in EHH plot) and scale-footed birds (n=28, black). (E) H3K27ac ChIP-seq enrichment differed significantly between leg buds of racing homer and Indian fantail embryos in regions immediately 5’ of Tbx5 (blue shading). Foot and wing drawings modified after Levi (1986).

https://doi.org/10.7554/eLife.12115.007

The peak pFst region on scaffold 79 contained a 44-kb deletion (from 6.719–6.763 Mb) that was homozygous in 10, and heterozygous in 2 of the 15 feather-footed birds (Figure 2B,C; Figure 2—figure supplement 2A). Birds homozygous for the deletion showed elevated levels of haplotype homozygosity relative to scaled birds, a signature of positive selection on this region (Figure 2B). This deletion spans an element orthologous to a known human limb enhancer, hs1473 (Spielmann et al., 2012), which contains active chromatin marks (Cotney et al., 2012) and is bound by the hindlimb-specific transcription factor Pitx1 in the developing mouse hindlimb (Infante et al., 2013) (Figure 2C). The locus was homozygous for the deletion in 35 of 54 additional feather-footed birds from 21 breeds, but was never homozygous in 96 scale-footed birds from 56 breeds (Chi-square, p<0.0001; Figure 2—figure supplement 2B). The 3 feather-footed birds from our whole-genome panel that did not have this deletion (including the male founder of the aforementioned genetic cross implicating this same region) also showed allelic differentiation from scale-footed birds over this interval, suggesting that an additional feathered-foot allele may also exist at this locus (Figure 2—figure supplement 2A).

In contrast to the differentiation signal we observed between scale-footed and all feather-footed birds on scaffold 79, only the muffed birds (more heavily feathered) showed signatures of selection and shared similar haplotypes on scaffold 70 (higher pFst signal than the adjacent scaffold 95) (Figure 2D, Figure 2—figure supplements 2C, 3). Thus, both QTL analyses and whole-genome scans show that mutation on scaffold 79 alone is sufficient for the grouse phenotype (gr locus), and point to scaffold 70 as the major-effect locus for longer toe feathers in birds with muffs (Sl locus) (Figure 1F,G).

Expression of limb outgrowth and identity genes differs between scale-footed and feather-footed breeds

Next, we examined scaffolds 79 and 70 for candidate genes that might control the scale-to-feather transition. The highest pFst peak on scaffold 79 – corresponding to the major-effect QTL on LG 11 for the proportion of tarsometatarsal feathering – was approximately 200 kb upstream of Pitx1, a gene that encodes a homeobox-containing transcription factor that is normally expressed in the vertebrate hindlimb but not the forelimb (Figure 2B). The highest pFst peak on scaffold 70 – corresponding to the major-effect QTL on LG 20 for toe feather length – was 40 kb upstream of Tbx5, a gene that encodes a T-box transcription factor that is normally expressed in the vertebrate forelimb but not the hindlimb (Figure 2D). These regions were especially intriguing because these two genes encode key transcriptional regulators of forelimb (Tbx5) and hindlimb (Pitx1) identity and development (Logan et al., 1998; Logan and Tabin, 1999; Rodriguez-Esteban et al., 1999; Szeto et al., 1999; Takeuchi et al., 1999). For example, misexpression of Pitx1 in the embryonic chick forelimb blocks feather development (Logan and Tabin, 1999), while misexpression of Tbx5 in the early hindlimb field of embryonic chickens is sufficient to induce a partial wing-like transformation, including feather formation on the feet (Takeuchi et al., 1999). In mouse, Pitx1 (but not Tbx5) plays a role in determining limb-type identity (Szeto et al., 1999; Minguillon et al., 2005; DeLaurier et al., 2006), suggesting that the roles of Tbx5 in limb outgrowth and identity have diversified during amniote evolution (Horton et al., 2008).

We did not identify any fixed non-synonymous coding changes in Pitx1 or Tbx5 between scale-footed and feather-footed breeds of pigeon. However, we found striking differences in embryonic hindlimb expression of these genes among three different breeds – racing homer (scale-footed), Indian fantail (small-muffed), and English trumpeter (large-muffed) – at Hamburger-Hamilton stage 25 (HH25; Hamburger and Hamilton, 1951). Pitx1 expression was reduced in both muffed breeds (expression relative to racing homer: Indian fantail 0.75 ± 0.06, p=0.0007; English trumpeter 0.40 ± 0.05, p=0.0007; n = 6 each) and was more severely reduced in the large-muffed English trumpeter (p=0.002) (Figure 3A). Conversely, Tbx5, the forelimb-specific transcription factor, was ectopically expressed in the hindlimbs of both muffed breeds (hindlimb expression relative to racing homer forelimb: racing homer 0.001 ± 0.0004; Indian fantail 0.01 ± 0.008, p=0.0007; English trumpeter 0.14 ± 0.05, p=0.0007; n = 6 each), and was higher in the large-muffed English trumpeter (p=0.002) (Figure 3B). Forelimb expression of Tbx5 was indistinguishable among the three breeds, indicating that upregulation of Tbx5 in feather-footed breeds is restricted to the hindlimb (Figure 3—figure supplement 1A).

Figure 3 with 3 supplements see all
Limb-type gene expression varies among feathered and scaled pigeons and chickens.

(A,B,F,G) qRT-PCR analyses of Pitx1 and Tbx5 expression in HH25 hindlimbs of pigeon (A,B) and chicken (F,G). Boxes span 1st to 3rd quartiles, bars extend to minimum and maximum observed values if within 1.5 times the interquartile range of the box, circles indicate values outside of this range, black line indicates median. **=p<0.01, ***=p<0.001. (C-E, H-J) RNA in situ hybridization for Tbx5 expression in HH25 embryos of racing homer (C), Indian fantail (D), and English trumpeter (E) pigeons; and white leghorn (H), Cochin (I), and silkie (J) chickens. Arrowheads indicate sites of ectopic Tbx5 expression. Scale bar = 1 mm.

https://doi.org/10.7554/eLife.12115.011
Figure 3—source data 1

Source data from quantitative RT-PCR experiments.

Abbreviations: fl, embryonic forelimb; hl, embryonic hindlimb.

https://doi.org/10.7554/eLife.12115.012

We examined expression of additional genes within the two candidate regions at HH25, and found that the Tbx5 paralog Tbx3 was also differentially expressed in both feather-footed pigeon breeds relative to racing homer (Figure 3—figure supplement 1B). This could be due to the fact that Tbx3 is a target of Tbx5 (Mori et al., 2006; Postma et al., 2008), and additional experiments confirm that cis-regulatory changes do not drive this expression difference (see below). The hindlimb-specific transcription factor Tbx4 is not contained in the candidate regions defined by our QTL mapping and genome-wide association studies, but this gene is a downstream transcriptional target of Pitx1 (Logan and Tabin, 1999; Takeuchi et al., 1999; Duboc and Logan, 2011). We therefore compared expression levels of Tbx4 among scaled and feathered breeds at HH25, but found no significant differences at this stage (racing homer 1.00 ± 0.21; Indian fantail 1.00 ± 0.17, p=0.85; English trumpeter 1.05 ± 0.19, p=0.66; n = 6 each) (Figure 3—figure supplement 1C). Thus, the embryonic hindlimbs of muffed pigeons show quantitative expression changes in transcription factors with reciprocal limb expression domains, including the striking downregulation of a key hindlimb identity gene (Pitx1), and the novel expression of a key forelimb-specific gene (Tbx5).

We next analyzed the patterns of Tbx5 and Pitx1 expression in scaled and muffed pigeon embryos at HH25. Interestingly, ectopic hindlimb expression of Tbx5 in muffed embryos was markedly different than its normal forelimb pattern in wild-type pigeons and other vertebrates (Gibson-Brown et al., 1998a; 1998b; Logan et al., 1998; Tamura et al., 1999; Ruvinsky et al., 2000). Tbx5 is typically expressed throughout the mesoderm of the forelimb, but ectopic Tbx5 expression was largely localized to the mesoderm of the proximal and posterior-dorsal hindlimb of the small-muffed Indian fantail (Figure 3D, Figure 3—figure supplement 3B). This domain was further expanded in the large-muffed English trumpeter (Figure 3E, Figure 3—figure supplement 3C), consistent with the quantitative differences in expression between the two breeds (Figure 3B). This domain shows a striking correlation with regions of epidermal transformation, as foot feathers are usually longest and most numerous on the posterior digits (Darwin, 1868; Levi, 1986) (Figure 3—figure supplement 3D–N). In contrast, and despite quantitative expression differences among breeds, Pitx1 had a qualitatively similar expression domain in embryos of scale-footed and feather-footed breeds at this stage (Figure 3—figure supplement 2). Therefore, consistent with the critical role of mesoderm in determining ectodermal fate (Hughes et al., 2011), regionalized ectopic expression of Tbx5 is correlated with enhanced local transformation of epidermal appendages.

Cis-regulatory changes contribute to expression differences in Pitx1 and Tbx5 in muffed pigeons

If cis-acting regulatory mutations are responsible for the differences in Pitx1, Tbx5, and/or Tbx3 expression between embryos of scale-footed and feather-footed pigeons, then differential expression of scaled-foot and feathered-foot alleles should persist in a common trans-acting cellular environment. To test this prediction, we generated F1 hybrid pigeon embryos by crossing an Old Dutch Capuchine (scale-footed) to a fairy swallow (muffed), and measured parent-of-origin allele expression in the hybrid embryonic hindlimb at HH25 (Figure 4A) (Domyan et al., 2014). Consistent with expression differences we observed among breeds (Figure 3A,B), expression of the feathered-foot allele of Pitx1 was approximately 20% lower than the scaled-foot allele (expression of feathered-foot relative to scaled-foot allele: 0.807 ± 0.039, p=0.003, n = 6 embryos), and expression of the feathered-foot allele of Tbx5 was nearly 1600% higher than the scaled-foot allele (relative expression of feathered-foot allele: 15.75 ± 4.69, p=0.002, n = 7 embryos) (Figure 4B,C). In contrast, expression levels of feathered-foot and scaled-foot alleles of Tbx3 were indistinguishable (relative expression of feathered-foot allele: 0.99 ± 0.05, p=0.68, n = 7 embryos) (Figure 4D). These results directly show that cis-acting genetic changes alter expression of feathered-foot alleles of Pitx1 and Tbx5, but not Tbx3, in the embryonic pigeon hindlimb.

Allele-specific expression (ASE) assays demonstrate cis-regulatory changes in Pitx1 and Tbx5.

(A) Schematic of ASE assay using Tbx5 expression as an example. Differences in Tbx5 expression between scale-footed and feather-footed breeds could be due to trans- and/or cis-acting mutations. If expression differences between parental breeds are due to trans changes only (stars), then expression of the two Tbx5 alleles in hybrid embryos will be the same (top right). In contrast, if cis-regulatory changes underlie differences in Tbx5 expression between the parental breeds, then expression of the two Tbx5 alleles in hybrid embryos will be different (bottom right). (B-D) ASE assay in hybrid hindlimb buds indicate cis-regulatory divergence between scale-footed (Old Dutch Capuchine) and muffed (fairy swallow) pigeon breeds in Pitx1 (B) and Tbx5 (C), but not in Tbx3 (D). Dashed blue line indicates null hypothesis of equal expression of alleles. (E) ASE assay in hybrid hindlimb buds indicate cis-regulatory divergence in Tbx5 between feather-footed (silkie) and scale-footed (white leghorn) chicken breeds. Boxes in (B-E) span 1st to 3rd quartiles, bars extend to minimum and maximum observed values if within 1.5 times the interquartile range of the box, circles indicate values outside of this range, black line indicates median. **p≤0.01 ***p≤0.005.

https://doi.org/10.7554/eLife.12115.016
Figure 4—source data 1

Source data from pyrosequencing ASE experiments.

Abbreviations: fl, embryonic forelimb; hl, embryonic hindlimb.

https://doi.org/10.7554/eLife.12115.017

Patterns of open chromatin differ between scale-footed and feather-footed pigeon breeds

We next performed genome-wide comparisons of differential enhancer activity in embryonic limbs of racing homers and Indian fantails, using H3K27ac as a marker for open chromatin (Creyghton et al., 2010). Strikingly, two regions of significantly different enrichment in racing homer hindlimbs relative to forelimbs were within the 44-kb genomic region that is deleted in the Indian fantail, and one of these regions (log10 likelihood ratio = 4.93) is adjacent to the hs1473 limb enhancer (Figure 2C). Furthermore, the most significant differentially enriched region in Indian fantail hindlimbs relative to racing homer hindlimbs was directly upstream of Tbx5 (log10 likelihood ratio = 10.3) (Figure 2E). The overlapping patterns of enrichment in Indian fantail hindlimbs and wild-type forelimbs suggest that ectopic hindlimb expression of Tbx5 could be due to de-repression of forelimb-specific enhancers. In summary, the differential expression of Pitx1 and Tbx5 among pigeon breeds (Figure 3) and between alleles (Figure 4) is also reflected by differential chromatin activation at these genes.

Muffed pigeon breeds incur musculoskeletal patterning changes

In mouse and chicken embryos, experimental manipulation of Pitx1 and Tbx5 expression can result in muscular and skeletal abnormalities. Experiments in both chick and mouse consistently demonstrate that Pitx1 plays a necessary (but not sufficient) role in determining hindlimb-type morphology of the skeleton, muscles, and tendons (Logan and Tabin, 1999; Takeuchi et al., 1999; Minguillon et al., 2005; DeLaurier et al., 2006; Duboc and Logan, 2011). Experimentally induced ectopic expression of Tbx5 in the hindlimbs of chick embryos can also induce muscular and skeletal anomalies, although Tbx5 does not directly control limb skeletal patterning or determine forelimb-type morphology in mice (Rodriguez-Esteban et al., 1999; Takeuchi et al., 1999; Minguillon et al., 2005Hasson et al., 2007) However, normal patterning of limb muscles and tendons is dependent on Tbx5 and Tbx4 in mice (Hasson et al., 2010). These apparent discrepancies between mammalian and avian systems point to subtle differences in limb development in different lineages (Horton et al., 2008).

Given the dramatic musculoskeletal defects observed in other organisms with experimentally altered Pitx1 and Tbx5 expression, we compared the hindlimb morphology of adult feral pigeons (scale-footed, n=2) to that of the English trumpeter (muffed, n=2) and the Pomeranian pouter (muffed, n=1). We found consistent soft-tissue patterning defects in both feather-footed pigeon breeds: the fibularis longus (FL) tendon inserts on the dorsal rather than ventral surface of the proximal tarsometatarsus, the flexor perforans et perforatus (FPP3) muscle adopts a longitudinal rather than pennate fiber orientation, and a slip of the FPP3 fuses with the FL tendon (Figure 5A). These changes are aberrations of normal patterning, although they are not necessarily clear transformations to a more forelimb-like configuration. We also found that the fibula, which is normally splint-like and shorter than the tibia in pigeons, was enlarged (Figure 5D,E) and two phalanges of digit 4 were fused in feather-footed breeds (not shown). These are not necessarily limb-type transformations, either. However, experimental ectopic expression of Tbx5 in the hindlimbs of chick embryos produces an enlargement of the fibula reminiscent of extreme pigeon phenotypes, and Takeuchi et al., 1999 compared this morphology to a forelimb-like condition (the fibula “makes a joint at its distal end like a normal ulna [the corresponding postaxial zeugopod bone of the forelimb],” p. 810). Notably, all of the modified structures of ptilopodous pigeons develop in the posterior (lateral in the adult) and dorsal hindlimb, which are the primary sites of ectopic Tbx5 expression. Thus, the morphological changes to the hindlimbs of feather-footed pigeon breeds are considerably more than skin deep.

Muffed pigeons have re-patterned hindlimb musculoskeletal system and wing-like foot feathers.

(A,B,C) Gross muscle morphology of scale-footed (feral) and muffed (Pomeranian pouter) left hindlimbs, dorsal view. (A) Skin and superficial muscles have been removed to reveal re-patterning of the fibularis longus (FL, red). Dashed black line, approximate position of ankle joint axis. (B,C) The FPP3 is a pinnate muscle in scale-footed pigeons (B), but a slip of fibers fuses with the adjacent FL in muffed pigeons (arrowhead in C). (D,E) X-ray computed tomography images of scale-footed (feral, right leg) and muffed (English trumpeter, left leg) hindlimbs. Arrowheads mark the proximal and distal ends of the fibula. The wild-type pigeon fibula (D) is short and splint-like. In the muffed bird (E), the fibula extends from the knee to the ankle. We observed distal elongation of the fibula in another muffed breed (fairy swallow) but the fibula did not completely extend to the ankle (not shown). t, tibia; tmt, tarsometatarsus. (F) Toe and wing (flight) feathers of a muffed pigeon (English trumpeter), highlighting vane width asymmetries. Blue bar, inner vane; red bar, outer vane. Scale bar = 2 cm.

https://doi.org/10.7554/eLife.12115.018

Tbx5 is ectopically expressed in the hindlimb buds of feather-footed chickens

Other bird species, including domestic chickens, independently evolved foot feathers. Similar to pigeons, Tbx5 was ectopically expressed at HH25 in hindlimb buds of two feather-footed chicken breeds, the Cochin and the silkie (hindlimb expression relative to white leghorn forelimb: white leghorn 0.005 ± 0.004; Cochin 0.066 ± 0.050, p=0.002; silkie 0.034 ± 0.014, p=0.009; n = 6 white leghorns, 6 Cochins, 4 silkies) (Figure 3G). Ectopic Tbx5 expression in feathered-foot chicken embryos had a similar domain to that of feathered-foot pigeon embryos at HH25 (Figure 3I,J) and, as in pigeons, cis-acting changes contributed to this expression (expression of feathered-foot allele relative to scaled-foot allele in HH25 silkie x white leghorn F1 hybrid hindlimbs: 1.80 ± 0.41, p=2.55 x 10–5, n = 11 hybrid embryos) (Figure 4E). Hence, Tbx5-related developmental mechanisms may, in part, underlie the evolution of foot feathering in two species that last shared a common ancestor over 80 million years ago (Claramunt and Cracraft, 2015).

Classical genetic studies implicate at least two loci in heavy foot feathering in chickens (Punnett and Bailey, 1918; Lambert and Knox, 1929; Warren, 1948; Somes, 1992), although the molecular genetic origins of the trait remain unknown. Previously, a chromosome region containing Pitx1 was implicated in foot feathering in silkie chickens (Dorshorst et al., 2010). However, we did not detect statistically significant changes in Pitx1 expression between scaled-foot (white leghorn) and feathered-foot (silkie and Cochin) chicken embryos at HH25 (expression relative to white leghorn: Cochin 0.92 ± 0.24, p=0.93; silkie 0.71 ± 0.18, p=0.18; n = 6 each) (Figure 3F). This apparent conflict could be because the causative gene in silkies is not actually Pitx1 but rather a gene closely linked to it, or because Pitx1 expression differences are more pronounced and consistent at developmental stages that we did not assay. Furthermore, different populations of breeds such as silkies appear to have different constellations of ptilopody loci and alleles, and it is possible that we used strains that do not have Pitx1 mutations (Wexelsen, 1934; Somes, 1992). Also in contrast to our results in feather-footed pigeons, Tbx3 was not upregulated in ptilopodous chicken breeds (white leghorn 1 ± 0.17, silkie 0.48 ± 0.18, p=0.004; Cochin 0.86 ± 0.44, p=0.40; silkie vs. Cochin p=0.07; n = 6 samples each) (Figure 3—figure supplement 1D). In all, these results suggest that both shared and distinct mechanisms regulate foot feathering among avian species.

Discussion

Genetic architecture of ptilopody in pigeons

Extensive classical breeding experiments in pigeons demonstrate that complex derived traits can often be parsed into component parts (Sell, 1994, 2012). Thus, while traits are not always simple, they are often genetically tractable when using an informed breeding strategy (Domyan et al., 2014). Equivalent insights about the genetic architecture of phenotypic divergence between wild vertebrate species are often considerably more difficult to acquire. With pigeons, however, we have documentation for specific breed selection criteria and direct evidence for the resulting genetic architecture of derived traits (Levi, 1965; 1986; Sell, 1994; National Pigeon Association, 2010; Sell, 2012). This information offers a crucial advantage because it informs how we design genetic crosses and choose breeds for whole-genome resequencing to identify causal genes and mutations. Thus, we can combine classical breeding strategies and genomics to identify the molecular basis of both simple and oligogenic traits, as well as dissect different components of a complex phenotype, and define functional interactions among genes (Shapiro and Domyan, 2013; Domyan et al., 2014; Vickrey et al., 2015).

Classical studies in pigeon suggest two major-effect loci – grouse (gr) and Slipper (Sl) – are responsible for most of the variation in foot feathering (Doncaster, 1912; Wexelsen, 1934; Hollander, 1937; Levi, 1986). Through a combination of genetic, genomic, and developmental approaches, our data implicate regulatory mutations in the limb outgrowth and identity genes Pitx1 and Tbx5 as the molecular identities of the gr and Sl locus, respectively (Figure 6). Each locus has significant and separable effects on qualitative and quantitative variation in epidermal appendages: derived alleles of Pitx1 increase the extent of foot feathering, while a derived allele of Tbx5 is associated with the more elaborate muffed phenotype (Figure 1). Notably, these feathers are most numerous on the central forward-facing toe (digit 3; Table 2; Figure 3—figure supplement 3), just as forelimb feathers are most numerous on the central forelimb digit in birds and their dinosaurian relatives (Yalden et al., 1985; Gishlick et al., 2001; Hieronymus, 2015). Further, we also find that muscular and skeletal morphology are altered in muffed pigeons.

Model describing link between Pitx1 and Tbx5 expression levels and foot epidermal appendage morphology.

Darker colors indicate higher expression levels. Decreased expression of Pitx1 and ectopic expression of Tbx5 are associated with foot feathering (and other morphological transformations) in domestic pigeons.

https://doi.org/10.7554/eLife.12115.019

Collectively, these findings point to a partial alteration of the identity of the developing hindlimb, rather than localized changes to individual epidermal placodes. These alterations do not represent a complete transformation of limb type, as the hindlimbs of feather-footed pigeons are still readily recognizable as legs. This suggests that limb-type identity is not a simple binary choice between two global fates. For example, feather-footed pigeons have a radical transformation of the distal hindlimb dermis, yet changes to other hindlimb mesoderm derivatives (muscle, skeleton) are subtler and largely restricted to lateral structures in the adult. Therefore, we propose that different aspects of fore- and hindlimb morphology could have different dosage- and/or stage-dependent requirements for exposure to identity cues. Our ongoing analyses of musculoskeletal phenotypes in our F2 cross, which includes individuals with different combinations of feathered-foot alleles of Pitx1 and Tbx5, will allow us to understand the separate and epistatic effects of these loci on musculoskeletal anatomy. We note that, although genetic manipulations indicate that Tbx5 does not specify forelimb identity in mice, the divergence time between mammals and birds is deep (>300 million years) and subtly different roles for this transcription factor in limb outgrowth and identity might have evolved in these lineages (Minguillon et al., 2005; Horton et al., 2008).

Roles of Pitx1 and Tbx5 in diversity and disease

Another important avenue of future research will be to determine the downstream molecular consequences of Pitx1 and Tbx5 misregulation, and how this ultimately results in the transformation of scaled into feathered epidermis. Mutations at these genes can cause congenital limb deformities in humans, including clubfoot and Liebenberg syndrome (Pitx1; Gurnett et al., 2008; Spielmann et al., 2012), and Holt-Oram syndrome (Tbx5; Basson et al., 1997). Notably, haploinsufficiency causes these human syndromes, and clubfoot is partially penetrant in Pitx1+/- mice (Alvarado et al., 2011), collectively pointing to an exquisite sensitivity of limb morphology to levels of Pitx1 and Tbx5 gene products. Pitx1 is also involved repeatedly in the evolution of adaptive pelvic fin loss in stickleback fish (Cresko et al., 2004; Shapiro et al., 2004; 2006; Coyle et al., 2007; Chan et al., 2010; Shikano et al., 2013). Threespine sticklebacks that are homozygous for a Pitx1 pelvic enhancer deletion have severely reduced or absent pelvises, but heterozygotes also have smaller pelvises, thereby allowing natural selection to act on fish carrying one mutant allele (Shapiro et al., 2004; Chan et al., 2010). Similarly, we observed increased foot feathering in pigeons with just one derived allele of Pitx1 or Tbx5 (Figure 1). Ancient pigeon breeders could have rapidly selected for ptilopodous pigeon lines starting with birds that were heterozygous for mutations at either locus, and later generated the extreme muffed phenotype by hybridization. Together, studies of diversity and disease indicate that modest changes to the amount and location of Pitx1 and Tbx5 gene expression can cause dramatic alterations to limb development and morphology.

In addition to implicating cis- acting mutations in Pitx1 and Tbx5 driving transformation of scales into feathers in domestic pigeon, our results suggest that additional, as yet unidentified, mutations contribute to the muff phenotype. Although all feathered-foot embryos examined in our gene expression experiments contained derived Pitx1 and Tbx5 haplotypes, misregulation of each gene was more severe in the large-muffed English trumpeter than in the small-muffed Indian fantail. The English trumpeter may therefore contain additional cis-acting mutations at one or both loci, and/or mutations in upstream regulators of Pitx1 and Tbx5. Additional studies will be required to discriminate between these possibilities. Our findings suggest a quantitative link between transcription factor abundance and skin appendage fate and morphology, thereby highlighting foot-feathering in pigeons as a model for studying the regulatory interactions that govern expression of these two important determinants of limb morphology.

Evolution of epidermal appendage distribution

How might pigeons help us understand the evolution of epidermal appendage distribution and limb morphology in other species? Our findings suggest the mechanistic basis for the development of feathered feet in two distantly related domestic bird species is due to a partial transformation of limb identity, through cis-acting regulatory mutations in limb-type specific transcription factors. Most modern wild birds have a scaled metatarsus and toes, although some species (e.g., ptarmigan, snowy owl, and golden eagle) have extensive foot feathering. However, recent paleontological evidence suggests that feathers – not scales – might be the ancestral hindlimb skin appendages in birds and some of their close non-avian dinosaur relatives (Hu et al., 2009; Zheng et al., 2013). Thus, the epidermis of feather-footed modern birds might actually represent a reversion to the ancestral avialan skin condition. In some cases, the large, asymmetric-vaned, pennaceous metatarsal feathers of basal birds and their non-avian dinosaur relatives are so extensive that their hindlimbs have been interpreted as 'hind wings', although they clearly retain hindlimb skeletal identity (Zheng et al., 2013). Furthermore, these long foot feathers are directed laterally from the foot, and they display vane width asymmetries reminiscent of flight feathers; we find a similar morphology in muffed pigeons (Figure 5F, Figure 3—figure supplement 3). Perhaps not coincidentally, Darwin, 1868 noted of the muffed English trumpeter pigeon, “Their feet are so heavily feathered, that they almost appear like little wings” (p. 155).

Building on classical breeding experiments in both pigeons and chickens, we find that a relatively small number of genetic changes account for a large proportion of the variation in epidermal appendage morphology and distribution. Thus, major determinants of dramatic phenotypic variation can be mechanistically simple and therefore potentially evolve rapidly. In pigeons, these mechanisms can generate wing-like feathers on a hindlimb that is not used for powered flight or gliding. This, in turn, suggests that wing-like foot and leg feathers in other species, such as non-avian dinosaurs, might result from developmental constraints on the morphology of large limb feathers, rather than from functional adaptations for flight (Gould and Lewontin, 1979; Foth et al., 2014).

Materials and methods

Animal husbandry and phenotyping of F2 offspring

Request a detailed protocol

Animals were housed in accordance with the University of Utah Institutional Animal Care and Use Committees of University of Utah (protocols 10–05007 and 13–04012). 130 F2 offspring were generated by mating a male Pomeranian pouter to two female Scandaroons, and DNA samples extracted (DNeasy Blood and Tissue Kit, Qiagen, Valencia, CA). 114 F2 offspring survived to 6 months of age, at which time they were euthanized and phenotypic measurements taken. Proportion of the tarsus was measured by dividing the length of the dorsal tarsus that was feathered by the total length of the tarsus (measured from the tibia-tarsometatarsus joint to the distal aspect of tarsometatarsal-phalangeal joint of digit 3), and averaged between the two tarsi. Toe feathers were counted on each toe, and summed across all 8 toes. The length of each of the longest three toe feathers on digit 3 (the central forward-directed toe), which bore the longest toe feathers on each foot, was measured to the nearest 1 mm and averaged for each bird.

Whole-genome genotyping by sequencing (GBS)

Request a detailed protocol

For genotyping, we used a previously published approach (Elshire et al., 2011) with minor modifications. Briefly, for each founder parent and 130 F2 offspring, 50 ng of DNA was digested with ApeKI, ligated to barcoded adapters, and then 10 ng of each barcoded sample was pooled in batches of 26 individuals and purified (Qiagen PCR Purification Kit). DNA fragments 550–650 bp in size were selected using Pippin Prep (Sage Science, Beverly, MA), and amplified by 10–12 cycles of PCR using custom indexed primers. Libraries were purified with Ampure beads (Sigma-Aldrich, St. Louis, MO) and sequenced using 100- or 125-bp, paired-end sequencing on the Illumina HiSeq2000 platform at the University of Utah Genomics Core Facility. Reads were trimmed to 90 bp, filtered for quality, and de-multiplexed using Stacks (Catchen et al., 2011). Reads were mapped to the pigeon reference genome (Shapiro and Domyan, 2013) using Bowtie2 (Langmead and Salzberg, 2012), filtering for MAPQ < 20. The average number of mapped reads among F2 individuals was 3,397,598, with a mean depth of 6.3x. Genotypes were called using Stacks (Catchen et al., 2011), with a minimum read-depth cutoff of 5. Markers that were genotyped in ≥ 70 of the 130 F2 individuals were retained.

Genetic map construction and QTL mapping

Request a detailed protocol

Genetic map construction and QTL mapping was performed using R/qtl (www.rqtl.org) (Broman et al., 2003). Markers showing segregation distortion (Chi-square, p<0.05) were removed. 3803 markers were assembled into linkage groups using the parameters (max.rf = 0.15, min.lod = 6). Linkage groups were numbered in descending order, based on the number of markers. Linkage group 11 and 20 QTL were initially mapped using the scanone function using Haley-Knott regression. Probable false-homozygote genotyping errors resulting from the low read-depth cutoff used (5x), identified as closely-spaced double-crossover events, were manually corrected on these linkage groups. Subsequently, the stepwiseqtl function was used to identify additional QTL, and the fitqtl function used to account for the effect of one linkage group while calculating the LOD scores and percent variance explained (PVE) of the other. Significance thresholds of α = 0.05 were calculated with 1000 permutations of each phenotype across all linkage groups. The peak markers for each phenotype were used to test for the effect of each QTL.

Genomic analyses

Request a detailed protocol

BAM files generated previously for a whole-genome resequencing panel (Shapiro and Domyan, 2013) were combined with BAM files for two new Pomeranian pouter whole-genome sequences to call genomic variants (SNVs and small indels) using the Genome Analysis Toolkit (Unified Genotyper and LeftAlignAndTrimVariants functions, default settings) (McKenna et al., 2010). We removed variant sites that were called in two or fewer genomes (i.e., all other genomes were no-calls) or that had variant alleles on only two or fewer chromosomes, as these categories of low-frequency variants would be uninformative to our analyses. The resulting variant call format (VCF) file was used for subsequent analyses.

Individual birds from different breeds were binned into the following phenotypic classifications:

Groused: Berlin long-faced tumbler, Lahore, Oriental frill, Shaksharli.

Muffed: English long-faced muffed tumbler, English pouter, English trumpeter, frillback, ice pigeon, Indian fantail (2 individuals), Pomeranian pouter (2 individuals), Saxon monk, Saxon pouter. The English pouter is an unusual breed that is sometimes classified by breeders as slipper only. Its foot feathering is far more extensive than groused breeds, which led us to include it in the muffed group for the purposes of the genomic analyses.

Scale-footed: African owl, archangel, Birmingham roller, carneau, Chinese owl, cumulet, Egyptian swift, English carrier, fantail, feral (2 individuals), Iranian tumbler, Jacobin, king, Lebanon, Marchenero pouter, mookee, Oriental roller, parlor roller, runt, Scandaroon, Spanish barb, starling, Syrian dewlap, Thai laugher.

pFst, a modified likelihood ratio test that accounts for genotype uncertainty, extended haplotype homozygosity (EHH), and haplotype network analyses were implemented using the GPAT++ software library (Kronenberg et al., 2014; see https://github.com/vcflib/vcflib for software updates).

Genotyping assays

Request a detailed protocol

Primers for genotyping the scaffold 79 deletion are listed in Supplementary file 1. Breeds used for association testing were as follows:

Feather-footed (21 breeds total): Berlin long-faced tumbler, Berlin short-faced tumbler, Bokhara trumpeter, classic Oriental frill, crested Saxon field color, English trumpeter, fairy swallow, frillback, German double-crested trumpeter, ice pigeon, Indian fantail, Lahore, Mindian fantail, Oriental frill, Persian roller, Pomeranian pouter, Russian tumbler, saint, Schmalkaldner moorhead, Uzbeck tumbler, West of England.

Scale-footed (56 breeds total): African owl, Altenburg trumpeter, American flying tumbler, American giant homer, American mini crest, American show racer, archangel, Bohemian pouter, Brunner cropper, Budapest tumbler, Cauchois, Chinese owl, cumulet, Danzig highflier, domestic show flight, dragoon, English baldhead long-faced clean-legged tumbler, English carrier, English magpie, English short-faced tumbler, exhibition homer, fantail, Franconian trumpeter, French mondaine, giant runt, Holle cropper, horseman pouter, Italian owl, Jacobin, Jiennesse pouter, king, Lebanon, Spanish little friar tumbler, Maltese, medium-faced crested helmet, Modena, mookee, Norwich cropper, nun, Old Dutch Capuchine, Old German owl, Oriental roller, parlor roller, Portuguese tumbler, Scandaroon, showtype racing homer, Spanish barb, starling, Syrian Baghdad, Texas pioneer, Thai laugher, Thuringer clean leg, Vienna medium-faced tumbler, Voorburg shield cropper, zitterhals.

High-throughput chromatin immunoprecipitation and sequencing (ChIP-seq)

Request a detailed protocol

Forelimb and hindlimb buds from HH25 racing homer and Indian fantail embryos were collected and placed in 1% formaldehyde for 20 min at room temperature, then washed 3x in ice-cold PBS and stored at -80°C until chromatin extraction. ChIP was performed on 200 micrograms of chromatin isolated from embryonic pigeon limbs. Control libraries were prepared using 100 ng of input chromatin. A total of 16 libraries were created (8 ChIP and 8 input controls for each breed and limb combination). A validated monoclonal antibody against H3K27ac (Millipore #05–1334, Billerica, MA) was used to perform ChIP, and sequencing libraries were prepared using NEBNext Ultra DNA Library Prep Kit for Illumina with NEBNext Multiplex Oligos for Illumina (Index Primers Set 1; New England BioLabs, Ipswich, MA). All libraries were size selected using SPRI magnetic beads to eliminate adapter dimers. All 8 ChIP libraries showed enrichment for a positive control site relative to input libraries (tested by qPCR). Single-end, 50-bp read sequencing was performed on Illumina HiSeq2000 platform at the University of Utah Genomics Core Facility.

Fold-enrichment plots were generated using MACS (Zhang et al., 2008; Feng et al., 2012) and visualized in IGV (Robinson et al., 2011; Thorvaldsdóttir et al., 2013). Regions of differential enrichment between racing homer and Indian fantail hindlimbs were identified using function bdgdiff in MACS2 (https://pypi.python.org/pypi/MACS2/2.0.10.20130522). Regions with a log10 likelihood ratio ≥ 3 were considered to have differential enrichment between the two groups.

RNA isolation and cDNA synthesis

Request a detailed protocol

To assay gene expression, limb buds from HH25 embryos were harvested and placed in RNAlater (Qiagen, Valencia, CA) at 4°C overnight. Total RNA was extracted, cleaned and DNase-treated (Qiagen RNeasy Kit). mRNA was reverse-transcribed to cDNA using oligo-dT and M-MLV RT (Invitrogen, Carlsbad, CA) according to the manufacturer’s protocol.

qRT-PCR analyses

Request a detailed protocol

cDNA was amplified using intron-spanning primers for each target using a CFX96 qPCR instrument and iTaq Universal Sybr Green Supermix (Bio-Rad, Hercules, CA). Results were compared by Mann-Whitney U test. Two technical replicates of each sample were performed, and the mean value determined. Differences were considered statistically significant if p<(0.05 / # genes assayed) to control for multiple-testing. Each experiment was performed three times, and the results presented are from one representative experiment. Primers used for each assay are listed in Supplementary file 1.

Allele-specific expression assay

Request a detailed protocol

SNPs in Pitx1 and Tbx5 transcripts were identified by Sanger sequencing in the parents of a cross between an Old Dutch Capuchine (scale-footed) and a fairy swallow (muffed) that was homozygous for the 44-kb deletion upstream of Pitx1, and PyroMark Custom Assays (Qiagen) for each SNP were designed using the manufacturer’s software. Pyrosequencing was performed on cDNA and gDNA derived from HH25 limb buds using a PyroMark Q24 instrument (Qiagen). The signal intensity ratio of feathered allele to scaled allele from cDNA samples was normalized to ratios obtained from gDNA samples from the same embryos to control for allele-specific amplification bias. Normalized ratios were analyzed by Mann-Whitney U test, and considered significant if p<(0.05 / # genes assayed) to control for multiple-testing. Each experiment was performed twice, and the results presented are representative. Primers used for each assay are listed in Supplementary file 1.

Whole-mount in situ hybridization

Request a detailed protocol

Linear templates for probe synthesis were amplified from cDNA by PCR using primers listed in Supplementary file 1. Binding sites for T3 and T7 polymerase were incorporated into the forward and reverse primers to facilitate subsequent transcription of sense and antisense probe, respectively.

Embryos used for RNA in situ hybridization were dissected from eggs, and fixed overnight in 4% paraformaldehyde at 4°C on a shaking table, then dehydrated into 100% MeOH and stored at -20°C. RNA in situ hybridization was performed as described (Abler et al., 2011). Hybridization with sense probe was performed as negative control.

Data availability

The following data sets were generated
The following previously published data sets were used

References

    1. Broman KW
    2. Wu H
    3. Sen S
    4. Churchill GA
    (2003)
    R/qtl: QTL mapping in experimental crosses
    Bioinformatics 19:889–890.
    1. Crowe R
    2. Henrique D
    3. Ish-Horowicz D
    4. Niswander L
    (1998)
    A new role for Notch and Delta in cell fate decisions: patterning the feather array
    Development 125:767–775.
    1. Danforth CH
    (1919)
    An Hereditary Complex in the Domestic Fowl
    Genetics 4:587–596.
  1. Book
    1. Darwin CR
    (1868)
    The Variation of Animals and Plants Under Domestication
    London: John Murray.
    1. Gibson-Brown JJ
    2. Agulnik SI
    3. Silver LM
    4. Niswander L
    5. Papaioannou VE
    (1998a)
    Involvement of t-box genes Tbx2-Tbx5 in vertebrate limb specification and development
    Development 125:2499–2509.
    1. Gishlick AD
    (2001)
    New Perspectives on the Origin and Early Evolution of Birds: Proceedings of the International Symposium in Honor of John H. Ostrom
    The function of the manus and forelimb of Deinonychus antirrhopus and its importance for the origin of avian flight , New Perspectives on the Origin and Early Evolution of Birds: Proceedings of the International Symposium in Honor of John H. Ostrom, New Haven, Peabody Museum of Natural History, Yale University.
    1. Gould SJ
    2. Lewontin RC
    (1979)
    The spandrels of San Marco and the Panglossian paradigm: a critique of the adaptationist programme
    Proceedings of the Royal Society of London. Series B, Biological Sciences of the United States of America 205:581–598.
    1. Hamburger V
    2. Hamilton HL
    (1951)
    A series of normal stages in the development of the chick embryo
    Journal of Morphology 88:49–92.
    1. Hollander WF
    (1937)
    Department of Genetics
    p 67, Hereditary Interrelationships of Certain Factors in Pigeons, Department of Genetics, University of Wisconsin.
    1. Lambert WV
    2. Knox CW
    (1929)
    The inheritance of of leg-feathering in the chicken
    Iowa Acad Sci Proc 36:383–384.
  2. Book
    1. Levi WM
    (1965)
    Encyclopedia of Pigeon Breeds
    Sumter: Levi Publishing Co., Inc.
  3. Book
    1. Levi WM
    (1986)
    The Pigeon
    Sumter: Levi Publishing Co., Inc.
    1. Logan M
    2. Simon HG
    3. Tabin C
    (1998)
    Differential regulation of T-box and homeobox transcription factors suggests roles in controlling chick limb-type identity
    Development 125:2825–2835.
    1. Logan M
    2. Tabin CJ
    (1999)
    Role of Pitx1 upstream of Tbx4 in specification of hindlimb identity
    Science 283:1736–1739.
  4. Book
    1. Lucas AM
    2. Stettenheim PR
    (1972)
    Avian Anatomy - Integument
    Washington, DC.
  5. Book
    1. National Pigeon Association
    (2010)
    National Pigeon Association Book of Standards
    Goodlettsville: Purebred Pigeon Publishing.
    1. Ruvinsky I
    2. Oates AC
    3. Silver LM
    4. Ho RK
    (2000)
    The evolution of paired appendages in vertebrates: T-box genes in the zebrafish
    Development Genes and Evolution 210:82–91.
  6. Book
    1. Sell A
    (1994)
    Breeding and Inheritance in Pigeons
    Germany: Schober Verlags-GmbH, Hengersberg.
  7. Book
    1. Sell A
    (2012)
    Pigeon Genetics: Applied Genetics in the Domestic Pigeon
    Germany: Verlag Karin und A. Sell, Achim.
    1. Shapiro MD
    2. Bell MA
    3. Kingsley DM
    (2006) Parallel genetic origins of pelvic reduction in vertebrates
    Proceedings of the National Academy of Sciences of the United States of America 103:13753–13758.
    https://doi.org/10.1073/pnas.0604706103
    1. Somes RG
    (1992)
    Identifying the ptilopody (feathered shank) loci of the chicken
    The Journal of Heredity 83:230–234.
    1. Szeto DP
    2. Rodriguez-Esteban C
    3. Ryan AK
    4. O'Connell SM
    5. Liu F
    6. Kioussi C
    7. Gleiberman AS
    8. Izpisúa-Belmonte JC
    9. Rosenfeld MG
    (1999)
    Role of the Bicoid-related homeodomain factor Pitx1 in specifying hindlimb morphogenesis and pituitary development
    Genes & Development 13:484–494.
    1. Tamura K
    2. Yonei-Tamura S
    3. Izpisúa Belmonte JC
    (1999)
    Differential expression of Tbx4 and Tbx5 in Zebrafish fin buds
    Mechanisms of Development 87:181–184.
    1. Warren DC
    (1948)
    Linkage relations of autosomal factors in the fowl
    Genetics 34:333–350.
    1. Wexelsen H
    (1934)
    Types of leg feathering in pigeons
    Hereditas 18:192–198.
  8. Book
    1. Wickham H
    (2009)
    ggplot2: Elegant Graphics for Data Analysis
    New York: Springer.
  9. Book
    1. Yalden DW
    (1985)
    Forelimb function in Archaeopteryx
    In: Hecht M. K, Ostrom J. H, Viohl G, Wellnhofer P, editors. The Beginnings of Birds: Proceedings of the International Archaeopteryx Conference, 1984. Eichstätt, Germany: Freunde des Jura-Museums Eichstätt. pp. 91–97.

Article and author information

Author details

  1. Eric T Domyan

    Department of Biology, University of Utah, Salt Lake City, United States
    Present address
    Department of Biology, Utah Valley University, Orem, United States
    Contribution
    ETD, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  2. Zev Kronenberg

    Department of Human Genetics, University of Utah, Salt Lake City, United States
    Present address
    Department of Genome Sciences, University of Washington, Seattle, United States
    Contribution
    ZK, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  3. Carlos R Infante

    Department of Genetics, University of Georgia, Athens, United States
    Present address
    Department of Molecular and Cellular Biology, University of Arizona, Tucson, United States
    Contribution
    CRI, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  4. Anna I Vickrey

    Department of Biology, University of Utah, Salt Lake City, United States
    Contribution
    AIV, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  5. Sydney A Stringham

    Department of Biology, University of Utah, Salt Lake City, United States
    Contribution
    SAS, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  6. Rebecca Bruders

    Department of Biology, University of Utah, Salt Lake City, United States
    Contribution
    RB, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  7. Michael W Guernsey

    Department of Biology, University of Utah, Salt Lake City, United States
    Present address
    Department of Genetics, Stanford University, Stanford, United States
    Contribution
    MWG, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  8. Sungdae Park

    Department of Genetics, University of Georgia, Athens, United States
    Contribution
    SP, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  9. Jason Payne

    Poultry Science Department, University of Georgia, Athens, United States
    Contribution
    JP, Acquisition of data, Drafting or revising the article, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  10. Robert B Beckstead

    Poultry Science Department, University of Georgia, Athens, United States
    Contribution
    RBB, Conception and design, Drafting or revising the article, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  11. Gabrielle Kardon

    Department of Human Genetics, University of Utah, Salt Lake City, United States
    Contribution
    GK, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  12. Douglas B Menke

    Department of Genetics, University of Georgia, Athens, United States
    Contribution
    DBM, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7109-1451
  13. Mark Yandell

    1. Department of Human Genetics, University of Utah, Salt Lake City, United States
    2. Utah Center for Genetic Discovery, University of Utah, Salt Lake City, United States
    Contribution
    MY, Conception and design, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  14. Michael D Shapiro

    Department of Biology, University of Utah, Salt Lake City, United States
    Contribution
    MDS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    shapiro@biology.utah.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2900-4331

Funding

National Institutes of Health (F32GM103077 Fellowship)

  • Eric T Domyan

National Institutes of Health (T32HD07491 Fellowship)

  • Eric T Domyan

National Institutes of Health (T32GM007464 Fellowships)

  • Zev Kronenberg
  • Sydney A Stringham

National Science Foundation (EDEN RCN Internship IOS0955517)

  • Anna I Vickrey

National Science Foundation (Graduate Research Fellowship)

  • Rebecca Bruders

National Institutes of Health (R01HD053728)

  • Gabrielle Kardon

National Science Foundation (CAREER IOS1149453)

  • Douglas B Menke

National Institutes of Health (R01GM104390)

  • Mark Yandell

National Institutes of Health (R01GM115996)

  • Michael D Shapiro

Burroughs Wellcome Fund (Career Award in the Biomedical Sciences)

  • Michael D Shapiro

National Science Foundation (CAREER DEB1149160)

  • Michael D Shapiro

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Della Fixsen, Cassandra Garner, Kamala Ganesh, Patrick Miller, the University of Utah Genomics Core, and the University of Utah Small Animal Imaging Core for technical assistance, and members of the Utah Pigeon Club and National Pigeon Association for generously providing samples. We also thank Elena Boer for comments and discussion on the manuscript. Wild-type pigeon image in Figure 5D is courtesy of Dr. M Scott Echols, Grey Parrot Anatomy Project, University of Utah and the Medical Center for Birds (Oakley, California). This work was supported by the National Science Foundation (CAREER DEB1149160 to MDS, CAREER IOS1149453 to DBM, IOS0955517 EDEN internship to AIV), a Burroughs Wellcome Career Award in the Biomedical Sciences (MDS), and the National Institutes of Health (R01GM115996 to MDS, F32GM103077 and T32HD07491 fellowships to ETD, T32GM007464 fellowships to ZK and SAS, R01HD053728 to GK, R01GM104390 to MY). We acknowledge a computer time allocation from the Center for High Performance Computing at the University of Utah. Illumina shotgun reads for two Pomeranian pouter pigeons are deposited in NCBI BioProject ID PRJNA284526. H3K27ac ChIP-seq data are deposited in the Gene Expression Omnibus (GEO accession number GSE67875). Correspondence should be addressed to MDS (shapiro@biology.utah.edu).

Ethics

Animal experimentation: This study was performed in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled and housed according to approved University of Utah institutional animal care and use committee (IACUC) protocols 10-05007 and 13-04012.

Copyright

© 2016, Domyan et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 5,975
    views
  • 925
    downloads
  • 68
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Eric T Domyan
  2. Zev Kronenberg
  3. Carlos R Infante
  4. Anna I Vickrey
  5. Sydney A Stringham
  6. Rebecca Bruders
  7. Michael W Guernsey
  8. Sungdae Park
  9. Jason Payne
  10. Robert B Beckstead
  11. Gabrielle Kardon
  12. Douglas B Menke
  13. Mark Yandell
  14. Michael D Shapiro
(2016)
Molecular shifts in limb identity underlie development of feathered feet in two domestic avian species
eLife 5:e12115.
https://doi.org/10.7554/eLife.12115

Share this article

https://doi.org/10.7554/eLife.12115

Further reading

    1. Developmental Biology
    2. Stem Cells and Regenerative Medicine
    Thi Thom Mac, Teddy Fauquier ... Thierry Brue
    Research Article

    Deficient Anterior pituitary with common Variable Immune Deficiency (DAVID) syndrome results from NFKB2 heterozygous mutations, causing adrenocorticotropic hormone deficiency (ACTHD) and primary hypogammaglobulinemia. While NFKB signaling plays a crucial role in the immune system, its connection to endocrine symptoms is unclear. We established a human disease model to investigate the role of NFKB2 in pituitary development by creating pituitary organoids from CRISPR/Cas9-edited human induced pluripotent stem cells (hiPSCs). Introducing homozygous TBX19K146R/K146R missense pathogenic variant in hiPSC, an allele found in congenital isolated ACTHD, led to a strong reduction of corticotrophs number in pituitary organoids. Then, we characterized the development of organoids harboring NFKB2D865G/D865G mutations found in DAVID patients. NFKB2D865G/D865G mutation acted at different levels of development with mutant organoids displaying changes in the expression of genes involved on pituitary progenitor generation (HESX1, PITX1, LHX3), hypothalamic secreted factors (BMP4, FGF8, FGF10), epithelial-to-mesenchymal transition, lineage precursors development (TBX19, POU1F1) and corticotrophs terminal differentiation (PCSK1, POMC), and showed drastic reduction in the number of corticotrophs. Our results provide strong evidence for the direct role of NFKB2 mutations in the endocrine phenotype observed in patients leading to a new classification of a NFKB2 variant of previously unknown clinical significance as pathogenic in pituitary development.

    1. Developmental Biology
    2. Genetics and Genomics
    Debashish U Menon, Prabuddha Chakraborty ... Terry Magnuson
    Research Article

    We present evidence implicating the BAF (BRG1/BRM Associated Factor) chromatin remodeler in meiotic sex chromosome inactivation (MSCI). By immunofluorescence (IF), the putative BAF DNA binding subunit, ARID1A (AT-rich Interaction Domain 1 a), appeared enriched on the male sex chromosomes during diplonema of meiosis I. Germ cells showing a Cre-induced loss of ARID1A arrested in pachynema and failed to repress sex-linked genes, indicating a defective MSCI. Mutant sex chromosomes displayed an abnormal presence of elongating RNA polymerase II coupled with an overall increase in chromatin accessibility detectable by ATAC-seq. We identified a role for ARID1A in promoting the preferential enrichment of the histone variant, H3.3, on the sex chromosomes, a known hallmark of MSCI. Without ARID1A, the sex chromosomes appeared depleted of H3.3 at levels resembling autosomes. Higher resolution analyses by CUT&RUN revealed shifts in sex-linked H3.3 associations from discrete intergenic sites and broader gene-body domains to promoters in response to the loss of ARID1A. Several sex-linked sites displayed ectopic H3.3 occupancy that did not co-localize with DMC1 (DNA meiotic recombinase 1). This observation suggests a requirement for ARID1A in DMC1 localization to the asynapsed sex chromatids. We conclude that ARID1A-directed H3.3 localization influences meiotic sex chromosome gene regulation and DNA repair.