Integration of chicken genomic resources to enable whole-genome sequencing

Modification of Keratin Integrations and the Associated Morphogenesis in Frizzling Chicken Feathers

The morphological and compositional complexities of keratinized components make feathers ingenious skin appendages adapted to diverse ecological needs. Frizzling feathers, characterized by their distinct curling phenotypes, offer a unique model to explore the intricate morphogenesis in developing a keratin-based bioarchitecture over a wide range of morphospace. Here, we investigated the heterogeneous allocation of α- and β-keratins in flight feather shafts of homozygous and heterozygous frizzle chickens by analyzing the medulla-cortex integrations using quantitative morphology characterizations across scales. Our results reveal the intriguing construction of the frizzling feather shaft through the modified medulla development, leading to a perturbed balance of the internal biomechanics and, therefore, introducing the inherent natural frizzling compared to those from wild-type chickens. We elucidate how the localized developmental suppression of the α-keratin in the medulla interferes with the growth of the hierarchical keratin organization by changing the internal stress in the frizzling feather shaft. This research not only offers insights into the morphogenetic origin of the inherent bending of frizzling feathers but also facilitates our in-depth understanding of the developmental strategies toward the diverse integuments adapted for ecological needs.

The chicken frizzle feather is due to an α-keratin (KRT75) mutation that causes a defective rachis

Feathers have complex forms and are an excellent model to study the development and evolution of morphologies. Existing chicken feather mutants are especially useful for identifying genetic determinants of feather formation. This study focused on the gene F, underlying the frizzle feather trait that has a characteristic curled feather rachis and barbs in domestic chickens. Our developmental biology studies identified defects in feather medulla formation, and physical studies revealed that the frizzle feather curls in a stepwise manner. The frizzle gene is transmitted in an autosomal incomplete dominant mode. A whole-genome linkage scan of five pedigrees with 2678 SNPs revealed association of the frizzle locus with a keratin gene-enriched region within the linkage group E22C19W28_E50C23. Sequence analyses of the keratin gene cluster identified a 69 bp in-frame deletion in a conserved region of KRT75, an α-keratin gene. Retroviral-mediated expression of the mutated F cDNA in the wild-type rectrix qualitatively changed the bending of the rachis with some features of frizzle feathers including irregular kinks, severe bending near their distal ends, and substantially higher variations among samples in comparison to normal feathers. These results confirmed KRT75 as the F gene. This study demonstrates the potential of our approach for identifying genetic determinants of feather forms.

With the availability of a sequenced chicken genome, the reservoir of variant plumage genes found in domestic chickens can provide insight into the molecular mechanisms underlying the diversity of feather forms. In this paper, we identify the molecular basis of the distinctive frizzle (F) feather phenotype that is caused by a single autosomal incomplete dominant gene in which heterozygous individuals show less severe phenotypes than homozygous individuals. Feathers in frizzle chickens curve backward. We used computer-assisted analysis to establish that the rachis of the frizzle feather was irregularly kinked and more severely bent than normal. Moreover, microscopic evaluation of regenerating feathers found reduced proliferating cells that give rise to the frizzle rachis. Analysis of a pedigree of frizzle chickens showed that the phenotype is linked to two single-nucleotide polymorphisms in a cluster of keratin genes within the linkage group E22C19W28_E50C23. Sequencing of the gene cluster identified a 69-base pair in-frame deletion of the protein coding sequence of the α-keratin-75 gene. Forced expression of the mutated gene in normal chickens produced a twisted rachis. Although chicken feathers are primarily composed of beta-keratins, our findings indicate that alpha-keratins have an important role in establishing the structure of feathers.

Precision of distances and ordering of microsatellite markers in consensus linkage maps of chromosomes 1, 3 and 4 from two reciprocal chicken populations using bootstrap sampling

Some factors complicate comparisons between linkage maps from different studies. This problem can be resolved if measures of precision, such as confidence intervals and frequency distributions, are associated with markers. We examined the precision of distances and ordering of microsatellite markers in the consensus linkage maps of chromosomes 1, 3 and 4 from two F(2) reciprocal Brazilian chicken populations, using bootstrap sampling. Single and consensus maps were constructed. The consensus map was compared with the International Consensus Linkage Map and with the whole genome sequence. Some loci showed segregation distortion and missing data, but this did not affect the analyses negatively. Several inversions and position shifts were detected, based on 95% confidence intervals and frequency distributions of loci. Some discrepancies in distances between loci and in ordering were due to chance, whereas others could be attributed to other effects, including reciprocal crosses, sampling error of the founder animals from the two populations, F(2) population structure, number of and distance between microsatellite markers, number of informative meioses, loci segregation patterns, and sex. In the Brazilian consensus GGA1, locus LEI1038 was in a position closer to the true genome sequence than in the International Consensus Map, whereas for GGA3 and GGA4, no such differences were found. Extending these analyses to the remaining chromosomes should facilitate comparisons and the integration of several available genetic maps, allowing meta-analyses for map construction and quantitative trait loci (QTL) mapping. The precision of the estimates of QTL positions and their effects would be increased with such information.

