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

Expressed sequence tags (ESTs) provide a rapid and reliable method for gene discovery as well as a resource for the large-scale analysis of gene expression of known and unknown genes. Here we describe a normalized cDNA library developed from a 10-day-old White Leghorn chicken whole embryo. The utility of the library was evaluated by partial sequencing of 99 randomly selected insert-containing clones and the analysis of EST-targeted genomic regions for single nucleotide polymorphisms (SNPs) in the East Lansing chicken reference DNA mapping panel. Using stringent match criteria of percent identity of 80 or higher across a length of 50 or more bases, 46 ESTs matched database sequences including previously reported Gallus gallus genes. Thirty-seven of the 50 primer pairs developed from 50 unique ESTs amplified a single fragment. The size of the 37 amplicons ranged from 276 to 693 bp for a total of 17,508 and an average of 473. About 70% of the SNPs detected were either G-->A or C-->T transition. The number of SNPs detected within the amplicons from EST-targeted genomic regions ranged from 0 to 4 for a total of 65 and a frequency of about 1 every 470 bases. About 35% of the amplicons contained only 1 SNP, while 19% had 4 SNPs. Using the SNPs that were informative in the East Lansing reference panel, 17 ESTs were mapped on the East Lansing chicken genetic map. The ESTs described, as well as the nucleotide variants identified within the EST-targeted genomic regions, represent significant resources for genome analysis in the chicken.     

Genetic relationships between Japanese native and commercial breeds using 70 chicken autosomal SNP genotypes by the DigiTag2 assay

Recently, single nucleotide polymorphisms (SNPs) have been used to identify genes or genomic regions responsible for economic traits, including genetic diseases in domestic animals, and to examine genetic diversity of populations. In this study, we genotyped 70 chicken autosomal SNPs using DigiTag2 assay to understand the genetic structure of the Japanese native chicken breeds Satsumadori and Ingie, and the relationship of these breeds with other established breeds, Rhode Island Red (RIR), commercial broiler and layer. Five breeds, each consisting of approximately 20 chickens, were subjected to the assay, revealing the following: Average expected heterozygosities of broiler, Satsumadori, RIR, layer and Ingie were 0.265, 0.254, 0.244, 0.179 and 0.176, respectively. Phylogenetic analysis using the concatenated 70 autosomal SNP genotypes distinguished all chickens and formed clusters of chickens belonging to the respective breeds. In addition, the 2-D scatter plot of the first two principal components was consistent with the phylogenic tree. Taken together with the pairwise F(st) distances, broiler and RIR were closely positioned near each other, while Ingie was positioned far from the other breeds. Structure analysis revealed that the probable number of genetic clusters (K) was six and four with maximum likelihood and ΔK values, respectively. The clustering with maximum likelihood revealed that, in addition to the clustering of the other five breeds, the Satsumadori was subdivided into two genetic clusters. The clustering with ΔK value indicated that the broiler and Rhode Island Red were assigned to the same genetic cluster.     

SNP-specific extraction of haplotype-resolved targeted genomic regions

The availability of genotyping platforms for comprehensive genetic analysis of complex traits has resulted in a plethora of studies reporting the association of specific single-nucleotide polymorphisms (SNPs) with common diseases or drug responses. However, detailed genetic analysis of these associated regions that would correlate particular polymorphisms to phenotypes has lagged. This is primarily due to the lack of technologies that provide additional sequence information about genomic regions surrounding specific SNPs, preferably in haploid form. Enrichment methods for resequencing should have the specificity to provide DNA linked to SNPs of interest with sufficient quality to be used in a cost-effective and high-throughput manner. We describe a simple, automated method of targeting specific sequences of genomic DNA that can directly be used in downstream applications. The method isolates haploid chromosomal regions flanking targeted SNPs by hybridizing and enzymatically elongating oligonucleotides with biotinylated nucleotides based on their selective binding to unique sequence elements that differentiate one allele from any other differing sequence. The targeted genomic region is captured by streptavidin-coated magnetic particles and analyzed by standard genotyping, sequencing or microarray analysis. We applied this technology to determine contiguous molecular haplotypes across a ~150 kb genomic region of the major histocompatibility complex.     

Development and characterization of expressed sequence tags for the turkey (Meleagris gallopavo) genome and comparative sequence analysis with other birds

