Association and in Silico Assignment of Sequences from Turkey BACs

Bacterial artificial chromosomes (BACs) provide an important resource in genetic mapping. An initial set of BACs corresponding to microsatellite markers in the turkey (Meleagris gallopavo) was isolated from the CHORI-260 turkey BAC library. The selected markers were distributed on both macro- and microchromosomes and included a genetically unlinked marker. End sequences were obtained for a subset of the recovered BACs and compared to the chicken whole genome sequence. Close association of the turkey BAC-end sequences and original marker sequences was generally conserved in the chicken genome. Gene content of the turkey BACs is predicted from the comparative sequence alignments.     

Cross-species overgo hybridization and comparative physical mapping within avian genomes

The chicken genome sequence facilitates comparative genomics within other avian species. We performed cross-species hybridizations using overgo probes designed from chicken genomic and zebra finch expressed sequence tags (ESTs) to turkey and zebra finch BAC libraries. As a result, 3772 turkey BACs were assigned to 336 markers or genes, and 1662 zebra finch BACs were assigned to 164 genes. As expected, cross-hybridization was more successful with overgos within coding sequences than within untranslated region, intron or flanking sequences and between chicken and turkey, when compared with chicken-zebra finch or zebra finch-turkey cross-hybridization. These data contribute to the comparative alignment of avian genome maps using a 'one sequence, multiple genomes' strategy.     

An integrated and comparative genetic map of the turkey genome

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

Using the chicken genome sequence in the development and mapping of genetic markers in the turkey (Meleagris gallopavo)

The efficacy of employing the chicken genome sequence in developing genetic markers and in mapping the turkey genome was studied. Eighty previously uncharacterized microsatellite markers were identified for the turkey using BLAST alignment to the chicken genome. The chicken sequence was then used to develop primers for polymerase chain reaction where the turkey sequence was either unavailable or insufficient. A total of 78 primer sets were tested for amplification and polymorphism in the turkey, and informative markers were genetically mapped. Sixty-five (83%) amplified turkey genomic DNA, and 33 (42%) were polymorphic in the University of Minnesota/Nicholas Turkey Breeding Farms mapping families. All but one marker genetically mapped to the position predicted from the chicken genome sequence. These results demonstrate the usefulness of the chicken sequence for the development of genomic resources in other avian species.     

A comparative genetic map of the turkey genome

Genetic markers (microsatellites and SNPs) were used to create and compare maps of the turkey and chicken genomes. A physical map of the chicken genome was built by comparing sequences of turkey markers with the chicken whole-genome sequence by BLAST analysis. A genetic linkage map of the turkey genome (Meleagris gallopavo) was developed by segregation analysis of genetic markers within the University of Minnesota/Nicholas Turkey Breeding Farms (UMN/NTBF) resource population. This linkage map of the turkey genome includes 314 loci arranged into 29 linkage groups. An additional 40 markers are tentatively placed within linkage groups based on two-point LOD scores and 16 markers remain unlinked. Total map distance contained within linkage groups is 2,011 cM with the longest linkage group (47 loci) measuring 413.3 cM. Average marker interval over the 29 linkage groups was 6.4 cM. All but one turkey linkage group could be aligned with the physical map of the chicken genome. The present genetic map of the turkey provides a comparative framework for future genomic studies.     

Assignment of non-informative turkey genetic markers through comparative approaches

Molecular markers such as microsatellites, provide genetic signposts for navigating genomes. In general, genetic markers that are monomorphic or non-informative in mapping populations typically remain unmapped and as such are less likely to be included in future studies. The use of hybrid cell panels and in silico mapping via whole genome sequences allow for positional mapping of non-segregating markers. This study utilizes the INRA ChickRH6 whole-genome radiation hybrid panel and chicken whole-genome shotgun sequence to map microsatellite markers from the turkey (Meleagris gallopavo). Thirty-three of the 41 markers typed on the RH panel had significant linkage to at least one other marker and 83 of 100 sequences returned significant BLAST similarities. Positioning of these markers provides additional sequence tagged sites in the turkey genome and increases the potential use of these markers for future genetic studies.     

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.     

