Litter Commensal Bacteria Can Limit the Horizontal Gene Transfer of Antimicrobial Resistance to Salmonella in Chickens.

Fostering a "balanced" gut microbiome through the administration of beneficial microbes that can competitively exclude pathogens has gained a lot of attention and use in human and animal medicine. However, little is known about how microbes affect the horizontal gene transfer of antimicrobial resistance (AMR). To shed more light on this question, we challenged neonatal broiler chicks raised on reused broiler chicken litter-a complex environment made up of decomposing pine shavings, feces, uric acid, feathers, and feed-with Salmonella enterica serovar Heidelberg (S. Heidelberg), a model pathogen. Neonatal chicks challenged with S. Heidelberg and raised on reused litter were more resistant to S. Heidelberg cecal colonization than chicks grown on fresh litter. Furthermore, chicks grown on reused litter were at a lower risk of colonization with S. Heidelberg strains that encoded AMR on IncI1 plasmids. We used 16S rRNA gene sequencing and shotgun metagenomics to show that the major difference between chicks grown on fresh litter and those grown on reused litter was the microbiome harbored in the litter and ceca. The microbiome of reused litter samples was more uniform and enriched in functional pathways related to the biosynthesis of organic and antimicrobial molecules than that in fresh litter samples. We found that Escherichia coli was the main reservoir of plasmids encoding AMR and that the IncI1 plasmid was maintained at a significantly lower copy per cell in reused litter compared to fresh litter. These findings support the notion that commensal bacteria play an integral role in the horizontal transfer of plasmids encoding AMR to pathogens like Salmonella. IMPORTANCE Antimicrobial resistance spread is a worldwide health challenge, stemming in large part from the ability of microorganisms to share their genetic material through horizontal gene transfer. To address this issue, many countries and international organizations have adopted a One Health approach to curtail the proliferation of antimicrobial-resistant bacteria. This includes the removal and reduction of antibiotics used in food animal production and the development of alternatives to antibiotics. However, there is still a significant knowledge gap in our understanding of how resistance spreads in the absence of antibiotic selection and the role commensal bacteria play in reducing antibiotic resistance transfer. In this study, we show that commensal bacteria play a key role in reducing the horizontal gene transfer of antibiotic resistance to Salmonella, provide the identity of the bacterial species that potentially perform this function in broiler chickens, and also postulate the mechanism involved.

Grid search approach to discriminate between old and recent inbreeding using phenotypic, pedigree and genomic information.

BACKGROUND: Although inbreeding caused by the mating of animals related through a recent common ancestor is expected to have more harmful effects on phenotypes than ancient inbreeding (old inbreeding), estimating these effects requires a clear definition of recent (new) and ancient (old) inbreeding. Several methods have been proposed to classify inbreeding using pedigree and genomic data. Unfortunately, these methods are largely based on heuristic criteria such as the number of generations from a common ancestor or length of runs of homozygosity (ROH) segments. To mitigate these deficiencies, this study aimed to develop a method to classify pedigree and genomic inbreeding into recent and ancient classes based on a grid search algorithm driven by the assumption that new inbreeding tends to have a more pronounced detrimental effect on traits. The proposed method was tested using a cattle population characterized by a deep pedigree. RESULTS: Effects of recent and ancient inbreeding were assessed on four growth traits (birth, weaning and yearling weights and average daily gain). Thresholds to classify inbreeding into recent and ancient classes were trait-specific and varied across traits and sources of information. Using pedigree information, inbreeding generated in the last 10 to 11 generations was considered as recent. When genomic information (ROH) was used, thresholds ranged between four to seven generations, indicating, in part, the ability of ROH segments to characterize the harmful effects of inbreeding in shorter periods of time. Nevertheless, using the proposed classification method, the discrimination between new and old inbreeding was less robust when ROH segments were used compared to pedigree. Using several model comparison criteria, the proposed approach was generally better than existing methods. Recent inbreeding appeared to be more harmful across the growth traits analyzed. However, both new and old inbreeding were found to be associated with decreased yearling weight and average daily gain. CONCLUSIONS: The proposed method provided a more objective quantitative approach for the classification of inbreeding. The proposed method detected a clear divergence in the effects of old and recent inbreeding using pedigree data and it was superior to existing methods for all analyzed traits. Using ROH data, the discrimination between old and recent inbreeding was less clear and the proposed method was superior to existing approaches for two out of the four analyzed traits. Deleterious effects of recent inbreeding were detected sooner (fewer generations) using genomic information than pedigree. Difference in the results using genomic and pedigree information could be due to the dissimilarity in the number of generations to a common ancestor. Additionally, the uncertainty associated with the identification of ROH segments and associated inbreeding could have an effect on the results. Potential biases in the estimation of inbreeding effects may occur when new and old inbreeding are discriminated based on arbitrary thresholds. To minimize the impact of inbreeding, mating designs should take the different inbreeding origins into consideration.

