Turkey Genome Research Articles

A new haplotype-resolved turkey genome to enable turkey genetics and genomics research.

The domesticated turkey (Meleagris gallopavo) is a species of significant agricultural importance and is the second largest contributor, behind broiler chickens, to world poultry meat production. The previous genome is of draft quality and partly based on the chicken (Gallus gallus) genome. A high-quality reference genome of M. gallopavo is essential for turkey genomics and genetics research and the breeding industry. By adopting the trio-binning approach, we were able to assemble a high-quality chromosome-level F1 assembly and 2 parental haplotype assemblies, leveraging long-read technologies and genome-wide chromatin interaction data (Hi-C). From a total of 40 chromosomes (2n = 80), we captured 35 chromosomes in a single scaffold, showing much improved genome completeness and continuity compared to the old assembly build. The 3 assemblies are of higher quality than the previous draft quality assembly and comparable to the chicken assemblies (GRCg7) shown by the largest contig N50 (26.6 Mb) and comparable BUSCO gene set completeness scores (96-97%). Comparative analyses confirm a previously identified large inversion of around 19 Mbp on the Z chromosome not found in other Galliformes. Structural variation between the parent haplotypes was identified, which poses potential new target genes for breeding. We contribute a new high-quality turkey genome at the chromosome level, benefiting turkey genetics and other avian genomics research as well as the turkey breeding industry.

Investigating inbreeding in the turkey (Meleagris gallopavo) genome

The detrimental effects of increased homozygosity due to inbreeding have prompted the development of methods to reduce inbreeding. The detection of runs of homozygosity (ROH), or contiguous stretches of homozygous marker genotypes, can be used to describe and quantify the level of inbreeding in an individual. The estimation of inbreeding coefficients can be calculated based on pedigree information, ROH, or the genomic relationship matrix. The aim of this study was to detect and describe ROH in the turkey genome and compare estimates of pedigree-based inbreeding coefficients (FPED) with genomic-based inbreeding coefficients estimated from ROH (FROH) and the genomic relationship matrix (FGRM). A total of 2,616,890 pedigree records were available. Of these records, 6,371 genotyped animals from three purebred turkey (Meleagris gallopavo) lines between 2013 and 2019 were available, and these were obtained using a dense single nucleotide polymorphism array (56,452 SNPs). The overall mean length of detected ROH was 2.87 ± 0.29 Mb with a mean number of 84.87 ± 8.79 ROH per animal. Short ROH with lengths of 1 to 2 Mb long were the most abundant throughout the genome. Mean ROH coverage differed greatly between chromosomes and lines. Considering inbreeding coefficient means across all lines, genomic derived inbreeding coefficients (FROH = 0.27; FGRM = 0.32) were higher than coefficients estimated from pedigree records (FPED = 0.14). Correlations between FROH and FPED, FROH and FGRM, and FPED and FGRM ranged between 0.19 to 0.31, 0.68 to 0.73, and 0.17 to 0.30, respectively. Additionally, correlations between FROH from different lengths and FPED substantially increased with ROH length from -0.06 to 0.33. Results of the current research, including the distribution of ROH throughout the genome and ROH-derived inbreeding estimates, can provide a more comprehensive description of inbreeding in the turkey genome. This knowledge can be used to evaluate genetic diversity, a requirement for genetic improvement, and develop methods to minimize inbreeding in turkey breeding programs.

PSIII-2 Assessment of runs of homozygosity and estimates of inbreeding in three purebred turkey (Meleagris gallopavo) lines

