First report on chicken genes and chromosomes 2000

Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution

We present here a draft genome sequence of the red jungle fowl, Gallus gallus. Because the chicken is a modern descendant of the dinosaurs and the first non-mammalian amniote to have its genome sequenced, the draft sequence of its genome--composed of approximately one billion base pairs of sequence and an estimated 20,000-23,000 genes--provides a new perspective on vertebrate genome evolution, while also improving the annotation of mammalian genomes. For example, the evolutionary distance between chicken and human provides high specificity in detecting functional elements, both non-coding and coding. Notably, many conserved non-coding sequences are far from genes and cannot be assigned to defined functional classes. In coding regions the evolutionary dynamics of protein domains and orthologous groups illustrate processes that distinguish the lineages leading to birds and mammals. The distinctive properties of avian microchromosomes, together with the inferred patterns of conserved synteny, provide additional insights into vertebrate chromosome architecture.

Wide genome comparisons reveal the origins of the human X chromosome

The eutherian X chromosome has one of the most conserved gene arrangements in mammals. Although earlier comparisons with distantly related mammalian groups pointed towards separate origins for the short and long arms, much deeper comparisons are now possible using draft sequences of the chicken genome, in combination with genome sequences from pufferfish and zebrafish. This enables surprising new insights into the origins of the mammalian X chromosome.

Molecular cytogenetic definition of the chicken genome: the first complete avian karyotype

Chicken genome mapping is important for a range of scientific disciplines. The ability to distinguish chromosomes of the chicken and other birds is thus a priority. Here we describe the molecular cytogenetic characterization of each chicken chromosome using chromosome painting and mapping of individual clones by FISH. Where possible, we have assigned the chromosomes to known linkage groups. We propose, on the basis of size, that the NOR chromosome is approximately the size of chromosome 22; however, we suggest that its original assignment of 16 should be retained. We also suggest a definitive chromosome classification system and propose that the probes developed here will find wide utility in the fields of developmental biology, DT40 studies, agriculture, vertebrate genome organization, and comparative mapping of avian species.

Assignment of fourteen microsatellite markers to the chicken linkage map