Alteration of dietary cysteine affects activities of genes of the transsulfuration and glutathione pathways, and development of skin tissues and feather follicles in chickens.

The dietary requirement for cysteine is not determined in poultry since it is not an essential amino acid. The cysteine need is expected to be met through the transsulfuration pathway where homocysteine, a precursor of methionine, is converted to cysteine. Cysteine is a major component of plumage, and the degree to which cysteine is involved in plumage and other keratized proteins are unknown. We randomly assigned chicks to control and treatment (deficient in cysteine) diets for 49 d. The thickness of the skin layers, feather follicle length, and thickness were measured at days 10, 24, 34, and 49. We also measured the hepatic mRNA expressions of cystathionine beta synthase (CBS), cystathionine γ-lyase (CTL), cysteine dioxygenase (CDO), and glutathione synthetase (GSS). Chickens fed the treatment diet had reduced epidermis thickness and shorter feather follicles compared with the controls. The chicken fed the treatment diet also had increased mRNA expression of CBS and CTL indicating a disruption of the transsulfuration pathway. The treatment chickens also had a decreased hepatic CDO and increased GSS mRNA expressions which are in concordance with the homeostatic regulation of cysteine. Compromised cysteine metabolism could affect thermoregulation and subsequently affect feed efficiency and welfare of the birds.

A Machine Vision-Based Method for Monitoring Broiler Chicken Floor Distribution.

The proper spatial distribution of chickens is an indication of a healthy flock. Routine inspections of broiler chicken floor distribution are done manually in commercial grow-out houses every day, which is labor intensive and time consuming. This task requires an efficient and automatic system that can monitor the chicken's floor distributions. In the current study, a machine vision-based method was developed and tested in an experimental broiler house. For the new method to recognize bird distribution in the images, the pen floor was virtually defined/divided into drinking, feeding, and rest/exercise zones. As broiler chickens grew, the images collected each day were analyzed separately to avoid biases caused by changes of body weight/size over time. About 7000 chicken areas/profiles were extracted from images collected from 18 to 35 days of age to build a BP neural network model for floor distribution analysis, and another 200 images were used to validate the model. The results showed that the identification accuracies of bird distribution in the drinking and feeding zones were 0.9419 and 0.9544, respectively. The correlation coefficient (R), mean square error (MSE), and mean absolute error (MAE) of the BP model were 0.996, 0.038, and 0.178, respectively, in our analysis of broiler distribution. Missed detections were mainly caused by interference with the equipment (e.g., the feeder hanging chain and water line); studies are ongoing to address these issues. This study provides the basis for devising a real-time evaluation tool to detect broiler chicken floor distribution and behavior in commercial facilities.

Transcriptional profiling and pathway analysis reveal differences in pituitary gland function, morphology, and vascularization in chickens genetically selected for high or low body weight.