Abstract Runs of homozygosity (ROH) are continuous stretches of homozygous genotypes in an individual that have been passed down from a common ancestor and can be used to accurately characterize genomic data. These ROH are correlated with other measures of inbreeding and have been applied to quantify individual autozygosity. The aim of this study was to detect and describe ROH in the turkey genome and estimate and compare measures of pedigree-based inbreeding coefficients (FPED) and genomic-based inbreeding coefficients estimated from ROH (FROH) and the genomic relationship matrix (FGRM). Pedigree records (n = 2,616,890) and genotypic records (n = 6,371) were available from three purebred turkey (Meleagris gallopavo) lines. Genotypic records were collected between 2013 and 2019 and were obtained using a dense single nucleotide polymorphism array (56,452 SNP). The overall mean length of detected ROH per animal was 2.87±0.29 Mb and mean number of ROH per animal was 84.87±8.79. Short ROH with lengths of 1–2 Mb long were the most abundant throughout the genome, accounting for approximately 45% of the identified segments. Mean ROH coverage differed greatly between chromosomes and lines. Across all lines, genomic derived inbreeding coefficients (FROH=0.27; FGRM=0.32) were higher than coefficients estimated from pedigree records (FPED=0.14). Ranges of correlations between FROH and FPED (0.19–0.31), FROH and FGRM (0.68–0.73), and FPED and FGRM (0.17–0.30) were estimated. Results from the current research provide a fundamental description of inbreeding in the turkey genome which is critical considering the growing concerns of the detrimental effects of increased inbreeding on fitness and health in livestock production. The knowledge gained from this study may subsequently be used to evaluate and maintain genetic diversity that is necessary for genetic improvement and minimizing inbreeding in turkey breeding programs.

PSIII-16 Genome-wide association mapping and functional analysis of body weight, feed intake and walking ability in turkeys

Abstract The underlying genetic mechanisms affecting turkey growth traits have not been widely investigated. Over the last few years, genome-wide association studies (GWAS) became the de facto approach to identify candidate regions associated with complex phenotypes and diseases in livestock. In the present study, we performed GWAS to identify regions associated with growth traits, feed intake and walking ability in a breeding turkey line. This was followed by studying the functional evidence that may support the impact of those regions on the economic traits in turkeys. A total of 31,950 phenotypic records for body weight, feed intake and walking ability with genomic (56,393 SNP) data were provided by Hybrid Turkeys, Kitchener, Canada. The analysis was carried out using a mixed linear model with hatch-week-year and sex fitted as fixed effects and the accumulated effect of all markers captured by the genomic relationship matrix fitted as random polygenic effects. Significant markers were observed on several chromosomes across the turkey genome. For example, COL8A1 and RBPMS2 genes were identified on chromosome 1 and 12, respectively, and associated with body weight. Furthermore, a gene-set enrichment analysis was performed for each trait. These functional and positional analyses uncovered a number of gene ontology functional terms, Reactome pathways and Medical Subject Headings that showed significant enrichment of genes associated with the studied traits. Many of the observed gene ontology functional terms (e.g., skeletal muscle tissue growth, regulation of digestive system process, and adult walking behavior) are known to be related to body growth, feed intake and walking ability. The results of this study revealed novel candidate genomic regions and candidate genes that could be managed within a turkey breeding program.

Discovering lethal alleles across the turkey genome using a transmission ratio distortion approach.

SummaryDeviation from Mendelian inheritance expectations (transmission ratio distortion, TRD) has been observed in several species, including the mouse and humans. In this study, TRD was characterized in the turkey genome using both allelic (specific‐ and unspecific‐parent TRD) and genotypic (additive‐ and dominance‐TRD) parameterizations within a Bayesian framework. In this study, we evaluated TRD for 23 243 genotyped Turkeys across 56 393 autosomal SNPs. The analyses included 500 sires, 2013 dams and 11 047 offspring (trios). Three different haplotype sliding windows of 4, 10 and 20 SNPs were used across the autosomal chromosomes. Based on the genotypic parameterizations, 14 haplotypes showed additive and dominance TRD effects highlighting regions with a recessive TRD pattern. In contrast, the allelic model uncovered 12 haplotype alleles with the allelic TRD pattern which showed an underrepresentation of heterozygous offspring in addition to the absence of homozygous animals. For regions with the allelic pattern, only one particular region showed a parent‐specific TRD where the penetrance was high via the dam, but low via the sire. The gene set analysis uncovered several gene ontology functional terms, Reactome pathways and several Medical Subject Headings that showed significant enrichment of genes associated with TRD. Many of these gene ontology functional terms (e.g. mitotic spindle assembly checkpoint, DRM complex and Aneuploidy), Reactome pathways (e.g. Mismatch repair) and Medical Subject Headings (e.g. Adenosine monophosphate) are known to be related to fertility, embryo development and lethality. The results of this study revealed potential novel candidate lethal haplotypes, functional terms and pathways that may enhance breeding programs in Turkeys through reducing mortality and improving reproduction rate.

