Genetic architecture of white striping in turkeys (Meleagris gallopavo)

White striping (WS) is a myopathy of growing concern to the turkey industry. It is rising in prevalence and has negative consequences for consumer acceptance and the functional properties of turkey meat. The objective of this study was to conduct a genome-wide association study (GWAS) and functional analysis on WS severity. Phenotypic data consisted of white striping scored on turkey breast fillets (N = 8422) by trained observers on a 0-3 scale (none to severe). Of the phenotyped birds, 4667 genotypic records were available using a proprietary 65 K single nucleotide polymorphism (SNP) chip. The SNP effects were estimated using a linear mixed model with a 30-SNP sliding window approach used to express the percentage genetic variance explained. Positional candidate genes were those located within 50 kb of the top 1% of SNP windows explaining the most genetic variance. Of the 95 positional candidate genes, seven were further classified as functional candidate genes because of their association with both a significant gene ontology and molecular function term. The results of the GWAS emphasize the polygenic nature of the trait with no specific genomic region contributing a large portion to the overall genetic variance. Significant pathways relating to growth, muscle development, collagen formation, circulatory system development, cell response to stimulus, and cytokine production were identified. These results help to support published biological associations between WS and hypoxia and oxidative stress and provide information that may be useful for future-omics studies in understanding the biological associations with WS development in turkeys.

Detection of lymphoproliferative disease virus in Iowa Wild Turkeys (Meleagris gallopavo): Comparison of two sections of the proviral genome

An accurate diagnostic test is an essential aspect of successfully monitoring and managing wildlife diseases. Lymphoproliferative Disease Virus (LPDV) is an avian retrovirus that was first identified in domestic turkeys in Europe and was first reported in a Wild Turkey (Meleagris gallopavo) in the United States in 2009. It has since been found to be widely distributed throughout North America. The majority of studies have utilized bone marrow and PCR primers targeting a 413-nucleotide sequence of the gag gene of the provirus to detect infection. While prior studies have evaluated the viability of other tissues for LPDV detection (whole blood, spleen, liver, cloacal swabs) none to date have studied differences in detection rates when utilizing different genomic regions of the provirus. This study examined the effectiveness of another section of the provirus, a 335-nucleotide sequence starting in the U3 region of the LTR (Long Terminal Repeat) and extending into the Matrix of the gag region (henceforth LTR), for detecting LPDV. Bone marrow samples from hunter-harvested Wild Turkeys (n = 925) were tested for LPDV with the gag gene and a subset (n = 417) including both those testing positive and those where LPDV was not detected was re-tested with LTR. The positive percent agreement (PPA) was 97.1% (68 of 70 gag positive samples tested positive with LTR) while the negative percent agreement (NPA) was only 68.0% (236 of 347 gag negative samples tested negative with LTR). Cohen's Kappa (κ = 0.402, Z = 10.26, p<0.0001) and the McNemar test (OR = 55.5, p<0.0001) indicated weak agreement between the two gene regions. We found that in Iowa Wild Turkeys use of the LTR region identified LPDV in many samples in which we failed to detect LPDV using the gag region and that LTR may be more appropriate for LPDV surveillance and monitoring. However, neither region of the provirus resulted in perfect detection and additional work is necessary to determine if LTR is more reliable in other geographic regions where LPDV occurs.

Targeted detection and molecular epidemiology of turkey coronavirus spike gene variants in turkeys and chickens

Turkey coronavirus (TCoV) is a member of the Avian coronavirus species with infectious bronchitis virus (IBV), which is considered to be the source of TCoV. These 2 viruses are highly similar in all regions of their genomes, except for the spike gene, which is necessary for virus attachment. Although TCoV causes severe enteric disease in turkey poults, it does not cause clinical disease in chickens. However, considering that TCoV can infect chickens, it is important to distinguish TCoV from IBV in chickens. This is particularly true for chickens that are housed near turkeys and thus might be infected with TCoV and serve as a silent source of TCoV for turkeys. We developed and validated a real-time PCR assay to detect the spike gene of TCoV and sequenced a portion of this gene to evaluate the molecular epidemiology of TCoV infections associated with a commercial turkey premises in the United States in 2020-2021. We identified natural infections of TCoV in chickens, and based on the molecular epidemiology of the viruses detected, these chickens may have served as a source of infection for the commercial turkey premises located nearby.

Listeriosis with viral coinfections in 8 gray foxes, 8 wild turkeys, and 2 young cervids in the southeastern United States

Listeria monocytogenes is a bacterium that can cause disease in many species, including humans, livestock, and wildlife. Increased interactions via shared habitats may promote pathogen transmission among these groups. Our objectives were to evaluate the Southeastern Cooperative Wildlife Disease Study diagnostic data to characterize and compare L. monocytogenes-induced lesions and comorbidities in gray foxes and wild turkeys, and to describe cases of listeriosis in 2 cervids. From 1991-2020, 8 gray foxes, 8 wild turkeys, a neonatal elk, and a white-tailed deer fawn from several eastern states in the United States were diagnosed with listeriosis. All 8 foxes had hepatitis and/or hepatic necrosis with intralesional gram-positive bacilli, and concurrent canine distemper virus (CDV) infection; 2 of the foxes had been vaccinated recently for CDV. L. monocytogenes was cultured from the liver (6 of 8) or lung (2 of 8) of foxes. Lesions in wild turkeys included hepatocellular necrosis (3 of 8), heterophilic hepatitis (1 of 8), heterophilic granulomas (1 of 8), intrasinusoidal gram-positive bacilli without hepatic lesions (1 of 8), granulomatous dermatitis (1 of 8), and/or granulomatous myocarditis (2 of 8). Lymphoproliferative disease viral DNA was detected in 5 of 6 turkeys tested; reticuloendotheliosis viral DNA was detected in 2 of 3 turkeys tested. Both cervids had systemic listeriosis, with L. monocytogenes isolated from liver. Immunohistochemistry for Listeria spp. on select cases revealed immunolabeling in affected organs. Listeriosis was thus established as a cause of morbidity and mortality in 3 wildlife species, which often suffered from concurrent infections and likely immunosuppression.

Genome-wide association study reveals candidate genes relevant to body weight in female turkeys (Meleagris gallopavo)

The underlying genetic mechanisms affecting turkey growth traits have not been widely investigated. Genome-wide association studies (GWAS) is a powerful 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 18-week body weight in a turkey population. The data included body weight observations for 24,989 female turkeys genotyped based on a 65K SNP panel. The analysis was carried out using a univariate mixed linear model with hatch-week-year and the 2 top principal components fitted as fixed effects and the accumulated polygenic effect of all markers captured by the genomic relationship matrix as random. Thirty-three significant markers were observed on 1, 2, 3, 4, 7 and 12 chromosomes, while 26 showed strong linkage disequilibrium extending up to 410 kb. These significant markers were mapped to 37 genes, of which 13 were novel. Interestingly, many of the investigated genes are known to be involved in growth and body weight. For instance, genes AKR1D1, PARP12, BOC, NCOA1, ADCY3 and CHCHD7 regulate growth, body weight, metabolism, digestion, bile acid biosynthetic and development of muscle cells. In summary, the results of our study revealed novel candidate genomic regions and candidate genes that could be managed within a turkey breeding program and adapted in fine mapping of quantitative trait loci to enhance genetic improvement in this species.

Accuracy of genomic selection for reducing susceptibility to pendulous crop in turkey (Meleagris gallopavo)