BACKGROUND: Though intensive genetic selection has led to extraordinary advances in growth rate and feed efficiency in production of meat-type chickens, endocrine processes controlling these traits are still poorly understood. The anterior pituitary gland is a central component of the neuroendocrine system and plays a key role in regulating important physiological processes that directly impact broiler production efficiency, though how differences in pituitary gland function contribute to various growth and body composition phenotypes is not fully understood. RESULTS: Global anterior pituitary gene expression was evaluated on post-hatch weeks 1, 3, 5, and 7 in male broiler chickens selected for high (HG) or low (LG) growth. Differentially expressed genes (DEGs) were analyzed with gene ontology categorization, self-organizing maps, gene interaction network determination, and upstream regulator identification to uncover novel pituitary genes and pathways contributing to differences in growth and body composition. A total of 263 genes were differentially expressed between HG and LG anterior pituitary glands (P ≤ 0.05 for genetic line-by-age interaction or main effect of line; ≥1.6-fold difference between lines), including genes encoding four anterior pituitary hormones. Genes involved in signal transduction, transcriptional regulation, and vesicle-mediated transport were differentially expressed and are predicted to influence expression and secretion of pituitary hormones. DEGs involved in immune regulation provide evidence that inflammation and response to cellular stressors may compromise pituitary function in LG birds, affecting their ability to adequately produce pituitary hormones. Many DEGs were also predicted to function in processes that regulate organ morphology and angiogenesis, suggesting pituitary gland structure differs between the divergently selected lines. CONCLUSIONS: The large number of DEGs within the anterior pituitary gland of birds selected for high or low body weight highlights the importance of this gland in regulating economically important traits such as growth and body composition in broiler chickens. Intracellular signaling, transcriptional regulation, and membrane trafficking are important cellular processes contributing to proper hormone production and secretion. The data also suggest that pituitary function is intimately tied to structure, and organization of the gland could influence hypothalamic and systemic metabolic inputs and delivery of hormones regulating growth and metabolism into peripheral circulation.

Increasing accuracy of genomic selection in presence of high density marker panels through the prioritization of relevant polymorphisms.

BACKGROUND: It becomes clear that the increase in the density of marker panels and even the use of sequence data didn't result in any meaningful increase in the accuracy of genomic selection (GS) using either regression (RM) or variance component (VC) approaches. This is in part due to the limitations of current methods. Association model are well over-parameterized and suffer from severe co-linearity and lack of statistical power. Even when the variant effects are not directly estimated using VC based approaches, the genomic relationships didn't improve after the marker density exceeded a certain threshold. SNP prioritization-based fixation index (FST) scores were used to track the majority of significant QTL and to reduce the dimensionality of the association model. RESULTS: Two populations with average LD between adjacent markers of 0.3 (P1) and 0.7 (P2) were simulated. In both populations, the genomic data consisted of 400 K SNP markers distributed on 10 chromosomes. The density of simulated genomic data mimics roughly 1.2 million SNP markers in the bovine genome. The genomic relationship matrix (G) was calculated for each set of selected SNPs based on their FST score and similar numbers of SNPs were selected randomly for comparison. Using all 400 K SNPs, 46% of the off-diagonal elements (OD) were between - 0.01 and 0.01. The same portion was 31, 23 and 16% when 80 K, 40 K and 20 K SNPs were selected based on FST scores. For randomly selected 20 K SNP subsets, around 33% of the OD fell within the same range. Genomic similarity computed using SNPs selected based on FST scores was always higher than using the same number of SNPs selected randomly. Maximum accuracies of 0.741 and 0.828 were achieved when 20 and 10 K SNPs were selected based on FST scores in P1 and P2, respectively. CONCLUSIONS: Genomic similarity could be maximized by the decrease in the number of selected SNPs, but it also leads to a decrease in the percentage of genetic variation explained by the selected markers. Finding the balance between these two parameters could optimize the accuracy of GS in the presence of high density marker panels.

Inhibition of the Transsulfuration Pathway Affects Growth and Feather Follicle Development in Meat-Type Chickens.