Transcriptome analysis of turkey (Meleagris gallopavo) reproductive tract revealed key pathways regulating spermatogenesis and post-testicular sperm maturation

The application of transcriptomics to the study of the reproductive tract in male turkeys can significantly increase our current knowledge regarding the specifics of bird reproduction. To characterize the complex transcriptomic changes that occur in the testis, epididymis, and ductus deferens, deep sequencing of male turkey RNA samples (n = 6) was performed, using Illumina RNA-Seq. The obtained sequence reads were mapped to the turkey genome, and relative expression values were calculated to analyze differentially expressed genes (DEGs). Statistical analysis revealed 1,682; 2,150; and 340 DEGs in testis/epididymis, testis/ductus deferens, and epididymis/ductus deferens comparisons, respectively. The expression of selected genes was validated using quantitative real-time reverse transcriptase-polymerase chain reaction. Bioinformatics analysis revealed several potential candidate genes involved in spermatogenesis, spermiogenesis and flagellum formation in the testis, and in post-testicular sperm maturation in the epididymis and ductus deferens. In the testis, genes were linked with the mitotic proliferation of spermatogonia and the meiotic division of spermatocytes. Histone ubiquitination and protamine phosphorylation were shown to be regulatory mechanisms for nuclear condensation during spermiogenesis. The characterization of testicular transcripts allowed a better understanding of acrosome formation and development and flagellum formation, including axoneme structures and functions. Spermatozoa motility during post-testicular maturation was linked to the development of flagellar actin filaments and biochemical processes, including Ca2+ influx and protein phosphorylation/dephosphorylation. Spermatozoa quality appeared to be controlled by apoptosis and antioxidant systems in the epididymis and ductus deferens. Finally, genes associated with reproductive system development and morphogenesis were identified. To the best of our knowledge, this is the first genome-wide functional investigation of genes associated with tissue-specific processes in turkey reproductive tract. A catalog of genes worthy of further studies to understand the avian reproductive physiology and regulation was provided.

De novo annotation of the transcriptome of the Northern Wheatear (Oenanthe oenanthe).

We have sequenced a partial transcriptome of the Northern Wheatear (Oenanthe oenanthe), a species with one of the longest migrations on Earth. The transcriptome was constructed de novo using RNA-Seq sequence data from the pooled mRNA of six different tissues: brain, muscle, intestine, liver, adipose tissue and skin. The samples came from nine captive-bred wheatears collected at three different stages of the endogenous autumn migratory period: (1) lean birds prior the onset of migration, (2) during the fattening stage and (3) individuals at their migratory body mass plateau, when they have almost doubled their lean body mass. The sample structure used to build up the transcriptome of the Northern Wheatears concerning tissue composition and time guarantees the future survey of the regulatory genes involved in the development of the migratory phenotype. Through the pre-migratory period, birds accomplish outstanding physical and behavioural changes that involve all organ systems. Nevertheless, the molecular mechanisms through which birds synchronize and control hyperphagia, fattening, restlessness increase, immunity boosting and tuning the muscles for such endurance flight are still largely unknown. The use of RNA-Seq has emerged as a powerful tool to analyse complex traits on a broad scale, and we believe it can help to characterize the migratory phenotype of wheatears at an unprecedented level. The primary challenge to conduct quantitative transcriptomic studies in non-model species is the availability of a reference transcriptome, which we have constructed and described in this paper. The cDNA was sequenced by pyrosequencing using the Genome Sequencer Roche GS FLX System; with single paired-end reads of about 400 bp. We estimate the total number of genes at 15,640, of which 67% could be annotated using Turkey and Zebra Finch genomes, or protein sequence information from SwissProt and NCBI databases. With our study, we have made a first step towards understanding the migratory phenotype regarding gene expression of a species that has become a model to study birds long-distance migrations.