Integration of microsatellite-based genetic maps for the turkey (Meleagris gallopavo).

Integration of turkey genetic maps and their associated markers is essential to increase marker density in support of map-based genetic studies. The objectives of this study were to integrate 2 microsatellite-based turkey genetic maps--the Roslin map and the University of Minnesota (UMN) map--by genotyping markers from the Roslin study on the mapping families of the UMN study. A total of 279 markers was tested, and 240 were subsequently screened for polymorphisms in the UMN/Nicholas Turkey Breeding Farms (NTBF) mapping families. Of the 240 markers, 89 were genetically informative and were used for genotyping the F2 offspring. Significant genetic linkages (log of odds > 3.0) were found for 84 markers from the Roslin study. BLASTn comparison of marker sequences with the draft assembly of the chicken genome found 263 significant matches. The combination of genetic and in silico mapping allowed for the alignment of all linkage groups of the Roslin map with those of the UMN map. With the addition of the markers from the Roslin map, 438 markers are now genetically linked in the UMN/NTBF families, and more than 1700 turkey sequences have now been assigned to likely positions in the chicken-genome sequence.     

Soybean genomic survey: BAC-end sequences near RFLP and SSR markers.

We are building a framework physical infrastructure across the soybean genome by using SSR (simple sequence repeat) and RFLP (restriction fragment length polymorphism) markers to identify BACs (bacterial artificial chromosomes) from two soybean BAC libraries. The libraries were prepared from two genotypes, each digested with a different restriction enzyme. The BACs identified by each marker were grouped into contigs. We have obtained BAC- end sequence from BACs within each contig. The sequences were analyzed by the University of Minnesota Center for Computational Genomics and Bioinformatics using BLAST algorithms to search nucleotide and protein databases. The SSR-identified BACs had a higher percentage of significant BLAST hits than did the RFLP-identified BACs. This difference was due to a higher percentage of hits to repetitive-type sequences for the SSR-identified BACs that was offset in part, however, by a somewhat larger proportion of RFLP-identified significant hits with similarity to experimentally defined genes and soybean ESTs (expressed sequence tags). These genes represented a wide range of metabolic functions. In these analyses, only repetitive sequences from SSR-identified contigs appeared to be clustered. The BAC-end sequences also allowed us to identify microsynteny between soybean and the model plants Arabidopsis thaliana and Medicago truncatula. This map-based approach to genome sampling provides a means of assaying soybean genome structure and organization.     

In Silco Mapping of ESTs from the Turkey (Meleagris Gallopavo)

Sequence similarity was used to predict the position of expressed sequence tags (ESTs) in the genome of the turkey (Meleagris gallopavo). Turkey EST sequences were compared with the draft assembly of the chicken whole-genome sequence and the chicken EST database by BLASTN. Among the 877 ESTs examined, 788 had significant matches in the chicken genome sequence. Position of orthologous sequences in the chicken genome and the predicted position of the EST loci in the turkey genome are presented. Genetic assignments suggest a high level of accuracy for the COMPASS predictions.     

In silco mapping of ESTs from the turkey (Meleagris gallopavo)

Sequence similarity was used to predict the position of expressed sequence tags (ESTs) in the genome of the turkey (Meleagris gallopavo). Turkey EST sequences were compared with the draft assembly of the chicken whole-genome sequence and the chicken EST database by BLASTN. Among the 877 ESTs examined, 788 had significant matches in the chicken genome sequence. Position of orthologous sequences in the chicken genome and the predicted position of the EST loci in the turkey genome are presented Genetic assignments suggest a high level of accuracy for the COMPASS predictions.     

Survey of a cDNA library from the turkey (Meleagris gallopavo).