Cysteine is a nonessential amino acid in poultry nutrition. Poultry diets are deficient in cysteine, but the bird's cysteine need is met through the transsulfuration pathway (TSP) where homocysteine is converted to cysteine: a process catalyzed by cystathionine ß-synthase (CBS) and cystathionine γ-lyase (CTH). Cysteine is also a major component of keratinized protein found in feathers, but the extent to which cysteine is involved in feather and skin development in poultry is unknown. We randomly assigned chicks to control and treatment (control diet plus 100 mg/kg body weight of propargylglycine which is an inhibitor of CTH) diets. The thickness of skin layers, primary feather follicle parameters, growth, and mRNA expression of CBS and CTH were measured. Inhibition of TSP corresponded with the upregulation of liver mRNA of both CBS and CTH and reduction in growth from 35 to 40 days of age. The epidermis thickness, feather follicle length, and diameter were reduced from 10 to 40 days of age. Incorporation of cysteine into keratinized protein may be more sensitive to the level of available cysteine than into nonkeratinized proteins. Thus, disruption of the TSP could affect the thermoregulatory ability of the bird.

Cellular antioxidant enzyme activity and biomarkers for oxidative stress are affected by heat stress.

Heat stress (HS) causes oxidative stress and cellular changes in an attempt to detoxify the harmful effects of reactive oxygen species (ROS). However, how ROS affect different organs in chickens under acute and chronic HS is relatively unknown. We investigated the cellular enzyme activity and biomarker changes in the liver and Pectoralis (P) major muscle in broiler chickens subjected to both acute and chronic HS. Forty-eight broiler chickens at 14 days old were randomly assigned to either 25 °C (control) or 35 °C (heat-stressed) for 12 days. Five birds per treatment at 1 and 12 days post-HS were euthanized, and the liver and P. major muscle were sampled. Superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GPx), glutathione (GSH), glutathione reductase (GR), glutathione S-transferase (GST) activity as well as 8-hydroxy-2'-deoxyguanosine (8-OHdG), advanced glycation end product (AGE), malondialdehyde (MDA), and protein carbonyl (PCO) were analyzed as biomarkers for DNA, carbohydrate, lipid, and protein oxidation, respectively. The SOD, CAT, and GSH-GPx activity levels in the liver and the P. major muscle changed under HS; however, some of the changes were tissue-specific or dependent on the duration of the HS. There were increased liver 8-OHdG during chronic HS and also increased liver AGE levels during both acute and chronic HS indicating significant carbohydrate and DNA oxidations. In the P. major muscle, we observed significant increases in lipid peroxidation and protein oxidation which may reflect that this tissue is less resilient to oxidative damage under heat stress. We show that heat stress caused tissue-specific changes to levels of oxidation biomarkers in chicken.

Expression of genes that encode cellular oxidant/antioxidant systems are affected by heat stress.

Heat stress causes critical molecular dysfunction that affects productivity in chickens. Thus, the purpose of this study was to evaluate the effect of heat stress (HS) on the expression of select genes in the oxidation/antioxidation machinery in the liver of chickens. Chickens at 14 days of age were randomly assigned to two treatment groups and kept under either a constant normal temperature (25 °C) or high temperature (35 °C) in individual cages for 12 days. mRNA expression of Nrf2, oxidants NADPH(NOX): [NOX1, NOX2, NOX3, NOX4, NOX5 and DUOX2], and antioxidants [SOD1, CAT, GR, GPx1, NQO1] in the liver were analyzed at 1 and 12 days post-HS. We show that, HS changes the mRNA expression of oxidants thereby increasing cellular reactive oxygen species (ROS). Additionally, persistent HS up-regulates SOD which converts superoxides to hydrogen peroxide. We further demonstrated the dynamic relationship between catalase, GSH peroxidase (GPx) and NADPH under both acute and chronic heat stress. The pentose phosphate pathway could be important under HS since it generates NADPH which serves as a cofactor for GPx. Also, methionine, a precursor of cysteine has been shown to have reducing properties and thereby makes for an alternative fuel for redox processes. Genes in the ROS and antioxidant generation pathways may provide insight into nutritional intervention strategies, especially the use of methionine and/or cysteine when birds are suffering from heat stress.

Effect of heat stress on protein utilization and nutrient transporters in meat-type chickens.