Few mitochondrial DNA sequences are inserted into the turkey (Meleagris gallopavo) nuclear genome: evolutionary analyses and informativity in the domestic lineage.

Mitochondrial DNA (mtDNA) insertions have been detected in the nuclear genome of many eukaryotes. These sequences are pseudogenes originated by horizontal transfer of mtDNA fragments into the nuclear genome, producing nuclear DNA sequences of mitochondrial origin (numt). In this study we determined the frequency and distribution of mtDNA-originated pseudogenes in the turkey (Meleagris gallopavo) nuclear genome. The turkey reference genome (Turkey_2.01) was aligned with the reference linearized mtDNA sequence using last. A total of 32 numt sequences (corresponding to 18 numt regions derived by unique insertional events) were identified in the turkey nuclear genome (size ranging from 66 to 1415bp; identity against the modern turkey mtDNA corresponding region ranging from 62% to 100%). Numts were distributed in nine chromosomes and in one scaffold. They derived from parts of 10 mtDNA protein-coding genes, ribosomal genes, the control region and 10 tRNA genes. Seven numt regions reported in the turkey genome were identified in orthologues positions in the Gallus gallus genome and therefore were present in the ancestral genome that in the Cretaceous originated the lineages of the modern crown Galliformes. Five recently integrated turkey numts were validated by PCR in 168 turkeys of six different domestic populations. None of the analysed numts were polymorphic (i.e. absence of the inserted sequence, as reported in numts of recent integration in other species), suggesting that the reticulate speciation model is not useful for explaining the origin of the domesticated turkey lineage.

Response of Turkey Muscle Satellite Cells to Thermal Challenge. II. Transcriptome Effects in Differentiating Cells.

Background: Exposure of poultry to extreme temperatures during the critical period of post-hatch growth can seriously affect muscle development and thus compromise subsequent meat quality. This study was designed to characterize transcriptional changes induced in turkey muscle satellite cells by thermal challenge during differentiation. Our goal is to better define how thermal stress alters breast muscle ultrastructure and subsequent development.Results: Skeletal muscle satellite cells previously isolated from the Pectoralis major muscle of 7-wk-old male turkeys (Meleagris gallopavo) from two breeding lines: the F-line (16 wk body weight-selected) and RBC2 (randombred control line) were used in this study. Cultured cells were induced to differentiate at 38°C (control) or thermal challenge temperatures of 33 or 43°C. After 48 h of differentiation, cells were harvested and total RNA was isolated for RNAseq analysis. Analysis of 39.9 Gb of sequence found 89% mapped to the turkey genome (UMD5.0, annotation 101) with average expression of 18,917 genes per library. In the cultured satellite cells, slow/cardiac muscle isoforms are generally present in greater abundance than fast skeletal isoforms. Statistically significant differences in gene expression were observed among treatments and between turkey lines, with a greater number of genes affected in the F-line cells following cold treatment whereas more differentially expressed (DE) genes were observed in the RBC2 cells following heat treatment. Many of the most significant pathways involved signaling, consistent with ongoing cellular differentiation. Regulation of Ca2+ homeostasis appears to be significantly affected by temperature treatment, particularly cold treatment.Conclusions: Satellite cell differentiation is directly influenced by temperature at the level of gene transcription with greater effects attributed to selection for fast growth. At lower temperature, muscle-associated genes in the satellite cells were among the genes with the greatest down regulation consistent with slower differentiation and smaller myotubes. Fewer expression differences were observed in the differentiating cells than previously observed for proliferating cells. This suggests the impact of temperature on satellite cells occurs primarily at early points in satellite cell activation.

