An integrated and comparative genetic map of the turkey genome

Genetic diversity in 10 populations of domestic Turkeys by using microsatellites markers

The domestic turkey is a native breed in danger of extinction due to the introduction of new breeds specializing in meat production and yield. Turkeys have lost some prominence in urban areas, and only certain breeds of turkeys are preserved in rural areas. Wild and domestic turkeys are different; rural or indigenous turkeys, with black plumage, were domesticated from Mexican turkeys and have been reproduced throughout Latin America. Some of them were taken to Europe in the 16th century and later arrived in North America, where they crossed with another wild species, from which the bronze turkey emerged: the ancestor of all commercial turkeys. The objective of the present work was to evaluate the genetic diversity in 10 populations of domestic turkeys worldwide by using breeds from Europe: Spain and Italy; America: Mexico, United States and Brazil; and the Near East: Iran and Egypt. A total of 522 blood samples of both sexes were collected from domestic turkey populations. Thirty-four microsatellites were used to obtain genetic parameters, and genetic diversity was evaluated. All microsatellites used were polymorphic, and a total of 427 alleles were detected across the 34 markers investigated. In this study, a mean number of 13.44 alleles was found. The four most diverse breeds were from the Andalusia, Mexico, United States, and wild populations, which had the highest mean heterozygosity expected (0.619, 0.612, 0.650, and 0.773) and heterozygosity observed (0.422, 0.521, 0.429, and 0.627), respectively. The MNT348 marker deviated from the HWE in all populations. Our study has shown that the populations close to the species origin are more diverse than those resulting from posterior expansions. Mexican birds were the most diverse, followed by the Spanish populations because Spain imported a large number of turkeys coming from America. Such information can be complementary to other genotypic data required to validate the evolutionary relationships among turkey populations.

Conserved MHC gene orthologs genetically map to the turkey MHC- B

The avian MHC-associated gene set includes orthologs to genes found throughout the human major histocompatibility complex (MHC), including some loci of the evolutionarily conserved class III region. In the turkey and other Galliformes, genes linked to the MHC have been identified because they are closely associated with class I or class II genes. This study was designed to evaluate additional class III genes for linkage to the avian MHC to further determine conservation of these loci in birds. BLAST searches were used to locate sequences in the turkey genome with similarity to genes shared between the MHC of Xenopus and humans. Primers were designed to target 25 genes, and putative orthologs were amplified by PCR and sequenced. Sequence polymorphisms were identified for 15 genes in turkey reference mapping families, and 8 genes showed significant genetic linkage to the turkey MHC-B locus. These new genetic markers and linkage relationships broaden our understanding of the composition of the avian MHC and expand the gene content for the turkey MHC-B.

Replicated high-density genetic maps of two great tit populations reveal fine-scale genomic departures from sex-equal recombination rates

Linking variation in quantitative traits to variation in the genome is an important, but challenging task in the study of life-history evolution. Linkage maps provide a valuable tool for the unravelling of such trait-gene associations. Moreover, they give insight into recombination landscapes and between-species karyotype evolution. Here we used genotype data, generated from a 10k single-nucleotide polymorphism (SNP) chip, of over 2000 individuals to produce high-density linkage maps of the great tit (Parus major), a passerine bird that serves as a model species for ecological and evolutionary questions. We created independent maps from two distinct populations: a captive F2-cross from The Netherlands (NL) and a wild population from the United Kingdom (UK). The two maps contained 6554 SNPs in 32 linkage groups, spanning 2010 cM and 1917 cM for the NL and UK populations, respectively, and were similar in size and marker order. Subtle levels of heterochiasmy within and between chromosomes were remarkably consistent between the populations, suggesting that the local departures from sex-equal recombination rates have evolved. This key and surprising result would have been impossible to detect if only one population was mapped. A comparison with zebra finch Taeniopygia guttata, chicken Gallus gallus and the green anole lizard Anolis carolinensis genomes provided further insight into the evolution of avian karyotypes.

Twisted signatures of GC-biased gene conversion embedded in an evolutionary stable karyotype