The aim of this study was to investigate the effect of heat stress (HS) on digestibility of protein and fat and the expression of nutrient transporters in broilers. Forty-eight male Cobb500 chicks were used in this study. At day 14, birds were randomly divided into two groups and kept under either constant normal temperature (25 °C) or high temperature (35 °C) in individual cages. Five birds per treatment at 1 and 12 days post-treatment were euthanized, and Pectoralis major (P. major) and ileum were sampled for gene expression analysis. At day 33, ileal contents were collected and used for digestibility analysis. The total consumption and retention of protein and fat were significantly lower in the HS group compared to the control group. Meanwhile, the retention of crude protein per BWG was significantly higher in the HS group compared to the control group. In P. major and ileum tissues at day 1, transporters FATP1 and SGLT1 were down-regulated in the HS group. Meanwhile, FABP1 and PepT1 were down-regulated only in the ileum of the HS group. The converse was shown in P. major. The nutrient transporter FABP1 at day 12 post-HS was down-regulated in the P. major and ileum, but GLUT1 and PepT2 were down-regulated only in the ileum, and PepT1 was down-regulated only in the P. major compared with the control group. These changes in nutrient transporters suggest that high ambient temperature might change the ileum and P. major lipids, glucose, and oligopeptide transporters.

Transcriptomic analysis to elucidate the molecular mechanisms that underlie feed efficiency in meat-type chickens.

Feed efficiency phenotypes defined by genotypes or gene markers are unknown. To date, there are only limited studies on global gene expression profiling on feed efficiency. The objective of this study was to identify genes and pathways associated with residual feed intake (RFI) through transcriptional profiling of duodenum at two different ages in a chicken population divergently selected for low (LRFI) or high (HRFI) RFI. The global gene expression differences in LRFI and HRFI were assessed by the Affymetrix GeneChip(®) Chicken Genome Array and RT-PCR using duodenal tissue on days 35 and 42. The Ingenuity Pathway Analysis program was used to identify canonical and gene network pathways associated with RFI. A global view of gene expression differences between LRFI and HRFI suggest that RFI can be explained by differences in cell division, growth, proliferation and apoptosis, protein synthesis, lipid metabolism, and molecular transport of cellular molecules. Chickens selected for improved RFI achieve efficiency by reducing feed intake with a nominal or no change in weight gain by either up-regulating CD36, PPARα, HMGCS2, GCG or down-regulating PCSK2, CALB1, SAT1, and SGK1 genes within the lipid metabolism, small molecule biochemistry, molecular transport, cell death, and protein synthesis molecular and cellular functions. Chickens selected for reduced RFI via reduced feed intake with no change in weight gain achieve feed efficiency for growth by the up-regulation of genes that reduce appetite with increased cellular oxidative stress, prolonged cell cycle, DNA damage, and apoptosis in addition to increased oxidation of dietary fat and efficient fatty acids transported from the intestines.

Phenotypic characterization and analysis of genetic diversity between commercial crossbred and indigenous chickens from three different agro-ecological zones using DArT-Seq technology.

Indigenous and were used to study genetic diversity and population structure analyses. Polymorphism information content (PIC) values ranged from 0.0 to 0.5, with 21,285 SNP markers (35%) being in the lowest PIC value range (0 to 0.15) while 13,511 (commercial chickens have developed unique adaptations to their environments, which may include nutrition, pathogens, and thermal stress. Besides, environmental pressures and artificial selection have generated significant genome-wide divergence in chickens, as those selection pressures contribute a considerable evolutionary force to phenotypic and genotypic differentiation. Herein, we determined genomic diversity of indigenous chickens from semi-deciduous rainforest (SDR), coastal savannah (CS) and Guinea savannah (GS) agro-ecological zones (AEZs) in Ghana and commercial crossbreds (CC) reared at the Kwame Nkrumah University of Science and Technology (KNUST). We generated SNP markers from 82 chickens (62 indigenous chicken ecotypes and 26 commercial crossbred ecotype) using DArT-Seq technology. A total of 85,396 SNP markers were generated and after filtering the data, 58,353 markers 21%) were in the highest PIC value range (0.45 to 0.50). The CC were more genetically diverse than the indigenous birds, with the highest expected heterozygosity value of 0.220. Between the commercial crossbreds population and the indigenous ecotypes, pairwise FST values were estimated to be 0.105 between CS, 0.096 between SDF, and 0.133 between GS. Furthermore, PCA analysis showed that the CC, SDF and GS chickens clustered together and are genetically distant from the commercial crossbred. We herein show that chickens from the AEZs studied can be considered as one population. However, due the abundance of agro-byproducts in the SDR compared to the CS and GS, chickens from the SDR AEZ had better growth compared to their counterparts. It is suggested that the genetic diversity within the local ecotypes could form the basis for genetic improvement.

Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens.