Twenty-one randomly selected clones from a turkey (Meleagris gallopavo) pituitary complementary DNA (cDNA) library were sequenced to develop expressed sequence tags (ESTs) for this economically important avian species whose genome is among the least understood. Primers specific for the ESTs were used to produce amplicons from the genomic DNA of turkey, chicken (Gallus gallus), guinea fowl (Numidia meleagris), pigeon (Columba domestica), and quail (Corturnix japonica). The amplicons were sequenced and analyzed for sequence variation within- and similarity among-species and with GenBank database sequences. The proportion of shared bases between the turkey sequence and the consensus sequence from each of the other species ranged from 72% to 93% between turkey and pigeon and quail and between turkey and chicken, respectively. The total number of single nucleotide polymorphisms (SNPs) observed ranged from 3 in quail to 18 in chicken out of 4898 and 5265 bases analyzed, respectively. The most frequent nucleotide variation observed was a C-->T transition. Linkage analysis of one such SNP in the backcross progeny of the East Lansing reference DNA panel, localized TUS0005, the chicken sequence derived from primers specific for turkey TUT2E EST, to chromosome 4. The ESTs reported, as well as the SNPs may provide a useful resource for ongoing efforts to develop high utility genome maps for the turkey and chicken. The primers described can also be used as a tool in future investigations directed at further understanding the biology of the guinea fowl, pigeon and quail and their relatedness to the turkey.     

Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo)

ABSTRACT: BACKGROUND: The turkey (Meleagris gallopavo) is an important agricultural species and the second largest contributor to the world's poultry meat production. Genetic improvement is attributed largely to selective breeding programs that rely on highly heritable phenotypic traits, such as body size and breast muscle development. Commercial breeding with small effective population sizes and epistasis can result in loss of genetic diversity, which in turn can lead to reduced individual fitness and reduced response to selection. The presence of genomic diversity in domestic livestock species therefore, is of great importance and a prerequisite for rapid and accurate genetic improvement of selected breeds in various environments, as well as to facilitate rapid adaptation to potential changes in breeding goals. Genomic selection requires a large number of genetic markers such as e.g. single nucleotide polymorphisms (SNPs) the most abundant source of genetic variation within the genome. RESULTS: Alignment of next generation sequencing data of 32 individual turkeys from different populations was used for the discovery of 5.49 million SNPs, which subsequently were used for the analysis of genetic diversity among the different populations. All of the commercial lines branched from a single node relative to the heritage varieties and the South Mexican turkey population. Heterozygosity of all individuals from the different turkey populations ranged from 0.17-2.73 SNPs/Kb, while heterozygosity of populations ranged from 0.73-1.64 SNPs/Kb. The average frequency of heterozygous SNPs in individual turkeys was 1.07 SNPs/Kb. Five genomic regions with very low nucleotide variation were identified in domestic turkeys that showed state of fixation towards alleles different than wild alleles. CONCLUSION: The turkey genome is much less diverse with a relatively low frequency of heterozygous SNPs as compared to other livestock species like chicken and pig. The whole genome SNP discovery study in turkey resulted in the detection of 5.49 million putative SNPs compared to the reference genome. All commercial lines appear to share a common origin. Presence of different alleles/haplotypes in the SM population highlights that specific haplotypes have been selected in the modern domesticated turkey.     

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.     

Cloning and sequence comparison of the mouse, human, and chicken engrailed genes reveal potential functional domains and regulatory regions.

We have isolated and characterized genomic DNA clones for the human and chicken homologues of the mouse En-1 and En-2 genes and determined the genomic structure and predicted protein sequences of both En genes in all three species. Comparison of these vertebrate En sequences with the Xenopus En-2 [Hemmati-Brivanlou et al., 1991) and invertebrate engrailed-like genes showed that the two previously identified highly conserved regions within the En protein ]reviewed in Joyner and Hanks, 1991] can be divided into five distinct subregions, designated EH1 to EH5. Sequences 5' and 3' to the predicted coding regions of the vertebrate En genes were also analyzed in an attempt to identify cis-acting DNA sequences important for the regulation of En gene expression. Considerable sequence similarity was found between the mouse and human homologues both within the putative 5' and 3' untranslated as well as 5' flanking regions. Between the mouse and Xenopus En-2 genes, shorter stretches of sequence similarity were found within the 3' untranslated region. The 5' untranslated regions of the mouse, chicken and Xenopus En-2 genes, however, showed no similarly conserved stretches. In a preliminary analysis of the expression pattern of the human En genes, En-2 protein and RNA were detected in the embryonic and adult cerebellum respectively and not in other tissues tested. These patterns are analogous to those seen in other vertebrates. Taken together these results further strengthen the suggestion that En gene function and regulation has been conserved throughout vertebrate evolution and, along with the five highly conserved regions within the En protein, raise an interesting question about the presence of conserved genetic pathways.     