Identification of CR1 retroposons in Arborophila rufipectus and their application to Phasianidae phylogeny.

Chicken repeat 1 (CR1), a member of non-LTR retroposon, is an important phylogenetic marker in avian systematics. In this study, we reported several characteristics of CR1 elements in a draft genome of Arborophila rufipectus (Sichuan partridge). According to the analyses of RepeatMasker, approximately 254966 CR1 elements were identified in A.rufipectus, covering 6.7% of the genome. Subsequently, we selected eighteen novel CR1 elements by comparing the chicken genome, turkey genome and assembled A.rufipectus scaffolds. Here, a combined data set comprising of 22 CR1 loci, mitochondrial genomes and eight unlinked introns was analysed to infer the evolutionary relationships of twelve Phasianidae species. The applicability of CR1 sequences for inferring avian phylogeny relative to mtDNA and intron sequences was investigated as well. Our results elucidated the position of A.rufipectus in Phasianidae with robust supports that it presented a sister clade to Arborophila ardens/Arborophila brunneopectus, and implied that genus Arborophila was in a basal phylogenetic position within Phasianidae and a phylogenetic affinity between Meleagris gallopavo and Pucrasia macrolopha. Therefore, this work not only resolved some of the confounding relationships among Phasianidae, but also suggested CR1 sequences could provide powerful complementary data for phylogeny reconstruction.

Cloning, Sequencing, and Expression of Selenoprotein Transcripts in the Turkey (Meleagris gallopavo).

The minimum Se requirement for male turkey poults is 0.3 μg Se/g – three times higher than requirements found in rodents – based on liver and gizzard glutathione peroxidase-4 (GPX4) and GPX1 activities. In addition, turkey liver GPX4 activity is 10-fold higher and GPX1 activity is 10-fold lower than in rats, and both GPX1 and GPX4 mRNA levels are dramatically down-regulated by Se deficiency. Currently, the sequences of all annotated turkey selenoprotein transcripts and proteins in the NCBI database are only “predicted.” Thus we initiated cloning and sequencing of the full turkey selenoprotein transcriptome to demonstrate expression of selenoprotein transcripts in the turkey, and to develop tools to investigate Se regulation of the full selenoproteome. Total RNA was isolated from six tissues of Se-adequate adult tom turkeys, and used to prepare reverse-transcription cDNA libraries. PCR primers were designed, based initially on chicken, rodent, porcine, bovine and human sequences and later on turkey shotgun cloning sequences. We report here the cloning of full transcript sequences for 9 selenoproteins, and 3'UTR portions for 15 additional selenoproteins, which include SECIS elements in 22 3'UTRs, and in-frame Sec (UGA) codons within coding regions of 19 selenoproteins, including 12 Sec codons in SEPP1. In addition, we sequenced the gap between two contigs from the shotgun cloning of the turkey genome, and found the missing sequence for the turkey Sec-tRNA. RTPCR was used to determine the relative transcript expression in 6 tissues. GPX3 expression was high in all tissues except kidney, GPX1 expression was high in kidney, SEPW1 expression was high in heart, gizzard and muscle, and SELU expression was high in liver. SEPP2, a selenoprotein not found in mammals, was highly expressed in liver but not in other tissues. In summary, transcripts for 24 selenoproteins are expressed in the turkey, not just predicted.

Lineage-specific loss of FGF17 within the avian orders Galliformes and Passeriformes