Pendulous crop (PC) in the turkey occurs when the crop distends from its normal position, thereby preventing the movement of feed and water from the crop down into the digestive system. This condition negatively impacts the turkey industry at both production and welfare levels. In this study, we estimated the genetic parameters for PC incidence and its genetic correlation with 5 production traits. Additionally, we evaluated the prediction accuracy and bias of breeding values for the selection candidates using pedigree (BLUP) or pedigree-genomic (ssGBLUP) relationships among the animals. A total of 245,783 turkey records were made available by Hybrid Turkeys, Kitchener, Canada. Of these, 6,545 were affected with PC. In addition, the data included 9,634 records for breast meat yield (BMY); 5,592 records for feed conversion ratio (FCR) and residual feed intake (RFI) in males; 170,844 records for body weight (BW) and walking score (WS) between 18 and 20 wk of age for males (71,012) and females (99,832), respectively. Among this population, 36,830 were genotyped using a 65K SNP Illumina Inc. chip. While all animals passed the quality control criteria, only 53,455 SNP markers were retained for subsequent analysis. Heritability for PC was estimated at 0.16 ± 0.00 and 0.17 ± 0.00 using BLUP and ssGBLUP, respectively. The incidence of PC was not genetically correlated with WS or FCR. Low unfavourable genetic correlations with BW (0.12 and 0.14), BMY (0.24 and 0.24) and RFI (-0.33 and -0.28) were obtained using BLUP and ssGBLUP, respectively. Using ssGBLUP showed higher prediction accuracy (0.51) for the breeding values for the selection candidates than the pedigree-based model (0.35). Whereas the bias of the prediction was slightly reduced with ssGBLUP (0.33 ± 0.05) than BLUP (0.30 ± 0.08), both models showed a regression coefficient lower than one, indicating inflation in the predictions. The results of this study suggest that PC is a heritable trait and selection for lower PC incidence rates is feasible. Although further investigation is necessary, selection for BW, BMY, and RFI may increase PC incidence. Incorporating genomic information would lead to higher accuracy in predicting the genetic merit for selection candidates.

Performance and processing yield comparisons of Large White male turkeys by genetic lines, sources, and seasonal rearing

Large White male turkey genetic lines (GL) comparison in performance and processing yields under the same conditions are rare in the literature. Two rearing experiments (EXP) were conducted to accomplish 2 objectives. The first objective was to test the effects of poult source and genetic lines on performance and processing yields. The second objective was to extract season and growth patterns when comparing both EXP common treatments. In EXP 1, male poults from 5 different sources were randomly assigned to 48 concrete: litter-covered floor pens. In EXP 2, male poults from 7 different genetic lines were randomly assigned to 48 concrete: litter-covered floor pens. For both EXP, the experimental design was a completely randomized block design with a one-factor arrangement. Both EXP were placed in the same house with the same management and nutrition in two separate seasons of the same year. Bird performance and carcass processing yield were analyzed in SAS 9.4 or JMP 15.1 in a mixed model. In EXP 1 no significant difference in BW or processing yield was observed. However, a similar GL from a commercial hatchery had an improved feed conversion ratio (FCR) over the same GL sourced directly from the genetic company hatchery. In EXP 2, statistical differences were observed in performance and breast meat yield depending on the GL. A season effect was observed when comparing the two EXP. Birds raised in the fall season had a 2 kg BW increase, on average, over their spring counterparts. This difference in BW can also be observed in a statistically higher breast meat yield by the birds raised in the fall over the ones raised in the spring. In conclusion, a comparison between GL resulted in effects due to genetic line, poult source, and rearing season on bird performance and carcass yield.

Transcriptome and Proteome Analysis Revealed Key Pathways Regulating Final Stage of Oocyte Maturation of the Turkey (Meleagris gallopavo)

In birds, the zona pellucida (ZP) matrix that surrounds the ovulated oocyte—called the inner perivitelline layer—is involved in sperm–zona interaction and successful fertilization. To identify the important genes and proteins connected with the final step of egg development, next-generation sequencing and two-dimensional electrophoresis, combined with mass spectrometry, were used for the analysis of mature oocytes at the F1 developmental stage. A total of 8161 genes and 228 proteins were annotated. Six subfamilies of genes, with codes ZP, ZP1–4, ZPD, and ZPAX, were identified, with the dominant expression of ZPD. The main expression site for ZP1 was the liver; however, granulosa cells may also participate in local ZP1 secretion. A ubiquitination system was identified in mature oocytes, where ZP1 was found to be the main ubiquitinated protein. Analysis of transcripts classified in estrogen receptor (ESR) signaling indicated the presence of ESR1 and ESR2, as well as a set of estrogen-dependent genes involved in both genomic and nongenomic mechanisms for the regulation of gene expression by estrogen. Oxidative phosphorylation was found to be a possible source of adenosine triphosphate, and the nuclear factor erythroid 2-related factor 2 signaling pathway could be involved in the response against oxidative stress. Oocyte–granulosa cell communication by tight, adherens, and gap junctions seems to be essential for the final step of oocyte maturation.

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

Deviation 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.

Transcriptome analysis of the hypothalamus and pituitary of turkey hens with low and high egg production

BACKGROUND

High egg producing hens (HEPH) show increased hypothalamic and pituitary gene expression related to hypothalamo-pituitary-gonadal (HPG) axis stimulation as well as increased in vitro responsiveness to gonadotropin releasing hormone (GnRH) stimulation in the pituitary when compared to low egg producing hens (LEPH). Transcriptome analysis was performed on hypothalamus and pituitary samples from LEPH and HEPH to identify novel regulators of HPG axis function.

RESULTS

In the hypothalamus and pituitary, 4644 differentially expressed genes (DEGs) were identified between LEPH and HEPH, with 2021 genes up-regulated in LEPH and 2623 genes up-regulated in HEPH. In LEPH, up-regulated genes showed enrichment of the hypothalamo-pituitary-thyroid (HPT) axis. Beta-estradiol was identified as an upstream regulator regardless of tissue. When LEPH and HEPH samples were compared, beta-estradiol was activated in HEPH in 3 of the 4 comparisons, which correlated to the number of beta-estradiol target genes up-regulated in HEPH. In in vitro pituitary cell cultures from LEPH and HEPH, thyroid hormone pretreatment negatively impacted gonadotropin subunit mRNA levels in cells from both LEPH and HEPH, with the effect being more prominent in HEPH cells. Additionally, the effect of estradiol pretreatment on gonadotropin subunit mRNA levels in HEPH cells was negative, whereas estradiol pretreatment increased gonadotropin subunit mRNA levels in LEPH cells.

CONCLUSIONS

Up-regulation of the HPT axis in LEPH and upstream beta-estradiol activation in HEPH may play a role in regulating HPG axis function, and ultimately ovulation rates. Thyroid hormone and estradiol pretreatment impacted gonadotropin mRNA levels following GnRH stimulation, with the inhibitory effects of thyroid hormone more detrimental in HEPH and estradiol stimulatory effects more prominent in LEPH. Responsiveness to thyroid hormone and estradiol may be due to desensitization to thyroid hormone and estradiol in LEPH and HEPH, respectively, due to up-regulation of the HPT axis in LEPH and of the HPG axis in HEPH. Further studies will be necessary to identify possible target gene desensitization mechanisms and elicit the regulatory role of the HPT axis and beta-estradiol on ovulation rates in turkey hens.

The hepatic transcriptome of the turkey poult (Meleagris gallopavo) is minimally altered by high inorganic dietary selenium

There is interest in supplementing animals and humans with selenium (Se) above Se-adequate levels, but the only good biomarker for toxicity is tissue Se. We targeted liver because turkeys fed 5 μg Se/g have hepatic Se concentrations 6-fold above Se-adequate (0.4 μg Se/g) levels without effects on growth or health. Our objectives were (i) to identify transcript biomarkers for high Se status, which in turn would (ii) suggest proteins and pathways used by animals to adapt to high Se. Turkey poults were fed 0, 0.025, 0.4, 0.75 and 1.0 μg Se/g diet in experiment 1, and fed 0.4, 2.0 and 5.0 μg Se/g in experiment 2, as selenite, and the full liver transcriptome determined by RNA-Seq. The major effect of Se-deficiency was to down-regulate expression of a subset of selenoprotein transcripts, with little significant effect on general transcript expression. In response to high Se intake (2 and 5 μg Se/g) relative to Se-adequate turkeys, there were only a limited number of significant differentially expressed transcripts, all with only relatively small fold-changes. No transcript showed a consistent pattern of altered expression in response to high Se intakes across the 1, 2 and 5 μg Se/g treatments, and there were no associated metabolic pathways and biological functions that were significant and consistently found with high Se supplementation. Gene set enrichment analysis also found no gene sets that were consistently altered by high-Se and supernutritional-Se. A comparison of differentially expressed transcript sets with high Se transcript sets identified in mice provided high Se (~3 μg Se/g) also failed to identify common differentially expressed transcript sets between these two species. Collectively, this study indicates that turkeys do not alter gene expression in the liver as a homeostatic mechanism to adapt to high Se.

Potential authentication of various meat-based products using simple and efficient DNA extraction method

BACKGROUND