Eimeria (E.) maxima parasite infects chickens' midgut disrupting the jejunal and ileal mucosa causing high morbidity and mortality. Heat stress (HS) is a seasonal stressor that impacts biological functions leading to poor performance. This study elucidates how HS, E. maxima infection, and their combination affect the ileum transcriptome. Two-hundred and forty 2-week-old males Ross708 chickens were randomly allocated into four treatment groups: thermoneutral-control (TNc), thermoneutral-infected (TNi), heat-stress control (HSc), and heat stress-infected (HSi), with 6 replicates each of 10 birds. Infected groups received 200x103 sporulated E. maxima oocysts/bird, and heat-treated groups were raised at 35°C. At 6-day post-treatment, ileums of five randomly selected chickens per group were sampled, RNA was extracted and sequenced. A total of 413, 3377, 1908, and 2304 DEGs were identified when applying the comparisons: TNc vs HSc, TNc vs TNi, HSi vs HSc, and TNi vs HSi, respectively, at cutoff ≥1.2-fold change (FDR: q<0.05). HSc vs TNc showed upregulation of lipid metabolic pathways and degradation/metabolism of multiple amino acids; and downregulation of most immune-related and protein synthesis pathways. TNc vs TNi displayed upregulation of most of immune-associated pathways and eukaryotic mRNA maturation pathways; and downregulation of fatty acid metabolism and multiple amino acid metabolism pathways including tryptophan. Comparing HSi versus HSc and TNi revealed that combining the two stressors restored the expression of some cellular functions, e.g., oxidative phosphorylation and protein synthesis; and downregulate immune response pathways associated with E. maxima infection. During E. maxima infection under HS the calcium signaling pathway was downregulated, including genes responsible for increasing the cytoplasmic calcium concentration; and tryptophan metabolism was upregulated, including genes that contribute to catabolizing tryptophan through serotonin and indole pathways; which might result in reducing the cytoplasmic pool of nutrients and calcium available for the parasite to scavenge and consequently might affect the parasite's reproductive ability.

Dissection of Koch's residual feed intake: implications for selection.

For 50 yr, residual feed intake (RFI) has remained a black box even though many researchers have touted it as a more biological estimate of efficiency of feed utilization than feed conversion ratio (FCR). We successfully dissected the efficiency of feed utilization by decomposing the components of RFI and ascertained the contributions of its components. Currently, a fixed effect model is used to predict RFI, which we term RFIF. We used a random effect model to predict RFIR, which allowed a separate estimation of RFI for maintenance (RFIM) and for growth (RFIG) and also ascertained their respective efficiencies. Judged by residual variance, R(2) and deviance information criterion, the random effect model was superior to the traditional fixed effect model used to generate RFIF. Under the traditional method, the h(2) of RFIF was 0.13 but h(2) of RFIR was 0.35. The heritability of RFIM and RFIG were moderate (~0.50), but the genetic correlation between them was highly negative (-0.95), suggesting that these 2 efficiencies contribute in an opposing way toward RFI. As a result, there should be caution in ascribing a biological basis to RFI. Under the current methodology, a biological basis can be ascribed to RFIM and RFIG. Selecting on RFIM will lead to smaller but efficient birds. The genetic gains in feed efficiency will be achieved by reductions in feed required for maintenance. The RFIG is not an efficiency parameter and should not be used as a sole criterion for selection. The ability of the current method to estimate efficiency values for metabolic BW and BW gain provides geneticists with additional parameters to use to discriminate between animals with similar RFIR. It also provides the flexibility to impose weights on RFIM and RFIG to meet a desired objective.