The genomic and developmental complexity of vertebrates is commonly attributed to two rounds of whole genome duplications which occurred at the base of the vertebrate radiation. These duplications led to the rise of several, multi-gene families of developmental proteins like the fibroblast growth factors (FGFs); a signaling protein family which functions at various stages of embryonic development. One of the major FGF assemblages arising from these duplications is the FGF8 subfamily, which includes FGF8, FGF17, and FGF18 in tetrapods. While FGF8 and FGF18 are found in all tetrapods and are critical for embryonic survival, genomic analyses suggest putative loss of FGF17 in various lineages ranging from frogs and fish, to the chicken. This study utilizes 27 avian genomes in conjunction with molecular analyses of chicken embryos to confirm the loss of FGF17 in chicken as a true, biological occurrence. FGF17 is also missing in the turkey, black grouse, Japanese quail and northern bobwhite genomes. These species, along with chicken, form a monophyletic clade in the order Galliformes. Four additional species, members of the clade Passeroidea, within the order Passeriformes, are also missing FGF17. Additionally, analysis of intact FGF17 in other avian lineages reveals that it is still under strong purifying selection, despite being seemingly dispensable. Thus, FGF17 likely represents a molecular spandrel arising from a genome duplication event and due to its high connectivity with FGF8/FGF18, and potential for interference with their function, is retained under strong purifying selection, despite itself not having a strong selective advantage.

In silico identification and characterisation of 17 polymorphic anonymous non-coding sequence markers (ANMs) for red grouse (Lagopus lagopus scotica)

Anonymous non-coding sequence markers (ANMs) are powerful neutral genetic markers with great utility in phylogeography, population genetics and population genomics. Developing ANMs has previously relied on sequencing random fragments of genomic DNA in the target species and then querying bioinformatics databases to identify unannotated, putatively neutral fragments. Here, we describe an alternative in silico approach that is based on identifying large unannotated genomic regions in model species to provide a priori neutral targets for candidate ANMs that are remote from exonic regions. We illustrate this approach by developing a set of 17 polymorphic ANMs for red grouse (Lagopus lagopus scotica) from c. 1 Mbp non-coding chromosome regions of chicken, turkey and zebrafinch genomes. This pipeline presents a powerful and efficient approach when appropriate model genomes are available for the target species of interest.

Genome-wide candidate regions for selective sweeps revealed through massive parallel sequencing of DNA across ten turkey populations.

BackgroundThe domestic turkey (Meleagris gallopavo) is an important agricultural species that is largely used as a meat-type bird. Characterizing genetic variation in populations of domesticated species and associating these variation patterns with the evolution, domestication, and selective breeding is critical for understanding the dynamics of genomic change in these species. Intense selective breeding and population bottlenecks are expected to leave signatures in the genome of domesticated species, such as unusually low nucleotide diversity or the presence of exceptionally extended haplotype homozygosity. These patterns of variation in selected populations are highly useful to not only understand the consequences of selective breeding and population dynamics, but also to provide insights into biological mechanisms that may affect physiological processes important to bring changes in phenotype of interest.ResultsWe observed 54 genomic regions in heritage and commercial turkey populations on 14 different chromosomes that showed statistically significant (P < 0.05) reduction in genomic variation indicating candidate selective sweeps. Areas with evidence of selective sweeps varied from 1.5 Mb to 13.8 Mb in length. Out of these 54 sweeps, 23 overlapped at least partially between two or more populations. Overlapping sweeps were found on 13 different chromosomes. The remaining 31 sweeps were population-specific and were observed on 12 different chromosomes, with 26 of these regions present only in commercial populations. Genes that are known to affect growth were enriched in the sweep regions.ConclusionThe turkey genome showed large sweep regions. The relatively high number of sweep regions in commercial turkey populations compared to heritage varieties and the enrichment of genes important to growth in these regions, suggest that these sweeps are the result of intense selection in these commercial lines, moving specific haplotypes towards fixation.Electronic supplementary materialThe online version of this article (doi:10.1186/s12863-014-0117-4) contains supplementary material, which is available to authorized users.

A Comparison of the Chicken and Turkey Proteomes and Phosphoproteomes in the Development of Poultry-Specific Immuno-Metabolism Kinome Peptide Arrays.