A large mapping population, with 874 F2 individuals, was generated by reciprocally intercrossing 2 chicken lines. A genetic map of 2,426.6 cM comprising 25 linkage groups was established based on 145 microsatellite markers. Chromosome locations were assigned for 14 previously unmapped markers. The marker ADL0132 was previously mapped to chromosome 9; however, here close linkage to the MCW0091 marker on chromosome 4 was found. With this exception, the derived linkage map was in excellent agreement with the chicken consensus map. A comparison with the chicken genome assembly (http://genome.ucsc.edu; February 2004) suggested a few minor errors in the assembly. A PCR-RFLP test was used to genotype a single nucleotide polymorphism in the melanocortin receptor 3 (MC3R) gene in the intercross, and pyrosequencing was used to map the genes for Hemopoetic Cell Kinase (HCK) and Bone Morphogenic Protein 7 (BMP7). The HCK and BMP7 genes on linkage group E32 showed significant linkage to MC3R on the distal end of linkage group E47W24, consequently joining the 2 linkage groups. A comparison between the linkage data in the current study and the physical location of markers as revealed in the chicken genome sequence assembly (February 2004) showed a 3-fold higher recombination rate on microchromosomes than on macrochromosomes.

Plasticity of human chromosome 3 during primate evolution

Comparative mapping of more than 100 region-specific clones from human chromosome 3 in Bornean and Sumatran orangutans, siamang gibbon, and Old and New World monkeys allowed us to reconstruct ancestral simian and hominoid chromosomes. A single paracentric inversion derives chromosome 1 of the Old World monkey Presbytis cristata from the simian ancestor. In the New World monkey Callithrix geoffroyi and siamang, the ancestor diverged on multiple chromosomes, through utilizing different breakpoints. One shared and two independent inversions derive Bornean orangutan 2 and human 3, implying that neither Bornean orangutans nor humans have conserved the ancestral chromosome form. The inversions, fissions, and translocations in the five species analyzed involve at least 14 different evolutionary breakpoints along the entire length of human 3; however, particular regions appear to be more susceptible to chromosome reshuffling. The ancestral pericentromeric region has promoted both large-scale and micro-rearrangements. Small segments homologous to human 3q11.2 and 3q21.2 were repositioned intrachromosomally independent of the surrounding markers in the orangutan lineage. Breakage and rearrangement of the human 3p12.3 region were associated with extensive intragenomic duplications at multiple orangutan and gibbon subtelomeric sites. We propose that new chromosomes and genomes arise through large-scale rearrangements of evolutionarily conserved genomic building blocks and additional duplication, amplification, and/or repositioning of inherently unstable smaller DNA segments contained within them.

Subcellular and molecular mechanisms regulating anti-Müllerian hormone gene expression in mammalian and nonmammalian species

Anti-Müllerian hormone (AMH) is best known for its role as an inhibitor of the development of female internal genitalia primordia during fetal life. In the testis, AMH is highly expressed by Sertoli cells of the testis from early fetal life to puberty, when it is downregulated by the action of testosterone, acting through the androgen receptor, and meiotic spermatocytes, probably acting through TNFalpha. Basal expression of AMH is induced by SOX9; GATA4, SF1, and WT1 enhance SOX9-activated expression. When the hypothalamic-pituitary axis is active and the negative effect of androgens and germ cells is absent, for example, in the fetal and neonatal periods or in disorders like androgen insensitivity, FSH upregulates AMH expression through a nonclassical cAMP-PKA pathway involving transcription factors AP2 and NFkappaB. The maintenance and hormonal regulation of AMH expression in late fetal and postnatal life requires distal AMH promoter sequences. In the ovary, granulosa cells express AMH from late fetal life at low levels; DAX1 and FOG2 seem to be responsible for negatively modulating AMH expression. Particular features are observed in AMH expression in nonmammalian species. In birds, AMH is expressed both in the male and female fetal gonads, and, like in reptiles, its expression is not preceded by that of SOX9.

A high-resolution radiation hybrid map of chicken chromosome 5 and comparison with human chromosomes

BACKGROUND

The resolution of radiation hybrid (RH) maps is intermediate between that of the genetic and BAC (Bacterial Artificial Chromosome) contig maps. Moreover, once framework RH maps of a genome have been constructed, a quick location of markers by simple PCR on the RH panel is possible. The chicken ChickRH6 panel recently produced was used here to construct a high resolution RH map of chicken GGA5. To confirm the validity of the map and to provide valuable comparative mapping information, both markers from the genetic map and a high number of ESTs (Expressed Sequence Tags) were used. Finally, this RH map was used for testing the accuracy of the chicken genome assembly for chromosome 5.

RESULTS

A total of 169 markers (21 microsatellites and 148 ESTs) were typed on the ChickRH6 RH panel, of which 134 were assigned to GGA5. The final map is composed of 73 framework markers extending over a 1315.6 cR distance. The remaining 61 markers were placed alongside the framework markers within confidence intervals.

CONCLUSION

The high resolution framework map obtained in this study has markers covering the entire chicken chromosome 5 and reveals the existence of a high number of rearrangements when compared to the human genome. Only two discrepancies were observed in relation to the sequence assembly recently reported for this chromosome.

A high-resolution radiation hybrid map of chicken chromosome 5 and comparison with human chromosomes

Background

The resolution of radiation hybrid (RH) maps is intermediate between that of the genetic and BAC (Bacterial Artificial Chromosome) contig maps. Moreover, once framework RH maps of a genome have been constructed, a quick location of markers by simple PCR on the RH panel is possible. The chicken ChickRH6 panel recently produced was used here to construct a high resolution RH map of chicken GGA5. To confirm the validity of the map and to provide valuable comparative mapping information, both markers from the genetic map and a high number of ESTs (Expressed Sequence Tags) were used. Finally, this RH map was used for testing the accuracy of the chicken genome assembly for chromosome 5.

Results

A total of 169 markers (21 microsatellites and 148 ESTs) were typed on the ChickRH6 RH panel, of which 134 were assigned to GGA5. The final map is composed of 73 framework markers extending over a 1315.6 cR distance. The remaining 61 markers were placed alongside the framework markers within confidence intervals.

Conclusion

The high resolution framework map obtained in this study has markers covering the entire chicken chromosome 5 and reveals the existence of a high number of rearrangements when compared to the human genome. Only two discrepancies were observed in relation to the sequence assembly recently reported for this chromosome.

Sexual differentiation of the zebra finch song system

The song system of zebra finches (Taeniopygia gutatta) is highly sexually dimorphic. Only males sing, and the brain regions and muscles controlling song are much larger in males than in females. Development of the song system is highly sensitive to steroid hormones. However, unlike similar sexually dimorphic systems in other animal models, masculinization of song system structure and function is most likely not induced by testosterone secreted from the testes. Instead, sex-specific development of the neural song system appears to be regulated by factors intrinsic to the brain, probably by the expression of sex chromosome gene(s) that influence the levels of estradiol synthesized in the brain and/or the responses of brain tissue to estradiol. However, the existing data are complex and in some cases contradictory. More work is required to identify the critical genes and their relationships with steroid hormones.

From biomedicine to natural history research: EST resources for ambystomatid salamanders

BACKGROUND

Establishing genomic resources for closely related species will provide comparative insights that are crucial for understanding diversity and variability at multiple levels of biological organization. We developed ESTs for Mexican axolotl (Ambystoma mexicanum) and Eastern tiger salamander (A. tigrinum tigrinum), species with deep and diverse research histories.

RESULTS

Approximately 40,000 quality cDNA sequences were isolated for these species from various tissues, including regenerating limb and tail. These sequences and an existing set of 16,030 cDNA sequences for A. mexicanum were processed to yield 35,413 and 20,599 high quality ESTs for A. mexicanum and A. t. tigrinum, respectively. Because the A. t. tigrinum ESTs were obtained primarily from a normalized library, an approximately equal number of contigs were obtained for each species, with 21,091 unique contigs identified overall. The 10,592 contigs that showed significant similarity to sequences from the human RefSeq database reflected a diverse array of molecular functions and biological processes, with many corresponding to genes expressed during spinal cord injury in rat and fin regeneration in zebrafish. To demonstrate the utility of these EST resources, we searched databases to identify probes for regeneration research, characterized intra- and interspecific nucleotide polymorphism, saturated a human - Ambystoma synteny group with marker loci, and extended PCR primer sets designed for A. mexicanum / A. t. tigrinum orthologues to a related tiger salamander species.

CONCLUSIONS

Our study highlights the value of developing resources in traditional model systems where the likelihood of information transfer to multiple, closely related taxa is high, thus simultaneously enabling both laboratory and natural history research.

Close linkage of genes encoding receptors for subgroups A and C of avian sarcoma/leucosis virus on chicken chromosome 28

Avian sarcoma and leucosis viruses (ASLV) are classified into six major subgroups (A to E and J) according to the properties of the viral envelope proteins and the usage of cellular receptors for virus entry. Subgroup A and B receptors are identified molecularly and their genomic positions TVA and TVB are mapped. The subgroup C receptor is unknown, its genomic locus TVC is reported to be genetically linked to TVA, which resides on chicken chromosome 28. In this study, we used two chicken inbred lines that carry different alleles coding for resistance (TVC(R) and sensitivity (TVC(S)) to infection by subgroup C viruses. A backross population of these lines was tested for susceptibility to subgroup C infection and genotyped for markers from chicken chromosome 28. We confirmed the close linkage between TVA and TVC loci. Further, we have described the position of TVC on chromosome 28 relative to markers from the consensus map of the chicken genome.

The repetitive landscape of the chicken genome

Cot-based cloning and sequencing (CBCS) is a powerful tool for isolating and characterizing the various repetitive components of any genome, combining the established principles of DNA reassociation kinetics with high-throughput sequencing. CBCS was used to generate sequence libraries representing the high, middle, and low-copy fractions of the chicken genome. Sequencing high-copy DNA of chicken to about 2.7 x coverage of its estimated sequence complexity led to the initial identification of several new repeat families, which were then used for a survey of the newly released first draft of the complete chicken genome. The analysis provided insight into the diversity and biology of known repeat structures such as CR1 and CNM, for which only limited sequence data had previously been available. Cot sequence data also resulted in the identification of four novel repeats (Birddawg, Hitchcock, Kronos, and Soprano), two new subfamilies of CR1 repeats, and many elements absent from the chicken genome assembly. Multiple autonomous elements were found for a novel Mariner-like transposon, Galluhop, in addition to nonautonomous deletion derivatives. Phylogenetic analysis of the high-copy repeats CR1, Galluhop, and Birddawg provided insight into two distinct genome dispersion strategies. This study also exemplifies the power of the CBCS method to create representative databases for the repetitive fractions of genomes for which only limited sequence data is available.

The repetitive landscape of the chicken genome

Cot-based cloning and sequencing (CBCS) is a powerful tool for isolating and characterizing the various repetitive components of any genome, combining the established principles of DNA reassociation kinetics with high-throughput sequencing. CBCS was used to generate sequence libraries representing the high, middle, and low-copy fractions of the chicken genome. Sequencing high-copy DNA of chicken to about 2.7 x coverage of its estimated sequence complexity led to the initial identification of several new repeat families, which were then used for a survey of the newly released first draft of the complete chicken genome. The analysis provided insight into the diversity and biology of known repeat structures such as CR1 and CNM, for which only limited sequence data had previously been available. Cot sequence data also resulted in the identification of four novel repeats (Birddawg, Hitchcock, Kronos, and Soprano), two new subfamilies of CR1 repeats, and many elements absent from the chicken genome assembly. Multiple autonomous elements were found for a novel Mariner-like transposon, Galluhop, in addition to nonautonomous deletion derivatives. Phylogenetic analysis of the high-copy repeats CR1, Galluhop, and Birddawg provided insight into two distinct genome dispersion strategies. This study also exemplifies the power of the CBCS method to create representative databases for the repetitive fractions of genomes for which only limited sequence data is available.

Detection and localization of quantitative trait loci affecting fatness in broilers

A cross between 2 genetically different outcross broiler dam lines, originating from the White Plymouth Rock breed, was used to produce a large 3-generation broiler population. This population was used to detect and localize QTL affecting fatness in chicken. Twenty full-sib birds in generation 1 and 456 full-sib birds in generation 2 were typed for microsatellite markers, and phenotypic observations were collected for 3 groups of generation 3 birds (approximately 1,800 birds per group). Body weight, abdominal fat weight, and percentage abdominal fat was recorded at the age of 7, 9, and 10 wk. To study the presence of QTL, an across-family weighted regression interval mapping approach was used in a full-sib QTL analysis. Genotypes from 410 markers mapped on 25 chromosomes were available. For the 3 traits, 26 QTL were found for 18 regions on 12 chromosomes. Two genomewise significant QTL (P < 0.05) were detected, one for percentage abdominal fat at the age of 10 wk on chicken chromosome 1 at 241 cM (MCW0058 to MCW0101) with a test statistic of 2.75 and the other for BW at the age of 10 wk on chicken chromosome 13 at 9 cM (MCW0322 to MCW0110) with a test statistic of 2.77. Significance levels were obtained using the permutation test. Multiple suggestive QTL were found on chromosomes 1, 2, 4, 13, 15, and 18, whereas chromosomes 3, 7, 10, 11, 14, and 27 had a single suggestive QTL.

A Comprehensive Screen for Chicken Proteins that Interact with Proteins Unique to Virulent Strains of Marek's Disease Virus

Genetic resistance to Marek's disease (MD) has been proposed as a method to augment current vaccinal control of MD. Although it is possible to identify QTL and candidate genes that are associated with MD resistance, it is necessary to integrate functional screens with linkage analysis to confirm the identity of true MD resistance genes. To help achieve this objective, a comprehensive 2-hybrid screen was conducted using genes unique to virulent Marek's disease virus (MDV) strains. Potential MDV-host protein interactions were tested by an in vitro binding assay to confirm the initial two-hybrid results. As a result, 7 new MDV-chicken protein interactions were identified and included the chicken proteins MHC class II beta (BLB) and invariant (Ii) chain (CD74), growth-related translationally controlled tumor protein (TPT1), complement component Clq-binding protein (C1QBP), retinoblastoma-binding protein 4 (RBBP4), and alpha-enolase (ENO1). Mapping of the encoding chicken genes suggests that BLB, the gene for MHC class II beta chain, is a positional candidate gene. In addition, the known functions of the chicken proteins suggest mechanisms that MDV might use to evade the chicken immune system and alter host gene regulation. Taken together, our results indicate that integrated genomic methods provide a powerful strategy to gain insights on complex biological processes and yield a manageable number of genes and pathways for further characterization.

A first-generation microsatellite linkage map of the Japanese quail

A linkage map of the Japanese quail (Coturnix japonica) genome was constructed based upon segregation analysis of 72 microsatellite loci in 433 F(2) progeny of 10 half-sib families obtained from a cross between two quail lines of different genetic origins. One line was selected for long duration of tonic immobility, a behavioural trait related to fearfulness, while the other was selected based on early egg production. Fifty-eight of the markers were resolved into 12 autosomal linkage groups and a Z chromosome-specific linkage group, while the remaining 14 markers were unlinked. The linkage groups range from 8 cM (two markers) to 206 cM (16 markers) and cover a total map distance of 576 cM with an average spacing of 10 cM between loci. Through comparative mapping with chicken (Gallus gallus) using orthologous markers, we were able to assign linkage groups CJA01, CJA02, CJA05, CJA06, CJA14 and CJA27 to chromosomes. This map, which is the first in quail based solely on microsatellites, is a major step towards the development of a quality molecular genetic map for this valuable species. It will provide an important framework for further genetic mapping and the identification of quantitative trait loci controlling egg production and fear-related behavioural traits in quail.

DNA sequence and analysis of human chromosome 9

Chromosome 9 is highly structurally polymorphic. It contains the largest autosomal block of heterochromatin, which is heteromorphic in 6-8% of humans, whereas pericentric inversions occur in more than 1% of the population. The finished euchromatic sequence of chromosome 9 comprises 109,044,351 base pairs and represents >99.6% of the region. Analysis of the sequence reveals many intra- and interchromosomal duplications, including segmental duplications adjacent to both the centromere and the large heterochromatic block. We have annotated 1,149 genes, including genes implicated in male-to-female sex reversal, cancer and neurodegenerative disease, and 426 pseudogenes. The chromosome contains the largest interferon gene cluster in the human genome. There is also a region of exceptionally high gene and G + C content including genes paralogous to those in the major histocompatibility complex. We have also detected recently duplicated genes that exhibit different rates of sequence divergence, presumably reflecting natural selection.

Male-biased mutation rate and divergence in autosomal, z-linked and w-linked introns of chicken and Turkey

To investigate mutation-rate variation between autosomes and sex chromosomes in the avian genome, we have analyzed divergence between chicken (Gallus gallus) and turkey (Meleagris galopavo) sequences from 33 autosomal, 28 Z-linked, and 14 W-linked introns with a total ungapped alignment length of approximately 43,000 bp. There are pronounced differences in the mean divergence among autosomes and sex chromosomes (autosomes [A] = 10.08%, Z chromosome = 10.99%, and W chromosome = 5.74%), and we use these data to estimate the male-to-female mutation-rate ratio (alpha(m)) from Z/A, Z/W, and A/W comparisons at 1.71, 2.37, and 2.52, respectively. Because the alpha(m) estimates of the three comparisons do not differ significantly, we find no statistical support for a specific reduction in the Z chromosome mutation rate (Z reduction estimated at 4.89%, P = 0.286). The idea of mutation-rate reduction in the sex chromosome hemizygous in one sex (i.e., X in mammals, Z in birds) has been suggested on the basis of theory on adaptive mutation-rate evolution. If it exists in birds, the effect would, thus, seem to be weak; a preliminary power analysis suggests that it is significantly less than 18%. Because divergence may vary within chromosomal classes as a result of variation in mutation and/or selection, we developed a novel double-bootstrapping method, bootstrapping both by introns and sites from concatenated alignments, to estimate confidence intervals for chromosomal class rates and for alpha(m). The narrowest interval for the alpha(m) estimate is 1.88 to 2.97 from the Z/W comparison. We also estimated alpha(m) using maximum likelihood on data from all three chromosome classes; this method yielded alpha(m) = 2.47 and approximate 95% confidence intervals of 2.27 to 2.68. Our data are broadly consistent with the idea that mutation-rate differences between chromosomal classes can be explained by the male mutation bias alone.

Molecular and cytogenetic organization of the 5S ribosomal DNA array in chicken (Gallus gallus)

The 5S ribosomal (r) RNA genes encode a small (approximately 120-bp) highly-conserved component of the large ribosomal subunit. The objective of the present research was to study the molecular and cytogenetic organization of the chicken 5S rDNA. A predominant 2.2-kb gene (5Salpha) consisting of a coding and intergenic spacer (IGS) region was identified in ten research and commercial populations. A variant gene repeat of 0.6kb (5Sbeta) was observed in some of the populations. Genetic linkage analysis and cytogenetic localization by fluorescence in-situ hybridization assigned the 5S rDNA to chromosome 9. The 5S rDNA array was determined to be 80.2 +/- 7.0 kb upon electrophoretic sizing following EcoRV digestion. Sequence analysis of 5Salpha IGS regions revealed considerable conservation between chicken subspecies (98.4% identity) as well as homology with vertebrate Pol III promoter and regulatory sequence motifs. Minor intraindividual sequence variation within 1000 bp of IGS was observed in four cloned Red Jungle Fowl (Gallus gallus gallus) 5Salpha repeats (95.5% identity in this region). Sequence comparisons between IGS regions of 5Salpha and 5Sbeta genes indicated two short continuous (>20bp) and many short non-continuous homologous regions as well as other conserved features such as promoter and termination motifs.

The telomeric protein Rap1 is conserved in vertebrates and is expressed from a bidirectional promoter positioned between the Rap1 and KARS genes

We have identified the chicken homolog of the mammalian telomere protein repression and activation protein 1 (Rap1). Although cRap1 has only 36% sequence identity to hRap1, it contains the same conserved BRCA1 C-terminal (BRCT), Myb and Rap C-terminus (RCT) domains. Two-hybrid analysis and immunolocalization experiments revealed that cRap1 interacts with the telomere-binding protein telomeric repeat binding factor (TRF)2 and localizes to telomeres. Thus, despite considerable sequence divergence, the identity and overall domain structure of telomere-associated proteins is conserved in vertebrates. Analysis of the cRap1 genomic locus revealed that the cRap1 gene lies immediately adjacent to the cKARS (lysyl-tRNA synthetase) gene with the two genes in a head-to-head orientation separated by only 57 nt. This same organization is conserved at the human Rap1-KARS locus. When 5' regions of the cRap1 and cKARS genes were tested for promoter activity, the promoters of both genes were found to lie in or near the intergenic spacer. The two promoters lack TATA boxes but appear to have downstream promoter elements (DPEs). Analysis of human Rap1 and KARS expressed sequence tags (ESTs) indicated that this localization of TATA-less promoters to the intergenic spacer is a conserved feature of the Rap1-KARS locus.

Preliminary linkage map of the turkey (Meleagris gallopavo) based on microsatellite markers

The turkey is an agriculturally important species for which, until now, there is no published genetic linkage map based on microsatellite markers--still the markers most used in the chicken and other farm animals. In order to increase the number of markers on a turkey genetic linkage map we decided to map new microsatellite sequences obtained from a GT-enriched turkey genomic library. In different chicken populations more than 35-55% of microsatellites are polymorphic. In the turkey populations tested here, 43% of all turkey primers tested were found to be polymorphic, in both commercial and wild type turkeys. Twenty linkage groups (including the Z chromosome) containing 74 markers have been established, along with 37 other unassigned markers. This map will lay the foundations for further genetic mapping and the identification of genes and quantitative trait loci in this economically important species. Genome comparisons, based on genetic maps, with related species such as the chicken would then also be possible. All primer information, polymerase chain reaction (PCR) conditions, allele sizes and genetic linkage maps can be viewed at http://roslin.thearkdb.org/. The DNA is also available on request through the Roslin Institute.

Comparative map between chicken chromosome 15 and human chromosomal region 12q24 and 22q11-q12

The physical and comparative map of GGA15 was improved by the construction of 9 BAC contigs around loci previously mapped on GGA15 by linkage analysis. In total, 240 BAC clones were isolated, covering 30-35% of GGA15, and 120 STS were developed (104 STS derived from BAC end sequences and 18 STS derived within genes). Seventeen chicken orthologues of human genes located on human Chr 22q11-q12 were directly mapped within BAC contigs of GGA15. Furthermore, the partial sequences of the chicken BAC clones were compared with sequences present in the EMBL/GenBank databases and revealed matches to 26 genes, ESTs, and genomic clones located on HSA22q11-q12 and HSA12q24. These results provide a better alignment of GGA15 with the corresponding regions in human and mouse, and improve our knowledge of the evolution and dynamics of the vertebrate genome.

Association of twelve candidate gene polymorphisms and response to challenge with Salmonella enteritidis in poultry

Breeding for disease resistance to Salmonella enteritidis (SE) could be an effective approach to control Salmonella in poultry. The candidate gene approach is a useful method to investigate genes that are involved in genetic resistance. In this study, 12 candidate genes that are involved in the pathogenesis of Salmonella infection were investigated using five different genetic groups of meat-type chicken. The genes were natural resistance associated macrophage protein 1 (SLC11A1, previously known as NRAMP1), inhibitor of apoptosis protein 1 (IAP1), prosaposin (PSAP), Caspase-1 (CASP1), inducible nitric oxide production (iNOS), interferon-gamma (IFNG), interleukin-2 (IL2), immunoglobulin light chain (IGL), ZOV3, and transforming growth factors B2, B3 and B4 (TGFB2, B3 and B4). In total, 117 birds of all groups were challenged with SE at the age of 3 weeks. In all birds at 7-day post-infection SE load in caecum content, spleen and liver were quantified. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assays were used to genotype all animals for each gene. Overall we found the most significant associations with caecum content, nine of 12 genes showed a significant association (SLC11A1, IAP1, PSAP, CASP1, iNOS, IL2, IGL, TGFB2 and TGFB4). For liver, five genes (SLC11A1, CASP1, IL2, IGL, and TGFB4) and for spleen, only one gene (TGFB3) showed a significant association with SE load. By showing associations of 12 PCR-RFLP assays with SE load after a pathogen challenge, this study confirmed the polygenic nature of disease resistance to SE.

A BAC-based physical map of the chicken genome

A genome-wide physical map constructed with bacterial artificial chromosomes (BACs) is an essential component in linking phenotypic traits to the responsible genetic variation in the genomes of plants and animals. We have constructed a physical map of the chicken genome from 57,091 BACs (7.9-fold haploid genome coverage) by restriction fingerprint analysis using high-resolution polyacrylamide gel electrophoresis. The physical map consists of 2331 overlapping BAC contigs and is estimated to span 1510 Mb in physical length. BAC contigs were verified manually and by screening the BACs with 367 DNA markers. A total of 361 of the contigs have been anchored to the existing chicken genetic map. This map represents the first genome-wide, BAC-based physical map of the chicken genome. It provides a powerful platform for many areas of chicken genomics, including targeted marker development, fine mapping of genes and QTL alleles, positional cloning, analysis of avian genome organization and evolution, chicken-mammalian comparative genomics, and large-scale genome sequencing.

Microsatellite Loci for Genetic Mapping in the Turkey (Meleagris gallopavo)

New microsatellite loci for the turkey (Meleagris gallopavo) were developed from two small insert DNA libraries. Polymorphism at these new loci was examined in domestic birds and two resource populations designed for genetic linkage mapping. The majority of loci (152 of 168) was polymorphic in domestic turkeys and informative in two mapping resource populations and thus will be useful for genetic linkage mapping.

The Bh (Black at Hatch) Gene that Causes Abnormal Feather Pigmentation Maps to Chromosome 1 of the Japanese Quail

Japanese quail embryos normally have longitudinal black and brown stripes formed by colored feather buds on their back whereas an autosomal dominant mutation, black at hatch (Bh), disrupts this pigmentation pattern by causing overall black and brown coating in heterozygotes and homozygotes, respectively. These phenotypes of the Bh mutant embryos suggest that the Bh locus plays an important role in the pigment pattern formation of plumage, but its genetic origin, including cloning of the responsible gene, has been insufficiently studied. In this study, we adapted genetically directed representational difference analysis with elimination of excessive clones (GDRDA-WEEC) to Bh quails and isolated two genetic markers linked to the Bh locus as DNA fragments. Cytogenetic study by fluorescence in situ hybridization (FISH) of the DNA fragments used as probes demonstrated that the marker loci were located in the same region on the long arm of chromosome 1. Close genetic linkage between the Bh and the marker loci, and the chromosomal location of the latter suggested that the Bh locus is located on the long-arm of chromosome 1 of the Japanese quail.

The properties of the single chicken MHC classical class II alpha chain ( B-LA) gene indicate an ancient origin for the DR/E-like isotype of class II molecules

In mammals, there are MHC class II molecules with distinctive sequence features, such as the classical isotypes DR, DQ and DP. These particular isotypes have not been reported in non-mammalian vertebrates. We have isolated the class II (B-L) alpha chain from outbred chickens as the basis for the cloning and sequencing of the cDNA. We found only one class II alpha chain transcript, which bears the major features of a classical class II alpha sequence, including the critical peptide-binding residues. The chicken sequence is more similar to human DR than to the DQ, DP, DO or DM isotypes, most significantly in the peptide-binding alpha(1) domain. The cDNA and genomic DNA sequences from chickens of diverse origins show few alleles, which differ in only four nucleotides and one amino acid. In contrast, significant restriction fragment length polymorphism is detected by Southern blot analysis of genomic DNA, suggesting considerable diversity around the gene. Analysis of a large back-cross family indicates that the class II alpha chain locus ( B-LA) is located roughly 5.6 cM from the MHC locus, which encodes the classical class II beta chains. Thus the chicken class II alpha chain gene is like the mammalian DR and E isotypes in three properties: the presence of the critical peptide-binding residues, the low level of polymorphism and sequence diversity, and the recombinational separation from the class II beta chain genes. These results indicate that the sequence features of this lineage are both functionally important and at least 300 million years old.

A first-generation map of the turkey genome

A primary linkage map of the domestic turkey (Meleagris gallopavo) was developed by segregation analysis of genetic markers within a backcross family. This reference family includes 84 offspring from one F1sire mated to two dams. Genomic DNA was digested using one of five restriction enzymes, and restriction fragment length polymorphisms were detected on Southern blots using probes prepared from 135 random clones isolated from a whole-embryo cDNA library. DNA sequence was subsequently determined for 114 of these cDNA clones. Sequence comparisons were done using BLAST searches of the GenBank database, and redundant sequences were eliminated. High similarity was found between 23% of the turkey sequences and mRNA sequences reported for the chicken. The current map, based on expressed genes, includes 138 loci, encompassing 113 loci arranged into 22 linkage groups and an additional 25 loci that remain unlinked. The average distance between linked markers is 6 cM and the longest linkage group (17 loci) measures 131 cM. The total map distance contained within linkage groups is 651 cM. The present map provides an important framework for future genome mapping in the turkey.Key words: genetic map, Meleagris gallopavo, expressed sequence tag, RFLP.

Melanocortin 1-receptor (MC1R) mutations are associated with plumage colour in chicken

The co-segregation of plumage colour and sequence polymorphism in the melanocortin 1-receptor gene (MC1R) was investigated using an intercross between the red junglefowl and White Leghorn chickens. The results provided compelling evidence that the Extended black (E) locus controlling plumage colour is equivalent to MC1R. E/MC1R was assigned to chromosome 11 with overwhelming statistical support. Sequence analysis indicated that the E92K substitution, causing a constitutively active receptor in the sombre mouse, is the most likely causative mutation for the Extended black allele carried by the White Leghorn founders in this intercross. The MC1R sequence associated with the recessive buttercup (ebc) allele indicated that this allele evolved from a dominant Extended black allele as it shared the E92K and M71T substitutions with some E alleles. It also carried a third missense mutation H215P which thus may interfere with the constitutive activation of the receptor caused by E92K (and possibly M71T).

The twofold difference in adult size between the red junglefowl and White Leghorn chickens is largely explained by a limited number of QTLs

A large intercross between the domestic White Leghorn chicken and the wild ancestor, the red junglefowl, has been used in a Quantitative Trait Loci (QTL) study of growth and egg production. The linkage map based on 105 marker loci was in good agreement with the chicken consensus map. The growth of the 851 F2 individuals was lower than both parental lines prior to 46 days of age and intermediate to the two parental lines thereafter. The QTL analysis of growth traits revealed 13 loci that showed genome-wide significance. The four major growth QTLs explained 50 and 80% of the difference in adult body weight between the founder populations for females and males, respectively. A major QTL for growth, located on chromosome 1 appears to have pleiotropic effects on feed consumption, egg production and behaviour. There was a strong positive correlation between adult body weight and average egg weight. However, three QTLs affecting average egg weight but not body weight were identified. An interesting observation was that the estimated effects for the four major growth QTLs all indicated a codominant inheritance.

Genetic markers and their application in poultry breeding

The current chicken genetic map contains at least 1,965 loci within 50 linkage groups, and it covers about 4,000 cM. About 235 of these loci have homology with known human or mammalian genes. The remaining loci are anonymous molecular DNA markers, including microsatellites, amplified fragment length polymorphism (AFLP), randomly amplified polymorphic DNA (RAPD), CR1 elements, and others. A third generation genetic map for human uses single nucleotide polymorphisms (SNP), which have allowed the mapping of complex traits by linkage disequilibrium. One advantage of SNP is that they are usually linked to the gene of interest, and association of the SNP with traits of economic importance can be analyzed using candidate gene approaches. With the tremendous advancements in characterizing chicken expressed sequence tags (EST), the identification of genetic polymorphisms such as SNP in chicken genes has become a reality. Our laboratory has undertaken an in silico analysis of the chicken EST at the University of Delaware by using a Phred/Phrap/Polyphred/Consed pipeline to identify candidate chicken SNP. Initial scanning of 23,427 chicken EST identified a total of 1,209 candidate SNP, with at least 182 non-synonymous SNP that result in an amino acid change observed. Validation of these candidate chicken SNP is ongoing. Placement of the SNP on the chicken genetic map will enhance marker density, thus allowing for mapping of complex traits through linkage analysis and linkage disequilibrium. Application of SNP to identify disease resistance genes in chickens is of special interest to our laboratory, especially in regards to Marek's disease and coccidiosis.

The dynamics of the genotype-phenotype association

The integrity of an organism is maintained by networks of interacting genes. Such networks predict that genetic variants affect phenotypes in a nonadditive fashion. That is, the effect of an allelic variation in one gene is dependent on the variations in other genes. We summarize the analyses of a series of genes in a White Leghorn strain that support the existence of such gene networks: 1) genes are pleiotropic, 2) genes affect trait correlations, 3) genes affect trait distributions in a nonadditive fashion, 4) genes interact with each other, and 5) genes are at linkage disequilibrium, even when located on different chromosomes. The latter observation indicated that certain gene combinations lead to reduced reproductive fitness. Each candidate genes we analyzed segregated for multiple alleles that affected production traits. This finding was surprising, even for a strain with a large effective population size. The shapes of trait distributions appear to be a better descriptor of gene effects than measures of central tendency. Despite this complexity, it is feasible to conduct DNA-based selection, starting from any of several different genes that affect a trait. Gene networks may be altered in many different ways to improve a particular phenotype, but networks may differ in their effects on other traits.