Symposium: experimental design for poultry production and genomics research.

This symposium dealt with the theoretical and practical aspects of choosing and evaluating experimental designs, and how experimental results may be related to poultry production through modeling. Additionally, recent advances in techniques for generating high-throughput genomic sequencing data, genomic breeding values, genomics selection, and genome-wide association studies have provided unique computational challenges to the poultry industry. Such challenges were presented and discussed.

Experiences with a single-step genome evaluation.

Genomic selection can be implemented based on the genomic relationship matrix (GBLUP) and can be combined with phenotypes from nongenotyped animals through the use of best linear unbiased prediction (BLUP). A common method to combine both sources of information involves multiple steps, but is difficult to use with complicated models and is nonoptimal. A simpler method, termed single-step GBLUP, or ssGBLUP, integrates the genomically derived relationships (G) with population-based pedigree relationships (A) into a combined relationship matrix (H) and allows for genomic selection in a single step. The ssGBLUP method is easy to implement and uses standard BLUP-based programs. Experiences with field data in chickens, pigs, and dairy indicate that ssGBLUP is more accurate yet much simpler than multi-step methods. The current limits of ssGBLUP are approximately 100,000 genotypes and 18 traits. Models involving 10 million animals have been run successfully. The inverse of H can also be used in existing programs for parameter estimationm, but a properly scaled G is needed for unbiased estimation. Also, as genomic predictions can be converted to SNP effects, ssGBLUP is useful for genomic-wide association studies. The single-step method for genomic selection translates the use of genomic information into standard BLUP, and variance-component estimation programs become a routine.

Misclassification in binary responses and effect on genome-wide association studies.

Misclassification of dependent variables is a major issue in many areas of science that can arise when indirect markers are used to classify subjects or continuous traits are treated as categorical. In human medicine, this can have significant impacts on diagnostic accuracy. In animal science applications, misclassification can negatively affect both the accuracy of selection and the ability to ascertain the biological mechanisms for traits of interest. When dealing with traits influenced by genetic factors, genomic markers, such as SNP, can provide direct measurements of the underlying mechanisms controlling phenotypic responses. Unfortunately, in the presence of misclassification in the discrete dependent variables, the robustness of the analysis and the validity of the results could be severely compromised. To quantify the impact of misclassification on genome-wide association studies for binary responses, a real databased simulation was carried out. The simulated data consisted of 2,400 animals genotyped for 50K SNP. A binary trait with heritability equal to 0.10 and prevalence of 20% was generated. A rate of 1, 5, and 10% misclassification was artificially introduced to the true binary responses. Using a latent-threshold model, 3 analyses were carried out for each misclassification rate using 1) the true data (M1), 2) the contaminated data and ignoring misclassification (M2), and 3) the contaminated data and accounting for misclassification (M3). The results indicate that ignoring misclassification, when it exists in the data such as in M2, will lead to major deterioration in the performance of the model. When misclassification was contemplated in the model (M2), the results indicated a strong capacity of the procedure in dealing with potential misclassification in the training set. In fact, a large portion of miscoded samples in the training set was identified and corrected. The results of this study suggest that the proposed method is adequate and effective for practical genome-wide association studies for binary response classification.

Genetic variability characterization of the Moroccan Houbara Bustard (Chlamydotis undulata undulata) inferred from pedigree analysis.

A Moroccan Houbara Bustard pedigree was analyzed to evaluate the genetic variability in captive breeding population using genealogical approaches. The whole Houbara breeding flock (WP) for the period 1993-2004 was made up of 531 birds comprising 346 females and 185 males. The reference population (RP) comprised 198 individuals ready for reproduction from 2000 to 2004 cohorts. The corresponding percentage of known ancestors was estimated as 98.23% for the parent generation, 41.19% for the grandparent generation and 7.00% for the great grandparents generation. The average generation interval for Houbara was computed as 4.64 years. Genetic variability loss per generation was ascertained using the effective population size (Ne), the founder genome equivalent (fge), the effective number of ancestors and founders (fa) and (fe), respectively, for the RP and across each cohort. The results showed no bottleneck events in the breed but some loss of genetic variability just after the initiation of the conservation program. However, the annual effective population size based on the realized increase in inbreeding (ΔF) was estimated to be 207 for the RP and 1,000 for the WP. With regard to conservation breeding schemes, the genealogical evidence presented here is very useful as it revealed the positive effect of migration on Houbara breeding. The mating strategies will assist in the future control and management of the genetic variability of this population.