The growth of halal food consumption worldwide has resulted in an increase in the request for halal authentication. DNA-based detection using powerful real-time polymerase chain reaction (PCR) technique has been shown to be highly specific and sensitive authentication tool. The efficient DNA extraction method in terms of quality and quantity is a backbone step to obtain successful real-time PCR assays. In this study, different DNA extraction methods using three lysis buffers were evaluated and developed to recommend a much more efficient method as well as achieve a successful detection using real-time PCR.

RESULTS

The lysis buffer 2 (LB2) has been shown to be the best lysis buffer for DNA extraction from both raw and processed meat samples comparing to other lysis buffers tested. Hence, the LB2 has been found to be ideal to detect meat and porcine DNAs by real-time PCR using pairs of porcine specific primers and universal primers which amplified at 119 bp fragment and 93 bp fragment, respectively. This assay allows detection as low as 0.0001 ng of DNA. Higher efficiency and sensitivity of real-time PCR via a simplified DNA extraction method using LB2 have been observed, as well as a reproducible and high correlation coefficient (R²  = 0.9979) based on the regression analysis of the standard curve have been obtained.

CONCLUSION

This study has established a fast, simple, inexpensive and efficient DNA extraction method that is feasible for raw and processed meat products. This extraction technique allows an accurate DNA detection by real-time PCR and can also be implemented to assist the halal authentication of various meat-based products available in the market. © 2019 Society of Chemical Industry.

Research Note: Effect of selection for body weight on the adipogenic conversion of turkey myogenic satellite cells by Syndecan-4 and its covalently attached N-glycosylation chains

Adult myoblasts, satellite cells, will proliferate, and differentiate into myotubes in vitro. However, changes in environmental and nutritional conditions will result in the satellite cells differentiating into adipocyte-like cells synthesizing lipids. Prior research has shown that levels of N-glycosylation and heparan sulfate can promote or prevent the adipogenic conversion of myogenic satellite cells. Syndecan-4, an N-glycosylated heparan sulfate proteoglycan, has been shown to play key roles in satellite cell proliferation and migration. The objective of the current study was to determine if syndecan-4, and syndecan-4 N-glycosylation and heparan sulfate chain levels altered the conversion of satellite cells to an adipogenic cell fate and if growth selection affected the response of the satellite cells. Over-expression of syndecan-4, syndecan-4 without N-glycosylated chains but with its heparan sulfate chains attached, syndecan-4 without heparan sulfate chains with its N-glycosylation chains, and syndecan-4 without N-glycosylation and heparan sulfate chains was measured for lipid accumulation in pectoralis major muscle satellite cells isolated from the Randombred Control line 2 (RBC2) and 16 wk body weight (F line) turkeys. The F line was selected from the RBC2 line for only 16 wk body weight. Results from this study demonstrated that wild type levels of syndecan-4 and its covalently attached N-glycosylation chains play a key role in regulating the conversion of pectoralis major muscle satellite cells to an adipogenic lineage while selection for body weight was not a major contributing factor in this conversion.

Research Note: The effect of selection for 16-week body weight on turkey serum metabolome

The phenotype of modern commercial turkeys is substantially different than that of unselected, heritage turkey lines. These phenotypic changes have arisen from alterations in the genome/transcriptome, as well as the influence of many external factors on growth performance including nutrition, environment, and management. To investigate the phenotypic changes resulting from genetic selection for increased body weight, The Ohio State University maintains 2 unique genetic turkey lines: the randombred control (RBC2) line, which is comprised of genetics from 1960 era commercial turkeys and has been maintained without conscious selection for any trait; and the F line, which was originally selected from the RBC2 line and has been selected for increased 16 wk body weight for over 50 generations. This study used broad-spectrum mass-spectrometry profiling techniques to identify and quantify differences in the metabolome of the serum of F and RBC2 turkey lines. Serum samples from both F and RBC2 turkeys were subject to quantitative time of flight liquid chromatography tandem mass spectrometry analyses. Principle component analyses showed distinct populations of metabolites in the F vs. RBC2 serum, suggesting that increased body weight is associated with the accumulation of several metabolites. Comparing the spectral features to online databases resulted in the selection of 104 features with potentially identifiable chemical structures. Of these 104 features, 25 were found at higher levels in the serum of the RBC2 line turkeys, while 79 were found at a greater abundance in the F line turkeys. A more detailed analysis of these 104 features allowed for the putative identification of 49 compounds, which were clustered into 6 functional groups: 1) energy metabolism; 2) vitamins; 3) hormones and signaling molecules; 4) lipid derivatives, fatty acid metabolites, and membrane components; 5) amino acid/protein metabolism; and 6) microbial metabolites. Further validation and experimentation is needed to confirm the identity of these metabolites and understand their biological relevance and association with selection for increased body weight.

Differential expression of cell-cell junction proteins in the testis, epididymis, and ductus deferens of domestic turkeys (Meleagris gallopavo) with white and yellow semen

Tight, adherens, and gap junctions are involved in the regulation of reproductive tissue function in male mammals. In birds, including domestic turkeys, intercellular interactions performed by junctional networks have not yet been studied. Furthermore, the cellular and molecular basis of yellow semen syndrome (YSS) in the turkey population remains poorly understood. Thus, the aim of the present study was 2-fold: first, to provide new information on the localization and expression of cell-cell junction proteins in the testis, epididymis, and ductus deferens of domestic turkeys and second, to compare expression of junctional protein genes between 2 turkey population, one that produces white normal semen (WNS) and the other that produces yellow abnormal semen. Expression of occludin, zonula occludens-1 (ZO-1), connexin 43 (Cx43), N- and E-cadherin, and β-catenin genes were investigated using 3 complementary techniques: quantitative real-time PCR, western blot, and immunohistochemistry. Compared to WNS testis, epididymis, and ductus deferens, YSS tissues exhibited downregulation of occludin and β-catenin mRNA (P < 0.05) and protein (P < 0.05 and P < 0.01, respectively) and upregulation of N- and E-cadherin mRNA (P < 0.001, P < 0.05, P < 0.01, respectively) and protein (P < 0.01, P < 0.05, and P < 0.05, respectively). In contrast, ZO-1 and Cx43 mRNA and protein were upregulated in YSS testis (P < 0.05 and P < 0.001, respectively) but not in epididymis and ductus deferens; both mRNAs and proteins were downregulated (P < 0.05) compared to the respective WNS epididymis and ductus deferens. Altered staining intensity of immunoreactive proteins in YSS vs. WNS reproductive tissue sections confirmed the gene expression results. The present study is the first to demonstrate altered levels of junctional protein gene expression in reproductive tissues of male YSS turkeys. These findings may suggest that subtle changes in junctional protein expression affect the microenvironment in which spermatozoa develop and mature and thus may have an impact on the appearance of yellow semen in domestic turkeys.

Transmission dynamics between infected waterfowl and terrestrial poultry: Differences between the transmission and tropism of H5N8 highly pathogenic avian influenza virus (clade 2.3.4.4a) among ducks, chickens and turkeys

H5N8 highly-pathogenic avian influenza viruses (HPAIVs, clade 2.3.4.4) have spread globally via migratory waterfowl. Pekin ducks infected with a UK virus (H5N8-2014) served as the donors of infection in three separate cohousing experiments to attempt onward transmission chains to sequentially introduced groups of contact ducks, chickens and turkeys. Efficient transmission occurred among ducks and turkeys up to the third contact stage, with all (100%) birds becoming infected. Introduction of an additional fourth contact group of ducks to the turkey transmission chain demonstrated retention of H5N8-2014's waterfowl-competent adaptation. However, onward transmission ceased in chickens at the second contact stage where only 13% became infected. Analysis of viral progeny at this contact stage revealed no emergent polymorphisms in the intra-species (duck) transmission chain, but both terrestrial species included changes in the polymerase and accessory genes. Typical HPAIV pathogenesis and mortality occurred in infected chickens and turkeys, contrasting with 5% mortality among ducks.

Differential Transcriptome Responses to Aflatoxin B₁ in the Cecal Tonsil of Susceptible and Resistant Turkeys