Telomeres in the chicken: genome stability and chromosome ends

Telomeres are the complex nucleoprotein structures at the termini of linear chromosomes. Telomeric DNA consists of a highly conserved hexanucleotide arranged in tandem repeats. Telomerase, a ribonucleoprotein of the reverse transcriptase family, specifies the sequence of telomeric DNA and maintains telomere array length. Numerous studies in model organisms established the significance of telomere structure and function in regulating genome stability, cellular aging, and oncogenesis. Our overall research objectives are to understand the organization of the telomere arrays in chicken in the context of the unusual organization and specialized features of this higher vertebrate genome (which include a compact genome, numerous microchromosomes, and high recombination rate) and to elucidate the role telomeres play in genome stability impacting cell function and life span. Recent studies found that the chicken genome contains three overlapping size classes of telomere arrays that differ in location and age-related stability: Class I 0.5 to 10 kb, Class II 10 to 40 kb, and Class III 40 kb to 2 Mb. Some notable features of chicken telomere biology are that the chicken genome contains ten times more telomeric DNA than the human genome and the Class III telomere arrays are the largest described for any vertebrate species. In vivo, chicken telomeres (Class II) shorten in an age-related fashion and telomerase activity is high in early stage embryos and developing organs but down-regulates during late embryogenesis or postnatally in most somatic tissues. In vitro, chicken cells down-regulate telomerase activity unless transformed. Knowledge of chicken telomere biology contributes information relevant to present and future biotechnology applications of chickens in vivo and chicken cells in vitro.