The use of species-specific peptide arrays for the study of animal kinomes has a proven track record of success. This technique has been used in a variety of species for the study of host–pathogen interactions and metabolism. Species-specific peptide arrays have been designed previously for use with chicken but a turkey array has never been attempted. In addition, arrays designed around individual cellular functions have been designed and utilized, but cross-function immuno-metabolic arrays have not been considered previously. Antecedent to designing separate chicken and turkey immuno-metabolic kinome peptide arrays, we show that while the chicken and turkey genomes are quite similar, the two species are much more distinct at the proteome and phosphoproteome levels. Despite a genome identity of approximately 90%, we observe that only 83% of chicken and turkey orthologous proteins display sequence matches between the two species. Further, less than 70% of kinase recognition target sequences are exact matches between chicken and turkey. Thus, our analysis shows that, at the proteome and kinome level, these two species must be considered separately in the design of novel peptide arrays. Our ultimate array design covers numerous immune and metabolic processes including innate and adaptive immunity, inflammatory responses, carbohydrate, protein, and fat metabolism, and response to hormones. We have shown the proteomic and phosphoproteomic diversity of chicken and turkey and have designed a valuable research tool for the study of immuno-metabolism within these two species.

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.

Next-generation sequencing strategies for characterizing the turkey genome

The turkey genome sequencing project was initiated in 2008 and has relied primarily on next-generation sequencing (NGS) technologies. Our first efforts used a synergistic combination of 2 NGS platforms (Roche/454 and Illumina GAII), detailed bacterial artificial chromosome (BAC) maps, and unique assembly tools to sequence and assemble the genome of the domesticated turkey, Meleagris gallopavo. Since the first release in 2010, efforts to improve the genome assembly, gene annotation, and genomic analyses continue. The initial assembly build (2.01) represented about 89% of the genome sequence with 17X coverage depth (931 Mb). Sequence contigs were assigned to 30 of the 40 chromosomes with approximately 10% of the assembled sequence corresponding to unassigned chromosomes (ChrUn). The sequence has been refined through both genome-wide and area-focused sequencing, including shotgun and paired-end sequencing, and targeted sequencing of chromosomal regions with low or incomplete coverage. These additional efforts have improved the sequence assembly resulting in 2 subsequent genome builds of higher genome coverage (25X/Build3.0 and 30X/Build4.0) with a current sequence totaling 1,010 Mb. Further, BAC with end sequences assigned to the Z/W and MG18 (MHC) chromosomes, ChrUn, or not placed in the previous build were isolated, deeply sequenced (Hi-Seq), and incorporated into the latest build (5.0). To aid in the annotation and to generate a gene expression atlas of major tissues, a comprehensive set of RNA samples was collected at various developmental stages of female and male turkeys. Transcriptome sequencing data (using Illumina Hi-Seq) will provide information to enhance the final assembly and ultimately improve sequence annotation. The most current sequence covers more than 95% of the turkey genome and should yield a much improved gene level of annotation, making it a valuable resource for studying genetic variations underlying economically important traits in poultry.

The Turkey Ig-Like Receptor Family: Identification, Expression and Function