The nearly-ubiquitous food and feed-borne mycotoxin aflatoxin B₁ (AFB₁) is carcinogenic and mutagenic, posing a food safety threat to humans and animals. One of the most susceptible animal species known and thus a good model for characterizing toxicological pathways, is the domesticated turkey (DT), a condition likely due, at least in part, to deficient hepatic AFB₁-detoxifying alpha-class glutathione S-transferases (GSTAs). Conversely, wild turkeys (Eastern wild, EW) are relatively resistant to the hepatotoxic, hepatocarcinogenic and immunosuppressive effects of AFB₁ owing to functional gene expression and presence of functional hepatic GSTAs. This study was designed to compare the responses in gene expression in the gastrointestinal tract between DT (susceptible phenotype) and EW (resistant phenotype) following dietary AFB₁ challenge (320 ppb for 14 days); specifically in cecal tonsil which functions in both nutrient absorption and gut immunity. RNAseq and gene expression analysis revealed significant differential gene expression in AFB₁-treated animals compared to control-fed domestic and wild birds and in within-treatment comparisons between bird types. Significantly upregulated expression of the primary hepatic AFB₁-activating P450 (CYP1A5) as well as transcriptional changes in tight junction proteins were observed in AFB₁-treated birds. Numerous pro-inflammatory cytokines, TGF-β and EGF were significantly down regulated by AFB₁ treatment in DT birds and pathway analysis suggested suppression of enteroendocrine cells. Conversely, AFB₁ treatment modified significantly fewer unique genes in EW birds; among these were genes involved in lipid synthesis and metabolism and immune response. This is the first investigation of the effects of AFB₁ on the turkey gastro-intestinal tract. Results suggest that in addition to the hepatic transcriptome, animal resistance to this mycotoxin occurs in organ systems outside the liver, specifically as a refractory gastrointestinal tract.

Comparative Response of the Hepatic Transcriptomes of Domesticated and Wild Turkey to Aflatoxin B₁

The food-borne mycotoxin aflatoxin B₁ (AFB₁) poses a significant risk to poultry, which are highly susceptible to its hepatotoxic effects. Domesticated turkeys (Meleagris gallopavo) are especially sensitive, whereas wild turkeys (M. g. silvestris) are more resistant. AFB₁ toxicity entails bioactivation by hepatic cytochrome P450s to the electrophilic exo-AFB₁-8,9-epoxide (AFBO). Domesticated turkeys lack functional hepatic GST-mediated detoxification of AFBO, and this is largely responsible for the differences in resistance between turkey types. This study was designed to characterize transcriptional changes induced in turkey livers by AFB₁, and to contrast the response of domesticated (susceptible) and wild (more resistant) birds. Gene expression responses to AFB₁ were examined using RNA-sequencing. Statistically significant differences in gene expression were observed among treatment groups and between turkey types. Expression analysis identified 4621 genes with significant differential expression (DE) in AFB₁-treated birds compared to controls. Characterization of DE transcripts revealed genes dis-regulated in response to toxic insult with significant association of Phase I and Phase II genes and others important in cellular regulation, modulation of apoptosis, and inflammatory responses. Constitutive expression of GSTA3 was significantly higher in wild birds and was significantly higher in AFB₁-treated birds when compared to controls for both genetic groups. This pattern was also observed by qRT-PCR in other wild and domesticated turkey strains. Results of this study emphasize the differential response of these genetically distinct birds, and identify genes and pathways that are differentially altered in aflatoxicosis.

Karyotype analysis and sex determination in Australian Brush-turkeys (Alectura lathami)

Sexual differentiation across taxa may be due to genetic sex determination (GSD) and/or temperature sex determination (TSD). In many mammals, males are heterogametic (XY); whereas females are homogametic (XX). In most birds, the opposite is the case with females being heterogametic (ZW) and males the homogametic sex (ZZ). Many reptile species lack sex chromosomes, and instead, sexual differentiation is influenced by temperature with specific temperatures promoting males or females varying across species possessing this form of sexual differentiation, although TSD has recently been shown to override GSD in Australian central beaded dragons (Pogona vitticeps). There has been speculation that Australian Brush-turkeys (Alectura lathami) exhibit TSD alone and/or in combination with GSD. Thus, we sought to determine if this species possesses sex chromosomes. Blood was collected from one sexually mature female and two sexually mature males residing at Sylvan Heights Bird Park (SHBP) and shipped for karyotype analysis. Karyotype analysis revealed that contrary to speculation, Australian Brush-turkeys possess the classic avian ZW/ZZ sex chromosomes. It remains a possibility that a biased primary sex ratio of Australian Brush-turkeys might be influenced by maternal condition prior to ovulation that result in her laying predominantly Z- or W-bearing eggs and/or sex-biased mortality due to higher sensitivity of one sex in environmental conditions. A better understanding of how maternal and extrinsic factors might differentially modulate ovulation of Z- or W-bearing eggs and hatching of developing chicks possessing ZW or ZZ sex chromosomes could be essential in conservation strategies used to save endangered members of Megapodiidae.

Satiety splits drinking behavior into bouts: Organization of drinking in turkeys

The regulation of the drinking behavior of animals is usually overlooked, and the traits associated with it are not well defined. We used a unique data set of measurements of individual drinking behavior in turkeys 1) to validate the system of data generation, 2) to develop a methodology to allow clustering of drinking events and splitting behavior into bouts, and 3) to develop traits related to drinking behavior and its regulation and investigate how these traits may be affected by bird genotype. Visits to drinkers were generated by an electronic, custom-made equipment that automatically measures the individual drinking behavior of a large number of turkeys from 3 different genetic lines. The overall reliability of the electronic system was estimated from video observations and resulted in a predictability of 98.8% and sensitivity of 98.6%. A novel method based on mixture distribution models allowed clustering of drinking events and splitting behavior into bouts by estimating the shortest interval between visits to the drinker that was considered to be part of a bout (bout criterion). The method predicted that after the end of a given bout the probability of the bird initiating the next bout was low but increased with time since the last bout. As a result, drinking bouts were not randomly distributed but were predicated on the physiological principle of satiety, suggesting that they constitute biologically appropriate units for expressing drinking behavior. The applied method resulted in bout criteria estimates of 665, 672, and 602 s for genetic lines A, B, and C, respectively. On the basis of this methodology, a number of drinking behavior traits, such as bout duration and frequency, and water intake per bout were identified that revealed differences ( < 0.01) in the drinking behavior between the turkey genetic lines. Similarly, time accumulation patterns of drinking behavior traits within a day differed ( < 0.01) within and between genetic lines, suggesting that variation in drinking behavior exists and birds use different behavioral strategies to meet their water intake requirements. Development of drinking behavior over time was similar between the lines, suggesting conservation of this behavioral organization. As well as providing ideas about the regulation of drinking behavior, the developed behavioral traits may be of practical relevance because water utilization, along with feed efficiency, is part of overall biological efficiency.

Comparative Quantitative Studies on the Microvasculature of the Heart of a Highly Selected Meat-Type and a Wild-Type Turkey Line

In this study the macroscopic and microscopic structure of the heart of a fast growing, meat-type turkey line (British United turkeys BUT Big 6) and a wild-type turkey line (Canadian Wild turkey) were compared. At 8 and 16 weeks of age, 10 birds of each genotype and sex were sampled. The body mass and heart mass of the meat-type turkey both increased at a faster rate than those of the wild-type turkey. However in both turkey lines, the relative heart mass decreased slightly with age, the decrease was statistically significant only in the male turkeys. Furthermore meat-type turkeys had a significantly (p < 0.01) lower relative heart mass and relative thickness of the left ventricle compared to the wild-type turkeys of the same age. The wild-type turkeys showed no significant change in the size of cardiomyocytes (cross sectional area and diameter) from 8 weeks to 16 weeks. In contrast, the size of cardiomyocytes increased significantly (p < 0.001) with age in the meat-type turkeys. The number of capillaries in the left ventricular wall increased significantly (p < 0.001) in wild-type turkeys from 2351 per mm² at the age of 8 weeks to 2843 per mm² at 16 weeks. However, in the meat-type turkeys there were no significant changes, capillary numbers being 2989 per mm² at age 8 weeks and 2915 per mm² at age 16 weeks. Correspondingly the area occupied by capillaries in the myocardium increased in wild-type turkeys from 8.59% at the age of 8 weeks to 9.15% at 16 weeks, whereas in meat-type turkeys this area decreased from 10.4% at 8 weeks to 9.95% at 16 weeks. Our results indicate a mismatch in development between body mass and heart mass and a compromised cardiac capillary density and architecture in the meat-type turkeys in comparison to the wild-type turkeys.

Hepatic Transcriptome Responses of Domesticated and Wild Turkey Embryos to Aflatoxin B₁