Molecular cloning and characterization of avian bombesin-like peptide receptors: new tools for investigating molecular basis for ligand selectivity

(1) Bombesin (BN), originally isolated from amphibians, is structurally related to a family of BN-like peptides found in mammals, which include gastrin-releasing peptide (GRP) and neuromedin B (NMB). These peptides have important effects on secretion, smooth muscle contraction, metabolism and behavior. Here we report cloning and characterization of two subtypes of BN-like peptide receptors in Aves. (2) The amino-acid sequence of chick GRP-R (chGRP-R) is highly identical with mammalian and amphibian GRP-R, and this receptor showed high affinity for GRP, BN and synthetic bombesin agonist, [D-Phe(6), beta-Ala(11), Phe(13), Nle(14)]bombesin(6-14) ([FAFNl]BN(6-14)). The chGRP-R gene was localized to chicken chromosome 1q23distal-q24proximal, where chick homologs of other human X-linked genes have also been mapped. (3) ChBRS-3.5, having sequence similarities to both mammalian bombesin-like peptide receptor subtype-3 and amphibian bombesin-like peptide receptor subtype-4, showed high affinity for [FAFNl]BN(6-14), moderate affinity for BN, but low affinity for both GRP and NMB. (4) Expression of both receptors was detected in brain, but only chGRP-R was expressed in gastrointestinal (GI) tissues. (5) When expressed in Chinese hamster ovary K1 cells, these receptors mediate intracellular calcium mobilization upon agonist stimulation. These results suggest that a novel BN peptide may occur in Aves as an endogenous ligand for chBRS-3.5. (6) The receptor sequences responsible for ligand selectivities were discussed and this knowledge about avian BN-like peptide receptors will help us to understand the molecular basis for agonist sensitivities of BN-like peptide receptors.