The chicken leukocyte receptor complex located on microchromosome 31 encodes the chicken Ig-like receptors (CHIR), a vastly expanded gene family which can be further divided into three subgroups: activating CHIR-A, bifunctional CHIR-AB and inhibitory CHIR-B. Here, we investigated the presence of CHIR homologues in other bird species. The available genome databases of turkey, duck and zebra finch were screened with different strategies including BLAST searches employing various CHIR sequences, and keyword searches. We could not identify CHIR homologues in the distantly related zebra finch and duck, however, several partial and complete sequences of CHIR homologues were identified on chromosome 3 of the turkey genome. They were designated as turkey Ig-like receptors (TILR). Using cDNA derived from turkey blood and spleen RNA, six full length TILR could be amplified and further divided according to the typical sequence features into one activating TILR-A, one inhibitory TILR-B and four bifunctional TILR-AB. Since the TILR-AB sequences all displayed the critical residues shown to be involved in binding to IgY, we next confirmed the IgY binding using a soluble TILR-AB1-huIg fusion protein. This fusion protein reacted with IgY derived from various gallinaceous birds, but not with IgY from other bird species. Finally, we tested various mab directed against CHIR for their crossreactivity with either turkey or duck leukocytes. Whereas no staining was detectable with duck cells, the CHIR-AB1 specific mab 8D12 and the CHIR-A2 specific mab 13E2 both reacted with a leukocyte subpopulation that was further identified as thrombocytes by double immunofluorescence employing B-cell, T-cell and thrombocyte specific reagents. In summary, although the turkey harbors similar LRC genes as the chicken, their distribution seems to be distinct with predominance on thrombocytes rather than lymphocytes.

The proteome of the calcified layer organic matrix of turkey (Meleagris gallopavo) eggshell

BackgroundChicken eggshell mineralization is a prominent model for biomineralization not only because of its importance for avian reproduction but also because of the commercial interest associated with eggshell quality. An analysis and comparison of the protein constituents of eggshells of several species would contribute to a better understanding of the shell mineralization process. The recent publication of the turkey genome sequence now provides a basis for the in-depth analysis of the turkey eggshell proteome.ResultsProteomic analysis of turkey acid-soluble and acid-insoluble organic eggshell matrix yielded 697 identified proteins/protein groups. However, intensity-based absolute quantification (iBAQ) results indicated that the 47 most abundant identified proteins already constituted 95% of the total turkey eggshell matrix proteome. Forty-four of these proteins were also identified in chicken eggshell matrix previously. Despite these similarities there were important and unexpected differences. While ovocleidin-116 and ovocalyxin-36 were major proteins constituting approximately 37% of the identified proteome, other members of the group of so-called eggshell-specific proteins were not identified. Thus ovocalyxin-21 and ovocalyxin-32 were missing among matrix proteins. Conversely, major turkey eggshell proteins were not detected in chicken, such as the bone protein periostin, the mammalian counterpart of which is involved in many aspects of bone metabolism and which represented 10-11% of the total identified proteome.ConclusionsEven members of the same avian family show important differences in eggshell matrix composition and more studies on the proteome and the transcriptome level will be necessary to identify a common toolkit of eggshell mineralization and to work out species differences among functional eggshell protein sets and their role in eggshell production.

Unexpected Diversity and Expression of Avian Endogenous Retroviruses

ABSTRACTEndogenous retroviruses (ERVs) were identified and characterized in three avian genomes to gain insight into early retroviral evolution. Using the computer program RetroTector to detect relatively intact ERVs, we identified 500 ERVs in the chicken genome, 150 in the turkey genome, and 1,200 in the zebra finch genome. Previous studies suggested that endogenous alpharetroviruses were present in chicken genomes. In this analysis, a small number of alpharetroviruses were seen in the chicken and turkey genomes; however, these were greatly outnumbered by beta-like, gamma-like, and alphabeta proviruses. While the avian ERVs belonged to the same major groups as mammalian ERVs, they were more heterogeneous. In particular, the beta-like viruses revealed an evolutionary continuum with the gradual acquisition and loss of betaretroviral markers and a transition from beta to alphabeta and then to alpharetroviruses. Thus, it appears that birds may resemble a melting pot for early ERV evolution. Many of the ERVs were integrated in clusters on chromosomes, often near centromeres. About 25% of the chicken ERVs were in or near cellular transcription units; this is nearly random. The majority of these integrations were in the sense orientation in introns. A higher-than-random number of integrations were >100 kb from the nearest gene. Deep-sequencing studies of chicken embryo fibroblasts revealed that about 20% of the 500 ERVs were transcribed and translated. A subset of these were also transcribed in vivo in chickens, showing tissue-specific patterns of expression.