The mycotoxin, aflatoxin B₁ (AFB₁) is a hepatotoxic, immunotoxic, and mutagenic contaminant of food and animal feeds. In poultry, AFB₁ can be maternally transferred to embryonated eggs, affecting development, viability and performance after hatch. Domesticated turkeys (Meleagris gallopavo) are especially sensitive to aflatoxicosis, while Eastern wild turkeys (M. g. silvestris) are likely more resistant. In ovo exposure provided a controlled AFB₁ challenge and comparison of domesticated and wild turkeys. Gene expression responses to AFB₁ in the embryonic hepatic transcriptome were examined using RNA-sequencing (RNA-seq). Eggs were injected with AFB₁ (1 μg) or sham control and dissected for liver tissue after 1 day or 5 days of exposure. Libraries from domesticated turkey (n = 24) and wild turkey (n = 15) produced 89.2 Gb of sequence. Approximately 670 M reads were mapped to a turkey gene set. Differential expression analysis identified 1535 significant genes with |log₂ fold change| ≥ 1.0 in at least one pair-wise comparison. AFB₁ effects were dependent on exposure time and turkey type, occurred more rapidly in domesticated turkeys, and led to notable up-regulation in cell cycle regulators, NRF2-mediated response genes and coagulation factors. Further investigation of NRF2-response genes may identify targets to improve poultry resistance.

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.

A next generation semiconductor based sequencing approach for the identification of meat species in DNA mixtures

The identification of the species of origin of meat and meat products is an important issue to prevent and detect frauds that might have economic, ethical and health implications. In this paper we evaluated the potential of the next generation semiconductor based sequencing technology (Ion Torrent Personal Genome Machine) for the identification of DNA from meat species (pig, horse, cattle, sheep, rabbit, chicken, turkey, pheasant, duck, goose and pigeon) as well as from human and rat in DNA mixtures through the sequencing of PCR products obtained from different couples of universal primers that amplify 12S and 16S rRNA mitochondrial DNA genes. Six libraries were produced including PCR products obtained separately from 13 species or from DNA mixtures containing DNA from all species or only avian or only mammalian species at equimolar concentration or at 1:10 or 1:50 ratios for pig and horse DNA. Sequencing obtained a total of 33,294,511 called nucleotides of which 29,109,688 with Q20 (87.43%) in a total of 215,944 reads. Different alignment algorithms were used to assign the species based on sequence data. Error rate calculated after confirmation of the obtained sequences by Sanger sequencing ranged from 0.0003 to 0.02 for the different species. Correlation about the number of reads per species between different libraries was high for mammalian species (0.97) and lower for avian species (0.70). PCR competition limited the efficiency of amplification and sequencing for avian species for some primer pairs. Detection of low level of pig and horse DNA was possible with reads obtained from different primer pairs. The sequencing of the products obtained from different universal PCR primers could be a useful strategy to overcome potential problems of amplification. Based on these results, the Ion Torrent technology can be applied for the identification of meat species in DNA mixtures.

Tuberal hypothalamic expression of the glial intermediate filaments, glial fibrillary acidic protein and vimentin across the turkey hen (Meleagris gallopavo) reproductive cycle: Further evidence for a role of glial structural plasticity in seasonal reproduction

Glia regulate the hypothalamic-pituitary-gonadal (HPG) axis in birds and mammals. This is accomplished mechanically by ensheathing gonadotrophin-releasing hormone I (GnRH) nerve terminals thereby blocking access to the pituitary blood supply, or chemically in a paracrine manner. Such regulation requires appropriate spatial associations between glia and nerve terminals. Female turkeys (Meleagris gallopavo) use day length as a primary breeding cue. Long days activate the HPG-axis until the hen enters a photorefractory state when previously stimulatory day lengths no longer support HPG-axis activity. Hens must then be exposed to short days before reactivation of the reproductive axis occurs. As adult hens have discrete inactive reproductive states in addition to a fertile state, they are useful for examining the glial contribution to reproductive function. We immunostained tuberal hypothalami from short and long-day photosensitive hens, plus long-day photorefractory hens to examine expression of two intermediate filaments that affect glial morphology: glial fibrillary acidic protein (GFAP) and vimentin. GFAP expression was drastically reduced in the central median eminence of long day photosensitive hens, especially within the internal zone. Vimentin expression was similar among groups. However, vimentin-immunoreactive fibers abutting the portal vasculature were significantly negatively correlated with GFAP expression in the median eminence, which is consistent with our hypothesis for a reciprocal relationship between GFAP and vimentin expression. It appears that up-regulation of GFAP expression in the central median eminence of turkey hens is associated with periods of reproductive quiescence and that photofractoriness is associated with the lack of a glial cytoskeletal response to long days.

Sequencing of the core MHC region of black grouse (Tetrao tetrix) and comparative genomics of the galliform MHC

BACKGROUND

The MHC, which is regarded as the most polymorphic region in the genomes of jawed vertebrates, plays a central role in the immune system by encoding various proteins involved in the immune response. The chicken MHC-B genomic region has a highly streamlined gene content compared to mammalian MHCs. Its core region includes genes encoding Class I and Class IIB molecules but is only ~92Kb in length. Sequences of other galliform MHCs show varying degrees of similarity as that of chicken. The black grouse (Tetrao tetrix) is a wild galliform bird species which is an important model in conservation genetics and ecology. We sequenced the black grouse core MHC-B region and combined this with available data from related species (chicken, turkey, gold pheasant and quail) to perform a comparative genomics study of the galliform MHC. This kind of analysis has previously been severely hampered by the lack of genomic information on avian MHC regions, and the galliformes is still the only bird lineage where such a comparison is possible.

RESULTS

In this study, we present the complete genomic sequence of the MHC-B locus of black grouse, which is 88,390 bp long and contains 19 genes. It shows the same simplicity as, and almost perfect synteny with, the corresponding genomic region of chicken. We also use 454-transcriptome sequencing to verify expression in 17 of the black grouse MHC-B genes. Multiple sequence inversions of the TAPBP gene and TAP1-TAP2 gene block identify the recombination breakpoints near the BF and BLB genes. Some of the genes in the galliform MHC-B region also seem to have been affected by selective forces, as inferred from deviating phylogenetic signals and elevated rates of non-synonymous nucleotide substitutions.

CONCLUSIONS

We conclude that there is large synteny between the MHC-B region of the black grouse and that of other galliform birds, but that some duplications and rearrangements have occurred within this lineage. The MHC-B sequence reported here will provide a valuable resource for future studies on the evolution of the avian MHC genes and on links between immunogenetics and ecology of black grouse.

Significant selective constraint at 4-fold degenerate sites in the avian genome and its consequence for detection of positive selection

A major conclusion from comparative genomics is that many sequences that do not code for proteins are conserved beyond neutral expectations, indicating that they evolve under the influence of purifying selection and are likely to have functional roles. Due to the degeneracy of the genetic code, synonymous sites within protein-coding genes have previously been seen as "silent" with respect to function and thereby invisible to selection. However, there are indications that synonymous sites of vertebrate genomes are also subject to selection and this is not necessarily because of potential codon bias. We used divergence in ancestral repeats as a neutral reference to estimate the constraint on 4-fold degenerate sites of avian genes in a whole-genome approach. In the pairwise comparison of chicken and zebra finch, constraint was estimated at 24-32%. Based on three-species alignments of chicken, turkey, and zebra finch, lineage-specific estimates of constraint were 43%, 29%, and 24%, respectively. The finding of significant constraint at 4-fold degenerate sites from data on interspecific divergence was replicated in an analysis of intraspecific diversity in the chicken genome. These observations corroborate recent data from mammalian genomes and call for a reappraisal of the use of synonymous substitution rates as neutral standards in molecular evolutionary analysis, for example, in the use of the well-known d(N)/d(S) ratio and in inferences on positive selection. We show by simulations that the rate of false positives in the detection of positively selected genes and sites increases several-fold at the levels of constraint at 4-fold degenerate sites found in this study.

Fibroblast growth factor 2 and protein kinase C alpha are involved in syndecan-4 cytoplasmic domain modulation of turkey myogenic satellite cell proliferation