The genomes of many vertebrates show a characteristic heterogeneous distribution of GC content, the so-called GC isochore structure. The origin of isochores has been explained via the mechanism of GC-biased gene conversion (gBGC). However, although the isochore structure is declining in many mammalian genomes, the heterogeneity in GC content is being reinforced in the avian genome. Despite this discrepancy, which remains unexplained, examinations of individual substitution frequencies in mammals and birds are both consistent with the gBGC model of isochore evolution. On the other hand, a negative correlation between substitution and recombination rate found in the chicken genome is inconsistent with the gBGC model. It should therefore be important to consider along with gBGC other consequences of recombination on the origin and fate of mutations, as well as to account for relationships between recombination rate and other genomic features. We therefore developed an analytical model to describe the substitution patterns found in the chicken genome, and further investigated the relationships between substitution patterns and several genomic features in a rigorous statistical framework. Our analysis indicates that GC content itself, either directly or indirectly via interrelations to other genomic features, has an impact on the substitution pattern. Further, we suggest that this phenomenon is particularly visible in avian genomes due to their unusually low rate of chromosomal evolution. Because of this, interrelations between GC content and other genomic features are being reinforced, and are as such more pronounced in avian genomes as compared with other vertebrate genomes with a less stable karyotype.

A comparative physical map reveals the pattern of chromosomal evolution between the turkey (Meleagris gallopavo) and chicken (Gallus gallus) genomes

Background

A robust bacterial artificial chromosome (BAC)-based physical map is essential for many aspects of genomics research, including an understanding of chromosome evolution, high-resolution genome mapping, marker-assisted breeding, positional cloning of genes, and quantitative trait analysis. To facilitate turkey genetics research and better understand avian genome evolution, a BAC-based integrated physical, genetic, and comparative map was developed for this important agricultural species.

Results

The turkey genome physical map was constructed based on 74,013 BAC fingerprints (11.9 × coverage) from two independent libraries, and it was integrated with the turkey genetic map and chicken genome sequence using over 41,400 BAC assignments identified by 3,499 overgo hybridization probes along with > 43,000 BAC end sequences. The physical-comparative map consists of 74 BAC contigs, with an average contig size of 13.6 Mb. All but four of the turkey chromosomes were spanned on this map by three or fewer contigs, with 14 chromosomes spanned by a single contig and nine chromosomes spanned by two contigs. This map predicts 20 to 27 major rearrangements distinguishing turkey and chicken chromosomes, despite up to 40 million years of separate evolution between the two species. These data elucidate the chromosomal evolutionary pattern within the Phasianidae that led to the modern turkey and chicken karyotypes. The predominant rearrangement mode involves intra-chromosomal inversions, and there is a clear bias for these to result in centromere locations at or near telomeres in turkey chromosomes, in comparison to interstitial centromeres in the orthologous chicken chromosomes.

Conclusion

The BAC-based turkey-chicken comparative map provides novel insights into the evolution of avian genomes, a framework for assembly of turkey whole genome shotgun sequencing data, and tools for enhanced genetic improvement of these important agricultural and model species.

Functional characterization of alpha-class glutathione s-transferases from the Turkey (meleagris gallopavo)

Six Alpha-class glutathione S-transferase (GST) subunits were cloned from domestic turkey livers, which are one of the most susceptible animals known to the carcinogenic mycotoxin aflatoxin B₁. In most animals, GST dysfunction is a risk factor for susceptibility toward AFB₁, and we have shown that turkeys lack GSTs with affinity toward the carcinogenic intermediate exo-aflatoxin B(1)-8-9-epoxide (AFBO). Conversely, mice are resistant to AFB₁ carcinogenesis, due to high constitutive expression of mGSTA3 that has high affinity toward AFBO. When expressed in Escherichia coli, all six tGSTA subunits possessed conjugating activities toward substrates 1-chloro-2,4-dinitrobenzene (CDNB), 1,2-dichloro-4-nitrobenzene (DCNB), ethacrynic acid (ECA), and cumene hydroperoxide (CHP) with tGSTA1.2 appearing most active. Interestingly, tGSTA1.1, which lacks one of the four Alpha-class signature motifs, possessed enzymatic activities toward all substrates. All had comparable activities toward AFBO conjugation, an activity absent in turkey liver cytosols. E. coli-expressed mGSTA3 conjugated AFBO with more than 3-fold greater activity than that of tGSTAs and had higher activity toward GST prototype substrates. Mouse hepatic cytosols had approximately 900-fold higher catalytic activity toward AFBO compared with those from turkey. There was no apparent amino acid profile in tGSTAs that might correspond to specificity toward AFBO, although tGSTA1.2, which had slightly higher AFBO-trapping ability, shared Tyr¹⁰⁸ with mGSTA3, a residue postulated to be critical for AFBO trapping activity in mammalian systems. The observation that recombinant tGSTAs detoxify AFBO, whereas their hepatic forms do not, implies that the hepatic forms of these enzymes are silenced by one or more regulatory mechanisms.