Genome characterization and analysis is an imperative step in identifying and selectively breeding for improved traits of agriculturally important species. Expressed sequence tags (ESTs) represent a transcribed portion of the genome and are an effective way to identify genes within a species. Downstream applications of EST projects include DNA microarray construction and interspecies comparisons. In this study, 694 ESTs were sequenced and analyzed from a library derived from a 24-day-old turkey embryo. The 437 unique sequences identified were divided into 76 assembled contigs and 361 singletons. The majority of significant comparative matches occurred between the turkey sequences and sequences reported from the chicken. Whole genome sequence from the chicken was used to identify potential exon-intron boundaries for selected turkey clones and intron-amplifying primers were developed for sequence analysis and single nucleotide polymorphism (SNP) discovery. Identified SNPs were genotyped for linkage analysis on two turkey reference populations. This study significantly increases the number of EST sequences available for the turkey.     

Integration of animal linkage and BAC contig maps using overgo hybridization

The alignment of genome linkage maps, defined primarily by segregation of sequence-tagged site (STS) markers, with BAC contig physical maps and full genome sequences requires high throughput mechanisms to identify BAC clones that contain specific STS. A powerful technique for this purpose is multi-dimensional hybridization of "overgo" probes. The probes are chosen from available STS sequence data by selecting unique probe sequences that have a common melting temperature. We have hybridized sets of 216 overgo probes in subset pools of 36 overgos at a time to filter-spotted chicken BAC clone arrays. A four-dimensional pooling strategy, including one degree of redundancy, has been employed. This requires 24 hybridizations to completely assign BACs for all 216 probes. Results to date are consistent with about a 10% failure rate in overgo probe design and a 15-20% false negative detection rate within a group of 216 markers. Three complete rounds of overgo hybridization, each to sets of about 39,000 BACs (either BAMHI or ECORI partial digest inserts) generated a total of 1853 BAC alignments for 517 mapped chicken genome STS markers. These data are publicly available, and they have been used in the assembly of a first generation BAC contig map of the chicken genome.     

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.     

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 F1 sire 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.     

Comparative analysis of microsatellite loci in chicken and turkey.

Cross-species amplification of 520 chicken microsatellite markers was tested by polymerase chain reaction with genomic DNA of the turkey (Meleagris gallopavo). Each primer pair was tested at six different combinations of annealing temperature and MgCl2 concentration. A total of 280 (54%) of the primer pairs produced amplification products. The majority of these products were similar, if not identical in size to those expected based on the fragment sizes of the corresponding chicken loci. Structure of the dinucleotide repeat and flanking sequences was examined for 13 turkey fragments (amplified with chicken primers) and 5 chicken fragments (amplified with turkey primers). Sequence analysis found a wide array of mutations between species in addition to differences in repeat length. To estimate the usefulness of the amplified loci for genetic mapping in the turkey, allelic polymorphism was determined for 57 of the 280 amplified loci. A total of 20 of 57 markers (35%) were polymorphic with an average of 1.4 alleles per locus. The results of this study suggest that approximately 20% of the chicken microsatellite markers will be useful for mapping the turkey genome.     

Construction, characterization, and preliminary BAC-end sequencing analysis of a bacterial artificial chromosome library of white clover (Trifolium repens L.).

White clover (Trifolium repens L.) is a forage legume widely used in combination with grass in pastures because of its ability to fix nitrogen. We have constructed a bacterial artificial chromosome (BAC) library of an advanced breeding line of white clover. The library contains 37 248 clones with an average insert size of approximately 85 kb, representing an approximate 3-fold coverage of the white clover genome based on an estimated genome size of 960 Mb. The BAC library was pooled and screened by polymerase chain reaction (PCR) amplification using both white clover microsatellites and PCR-based markers derived from Medicago truncatula, resulting in an average of 6 hits per marker; this supports the estimated 3-fold genome coverage in this allotetraploid species. PCR-based screening of 766 clones with a multiplex set of chloroplast primers showed that only 0.5% of BAC clones contained chloroplast-derived inserts. The library was further evaluated by sequencing both ends of 724 of the clover BACs. These were analysed with respect to their sequence content and their homology to the contents of a range of plant gene, expressed sequence tag, and repeat element databases. Forty-three microsatellites were discovered in the BAC-end sequences (BESs) and investigated as potential genetic markers in white clover. The BESs were also compared with the partially sequenced genome of the model legume M. truncatula with the specific intention of identifying putative comparative-tile BACs, which represent potential regions of microsynteny between the 2 species; 14 such BACs were discovered. The results suggest that a large-scale BAC-end sequencing strategy has the potential to anchor a significant proportion of the genome of white clover onto the gene-space sequence of M. truncatula.     