Syndecan-4 core protein is composed of extracellular, transmembrane, and cytoplasmic domains. The cytoplasmic domain functions in transmitting signals into the cell through the protein kinase C alpha (PKCα) pathway. The glycosaminoglycan (GAG) and N-linked glycosylated (N-glycosylated) chains attached to the extracellular domain influence cell proliferation. The current study investigated the function of syndecan-4 cytoplasmic domain in combination with GAG and N-glycosylated chains in turkey muscle cell proliferation, differentiation, fibroblast growth factor 2 (FGF2) responsiveness, and PKCα membrane localization. Syndecan-4 or syndecan-4 without the cytoplasmic domain and with or without the GAG and N-glycosylated chains were transfected or co-transfected with a small interfering RNA targeting syndecan-4 cytoplasmic domain into turkey muscle satellite cells. The overexpression of syndecan-4 mutants increased cell proliferation but did not change differentiation. Syndecan-4 mutants had increased cellular responsiveness to FGF2 during proliferation. Syndecan-4 increased PKCα cell membrane localization, whereas the syndecan-4 mutants decreased PKCα cell membrane localization compared to syndecan-4. However, compared to the cells without transfection, syndecan-4 mutants increased cell membrane localization of PKCα. These data indicated that the syndecan-4 cytoplasmic domain and the GAG and N-glycosylated chains are critical in syndecan-4 regulating satellite cell proliferation, responsiveness to FGF2, and PKCα cell membrane localization.

Dopamine and mesotocin neurotransmission during the transition from incubation to brooding in the turkey

We investigated the neuroendocrine changes involved in the transition from incubating eggs to brooding of the young in turkeys. Numbers of mesotocin (MT; the avian analog of mammalian oxytocin) immunoreactive (ir) neurons were higher in the nucleus paraventricularis magnocellularis (PVN) and nucleus supraopticus, pars ventralis (SOv) of late stage incubating hens compared to the layers. When incubating and laying hens were presented with poults, all incubating hens displayed brooding behavior. c-fos mRNA expression was found in several brain areas in brooding hens. The majority of c-fos mRNA expression by MT-ir neurons was observed in the PVN and SOv while the majority of c-fos mRNA expression in dopaminergic (DAergic) neurons was observed in the ventral part of the nucleus preopticus medialis (POM). Following intracerebroventricular injection of DA or oxytocin (OT) receptor antagonists, hens incubating eggs were introduced to poults. Over 80% of those injected with vehicle or the D1 DA receptor antagonist brooded poults, while over 80% of those receiving the D2 DA receptor antagonist or the OT receptor antagonist failed to brood the poults. The D2 DA/OT antagonist groups also displayed less c-fos mRNA in the dorsal part of POM and the medial part of the bed nucleus of the stria terminalis (BSTM) areas than did the D1 DA/vehicle groups. These data indicate that numerous brain areas are activated when incubating hens initially transition to poult brooding behavior. They also indicate that DAergic, through its D2 receptor, and MTergic systems may play a role in regulating brooding behaviors in birds.

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 SNP based linkage map of the turkey genome reveals multiple intrachromosomal rearrangements between the turkey and chicken genomes

BACKGROUND

The turkey (Meleagris gallopavo) is an important agricultural species that is the second largest contributor to the world's poultry meat production. The genomic resources of turkey provide turkey breeders with tools needed for the genetic improvement of commercial breeds of turkey for economically important traits. A linkage map of turkey is essential not only for the mapping of quantitative trait loci, but also as a framework to enable the assignment of sequence contigs to specific chromosomes. Comparative genomics with chicken provides insight into mechanisms of genome evolution and helps in identifying rare genomic events such as genomic rearrangements and duplications/deletions.

RESULTS

Eighteen full sib families, comprising 1008 (35 F1 and 973 F2) birds, were genotyped for 775 single nucleotide polymorphisms (SNPs). Of the 775 SNPs, 570 were informative and used to construct a linkage map in turkey. The final map contains 531 markers in 28 linkage groups. The total genetic distance covered by these linkage groups is 2,324 centimorgans (cM) with the largest linkage group (81 loci) measuring 326 cM. Average marker interval for all markers across the 28 linkage groups is 4.6 cM. Comparative mapping of turkey and chicken revealed two inter-, and 57 intrachromosomal rearrangements between these two species.

CONCLUSION

Our turkey genetic map of 531 markers reveals a genome length of 2,324 cM. Our linkage map provides an improvement of previously published maps because of the more even distribution of the markers and because the map is completely based on SNP markers enabling easier and faster genotyping assays than the microsatellitemarkers used in previous linkage maps. Turkey and chicken are shown to have a highly conserved genomic structure with a relatively low number of inter-, and intrachromosomal rearrangements.

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.

A mitogenomic perspective on the ancient, rapid radiation in the Galliformes with an emphasis on the Phasianidae

BACKGROUND

The Galliformes is a well-known and widely distributed Order in Aves. The phylogenetic relationships of galliform birds, especially the turkeys, grouse, chickens, quails, and pheasants, have been studied intensively, likely because of their close association with humans. Despite extensive studies, convergent morphological evolution and rapid radiation have resulted in conflicting hypotheses of phylogenetic relationships. Many internal nodes have remained ambiguous.

RESULTS

We analyzed the complete mitochondrial (mt) genomes from 34 galliform species, including 14 new mt genomes and 20 published mt genomes, and obtained a single, robust tree. Most of the internal branches were relatively short and the terminal branches long suggesting an ancient, rapid radiation. The Megapodiidae formed the sister group to all other galliforms, followed in sequence by the Cracidae, Odontophoridae and Numididae. The remaining clade included the Phasianidae, Tetraonidae and Meleagrididae. The genus Arborophila was the sister group of the remaining taxa followed by Polyplectron. This was followed by two major clades: ((((Gallus, Bambusicola) Francolinus) (Coturnix, Alectoris)) Pavo) and (((((((Chrysolophus, Phasianus) Lophura) Syrmaticus) Perdix) Pucrasia) (Meleagris, Bonasa)) ((Lophophorus, Tetraophasis) Tragopan))).

CONCLUSIONS

The traditional hypothesis of monophyletic lineages of pheasants, partridges, peafowls and tragopans was not supported in this study. Mitogenomic analyses recovered robust phylogenetic relationships and suggested that the Galliformes formed a model group for the study of morphological and behavioral evolution.

DNA fingerprints of poultry

Human minisatellite probes cross-hybridize to DNA of several species of poultry (chicken, duck, turkey and goose), and detect high levels of polymorphism. The resulting DNA fingerprints are individual specific, and allow the discrimination even between closely related birds. The pattern of poultry DNA fingerprints is different from that of humans and other animals, having a higher average proportion of large DNA fragments. Pedigree analysis revealed a low number of allelic pairs of variable DNA fragments, indicating that most of the alleles are unresolved in the DNA fingerprint or too small to be detected. The total number of detectable loci in broilers, using probe 33.6, was estimated as 62, of which 13 loci are on average scoreable and available for use. Poultry DNA fingerprints can be used for individual identification, linkage studies and as an aid in breeding programmes.

Inheritance of adenosine deaminase variants in chickens and turkeys

Blood cell lysates of chickens and turkeys were subjected to starch gel electrophoresis and the gels were stained for adenosine deaminase. Two zones were observed singly or together in the electrophoretic patterns of each lysate. Zones of chicken lysates were analogous in electrophoretic mobility to those of turkeys. An extra zone which appeared in patterns of a sample stored over one month was not detected in patterns of a second aliquot of stored sample treated with a reducing agent prior to electrophoresis. Family data involving 110 chicken progeny and 221 turkey progeny supported the hypothesis that these zones were controlled by two codominant alleles designated ADAA and ADAB. In the two Leghorn strains studied ADAB was much more frequent than ADAA, but the frequency distribution was reversed in the Small White turkey strain examined.

Development of reagents to study the turkey's immune response: Identification and molecular cloning of turkey CD4, CD8alpha and CD28

The cDNAs of three turkey CD markers, CD4, CD8alpha and CD28, were identified by screening a turkey cDNA library. The coding regions of the chicken and turkey genes are highly conserved, with 91.3-96.1% nucleotide (nt) and 84.2-95.5% amino acid (aa) identity. Identity was less conserved between avian CD markers and their mammalian homologues, ranging from 44.7 to 59.8% and 22.4 to 50.4% at the nt and aa levels, respectively. Anti-chicken CD8alpha and CD28 monoclonal antibodies were demonstrated to specifically cross-react with turkey CD8alpha and CD28, respectively.

The mitochondrial genome sequence and molecular phylogeny of the turkey, Meleagris gallopavo

