Expressed sequence tags for the chicken genome from a normalized 10-day-old White Leghorn whole embryo cDNA library: 1. DNA sequence characterization and linkage analysis

Association of resistance to avian coccidiosis with single nucleotide polymorphisms in the zyxin gene

Our previous genetic studies demonstrated that resistance to avian coccidiosis is linked with microsatellite markers LEI0071 and LEI0101 on chromosome 1. In this study, the associations between parameters of resistance to coccidiosis and single nucleotide polymorphisms (SNP) in 3 candidate genes located between LEI0071 and LEI0101 [zyxin, CD4, and tumor necrosis factor receptor super family 1A (TNFRSF1A)] were determined. The SNP were genotyped in 24 F(1) generation and 290 F(2) generation animals. No SNP were identified in the TNFRSF1A gene, whereas 10 were located in the zyxin gene and 4 in the CD4 gene. At various times following experimental infection of the F(2) generation with Eimeria maxima, BW, fecal oocyst shedding, and plasma levels of carotenoid, nitrite plus nitrate (NO(2)(-) + NO(3)(-)), and interferon-gamma (IFN-gamma) were measured as parameters of resistance. Single marker and haplotype-based tests were applied to determine the associations between the 14 SNP and the parameters of coccidiosis resistance. None of the CD4 SNP were correlated with disease resistance. However, by single marker association, several of the zyxin SNP were significantly associated with carotenoid or NO(2)(-) + NO(3)(-) concentrations. These were the SNP at nucleotide 149 associated with carotenoid at d 3 postinfection (PI), nucleotide 187 with carotenoid at d 6 and 9 PI, and nucleotide 159 with carotenoid between d 3 and 9 PI. In addition, the zyxin SNP at nucleotide 191 was significantly associated with increased levels of NO(2)(-) + NO(3)(-) at d 3 PI. By haplotype association, the zyxin SNP also were found to be highly associated with NO(2)(-) + NO(3)(-) at d 3 PI and increased IFN-gamma at d 6 PI. These results suggest that zyxin is a candidate gene potentially associated with increased resistance to experimental avian coccidiosis.

Characterization of expressed sequence tags from a gallus gallus pineal gland cDNA library

The pineal gland is the circadian oscillator in the chicken, regulating diverse functions ranging from egg laying to feeding. Here, we describe the isolation and characterization of expressed sequence tags (ESTs) isolated from a chicken pineal gland cDNA library. A total of 192 unique sequences were analysed and submitted to GenBank; 6% of the ESTs matched neither GenBank cDNA sequences nor the newly assembled chicken genomic DNA sequence, three ESTs aligned with sequences designated to be on the Z_random, while one matched a W chromosome sequence and could be useful in cataloguing functionally important genes on this sex chromosome. Additionally, single nucleotide polymorphisms (SNPs) were identified and validated in 10 ESTs that showed 98% or higher sequence similarity to known chicken genes. Here, we have described resources that may be useful in comparative and functional genomic analysis of genes expressed in an important organ, the pineal gland, in a model and agriculturally important organism.

Mitochondrial DNA sequence and haplotype variation analysis in the chicken (Gallus gallus)

Although it is known to be useful for certain genotype:phenotype assignments, our knowledge of the nature and extent of variation in the entire chicken (Gallus gallus) mitochondrial genome (mtGenome) is limited. Here, we used experimental and in silico tools to identify nucleotide variants in the mtGenome, including the coding and non-coding (D-loop) regions. The distribution of the experimentally identified mitochondrial DNA variants in meat- (broilers) and egg-type (White Leghorn) chickens was also assessed. A total of 113 single-nucleotide polymorphisms (SNPs) were identified. The in silico analysis revealed a total of 91 SNPs, with 70 in the coding region and 21 in the non-coding region. Of the 41 experimentally identified SNPs, 27 were in the D-loop. Together, the experimentally identified SNPs in the non-coding region formed 11 haplotypes, whereas the 14 SNPs in the coding region formed 6. Though, 9 of the D-loop region haplotypes were observed only in broilers, 3 of the 6 haplotypes from the coding region occurred at a significantly higher frequency in broilers. To our knowledge, this investigation represents the first whole-mtGenome scan for variation and an evaluation, though limited in sample size, of the haplotype distribution in meat- and egg-type populations, using the SNPs and haplotypes identified.

Expressed transcripts associated with high rates of egg production in chicken ovarian follicles