Candidate gene approach: potentional association of caspase-1, inhibitor of apoptosis protein-1, and prosaposin gene polymorphisms with response to Salmonella enteritidis challenge or vaccination in young chicks

Salmonella enteritidis (SE) contamination of poultry products is a major cause of foodborne disease worldwide. Caspase-1 and inhibitor of apoptosis protein-1 (IAP-1) were selected as candidate genes for chicken response to SE because their proteins play critical roles in the apoptotic pathway when intracellular bacteria interact with host cells. Prosaposin (PSAP) was selected as a positional candidate gene based on a previous quantitative trait loci (QTL) linkage study using the same population. The F1 offspring of outbred sires crossed with three diverse, highly inbred dam lines (two major histocompatibility complex-congenic Leghorn lines named G-B1 and G-B2, and one Fayoumi line) were used to define the phenotypes. The F1 birds were involved in either pathogenic SE challenge, in which spleen and cecum content bacterial load were quantified, or SE vaccination, in which plasma antibody level to SE vaccine was evaluated. A polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLP) assay was developed to identify single-nucleotide polymorphism (SNP) in the three genes. The F1 offspring of heterozygous sires for each gene were genotyped. The sire caspase-1 gene was significantly associated with cecum content bacterial load (P = 0.04) in the three combined dam line crosses, and with spleen bacterial load in the G-B1 cross (P=0.02). The sire caspase-1 gene was also significantly associated with antibody level to SE vaccine (P=0.03) in F1 males in the three combined dam line crosses. The sire IAP-1 gene was significantly associated with spleen bacterial load (P=0.04) in the three combined dam-line crosses, and interacted with dam-line genetics (P = 0.01) for cecum content bacterial load. The sire PSAP gene significantly interacted with sex for spleen bacterial load (P = 0.004). This study is the first to demonstrate the association of SNPs for caspase-1, IAP-1, and PSAP genes with SE vaccine and with pathogen challenge response in chickens.