The mitochondrial genome (mtGenome) has been little studied in the turkey (Meleagris gallopavo), a species for which there is no publicly available mtGenome sequence. Here, we used PCR-based methods with 19 pairs of primers designed from the chicken and other species to develop a complete turkey mtGenome sequence. The entire sequence (16,717 bp) of the turkey mtGenome was obtained, and it exhibited 85% similarity to the chicken mtGenome sequence. Thirteen genes and 24 RNAs (22 tRNAs and 2 rRNAs) were annotated. An mtGenome-based phylogenetic analysis indicated that the turkey is most closely related to the chicken, Gallus gallus, and quail, Corturnix japonica. Given the importance of the mtGenome, the present work adds to the growing genomic resources needed to define the genetic mechanisms that underlie some economically significant traits in the turkey.

Cardiovascular system diseases in turkeys

Genetic manipulation conducted for many years has resulted in breeding turkeys with very intensive growth rate and high percentage of muscle tissue. These features promote dysfunctions in the cardiovascular system. The aim of this article was to present current data on the etiology of cardiovascular system diseases in turkeys. In this paper the most recent data on dilated cardiomyopathy (round heart disease), spontaneous aortic rupture and perirenal hemorrhage syndrome are described.

The Effect of Genetic Increases in Egg Production and Age and Sex on Breast Muscle Development of Turkeys

Pectoralis major muscle morphology was studied in both sexes of a turkey line (E) selected long-term for increased egg production and its randombred control (RBC1) from 25 d of incubation through 20 wk posthatch. Pectoralis major muscle samples from 10 individuals from each line-sex-age subgroup were obtained in a manner to prevent contraction. The muscle samples were dehydrated, cleared, embedded in paraffin, sectioned, incubated, and rehydrated before staining with hematoxylin and eosin. Representative sections were given a score by 4 individuals based on breast muscle morphology. The scores ranged from 1 (little extracellular matrix and indistinct muscle fibers) to 5 (large extracellular space and distinct muscle fibers). Scores from 2 to 4 were intermediate to these extremes. The pectoralis major muscle morphology scores were highest at 25 d of incubation and declined greatly at 1 wk of age. The scores increased from 1 to 4 wk of age and remained constant through 20 wk of age. Males had higher scores than females. In the current study, there was no significant difference between the E and RBC1 lines. Based on the results of 3 experiments (2 published and the present one) using the E and RBC1 lines, it appears that genetic increases in egg production may be associated with a slight reduction in pectoralis major muscle morphology scores at 16 wk of age.

Dam line and sire line effects on turkey embryo survival and thyroid hormone concentrations at the plateau stage in oxygen consumption

Inheritance of embryo thyroid function was measured in lines of turkeys. Two lines that had been selected for either increased egg production (E) or increased 16-wk BW (F) and their respective randombred controls (i.e., RBC1 and RBC2) were examined. Reciprocal crosses of dams and sires from each selected line and its randombred control were made to estimate sire line and dam line effects. Orthogonal contrasts were used to determine if the differences found were due to the presence of additive, nonadditive, or maternal, sex-linked, or both, gene effects. With the data involved, sex-linkage and maternal effects could not be separated. Embryo survival was measured for all lines and their reciprocal crosses. Crossing the RBC1 sire and E dam also resulted in better embryo survival and lower death losses at pipping than for the other cross- or purelines. Reciprocal crosses of the F and RBC2 lines showed better total embryo survival, and they survived pipping better than the F or RBC2 purelines. Thyroxine (T(4)) and triiodothyronine (T(3)) concentrations differed between the reciprocal crosses at external pipping, but the effects were inconsistent for the 2 data sets. Reciprocal tests indicated that maternal, sex-linked, or both, effects were present for T(3) concentrations at internal pipping in the E and RBC1 lines and at external pipping for the F and RBC2 lines. Reciprocal effects were significant for T(4) at internal pipping for both data sets. The RBC1 sire embryos had significantly higher T(3):T(4) ratios than the E line sire embryos at internal and external pipping, and the pureline RBC1 embryos had consistently higher ratios than the pureline E embryos. The differences for the T(3):T(4) ratios between these 2 lines at internal pipping, external pipping, and hatch appeared to be consistently additive in nature, although significant nonadditive or heterotic effects were present for the ratio at external pipping. Similar effects on the T(3):T(4) ratio were observed for the F and RBC2 lines at external pipping.

Characterization of the turkey MHC chromosome through genetic and physical mapping

Previous studies in the chicken have identified a single microchromosome (GGA16) containing the ribosomal DNA (rDNA) and two genetically unlinked MHC regions, MHC-B and MHC-Y. Chicken DNA sequence from these loci was used to develop PCR primers for amplification of homologous fragments from the turkey (Meleagris gallopavo). PCR products were sequenced and overgo probes were designed to screen the CHORI 260 turkey BAC library. BAC clones corresponding to the turkey rDNA, MHC-B and MHC-Y were identified. BAC end and subclone sequencing confirmed identity and homology of the turkey BAC clones to the respective chicken loci. Based on subclone sequences, single-nucleotide polymorphisms (SNPs) segregating within the UMN/NTBF mapping population were identified and genotyped. Analysis of SNP genotypes found the B and Y to be genetically unlinked in the turkey. Silver staining of metaphase chromosomes identified a single pair of microchromosomes with nucleolar organizer regions (NORs). Physical locations of the rDNA and MHC loci were determined by fluorescence in situ hybridization (FISH) of the BAC clones to metaphase chromosomes. FISH clearly positioned the rDNA distal to the Y locus on the q-arm of the MHC chromosome and the MHC-B on the p-arm. An internal telomere array on the MHC chromosome separates the B and Y loci.

Tandem 41-bp repeats in chicken and Japanese quail genomes: FISH mapping and transcription analysis on lampbrush chromosomes

The chromosomal distribution of 41-bp repeats, known as CNM and PO41 repeats in the chicken genome and BglII repeats in the Japanese quail, was analyzed precisely using giant lampbrush chromosomes (LBC) from chicken, Japanese quail, and turkey growing oocytes. The PO41 repeat is conserved in all galliform species, whereas the other repeats are species specific. In chicken and quail, the centromere and subtelomere regions share homologous satellite sequences. RNA polymerase II transcribes the 41-bp repeats in both centromere and subtelomere regions. Ongoing transcription of these repeats was demonstrated by incorporation of BrUTP injected into oocytes at the lampbrush stage. RNA complementary to both strands of CNM and PO41 repeats is present on chicken LBC loops, whereas strand-specific G-rich transcripts are characteristic of BglII repeats in the Japanese quail. The RNA from 41-bp repeats does not undergo cotranscriptional U snRNP-dependent splicing. At the same time, the ribonucleoprotein matrix of transcription units with C-rich RNA of CNM and PO41 repeats was enriched with hnRNP protein K. Potential promoters for satellite transcription are discussed.

Differential Effects of Sex and Genetics on Behavior and Stress Response of Turkeys

Three lines of turkeys were tested for response in T-maze and open-field tests during the first 8 d after hatch, and behavior was observed after catching, moving, and transport. They were also compared for corticosterone (CORT) levels and heterophil:lymphocyte ratios (H:L) at 15 wk of age in response to an Escherichia coli challenge followed by transport stress. Large commercial-(COMM) line birds were faster and more active in the T-maze at d 2 than egg-line birds. Male COMM-line birds were faster than male egg-line birds when tested in an open field at d 8. Egg-line birds had more sleeping behavior after moving to a new floor pen as compared with both an intermediate-sized line (F line) and the COMM line. Transport stress increased CORT levels in all 3 lines, and the increase was greater in males compared with females. The egg line had higher basal CORT levels (P = 0.03) and higher levels after transport (P < 0.0001). The H:L ratios were affected by both transport stress and line but not by sex. The H:L ratio was lower in the egg line as compared with both the F line and the COMM line (P < 0.0001), with the COMM line having the greatest increase in response to transport. These data, combined with those from previous studies of these lines, suggest that differences in activity of fast-growing turkeys may be used to select birds that are less susceptible to inflammatory bacterial disease and that the H:L ratio may be more useful than serum CORT in evaluating the deleterious effects of stress.

The application of random regression models in the genetic analysis of monthly egg production in turkeys and a comparison with alternative longitudinal models