Defining the turkey MHC: identification of expressed class I- and class IIB-like genes independent of the MHC-B

The MHC of the turkey (Meleagris gallopavo) is divided into two genetically unlinked regions; the MHC-B and MHC-Y. Although previous studies found the turkey MHC-B to be highly similar to that of the chicken, little is known of the gene content and extent of the MHC-Y. This study describes two partially overlapping large-insert BAC clones that genetically and physically map to the turkey MHC chromosome (MGA18) but to a region that assorts independently of MHC-B. Within the sequence assembly, 14 genes were predicted including new class I- and class IIB-like loci. Additional unassembled sequences corresponded to multiple copies of the ribosomal RNA repeat unit (18S-5.8S-28S). Thus, this newly identified MHC region appears to represent a physical boundary of the turkey MHC-Y. High-resolution multi-color fluorescence in situ hybridization studies confirm rearrangement of MGA18 relative to the orthologous chicken chromosome (GGA16) in regard to chromosome architecture, but not gene order. The difference in centromere position between the species is indicative of multiple chromosome rearrangements or alternate events such as neocentromere formation/centromere inactivation in the evolution of the MHC chromosome. Comparative sequencing of commercial turkeys (six amplicons totaling 7.6 kb) identified 68 single nucleotide variants defining nine MHC-Y haplotypes. Sequences of the new class I- and class IIB-like genes are most similar to MHC-Y genes in the chicken. All three loci are expressed in the spleen. Differential transcription of the MHC-Y class IIB-like loci was evident as one class IIB-like locus was only expressed in some individuals.

A locus-wide approach to assessing variation in the avian MHC: the B-locus of the wild turkey

Studies of major histocompatibility complex (MHC) diversity in non-model vertebrates typically focus on structure and sequence variation in the antigen-presenting loci: the highly variable and polymorphic class I and class IIB genes. Although these studies provide estimates of the number of genes and alleles/locus, they often overlook variation in functionally related and co-inherited genes important in the immune response. This study utilizes the sequence of the MHC B-locus derived from a commercial turkey to investigate MHC variation in wild birds. Sequences were obtained for nine interspersed MHC amplicons (non-class I/II) from each of 40 birds representing 3 subspecies of wild turkey (Meleagris gallopavo). Analysis of aligned sequences identified 238 single-nucleotide variants approximately one-third of which had minor allele frequencies >0.2 in the sampled birds. PHASE analysis identified 70 prospective MHC haplotypes in the wild turkeys, whereas a combined analysis with commercial birds identified almost 100 haplotypes in the species. Denaturing gradient gel electrophoresis (DGGE) of the class IIB loci was used to test the efficacy of single-nucleotide polymorphism (SNP) haplotyping to capture locus-wide variation. Diversity in SNP haplotypes and haplotype sharing among individuals was directly reflected in the DGGE patterns. Utilization of a reference haplotype to sequence interspersed regions of the MHC has significant advantages over other methods of surveying diversity while identifying high-frequency SNPs for genotyping. SNP haplotyping provides a means to identify both divergent haplotypes and homozygous individuals for assessment of immunological variation in wild and domestic populations.

Comparative physical maps derived from BAC end sequences of tilapia (Oreochromis niloticus)

Background

The Nile tilapia is the second most important fish in aquaculture. It is an excellent laboratory model, and is closely related to the African lake cichlids famous for their rapid rates of speciation. A suite of genomic resources has been developed for this species, including genetic maps and ESTs. Here we analyze BAC end-sequences to develop comparative physical maps, and estimate the number of genome rearrangements, between tilapia and other model fish species.

Results

We obtained sequence from one or both ends of 106,259 tilapia BACs. BLAST analysis against the genome assemblies of stickleback, medaka and pufferfish allowed identification of homologies for approximately 25,000 BACs for each species. We calculate that rearrangement breakpoints between tilapia and these species occur about every 3 Mb across the genome. Analysis of 35,000 clones previously assembled into contigs by restriction fingerprints allowed identification of longer-range syntenies.

Conclusions

Our data suggest that chromosomal evolution in recent teleosts is dominated by alternate loss of gene duplicates, and by intra-chromosomal rearrangements (~one per million years). These physical maps are a useful resource for comparative positional cloning of traits in cichlid fishes. The paired BAC end sequences from these clones will be an important resource for scaffolding forthcoming shotgun sequence assemblies of the tilapia genome.

Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis

A synergistic combination of two next-generation sequencing platforms with a detailed comparative BAC physical contig map provided a cost-effective assembly of the genome sequence of the domestic turkey (Meleagris gallopavo). Heterozygosity of the sequenced source genome allowed discovery of more than 600,000 high quality single nucleotide variants. Despite this heterozygosity, the current genome assembly (∼1.1 Gb) includes 917 Mb of sequence assigned to specific turkey chromosomes. Annotation identified nearly 16,000 genes, with 15,093 recognized as protein coding and 611 as non-coding RNA genes. Comparative analysis of the turkey, chicken, and zebra finch genomes, and comparing avian to mammalian species, supports the characteristic stability of avian genomes and identifies genes unique to the avian lineage. Clear differences are seen in number and variety of genes of the avian immune system where expansions and novel genes are less frequent than examples of gene loss. The turkey genome sequence provides resources to further understand the evolution of vertebrate genomes and genetic variation underlying economically important quantitative traits in poultry. This integrated approach may be a model for providing both gene and chromosome level assemblies of other species with agricultural, ecological, and evolutionary interest.

Defining the turkey MHC: sequence and genes of the B locus

The MHC, the most polymorphic and gene dense region in the vertebrate genome, contains many loci essential to immunity. In mammals, this region spans approximately 4 Mb. Studies of avian species have found the MHC to be greatly reduced in size and gene content with an overall locus organization differing from that of mammals. The chicken MHC has been mapped to two distinct regions (MHC-B and -Y) of a single chromosome. MHC-B haplotypes possess tightly linked genes encoding the classical MHC molecules and few other disease resistance genes. Furthermore, chicken haplotypes possess a dominantly expressed class I and class II B locus that have a significant effect on the progression or regression of pathogenic disease. In this study, we present the MHC-B region of the turkey (Meleagris gallopavo) as a similarly constricted locus, with 34 genes identified within a 0.2-Mb region in near-perfect synteny with that of the chicken MHC-B. Notable differences between the two species are three BG and class II B loci in the turkey compared with one BG and two class II B loci in the chicken MHC-B. The relative size and high level of similarity of the turkey MHC in relation to that of the chicken suggest that similar associations with disease susceptibility and resistance may also be found in turkey.

Genome-wide genetic diversity of 'Nici', the DNA source for the CHORI-260 turkey BAC library and candidate for whole genome sequencing

Vertebrate whole genome sequence assembly can benefit from a priori knowledge of variability in the target genome, with researchers often selecting highly inbred individuals for sequencing. However, for most species highly inbred research lines are lacking, requiring the use of an outbred individual(s). Here we examined the source DNA [Nicholas inbred (Nici)] of the CHORI-260 turkey bacterial artificial chromosome (BAC) library through analysis of microsatellites and BAC sequences. Heterozygosity of Nici was compared with that of individuals from several breeder lines. Seventy-eight microsatellites were screened for polymorphism in a total of 43 birds, identifying an average individual heterozygosity of 0.39, with Nici at 0.35. Additional loci (total of 147) were examined on a subset of individuals to obtain better genome coverage. The mean heterozygosity for this subset was 0.33 with Nici at 0.31. Examination of approximately 200 kb of genome sequence identified SNPs in the order of one per 200 bp in Nici. These data suggest that the heterozygosity of Nici is comparable to other birds of selected breeder lines and that whole genome sequencing would result in an abundant resource of genome-wide polymorphisms.

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.

Characterization of a cardiac complementary deoxyribonucleic acid library from the Turkey (Meleagris gallopavo)

Although the domestic turkey (Meleagris gallopavo) is a valuable agricultural commodity, genetic studies on this species lag behind those of other agricultural species. In this study, we examined expressed sequence tags (EST) from a turkey cardiac cDNA library constructed from 4 birds representing 2 developmental stages. A collection of 3,937 EST sequences were sequenced and analyzed for gene annotation and sequence variation. Clustering of sequences resulted in 353 contigs and 874 singletons (1,227 putative transcripts). All EST sequences were compared by BLASTN to the chicken whole genome sequence and to Ensembl and National Center for Biotechnology Information databases. The majority of significant matches correspond to genes found in the chicken. Sequence polymorphisms were identified in 310 contigs, 64 where the minor allele was observed to be present in more than 1 sequence. This study gives species-specific insight into the cardiac transcriptome of turkeys and provides resources for future studies of cardiac function.