Isolation and characterization of six different chicken actin genes.

Genes representing six different actin isoforms were isolated from a chicken genomic library. Cloned actin cDNAs as well as tissue-specific mRNAs enriched in different actin species were used as hybridization probes to group individual actin genomic clones by their relative thermal stability. Restriction maps showed that these actin genes were derived from separate and nonoverlapping regions of genomic DNA. Of the six isolated genes, five included sequences from both the 5' and 3' ends of the actin-coding area. Amino acid sequence analysis from both the NH2- and COOH-terminal regions provided for the unequivocal identification of these genes. The striated isoforms were represented by the isolated alpha-skeletal, alpha-cardiac, and alpha-smooth muscle actin genes. The nonmuscle isoforms included the beta-cytoplasmic actin gene and an actin gene fragment which lacked the 5' coding and flanking sequence; presumably, this region of DNA was removed from this gene during construction of the genomic library. Unexpectedly, a third nonmuscle chicken actin gene was found which resembled the amphibian type 5 actin isoform (J. Vandekerckhove, W. W. Franke, and K. Weber, J. Mol. Biol., 152:413-426). This nonmuscle actin type has not been previously detected in warm-blooded vertebrates. We showed that interspersed, repeated DNA sequences closely flanked the alpha-skeletal, alpha-cardiac, beta-, and type 5-like actin genes. The repeated DNA sequences which surround the alpha-skeletal actin-coding regions were not related to repetitious DNA located on the other actin genes. Analysis of genomic DNA blots showed that the chicken actin multigene family was represented by 8 to 10 separate coding loci. The six isolated actin genes corresponded to 7 of 11 genomic EcoRI fragments. Only the alpha-smooth muscle actin gene was shown to be split by an EcoRI site. Thus, in the chicken genome each actin isoform appeared to be encoded by a single gene.     

Comparative DNA sequence analysis of genetic variation in the African grey parrot, Psittacus erythacus

Comparative genome analysis promises to provide an insight into avian species that have been very little studied. To test the feasibility of this approach, we investigated the use of heterologous primers to generate single nucleotide polymorphisms (SNP) in the African grey parrot, Psittacus erythacus, using primers specific for chicken and turkey DNA fragments. Three of the primers were specific for three expressed sequence tagged sites in the turkey and the fourth for a chicken proteoglycan core protein-like DNA sequence. A total of about 2200 bp of the parrot genome was evaluated for DNA sequence variation. Seven SNPs were identified and confirmed by Mendelian segregation. The frequency distribution of the most common nucleotide at each SNP locus in an unrelated group of parrots ranged from 0.84 to 0.97. The percent similarity of each parrot sequence to the reference sequence was inconsistent and ranged from zero to 100%. The primers as well as the nucleotide variants described represent valuable resources for genetic analysis in the parrot. This revised version was published online in July 2006 with corrections to the Cover Date.     

Genomic DNA sequence comparison between two inbred soybean cyst nematode biotypes facilitated by massively parallel 454 micro-bead sequencing

Heterodera glycines, the soybean cyst nematode (SCN), is a damaging agricultural pest that could be effectively managed if critical phenotypes, such as virulence and host range could be understood. While SCN is amenable to genetic analysis, lack of DNA sequence data prevents the use of such methods to study this pathogen. Fortunately, new methods of DNA sequencing that produced large amounts of data and permit whole genome comparative analyses have become available. In this study, 400 million bases of genomic DNA sequence were collected from two inbred biotypes of SCN using 454 micro-bead DNA sequencing. Comparisons to a BAC, sequenced by Sanger sequencing, showed that the micro-bead sequences could identify low and high copy number regions within the BAC. Potential single nucleotide polymorphisms (SNPs) between the two SCN biotypes were identified by comparing the two sets of sequences. Selected resequencing revealed that up to 84% of the SNPs were correct. We conclude that the quality of the micro-bead sequence data was sufficient for de novo SNP identification and should be applicable to organisms with similar genome sizes and complexities. The SNPs identified will be an important starting point in associating phenotypes with specific regions of the SCN genome.     