Genetic markers associated with antibody response kinetics in adult chickens

A linkage disequilibrium approach with microsatellites was employed to investigate QTL affecting immune response. Highly inbred males of two MHC-congenic Fayoumi chicken lines were mated with highly inbred G-B1 Leghorn hens. Adult F2 hens (n = 158) were injected twice with SRBC and fixed Brucella abortus (BA). Agglutinating antibody titers were measured. Secondary phase parameters of maximum titers (Ymax) and time (Tmax) needed to achieve Ymax were estimated from postsecondary titers by using a nonlinear regression model. A three-step genotype strategy (DNA pooling, selective genotyping, and whole population genotyping) was used to identify microsatellite markers that are associated with immune response to SRBC and BA. The linkage distances between adjacent markers in the F2 population were estimated by Crimap. The QTL affecting immune response to SRBC and BA were detected based on F statistic by interval mapping. A total of five significant QTL, as determined by a permutation test, were detected at the 5% chromosome-wise level on Chromosomes 3, 5, 6, and Z. Two (Chromosome 3 and 6) of five QTL were significant at the 1% chromosome-wise level. The variance explained by the QTL ranged from 6.46 to 7.50%. The results suggest that regions on Chromosomes 3, 5, 6, and Z contain QTL that affect antibody kinetics in the hen.