Duck (Anas platyrhynchos) linkage mapping by AFLP fingerprinting

Amplified fragment length polymorphism (AFLP) with multicolored fluorescent molecular markers was used to analyze duck (Anas platyrhynchos) genomic DNA and to construct the first AFLP genetic linkage map. These markers were developed and genotyped in 766 F2 individuals from six families from a cross between two different selected duck lines, brown Tsaiya and Pekin. Two hundred and ninety-six polymorphic bands (64% of all bands) were detected using 18 pairs of fluorescent TaqI/EcoRI primer combinations. Each primer set produced a range of 7 to 29 fragments in the reactions, and generated on average 16.4 polymorphic bands. The AFLP linkage map included 260 co-dominant markers distributed in 32 linkage groups. Twenty-one co-dominant markers were not linked with any other marker. Each linkage group contained three to 63 molecular markers and their size ranged between 19.0 cM and 171.9 cM. This AFLP linkage map provides important information for establishing a duck chromosome map, for mapping quantitative trait loci (QTL mapping) and for breeding applications.

Histocompatible chicken inbred lines: homogeneities in the major histocompatibility complex antigens of the GSP, GSN/1, PNP/DO and BM-C inbred lines assessed by hemagglutination, mixed lymphocyte reaction and skin transplantation

Chicken inbred lines of the GSP, GSN/1, PNP/DO and BM-C have been established by selection of a specific allele at the B blood group locus (MHC B-G region) and other polymorphic loci through pedigree mating. To extend the potential of these inbred lines as experimental animals in Aves, we assessed the antigenic homogeneities of the MHC antigens by three immunological methods. Antigenic variations of red blood cells (RBCs) were surveyed in the inbred lines and a random-bred line (NG) derived from the Nagoya breed by using ten kinds of intact antisera produced in the inbred line of chickens against RBCs of a red junglefowl and hybrids. In the hemagglutination test, no individual variations were found within the inbred line at all, while all the ten antisera detected highly heterogeneous reactions in individuals of the NG. The reciprocal one-way mixed lymphocyte reactions gave constantly higher stimulation responses (P<0.01) between individual pairs from the inbred lines having different B alleles compared to pairs within the inbred line, while lower stimulation was observed between pairs of the GSP and GSN/1 inbred lines both having the B(21) allele. In reciprocal skin transplantation, the transplanted skingrafts within the inbred line and between individuals from the GSP and GSN/1 inbred lines survived more than 100 days, while all the skingrafts showed signs of rejection within 7 days among the inbred lines having different B alleles. The results obtained by the three practical methods coincidentally indicated that the individuals in the respective four inbred lines were histocompatible, and further, that the GSP and GSN/1 individuals were histocompatible.

An integrated and comparative genetic map of the turkey genome

An integrated genetic linkage map was developed for the turkey (Meleagris gallopavo) that combines the genetic markers from the three previous mapping efforts. The UMN integrated map includes 613 loci arranged into 41 linkage groups. An additional 105 markers are tentatively placed within linkage groups based on two-point LOD scores and 19 markers remain unlinked. A total of 210 previously unmapped markers has been added to the UMN turkey genetic map. Markers from each of the 20 linkage groups identified in the Roslin map and the 22 linkage groups of the Nte map are incorporated into the new integrated map. Overall map distance contained within the 41 linkage groups is 3,365 cM (sex-averaged) with the largest linkage group (94 loci) measuring 533.1 cM. Average marker interval for the map was 7.86 cM. Sequences of markers included in the new map were compared to the chicken genome sequence by 'BLASTN'. Significant similarity scores were obtained for 95.6% of the turkey sequences encompassing an estimated 91% of the chicken genome. A physical map of the chicken genome based on positions of the turkey sequences was built and 36 of the 41 turkey linkage groups were aligned with the physical map, five linkage groups remain unassigned. Given the close similarities between the turkey and chicken genomes, the chicken genome sequence could serve as a scaffold for a genome sequencing effort in the turkey.

Avian genomics in the 21st century