Isolation of chicken cellular DNA sequences with homology to the region of viral oncogenes that encodes the tyrosine kinase domain.

A library of chicken genomic DNA was screened for sequences that could hybridize to a cloned DNA fragment containing the transforming gene (v-fps) of Fujinami sarcoma virus. In addition to c-fps, two unique chicken cellular DNA sequences were isolated that hybridized weakly to v-fps. These sequences hybridized with many other viral oncogenes encoding tyrosine kinases. Sequence analysis of the region where homology was detected revealed a region that is highly conserved among the tyrosine kinases both at the nucleotide and amino acid levels. Although we were unable to detect expression of either chicken cellular DNA sequence in a variety of avian tissues, the data suggest the existence of additional members of the tyrosine kinase gene family. Screening genomic libraries for sequences that hybridize weakly to functional regions of other genes may prove useful for the isolation and characterization of additional members of other gene families.     

Amplification of sequence tagged sites in five avian species using heterologous oligonucleotides

Short of a complete genomic DNA sequence, sequence tagged sites (STSs) have emerged as major genomic reagents for the genetic analysis of little-studied ecologically and agriculturally important organisms. Here, we report STS developed for the turkey (Meleagris gallopavo), guinea fowl (Numidea meleagris), Japanese quail (Coturnix coturnix) and pigeon using primers specific for reference DNA sequences of two chicken (Gallus gallus) genes, aggrecan (agc1) and type X collagen (col10). Additional STSs were also developed for turkey, quail and chicken using primers specific for the human apobec-1 gene. The total length of the STSs developed was 5990, 2522, 4127, 1539 and 6600 bp for the turkey, guinea fowl, Japanese quail, pigeon and chicken, respectively. Based on splice site consensus GT and AG sequences, four of the seven agc1-based chicken STS appear to contain introns. The human gene-based STSs showed no significant sequence identity with the reference GenBank sequences. Maximum likelihood, maximum parsimony and neighbour-joining analysis of an agc1-based STS that was common to all five species showed phylogenetic relationships consistent with those previously defined using mitochondria DNA sequences and nuclear gene restriction maps. Additionally, several putative single nucleotide polymorphisms (SNPs) were detected within the STSs, including eight in the turkey, two in the quail, and two in the chicken when multiple sequences were evaluated from each species. This report describes new STSs that are resources for genetic and physical mapping and genome analysis within and among avian species. These resources should further aid in our understanding of the biology of agriculturally important but little-studied guinea fowl and turkey. This revised version was published online in July 2006 with corrections to the Cover Date.     

Genomic analysis of genetic markers associated with inherited cardiomyopathy (round heart disease) in the turkey (Meleagris gallopavo)

In turkeys, spontaneous cardiomyopathy or round heart (RH) disease is characterised by dilated ventricles and cardiac muscle hypertrophy. Although the aetiology of RH is still unknown, the disease can have a significant economic impact on turkey producers. In an initial attempt to identify genomic regions associated with RH, we utilised the chicken genome sequence to target short DNA sequences (sequence-characterised amplified regions, SCARs) identified in previous studies that had significant differences in frequency distribution between RH+ and RH- turkeys. SCARs were comparatively aligned with the chicken whole-genome sequence to identify flanking regions for primer design. Primers from 32 alignments were tested and target sequences were successfully amplified for 30 loci (94%). Comparative re-sequencing identified putative SNPs in 20 of the 30 loci (67%). Genetically informative SNPs at 16 loci were genotyped in the UMN/NTBF turkey mapping population. As a result of this study, 34 markers were placed on the turkey/chicken comparative map and 15 markers were added to the turkey genetic linkage map. The position of these markers relative to cardiac-related genes is presented. In addition, analysis of genotypes at 109 microsatellite loci presumed to flank the SCAR sequences in the turkey genome identified four significant associations with RH.     

Assessing genetic diversity in indigenous Veneto chicken breeds using AFLP markers