Research Note: Analysis of body weight and walking ability in turkeys and the prediction of categorical responses across systematic effect classes using a linear threshold model.

Heavy selection for growth in turkeys has led to a decay in leg soundness and walking ability. In this study, different models and traits were used to investigate the genetic relationships between body weight (BW) and walking ability (WA) in a turkey population. The data consisted of BW and WA traits collected on 276,059 male birds. Body weight was measured at 12 and 20 wk and WA at 20 wk of age. For WA, birds were scored based on a 1 (bad) to 6 (good) grading system. Due to the small number of records with scores 5 and 6, birds with WA scores of 4, 5, and 6 were grouped together resulting in only 4 classes. Additionally, a binary classification of WA (scores 1 and 2 = Similarly, an estimate of the genetic correlation between WA and BW at 20 wk was -0.45, indicating a more pronounced class 1; scores 3, 4, 5, and 6 = class 2) was evaluated. The inheritability estimates of WA ranged between 0.25 and 0.27 depending on the number of classes. The Heritability of BW at 12 and 20 wk was 0.44 and 0.51, respectively. The genetic correlation between WA and BW at 12 wk was around -0.35, indicating that heavy birds tend to have poor WA. antagonistic relationship between BW and WA. The genetic correlation between BW at 12 and 20 wk was positive and high (0.80). The residual correlation between WA and BW at 12 and 20 wk of age was -0.07 and -0.02, respectively. The residual correlation between body weight traits was 0.57. Similar results were observed when a binary classification was adopted for WA. The probability of an individual with a given genetic merit expressing a certain class of WA was determined for different fixed effect designations. Predictive probabilities clearly showed that birds when hatched in the winter would have a small chance to exhibit good WA phenotypes.

Methionine sources at different dietary levels alters the growth and expression of genes related to homocysteine remethylation in the jejunum of broilers.

Sulfur amino acids are essential for the proper development of broilers and are required throughout the bird's life to perform important physiological functions. Studies that seek to understand the actions of sulfur amino acids in the body of birds are essential. The present study evaluated the influence of sulfur amino acid supplementation using DL-Methionine (DL-Met) and DL-Methionine hydroxy analogue (DL-HMTBA), on the performance and expression of genes related to methionine metabolism, in the jejunum of broilers. Four hundred and fifty male broilers (Cobb-700 slow feathering) were distributed in a completely randomized design, in a factorial scheme (2x3), with two sources of methionine (DL-Met and DL-HMTBA) and three levels of methionine (deficiency, requirement and excess). The mRNA expression of the MAT1, MTR, BHMT, MTRR, CBG and GSS genes, and performance data such as feed intake, weight gain, and feed conversion were evaluated. DL-HMTBA increased the expression of BHMT (p = 0.0072) and MTRR (p = 0.0003) in the jejunum of the birds. Methionine deficiency increased the expression of BHMT (p = 0.0805) and MTRR (p = 0.0018). Higher expression of GSS was observed in birds that were supplemented with DL-HMTBA (p = 0.0672). Analyzing our results, it is preferable to supplement sulfur amino acids with DL-Met at the requirement level. Birds fed with DL-HMTBA showed worse weight gain (p = 0.0117) and higher feed conversion (p = 0.0170); methionine deficiency resulted in higher feed intake (p = 0.0214), lower weight gain (p<0.0001) and consequently higher feed conversion (p<0.0001). Based on the information found in this work, it is recommended to supplement sulfur amino acids with DL-Met at the level of compliance with the requirement.