Random regression models (RR) have become a popular methodology for the genetic study of longitudinal data since the last decade. The first objective of the current study was to investigate the application of RR models for the genetic analysis of egg production in turkeys. Data collected from a heavy dam line were used to estimate genetic parameters with 2 RR models, one having second-order Legendre polynomials as regression over time (RR2) and another with third-order polynomials (RR3). The second objective was to benchmark the performance of RR models with more conventional methods, so genetic parameters were reestimated using a multitrait (MT) and a repeatability model. To assess the model efficiency of predicting missing values, a reduced data set was used, and for each model, the predicted values of the deleted records were compared with the true values. The RR models were further compared against each other by eliminating the last period and estimating the MS error of the predictions for both models. The repeatability model had the poorest performance in predicting missing values. Heritability estimates from RR2 and MT models were close, whereas the RR3 model estimates were different. Both RR models demonstrated better prediction ability than the MT model. However, when RR models were compared solely, the RR2 model resulted in the smallest MS error. The results indicated that the RR3 model overfitted the data, suggesting that the choice of the appropriate polynomial order requires careful consideration. The present study illustrated that the application of RR models for the genetic analysis of egg production in turkeys is not only feasible but also offers a high accuracy of prediction.

Ontogenesis of the Expression of Prolactin Receptor Messenger Ribonucleic Acid During Late Embryogenesis in Turkeys and Chickens

Changes in circulating levels of prolactin (PRL) and tissue content of PRL receptor (PRLR) messenger RNA (mRNA) in the liver, pancreas, kidney, and gonads (testis/ovary) were measured in turkey and chicken embryos from embryonic day (ED) 21 or ED15, respectively, to 1 d after hatch by real-time PCR. There were no differences between the sexes in chickens or turkeys. Both species had very similar patterns of PRL release and expression of PRLR mRNA, and no major differences were observed between turkey or chicken embryos. Plasma levels of PRL increased from low levels during the last week of embryonic development and were at significantly higher levels (about 4-fold) by 1 d after hatch. Similarly, in all tissues the content of PRLR mRNA was minimal at the outset and increased to reach maxima about the time of hatch. In both species, the highest levels of transcript were observed in the kidney followed by the gonad, liver, and pancreas. The tissue content of PRLR was correlated (0.6 to 0.8 dependent on the tissue) to circulating levels of PRL, which suggested that PRL may be associated with an increase in its receptor number around the time of hatch. Because levels of PRL and tissue content of PRLR mRNA increased around the time of hatch, this suggests that these tissues may be targets for PRL and may be involved in the physiologic changes occurring in embryos around the time of hatching.

DNA sequence and haplotype variation in two candidate genes for dilated cardiomyopathy in the turkey Meleagris gallopavo

Determining variation in genes is fundamental to understanding their function in the disease state. Cardiac troponin T (cTnT) and phospholamban (PLN) genes have been implicated in dilated cardiomyopathy (DCM) in human and model species. To investigate the role of these 2 candidate genes in DCM in the turkey Meleagris gallopavo, understanding sequence variants and map position distribution is necessary. To this end, a total of 1854 and 1771 bp of cTnT and PLN gene sequences, respectively, were scanned for single nucleotide polymorphisms (SNPs) in a randomly bred population. A total of 15 SNPs was identified in the cTnT and PLN genomic sequences. Nine haplotypes, 5 in cTnT and 4 in PLN, were identified. Observed heterozygosities (0.02–0.39) in the turkey population were low for both genes. Within each gene, 1 SNP corresponding to a restriction enzyme site was identified and used to develop a PCR–restriction fragment length polymorphism (RFLP) genotyping assay. The PLN gene was genetically mapped to turkey chromosome 2, equivalent to Gallus gallus chromosome 3, and cTnT mapped to a turkey microchromosome. Although limited because of the relatively small sample size of 55 birds, the data from this SNP analysis of PLN and cTnT provide a foundation from which to evaluate the function of cTnT and PLN in the turkey. Information about the distribution of the SNPs and haplotypes will facilitate future association and linkage studies.

Follicular Development and Expression of the Messenger Ribonucleic Acid for the Inhibin/Activin Subunits in Two Genetic Lines of Turkey Hens that Differ in Total Egg Production

The characterization of the follicular hierarchy and the expression of the mRNA for the inhibin/activin subunits was investigated in the follicles of 2 lines of turkey hens selected for over 40 generations for increased egg production (Egg line) or increased body weight (Growth line). The follicular hierarchies of 6 hens from the Egg and Growth lines were characterized in middle (45 wk of age) and late production (58 wk of age). Relative follicular weights for individual hierarchical follicles (>12 mm), pooled small yellow follicles (5 to 12 mm), and large white follicles (2 to 5 mm) were calculated. Total RNA was extracted for Northern blot analysis from individual granulosa cell layers of the F1 through F4 follicles, and from the combined granulosa and theca layers of small yellow follicles and large white follicles from an additional 6 hens from each genetic line. Egg line hens displayed a more distinct follicular size hierarchy than Growth line hens at 45 and 58 wk. Although total follicular weight relative to body size was greater at 45 and 58 wk of age for the Egg line hens than the Growth line hens, the total number of hierarchical follicles was greater in the Growth line hens at 45 and 58 wk of age. Expression of follistatin and the inhibin beta(B)-subunit was highest in nonhierarchical follicles, whereas the expression of the inhibin alpha- and beta(A)-subunits was highest in the hierarchical follicles. The inhibin alpha- and beta(A)-subunit mRNA expression pattern in the 4 largest follicles of the Growth line hens was not similar to the Egg line hens or characteristic of laying hens that have a high rate of egg production. The unusual inhibin subunit mRNA expression in the largest hierarchical follicles of the Growth line hens may account for their development of an abnormal follicular size hierarchy and for their poor egg production.

Molecular cloning and sequencing of MHC class II beta 1 domain of turkey reveals high sequence identity with chicken

We report the nucleotide sequences of turkey (Meleagris gallopavo) major histocompatibility complex (MHC) class II loci (beta 1 domain or exon 2 encoding the peptide-binding region). In the present investigation, three distinct sequences from the beta 1 domain of turkey MHC class II were isolated. A BLAST search and phylogenetic analysis revealed that turkey MHC sequences are most similar to chicken and peacock MHC. There was no strong evidence of recombination among the turkey MHC sequences or with other avian MHC, but diversity was high. The diversity in this peptide-binding region may be the result of point mutation and balancing selection or frequent gene conversion within turkey. However, more work and data are needed to understand the evolution of turkey and other avian MHC. Moreover, polymerase chain reaction-restriction fragment-length polymorphism analysis of exon 2 using the Hinf I restriction enzyme demonstrated three restriction patterns and a preliminary evidence of multiple beta loci in turkey. PCR-RFLP analysis of turkey MHC class II loci could be a promising method of MHC genotyping, when more sequences are available. Turkey MHC haplotypes identified earlier by RFLP analysis should be sequenced to standardize turkey MHC nomenclature and to develop DNA based method of haplotyping.

Comparison of the performance of 1966- versus 2003-type turkeys when fed representative 1966 and 2003 turkey diets: growth rate, livability, and feed conversion

Body weight, livability, and feed conversion of a randombred control turkey line (RBC2) started in 1966 at The Ohio State University was compared with that of modern commercial turkeys hatched in 2003 when fed representative 1966- and 2003-type diets from hatch (March 5, 2003) through 196 d of age. Each pen of modern turkeys consisted of 5 birds each of the Nicholas, British United Turkeys of America, and Hybrid strains. Eight groups (i.e., 2 strains (RBC2 vs. modern), 2 sexes, and 2 dietary regimens) were randomly assigned into each of 4 blocks of 8 litter floor pens (32 total) for growout. Using the BW performance of the 2 strains on the modern feed as the basis, the study showed that the 2003 turkeys were approximately twice as heavy as the 1966 RBC2 at the 4 slaughter ages and that tom weights have increased by 186, 208, 227, and 241 g/yr, and hen weights have increased by 164, 179, 186, and 205 g/yr at 112, 140, 168, and 196 d of age, respectively, over the past 37 yr. Cumulative feed conversion (kg of feed/kg of BW) was approximately 20% better in the 2003 tom turkey on the 2003 feed (2.638) than in the RBC2 tom on the 1966 feed (3.278) at 20 wk of age. Feed efficiency to 11 kg of BW in the 2003 toms (2.132 at 98 d of age) was approximately 50% better than in the RBC2 toms (4.208 at 196 d of age). The number of days to reach that weight was halved during this period of time. Growth performance during the different periods of the study appeared to be strongly affected by type of feed used and seasonal changes in ambient temperature. Overall livability was very good for all groups, but the mortality level of the RBC2 was consistently higher, although not significantly so, than for the modern birds.