Genetic variation in four indigenous chicken breeds from the Veneto region of Italy was assessed using amplified fragment length polymorphism (AFLP) markers. A total of 99 individuals were analysed using three AFLP primer combinations that produced 70 polymorphisms. Four indigenous Veneto chicken breeds (Ermellinata, Padovana, Pépoi and Robusta) and a reference broiler line were included in the analysis. Breed-specific markers were identified in each breed. The expected heterozygosity did not differ significantly among the indigenous Veneto chicken breeds and the broiler line. The coefficient of gene variation (Gst) value across loci indicated that almost half of the total variability was observed among breeds. Nei's standard genetic distance between pairs of breeds showed that the distance between the broiler line and the Pépoi breed was greater than the distances between the broiler line and the other three chicken breeds. Cluster analysis based on standard genetic distances between breeds indicated that the Padovana and Pépoi breeds were closely related. Factorial analysis based on a binary matrix of the AFLP data showed a clear distinction of all breeds.     

Infinium Assay for Large-scale SNP Genotyping Applications

Genotyping variants in the human genome has proven to be an efficient method to identify genetic associations with phenotypes. The distribution of variants within families or populations can facilitate identification of the genetic factors of disease. Illumina''s panel of genotyping BeadChips allows investigators to genotype thousands or millions of single nucleotide polymorphisms (SNPs) or to analyze other genomic variants, such as copy number, across a large number of DNA samples. These SNPs can be spread throughout the genome or targeted in specific regions in order to maximize potential discovery. The Infinium assay has been optimized to yield high-quality, accurate results quickly. With proper setup, a single technician can process from a few hundred to over a thousand DNA samples per week, depending on the type of array. This assay guides users through every step, starting with genomic DNA and ending with the scanning of the array. Using propriety reagents, samples are amplified, fragmented, precipitated, resuspended, hybridized to the chip, extended by a single base, stained, and scanned on either an iScan or Hi Scan high-resolution optical imaging system. One overnight step is required to amplify the DNA. The DNA is denatured and isothermally amplified by whole-genome amplification; therefore, no PCR is required. Samples are hybridized to the arrays during a second overnight step. By the third day, the samples are ready to be scanned and analyzed. Amplified DNA may be stockpiled in large quantities, allowing bead arrays to be processed every day of the week, thereby maximizing throughput.     

Detection of Genomic Variation by Selection of a 9 Mb DNA Region and High Throughput Sequencing

Detection of the rare polymorphisms and causative mutations of genetic diseases in a targeted genomic area has become a major goal in order to understand genomic and phenotypic variability. We have interrogated repeat-masked regions of 8.9 Mb on human chromosomes 21 (7.8 Mb) and 7 (1.1 Mb) from an individual from the International HapMap Project (NA12872). We have optimized a method of genomic selection for high throughput sequencing. Microarray-based selection and sequencing resulted in 260-fold enrichment, with 41% of reads mapping to the target region. 83% of SNPs in the targeted region had at least 4-fold sequence coverage and 54% at least 15-fold. When assaying HapMap SNPs in NA12872, our sequence genotypes are 91.3% concordant in regions with coverage≥4-fold, and 97.9% concordant in regions with coverage≥15-fold. About 81% of the SNPs recovered with both thresholds are listed in dbSNP. We observed that regions with low sequence coverage occur in close proximity to low-complexity DNA. Validation experiments using Sanger sequencing were performed for 46 SNPs with 15-20 fold coverage, with a confirmation rate of 96%, suggesting that DNA selection provides an accurate and cost-effective method for identifying rare genomic variants.     

Editor's choice: Functional classification of 15 million SNPs detected from diverse chicken populations

Next-generation sequencing has prompted a surge of discovery of millions of genetic variants from vertebrate genomes. Besides applications in genetic association and linkage studies, a fraction of these variants will have functional consequences. This study describes detection and characterization of 15 million SNPs from chicken genome with the goal to predict variants with potential functional implications (pfVars) from both coding and non-coding regions. The study reports: 183K amino acid-altering SNPs of which 48% predicted as evolutionary intolerant, 13K splicing variants, 51K likely to alter RNA secondary structures, 500K within most conserved elements and 3K from non-coding RNAs. Regions of local fixation within commercial broiler and layer lines were investigated as potential selective sweeps using genome-wide SNP data. Relationships with phenotypes, if any, of the pfVars were explored by overlaying the sweep regions with known QTLs. Based on this, the candidate genes and/or causal mutations for a number of important traits are discussed. Although the fixed variants within sweep regions were enriched with non-coding SNPs, some non-synonymous-intolerant mutations reached fixation, suggesting their possible adaptive advantage. The results presented in this study are expected to have important implications for future genomic research to identify candidate causal mutations and in poultry breeding.     

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.