Analyses of Beta-1 Syntrophin, Syndecan 2 and Gem GTPase as Candidates for Chicken Muscular Dystrophy

Despite intensive studies of muscular dystrophy of chicken, the responsible gene has not yet been identified. Our recent studies mapped the genetic locus for abnormal muscle (AM) of chicken with muscular dystrophy to chromosome 2q using the Kobe University (KU) resource family, and revealed the chromosome region where the AM gene is located has conserved synteny to human chromosome 8q11-24.3, where the beta-1 syntrophin (SNTB1), syndecan 2 (SDC2) and Gem GTPase (GEM) genes are located. It is reasonable to assume those genes might be candidates for the AM gene. In this study, we cloned and sequenced the chicken SNTB1, SDC2 and GEM genes, and identified sequence polymorphisms between parents of the resource family. The polymorphisms were genotyped to place these genes on the chicken linkage map. The AM gene of chromosome 2q was mapped 130 cM from the distal end, and closely linked to calbindin 1 (CALB1). SNTB1 and SDC2 genes were mapped 88.5 cM distal and 27.6 cM distal from the AM gene, while the GEM gene was mapped 18.5 cM distal from the AM gene and 9.1 cM proximal from SDC2. Orthologues of SNTB1, SDC2 and GEM were syntenic to human chromosome 8q. SNTB1, SDC2 and GEM did not correspond to the AM gene locus, suggesting it is unlikely they are related to chicken muscular dystrophy. However, this result also suggests that the genes located in the proximal region of the CALB1 gene on human chromosome 8q are possible candidates for this disease.