The purpose of this study was to characterize differentially expressed transcripts associated with varying rates of egg production in Taiwan country chickens. Ovarian follicles were isolated from two strains of chicken which showed low (B) or high (L2) rates of egg production, then processed for RNA extraction and cDNA library construction. Three thousand and eight forty clones were randomly selected from the cDNA library and amplified by PCR, then used in microarray analysis. Differentially expressed transcripts (P<0.05, log(2)> or = 1.75) were sequenced, and aligned using GenBank. This analysis revealed 20 non-redundant sequences which corresponded to known transcripts. Eight transcripts were expressed at a higher level in ovarian tissue prepared from chicken strain B, and 12 transcripts were expressed at a higher level in L2 birds. These differential patterns of expression were confirmed by semi-quantitative RT-PCR. We show that transcripts of cyclin B2 (cycB2), ferritin heavy polypeptide 1 (FTH1), Gag-Pol polyprotein, thymosin beta4 (TB4) and elongation factor 1 alpha1 (EEF1A1) were enriched in B strain ovarian follicles. In contrast, thioredoxin (TXN), acetyl-CoA dehydrogenase long chain (ACADL), inhibitor of growth family member 4 (ING4) and annexin II (ANXA2) were expressed in at higher levels in the L2 strain. We suggest that our approach may lead to the isolation of effective molecular markers that can be used in selection programs in Taiwan country chickens.

Strategies to Assess Structural Variation in the Chicken Genome and its Associations with Biodiversity and Biological Performance

A primary goal in the assessment of structural variation in the avian genome is to understand the relationship of this variation with biodiversity and with biological performance. To develop such knowledge, certain essential tools are needed. One set of tools includes the laboratory techniques used to assess molecular genetic variation. The current time is a transitional one for this field, in that the recently sequenced chicken genome will add significantly to the portfolio of existing methods used to identify molecular markers. To most efficiently discover marker-trait associations, the experimental mapping populations must be appropriately designed and the relevant statistical analyses applied. This paper reviews methods for assessment of molecular markers in poultry and their use in the characterization of avian biodiversity and in studies to identify marker associations with biological traits, including important considerations of population structure and statistical analysis.

Reduced variation on the chicken Z chromosome

Understanding the population genetic factors that shape genome variability is pivotal to the design and interpretation of studies using large-scale polymorphism data. We analyzed patterns of polymorphism and divergence at Z-linked and autosomal loci in the domestic chicken (Gallus gallus) to study the influence of mutation, effective population size, selection, and demography on levels of genetic diversity. A total of 14 autosomal introns (8316 bp) and 13 Z-linked introns (6856 bp) were sequenced in 50 chicken chromosomes from 10 highly divergent breeds. Genetic variation was significantly lower at Z-linked than at autosomal loci, with one segregating site every 39 bp at autosomal loci (theta(W) = 5.8 +/- 0.8 x 10(-3)) and one every 156 bp on the Z chromosome (theta(W) = 1.4 +/- 0.4 x 10(-3)). This difference may in part be due to a low male effective population size arising from skewed reproductive success among males, evident both in the wild ancestor-the red jungle fowl-and in poultry breeding. However, this effect cannot entirely explain the observed three- to fourfold reduction in Z chromosome diversity. Selection, in particular selective sweeps, may therefore have had an impact on reducing variation on the Z chromosome, a hypothesis supported by the observation of heterogeneity in diversity levels among loci on the Z chromosome and the lower recombination rate on Z than on autosomes. Selection on sex-linked genes may be particularly important in organisms with female heterogamety since the heritability of sex-linked sexually antagonistic alleles advantageous to males is improved when fathers pass a Z chromosome to their sons.

Detection of sequence polymorphisms in red junglefowl and White Leghorn ESTs

Over 16,000 high quality expressed sequence tags (ESTs) from red junglefowl (RJ) and White Leghorn (WL) brain and testis cDNA libraries were generated. Here, we have used this resource for detection of single nucleotide polymorphisms (SNPs), and also completed full-length sequencing of 46 pairs of clones, representing the same gene from both the RJ and WL libraries. From the main set of ESTs, which were assembled using Phrap, 746 putative SNPs were identified, of which 76% were transitions and 24% were transversions. A subset of SNPs was evaluated by sequence analysis of five RJ and five WL birds. Nine of 12 SNPs were verified in this limited sample, suggesting that a majority of the putative polymorphisms documented in this study represent real SNPs. During full-length sequencing of the 46 RJ/WL clones 100 SNPs were identified, which translated to a frequency of 1.90 SNPs/1000 bp. The number of transitions and transversions were 77% and 23%, respectively, and the proportion of non-synonymous vs. synonymous SNPs was 20% and 80%, respectively. Four large insertions/deletions were identified between the RJ and WL full-length sequences, and they appear to represent different splice variants.