The chicken has long been an important model organism for developmental biology, as well as a major source of protein with billions of birds used in meat and egg production each year. Chicken genomics has been transformed in recent years, with the characterisation of large EST collections and most recently with the assembly of the chicken genome sequence. As the first livestock genome to be fully sequenced it leads the way for others to follow--with zebra finch later this year. The genome sequence and the availability of three million genetic polymorphisms are expected to aid the identification of genes that control traits of importance in poultry. As the first bird genome to be sequenced it is a model for the remaining 9,600 species thought to exist today. Many of the features of avian biology and organisation of the chicken genome make it an ideal model organism for phylogenetics and embryology, along with applications in agriculture and medicine. The availability of new tools such as whole-genome gene expression arrays and SNP panels, coupled with information resources on the genes and proteins are likely to enhance this position.

Progress from chicken genetics to the chicken genome

The chicken has a proud history, both in genetic research and as a source of food. Here we attempt to provide an overview of past contributions of the chicken in both arenas and to link those contributions to the near future from a genetic perspective. Companion articles will discuss current poultry genetics research in greater detail. The chicken was the first animal species in which Mendelian inheritance was demonstrated. A century later, the chicken was the first among farm animals to have its genome sequenced. Between these firsts, the chicken remained a key organism used in genetic research. Breeding programs, based on sound genetic principles, facilitated the global emergence of the chicken meat and egg industries. Concomitantly, the chicken served as a model whose experimental populations and mutant stocks were used in basic and applied studies with broad application to other species, including humans. In this paper, we review some of these contributions, trace the path from the origin of molecular genetics to the sequence of the chicken genome, and discuss the merits of the chicken as a model organism for furthering our understanding of biology.

Integrated maps in quail (Coturnix japonica) confirm the high degree of synteny conservation with chicken (Gallus gallus) despite 35 million years of divergence

BACKGROUND

By comparing the quail genome with that of chicken, chromosome rearrangements that have occurred in these two galliform species over 35 million years of evolution can be detected. From a more practical point of view, the definition of conserved syntenies helps to predict the position of genes in quail, based on information taken from the chicken sequence, thus enhancing the utility of this species in biological studies through a better knowledge of its genome structure. A microsatellite and an Amplified Fragment Length Polymorphism (AFLP) genetic map were previously published for quail, as well as comparative cytogenetic data with chicken for macrochromosomes. Quail genomics will benefit from the extension and the integration of these maps.

RESULTS

The integrated linkage map presented here is based on segregation analysis of both anonymous markers and functional gene loci in 1,050 quail from three independent F2 populations. Ninety-two loci are resolved into 14 autosomal linkage groups and a Z chromosome-specific linkage group, aligned with the quail AFLP map. The size of linkage groups ranges from 7.8 cM to 274.8 cM. The total map distance covers 904.3 cM with an average spacing of 9.7 cM between loci. The coverage is not complete, as macrochromosome CJA08, the gonosome CJAW and 23 microchromosomes have no marker assigned yet. Significant sequence identities of quail markers with chicken enabled the alignment of the quail linkage groups on the chicken genome sequence assembly. This, together with interspecific Fluorescence In Situ Hybridization (FISH), revealed very high similarities in marker order between the two species for the eight macrochromosomes and the 14 microchromosomes studied.

CONCLUSION

Integrating the two microsatellite and the AFLP quail genetic maps greatly enhances the quality of the resulting information and will thus facilitate the identification of Quantitative Trait Loci (QTL). The alignment with the chicken chromosomes confirms the high conservation of gene order that was expected between the two species for macrochromosomes. By extending the comparative study to the microchromosomes, we suggest that a wealth of information can be mined in chicken, to be used for genome analyses in quail.

Development of a chicken 5 K microarray targeted towards immune function

BACKGROUND

The development of microarray resources for the chicken is an important step in being able to profile gene expression changes occurring in birds in response to different challenges and stimuli. The creation of an immune-related array is highly valuable in determining the host immune response in relation to infection with a wide variety of bacterial and viral diseases.

RESULTS

Here we report the development of chicken immune-related cDNA libraries and the subsequent construction of a microarray containing 5190 elements (in duplicate). Clones on the array originate from tissues known to contain high levels of cells related to the immune system, namely Bursa, Peyers patch, thymus and spleen. Represented on the array are genes that are known to cluster with existing chicken ESTs as well as genes that are unique to our libraries. Some of these genes have no known homologies and represent novel genes in the chicken collection. A series of reference genes (ie. genes of known immune function) are also present on the array. Functional annotation data is also provided for as many of the genes on the array as is possible.

CONCLUSION