Chromosomal localization of the chicken and mammalian orthologues of the orphan phosphatase PHOSPHO1 gene

PHOSPHO1 is a recently identified phosphatase expressed at high levels in the chicken growth plate and which may be involved in generating inorganic phosphate for skeletal matrix mineralization. Using a degenerate RT-PCR approach a fragment of human PHOSPHO1 was cloned. This enabled the identification of the human orthologue on HSA17q21, and the mouse orthologue on a region of MMU11 that exhibits conservation of synteny with HSA17q21. Chicken PHOSPHO1 was mapped by SSCP analysis to position 44 cM on GGA27, adjacent to the HOXB@ (44 cM) and COL1A1 (36 cM) loci. Comparison of genes on GGA27 with their orthologues on the preliminary draft of the human genome identifies regions of conserved synteny equivalent to 25 Mb on HSA17q21.2-23.3 and approximately 20 Mb on GGA27 in which the gene order appears to be conserved. Mapping of the PHOSPHO1 genes to regions of HSA17q21.3, MMU11 and GGA27 that exhibit conservation of synteny provides strong evidence that they are orthologous.

Mapping of quantitative trait loci for body weight at three, six, and nine weeks of age in a broiler layer cross

An F2 chicken population was established from a cross of a broiler sire-line and an egg laying (White Leghorn) line. There were two males and two females from both lines in the base population. The F1 progeny consisted of 8 males and 32 females. Over 500 F2 offspring from five hatches were reared to slaughter at a live weight of 2 kg at 9 wk of age. Body weights at 3, 6, and 9 wk were recorded. The DNA was extracted from blood samples, and genotypes for 101 microsatellite markers were determined. Data of 466 individuals from 30 families were available for analysis. Interval mapping QTL analyses were carried out. The QTL significant at the genome wide level that affected body weight at two ages were identified on chromosomes 1, 2, 4, 7, and 8 and a QTL on Chromosome 13 influenced body weight at all three ages. Genetic effects were generally additive, and the broiler allele increased body weight in all cases. The effects for significant individual QTL accounted for between 0.2 and 1.0 phenotypic standard deviations and the sum of the additive effects accounted for approximately 0.75 of the line difference in body weight at 6 wk of age. The largest single additive effect was on chromosome 4, and the effect of substituting one copy of the gene was an increase in weight of 249 g. Interactions of the QTL with sex or family were unimportant. There was no evidence for imprinting or of two or more QTL at the same location for any of the traits.

Quantitative trait loci affecting fatness in the chicken

An F2 chicken population of 442 individuals from 30 families, obtained by crossing a broiler line with a layer line, was used for detecting and mapping Quantitative Trait Loci (QTL) affecting abdominal fat weight, skin fat weight and fat distribution. Within-family regression analyses using 102 microsatellite markers in 27 linkage groups were carried out with genome-wide significance thresholds. The QTL for abdominal fat weight were found on chromosomes 3, 7, 15 and 28; abdominal fat weight adjusted for carcass weight on chromosomes 1, 5, 7 and 28; skin and subcutaneous fat on chromosomes 3, 7 and 13; skin fat weight adjusted for carcass weight on chromosomes 3 and 28; and skin fat weight adjusted for abdominal fat weight on chromosomes 5, 7 and 15. Interactions of the QTL with sex or family were unimportant and, for each trait, there was no evidence for imprinting or of multiple QTL on any chromosome. Significant dominance effects were obtained for all but one of the significant locations for QTL affecting the weight of abdominal fat, none for skin fat and one of the three QTL affecting fat distribution. The magnitude of each QTL ranged from 3.0 to 5.2% of the residual phenotypic variation or 0.2-0.8 phenotypic standard deviations. The largest additive QTL (on chromosome 7) accounted for more than 20% of the mean weight of abdominal fat. Significant positive and negative QTL were identified from both lines.

FHF-2 in the turkey (Meleagris gallopavo)

A cDNA clone homologous to the fibroblast growth factor homologous factor (FHF-2) was isolated and sequenced from the turkey (Meleagris gallopavo). The DNA sequence of the turkey was almost identical to that of the chicken (99% similarity) differing at only 8 of 770 nucleotides in the coding region resulting in a single amino acid difference between these poultry species. The 3'UTR of the turkey FHF-2 gene was 445 nucleotides in length and included an imperfect CT microsatellite (ms) repeat. The sequence of the 3'UTR was amplified from genomic DNA of the chicken and found to be highly conserved differing at only three nucleotides when compared to the turkey. Length of the CT repeat was indifferent in a sample of 52 turkeys (monomorphic) however, the number of CT repeats was greater in the turkey than in the chicken. No inter-individual polymorphism was detected in multiple sequences of the 3'UTR of the FHF-2 gene in the turkey. Based on comparison of the turkey and chicken sequences, the mutation rate for coding and associated non-coding (3'UTR) regions of FHF-2 are approximately equal.