Insertion events of CR1 retrotransposable elements elucidate the phylogenetic branching order in galliform birds

Using standard phylogenetic methods, it can be hard to resolve the order in which speciation events took place when new lineages evolved in the distant past and within a short time frame. As an example, phylogenies of galliform birds (including well-known species such as chicken, turkey, and quail) usually show low bootstrap support values at short internal branches, reflecting the rapid diversification of these birds in the Eocene. However, given the key role of chicken and related poultry species in agricultural, evolutionary, general biological and disease studies, it is important to know their internal relationships. Recently, insertion patterns of transposable elements such as long and short interspersed nuclear element markers have proved powerful in revealing branching orders of difficult phylogenies. Here we decipher the order of speciation events in a group of 27 galliform species based on insertion events of chicken repeat 1 (CR1) transposable elements. Forty-four CR1 marker loci were identified from the draft sequence of the chicken genome, and from turkey BAC clone sequence, and the presence or absence of markers across species was investigated via electrophoretic size separation of amplification products and subsequent confirmation by DNA sequencing. Thirty markers proved possible to type with electrophoresis of which 20 were phylogenetically informative. The distribution of these repeat elements supported a single homoplasy-free cladogram, which confirmed that megapodes, cracids, New World quail, and guinea fowl form outgroups to Phasianidae and that quails, pheasants, and partridges are each polyphyletic groups. Importantly, we show that chicken is an outgroup to turkey and quail, an observation which does not have significant support from previous DNA sequence- and DNA-DNA hybridization-based trees and has important implications for evolutionary studies based on sequence or karyotype data from galliforms. We discuss the potential and limitations of using a genome-based retrotransposon approach in resolving problematic phylogenies among birds.     

Reassignment of chicken W chromosome sequences to the Z chromosome by fluorescence in situ hybridization (FISH)

There is much interest in the gene content of the small heterochromatic W chromosome of the chicken, on the supposition that it may contain sex-determining genes. A considerable region in the chicken genome has been assigned to the W chromosome on the basis of its repetitive sequences. Using fluorescent in situ hybridization (FISH) we localized five Bacterial Artificial Chromosomes (BACs) onto female chicken metaphase spreads. We physically mapped these BACs to the Z chromosome. The chicken genome database, however, assigned all five BACs to the W chromosome. Our results demonstrate that the 17 genes on these BACs are Z-specific, and points to the inadequacy of assigning regions of the genome based exclusively on repetitive sequences.     

A genome-wide set of 106 microsatellite markers for the blue tit (Cyanistes caeruleus)

We have characterized a set of 106 microsatellite markers in 26-127 individual blue tits (Cyanistes caeruleus), and assigned their location on the zebra finch (Taeniopygia guttata) and on the chicken (Gallus gallus) genome on the basis of sequence homology. Thirty-one markers are newly designed from zebra finch EST (expressed sequence tags) sequences, 22 markers were developed by others from EST sequences using different methods and the remaining 53 loci were previously designed or modified passerine markers. The 106 microsatellite markers are distributed over 26 and 24 chromosomes in the zebra finch and in the chicken genome respectively and the number of alleles varies between 2 and 49. Eight loci deviate significantly from Hardy-Weinberg equilibrium and show a high frequency of null alleles, and three pairs of markers located in the same chromosome appear to be in linkage disequilibrium. With the exception of these few loci, the polymorphic microsatellite markers presented here provide a useful genome-wide resource for population and evolutionary genetic studies of the blue tit, in addition to their potential utility in other passerine birds.