Six new chicken immune cDNA libraries have been created and nearly 10,000 sequences submitted to GenBank [GenBank: AM063043-AM071350; AM071520-AM072286; AM075249-AM075607]. A 5 K immune-related array has been developed from these libraries. Individual clones and arrays are available from the ARK-Genomics resource centre.

Chromosomal assignment of chicken clone contigs by extending the consensus linkage map

The bacterial artificial clone-based physical map for chicken plays an important role in the integration of the consensus linkage map and the whole-genome shotgun sequence. It also provides a valuable resource for clone selection within applications such as fluorescent in situ hybridization and positional cloning. However, a substantial number of clone contigs have not yet been assigned to a chromosomal location or have an ambiguous chromosome assignment. In this study, 86 single nucleotide polymorphism markers derived from 86 clones were mapped on the genetic map. These markers added anchoring information for 56 clone contigs and 13 individual clones, covering a total of 57,145 clones.

The chicken genome and the developmental biologist

Recently the initial draft sequence of the chicken genome was released. The reasons for sequencing the chicken were to boost research and applications in agriculture and medicine, through its use as a model of vertebrate development. In addition, the sequence of the chicken would provide an important anchor species in the phylogenetic study of genome evolution. The chicken genome project has its roots in a decade of map building by genetic and physical mapping methods. Chicken genetic markers for map building have generally depended on labour intensive screening procedures. In recent years this has all changed with the availability of over 450,000 EST sequences, a draft sequence of the entire chicken genome and a map of over 1 million SNPs. Clearly, the future for the chicken genome and developmental biology is an exciting one. Through the integration of these resources, it will be possible to solve challenging scientific questions exploiting the power of a chicken model. In this paper we review progress in chicken genomics and discuss how the new tools and information on the chicken genome can help the developmental biologists now and in the future.

A radiation hybrid map of chicken Chromosome 4

The mapping resolution of the physical map for chicken Chromosome 4 (GGA4) was improved by a combination of radiation hybrid (RH) mapping and bacterial artificial chromosome (BAC) mapping. The ChickRH6 hybrid panel was used to construct an RH map of GGA4. Eleven microsatellites known to be located on GGA4 were included as anchors to the genetic linkage map for this chromosome. Based on the known conserved synteny between GGA4 and human Chromosomes 4 and X, sequences were identified for the orthologous chicken genes from these human chromosomes by BLAST analysis. These sequences were subsequently used for the development of STS markers to be typed on the RH panel. Using a logarithm of the odds (LOD) threshold of 5.0, nine linkage groups could be constructed which were aligned with the genetic linkage map of this chromosome. The resulting RH map consisted of the 11 microsatellite markers and 50 genes. To further increase the number of genes on the map and to provide additional anchor points for the physical BAC map of this chromosome, BAC clones were identified for 22 microsatellites and 99 genes. The combined RH and BAC mapping approach resulted in the mapping of 61 genes on GGA4 increasing the resolution of the chicken-human comparative map for this chromosome. This enhanced comparative mapping resolution enabled the identification of multiple rearrangements between GGA4 and human Chromosomes 4q and Xp.

A physical map of the chicken genome

Strategies for assembling large, complex genomes have evolved to include a combination of whole-genome shotgun sequencing and hierarchal map-assisted sequencing. Whole-genome maps of all types can aid genome assemblies, generally starting with low-resolution cytogenetic maps and ending with the highest resolution of sequence. Fingerprint clone maps are based upon complete restriction enzyme digests of clones representative of the target genome, and ultimately comprise a near-contiguous path of clones across the genome. Such clone-based maps are used to validate sequence assembly order, supply long-range linking information for assembled sequences, anchor sequences to the genetic map and provide templates for closing gaps. Fingerprint maps are also a critical resource for subsequent functional genomic studies, because they provide a redundant and ordered sampling of the genome with clones. In an accompanying paper we describe the draft genome sequence of the chicken, Gallus gallus, the first species sequenced that is both a model organism and a global food source. Here we present a clone-based physical map of the chicken genome at 20-fold coverage, containing 260 contigs of overlapping clones. This map represents approximately 91% of the chicken genome and enables identification of chicken clones aligned to positions in other sequenced genomes.

Genetic variants for chick biology research: from breeds to mutants

The availability of the draft sequence of the chicken genome will undoubtedly propel an already important vertebrate research model, the domestic chicken, to a new level. This review describes aspects of chicken natural history and cross-disciplinary biological value. The diversity of extant genetic variants available to researchers is reviewed along with institutional stock locations for North America. An overview of the problem of lack of long-term stability for these resources is presented.