A double-screening method to identify reliable candidate non-synonymous SNPs from chicken EST data

Discovery of non-synonymous single nucleotide polymorphisms (nsSNP), which cause amino acid substitutions, is important because they are more likely to alter protein function than synonymous SNPs (sSNP) or those SNPs that do not result in amino acid changes. By changing the coding sequences, nsSNP may play a role in heritable differences between individual organisms. In the chicken and many other vertebrates, the main obstacle for identifying nsSNP is that there is insufficient protein and mRNA sequence information for self-species referencing and thus, determination of the correct reading frame for expressed sequence tags (ESTs) is difficult. Therefore, in order to estimate the correct reading frame at nsSNP in chicken ESTs, a double-screening approach was designed using self- or cross-species protein referencing, in addition to the ESTScan coding region estimation programme. Starting with 23 427 chicken ESTs, 1210 potential SNPs were discovered using a phred/phrap/polyphred/consed pipeline process and among these, 108 candidate nsSNP were identified with the double screening method. A searchable SNP database (chicksnps) for the candidate chicken SNPs, including both nsSNPs and sSNPs is available at http://chicksnps.afs.udel.edu. The chicken SNP data described in this paper have been submitted to the data base SNP under National Center for Biotechnology Information assay ID ss4387050-ss4388259.

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.

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.

Expressed sequence tags for the chicken genome from a normalized 10-day-old white leghorn whole-embryo cDNA library. 3. DNA sequence analysis of genetic variation in commercial chicken populations

Single nucleotide polymorphisms (SNPs) have emerged as a major class of DNA markers with the advantage of permitting the development of high-density genetic maps adequate for quantitative trait loci (QTL) identification by linkage-disequilibrium analysis. Here we describe results of a relatively high-depth survey of chicken broiler and layer populations for SNPs in targeted genomic regions of chicken expressed sequence tag (EST) sites. The sequences scanned, representing the composite sequence of 12 amplified fragments for a total of 6489 bp, were randomly distributed, occurring on six different chromosomes or linkage groups in the chicken genome. Although one of the genomic DNA sequences did not match the reference cDNA sequence, another contained an intron that separated two putative exons. The number of SNPs observed within each of the 12 EST-targeted genomic regions ranged from 0 to 10 for a total of 44 and a frequency of 0.7%. About 70% of the polymorphisms were shared between layer and broiler populations. The average heterozygosity within the populations ranged from 0.15 to 0.48, with the layer populations showing the higher heterozygosity. SNPs and oligonucleotides described will provide a resource for genetic analysis in commercial chicken populations. The data appear to indicate that the relative frequency of SNPs in the targeted regions scanned is higher than the frequency reported for any of the other regions scanned to date in other eukaryotic genomes. Additionally, the results suggest that the use of DNA pools may offer an efficient approach to SNP detection in chickens, as has been shown in other vertebrates.

Expressed sequence tags for the chicken genome from a normalized, ten-day-old white leghorn whole embryo cDNA library. 2. Comparative DNA sequence analysis of guinea fowl, quail, and turkey genomes

Accelerated efforts to develop a high-utility chicken genome map have resulted in the development of resources that may be useful for genetic analysis in other economically important poultry species. Here we describe a total of 26 comparative genomic DNA sequences (CGS) for the guinea fowl, Japanese quail, and domestic turkey developed using 10 primer pairs specific for 10 previously reported, unique, chicken expressed sequence tags (EST). The total length of CGS developed for each of the three species was 4,193, 4,597, and 6,057 bp in quail, turkey, and guinea fowl, respectively. About 70% of the CGS showed significant sequence similarity to reference database sequences, including the reference chicken EST and other avian and nonavian genes. A majority of the between-species comparisons of the CGS from all but two primer pairs were significant and ranged from 81 to 99%. The percentage similarity of the CGS appears to be a function of phylogenetic relatedness and was generally higher for comparisons between the chicken, quail, and turkey and lower between the guinea fowl and chicken, quail, or turkey. Maximum likelihood estimation of the phylogenetic relationships using CGS from two primer pairs also showed a closer relationship, as expected, among chicken, quail, and turkey than between guinea fowl and either chicken, quail, or turkey. Within the guinea fowl, quail, and turkey CGS developed, the total number of single nucleotide polymorphisms detected was 28, 17, and 14, respectively. Together, these resources represent tools that will facilitate genetic analysis of species that have been studied very little and our understanding of their genomes and genome evolution.