Comparison of Bayesian models to estimate direct genomic values in multi-breed commercial beef cattle

Background

While several studies have examined the accuracy of direct genomic breeding values (DGV) within and across purebred cattle populations, the accuracy of DGV in crossbred or multi-breed cattle populations has been less well examined. Interest in the use of genomic tools for both selection and management has increased within the hybrid seedstock and commercial cattle sectors and research is needed to determine their efficacy. We predicted DGV for six traits using training populations of various sizes and alternative Bayesian models for a population of 3240 crossbred animals. Our objective was to compare alternate models with different assumptions regarding the distributions of single nucleotide polymorphism (SNP) effects to determine the optimal model for enhancing feasibility of multi-breed DGV prediction for the commercial beef industry.

Results

Realized accuracies ranged from 0.40 to 0.78. Randomly assigning 60 to 70% of animals to training (n ≈ 2000 records) yielded DGV accuracies with the smallest coefficients of variation. Mixture models (BayesB95, BayesCπ) and models that allow SNP effects to be sampled from distributions with unequal variances (BayesA, BayesB95) were advantageous for traits that appear or are known to be influenced by large-effect genes. For other traits, models differed little in prediction accuracy (~0.3 to 0.6%), suggesting that they are mainly controlled by small-effect loci.

Conclusions

The proportion (60 to 70%) of data allocated to training that optimized DGV accuracy and minimized the coefficient of variation of accuracy was similar to large dairy populations. Larger effects were estimated for some SNPs using BayesA and BayesB95 models because they allow unequal SNP variances. This substantially increased DGV accuracy for Warner-Bratzler Shear Force, for which large-effect quantitative trait loci (QTL) are known, while no loss in accuracy was observed for traits that appear to follow the infinitesimal model. Large decreases in accuracy (up to 0.07) occurred when SNPs that presumably tag large-effect QTL were over-regressed towards the mean in BayesC0 analyses. The DGV accuracies achieved here indicate that genomic selection has predictive utility in the commercial beef industry and that using models that reflect the genomic architecture of the trait can have predictive advantages in multi-breed populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0106-8) contains supplementary material, which is available to authorized users.     

The histone deacetylase inhibitor butyrate improves metabolism and reduces muscle atrophy during aging

Sarcopenia, the loss of skeletal muscle mass and function during aging, is a major contributor to disability and frailty in the elderly. Previous studies found a protective effect of reduced histone deacetylase activity in models of neurogenic muscle atrophy. Because loss of muscle mass during aging is associated with loss of motor neuron innervation, we investigated the potential for the histone deacetylase (HDAC) inhibitor butyrate to modulate age‐related muscle loss. Consistent with previous studies, we found significant loss of hindlimb muscle mass in 26‐month‐old C57Bl/6 female mice fed a control diet. Butyrate treatment starting at 16 months of age wholly or partially protected against muscle atrophy in hindlimb muscles. Butyrate increased muscle fiber cross‐sectional area and prevented intramuscular fat accumulation in the old mice. In addition to the protective effect on muscle mass, butyrate reduced fat mass and improved glucose metabolism in 26‐month‐old mice as determined by a glucose tolerance test. Furthermore, butyrate increased markers of mitochondrial biogenesis in skeletal muscle and whole‐body oxygen consumption without affecting activity. The increase in mass in butyrate‐treated mice was not due to reduced ubiquitin‐mediated proteasomal degradation. However, butyrate reduced markers of oxidative stress and apoptosis and altered antioxidant enzyme activity. Our data is the first to show a beneficial effect of butyrate on muscle mass during aging and suggests HDACs contribute to age‐related muscle atrophy and may be effective targets for intervention in sarcopenia and age‐related metabolic disease.     

The pyruvate, orthophosphate dikinase regulatory proteins of Arabidopsis are both bifunctional and interact with the catalytic and nucleotide-binding domains of pyruvate, orthophosphate dikinase

Pyruvate orthophosphate dikinase (PPDK) is a key enzyme in C(4) photosynthesis and is also found in C(3) plants. It is post-translationally modified by the PPDK regulatory protein (RP) that possesses both kinase and phosphotransferase activities. Phosphorylation and dephosphorylation of PPDK lead to inactivation and activation respectively. Arabidopsis thaliana contains two genes that encode chloroplastic (RP1) and cytosolic (RP2) isoforms of RP, and although RP1 has both kinase and phosphotransferase activities, to date RP2 has only been shown to act as a kinase. Here we demonstrate that RP2 is able to catalyse the dephosphorylation of PPDK, although at a slower rate than RP1 under the conditions of our assay. From yeast two-hybrid analysis we propose that RP1 binds to the central catalytic domain of PPDK, and that additional regions towards the carboxy and amino termini are required for a stable interaction between RP2 and PPDK. For 21 highly conserved amino acids in RP1, mutation of 15 of these reduced kinase and phosphotransferase activity, while mutation of six residues had no impact on either activity. We found no mutant in which only one activity was abolished. However, in some chimaeric fusions that comprised the amino and carboxy termini of RP1 and RP2 respectively, the kinase reaction was severely compromised but phosphotransferase activity remained unaffected. These findings are consistent with the findings that both RP1 and RP2 modulate reversibly the activity of PPDK, and possess one bifunctional active site or two separate sites in close proximity.     

Genetic impacts of Anacapa deer mice reintroductions following rat eradication

The Anacapa deer mouse is an endemic subspecies that inhabits Anacapa Island, part of Channel Islands National Park, California. We used mitochondrial DNA cytochrome c oxidase subunit II gene (COII) and 10 microsatellite loci to evaluate the levels of genetic differentiation and variation in ~1400 Anacapa deer mice sampled before and for 4 years after a black rat (Rattus rattus) eradication campaign that included trapping, captive holding and reintroduction of deer mice. Both mitochondrial and microsatellite analyses indicated significant differentiation between Anacapa deer mice and mainland mice, and genetic variability of mainland mice was significantly higher than Anacapa mice even prior to reintroduction. Bayesian cluster analysis and Principal Coordinates Analysis indicated that East, Middle and West Anacapa mice were genetically differentiated from each other, but translocation of mice among islands resulted in the East population becoming less distinct as a result of management. Levels of heterozygosity were similar before and after management. However, numerous private alleles in the founder populations were not observed after reintroduction and shifts in allele frequencies occurred, indicating that the reintroduced populations experienced substantial genetic drift. Surprisingly, two mitochondrial haplotypes observed in an earlier study of Anacapa deer mice were lost in the 20 years prior to the rat eradication program, leaving only a single haplotype in Anacapa deer mice. This study demonstrates how genetic monitoring can help to understand the re-establishment of endemic species after the eradication of invasive species and to evaluate the effectiveness of the management strategies employed.     

Characterization of chemical, radiochemical and optical properties of a dual-labeled MMP-9 targeting peptide

Optical imaging possesses similar sensitivity to nuclear imaging and has led to the emergence of multimodal approaches with dual-labeled nuclear/near-infrared (NIR) agents. The growing impact of (68)Ga (t(1/2)=68 min) labeled peptides on preclinical and clinical research offers a promising opportunity to merge the high spatial resolution of NIR imaging with the clinically-accepted positron emission tomography (PET). Previously, dual-labeled agents have been prepared with longer-lived radiometals and showed no detrimental effects on optical properties as a result of radiolabeling. In this study, we selected a peptide (M(2)) that targets MMP-2/9 and is dual-labeled with IRDye 800 CW and (68)Ga. Since (68)Ga chelation typically requires low pH (3.5-4) and elevated heating temperatures (95 °C), we sought to evaluate the impact of (68)Ga labeling on the optical properties of M(2). An efficient method for preparation of (68)Ga-M(2) was developed and reaction conditions were optimized. Stability studies in PBS, DTPA, and serum were performed and high levels of intact agent were evident under each condition. The addition of multiple reporters to a targeting agent adds further complexity to the characterization and validation and thus requires not only testing to ensure the agent is stable chemically and radiochemically, but also optically. Therefore, fluorescence properties were evaluated using a spectrofluorometer as well as by fluorescence detection via HPLC. It was determined that (68)Ga-labeling conditions did not impair the fluorescent properties of the agent. The agent was then used for in vivo imaging in a mouse model of heterotopic ossification (HO) with activated MMP-9 expression as an early biomarker which precedes mineralization. Although (68)Ga-complexation greatly reduced binding affinity of the peptide and negated tracer uptake on PET, NIR imaging showed consistent fluorescent signal that correlated to MMP-9 expression. This attests to the feasibility of using (68)Ga/NIR for dual-labeling of other peptides or small molecules for multimodality molecular imaging.     

Reagents for astatination of biomolecules. 5. Evaluation of hydrazone linkers in (211)At- and (125)I-labeled closo-decaborate(2-) conjugates of Fab' as a means of decreasing kidney retention

Evaluation of monoclonal antibody (mAb) fragments (e.g., Fab', Fab, or engineered fragments) as cancer-targeting reagents for therapy with the α-particle emitting radionuclide astatine-211 ((211)At) has been hampered by low in vivo stability of the label and a propensity of these proteins localize to kidneys. Fortunately, our group has shown that the low stability of the (211)At label, generally a meta- or para-[(211)At]astatobenzoyl conjugate, on mAb Fab' fragments can be dramatically improved by the use of closo-decaborate(2-) conjugates. However, the higher stability of radiolabeled mAb Fab' conjugates appears to result in retention of radioactivity in the kidneys. This investigation was conducted to evaluate whether the retention of radioactivity in kidney might be decreased by the use of an acid-cleavable hydrazone between the Fab' and the radiolabeled closo-decaborate(2-) moiety. Five conjugation reagents containing sulfhydryl-reactive maleimide groups, a hydrazone functionality, and a closo-decaborate(2-) moiety were prepared. In four of the five conjugation reagents, a discrete poly(ethylene glycol) (PEG) linker was used, and one substituent adjacent to the hydrazone was varied (phenyl, benzoate, anisole, or methyl) to provide varying acid sensitivity. In the initial studies, the five maleimido-closo-decaborate(2-) conjugation reagents were radioiodinated ((125)I or (131)I), then conjugated with an anti-PSMA Fab' (107-1A4 Fab'). Biodistributions of the five radioiodinated Fab' conjugates were obtained in nude mice at 1, 4, and 24 h post injection (pi). In contrast to closo-decaborate(2-) conjugated to 107-1A4 Fab' through a noncleavable linker, two conjugates containing either a benzoate or a methyl substituent on the hydrazone functionality displayed clearance rates from kidney, liver, and spleen that were similar to those obtained with directly radioiodinated Fab' (i.e., no conjugate). The maleimido-closo-decaborate(2-) conjugation reagent containing a benzoate substituent on the hydrazone was chosen for study with (211)At. That reagent was conjugated with 107-1A4 Fab', then labeled (separately) with (125)I and (211)At. The radiolabeled Fab' conjugates were coinjected into nude mice bearing LNCaP human tumor xenografts, and biodistribution data were obtained at 1, 4, and 24 h pi. Tumor targeting was achieved with both (125)I- and (211)At-labeled Fab', but the (211)At-labeled Fab' reached a higher concentration (25.56 ± 11.20 vs 11.97 ± 1.31%ID/g). Surprisingly, while the (125)I-labeled Fab' was cleared from kidney similar to earlier studies, the (211)At-labeled Fab'was not (i.e., kidney conc. for (125)I vs (211)At; 4 h, 13.14 ± 2.03 ID/g vs 42.28 ± 16.38%D/g; 24 h, 4.23 ± 1.57 ID/g vs 39.52 ± 15.87%ID/g). Since the Fab' conjugate is identical in both cases except for the radionuclide, it seems likely that the difference in tissue clearance seen is due to an effect that (211)At has on either the hydrazone cleavage or on the retention of a metabolite. Results from other studies in our laboratory suggest that the latter case is most likely. The hydrazone linkers tested do not provide the tissue clearance sought for (211)At, so additional hydrazones linkers will be evaluated. However, the results support the use of hydrazone linkers when Fab' conjugated with closo-decaborate(2-) reagents are radioiodinated.     

The carbon monoxide-releasing molecule, CORM-3 (RU(CO)(3) CL(glycinate)), targets respiration and oxidases in Campylobacter jejuni, generating hydrogen peroxide

Carbon monoxide (CO) is a classical respiratory inhibitor, but CO-releasing molecules (CO-RMs) have therapeutic value, increasing phagocytosis, and reducing sepsis-induced lethality. CORM-3, Ru(CO)(3) Cl(glycinate), a ruthenium-based carbonyl that liberates CO under physiological conditions, has previously been shown to inhibit bacterial growth and respiration, even at high concentrations of oxygen. Here, we report the effects of CORM-3 on the microaerophilic foodborne pathogen Campylobacter jejuni. Even at CO-RM (i.e., CO) concentrations that exceed dissolved oxygen levels, CORM-3 does not inhibit microaerobic growth. This insensitivity is not due to failure of CORM-3 to penetrate cells, as revealed by assay with extracellular myoglobin and by the ability of CO from externally added CORM-3 to bind intracellular membrane-associated respiratory oxidases. Even at almost 200 μ M oxygen, CORM-3 inhibits formate-dependent respiration and leads to generation of hydrogen peroxide. This work shows that CO-RMs have valuable properties as antimicrobial agents; however, growth inhibition does not always accompany inhibition of respiration, even when ambient oxygen concentrations are low.     

A LysM and SH3-domain containing region of the Listeria monocytogenes p60 protein stimulates accessory cells to promote activation of host NK cells

Listeria monocytogenes (Lm) infection induces rapid and robust activation of host natural killer (NK) cells. Here we define a region of the abundantly secreted Lm endopeptidase, p60, that potently but indirectly stimulates NK cell activation in vitro and in vivo. Lm expression of p60 resulted in increased IFNγ production by naïve NK cells co-cultured with treated dendritic cells (DCs). Moreover, recombinant p60 protein stimulated activation of naive NK cells when co-cultured with TLR or cytokine primed DCs in the absence of Lm. Intact p60 protein weakly digested bacterial peptidoglycan (PGN), but neither muropeptide recognition by RIP2 nor the catalytic activity of p60 was required for NK cell activation. Rather, the immune stimulating activity mapped to an N-terminal region of p60, termed L1S. Treatment of DCs with a recombinant L1S polypeptide stimulated them to activate naïve NK cells in a cell culture model. Further, L1S treatment activated NK cells in vivo and increased host resistance to infection with Francisella tularensis live vaccine strain (LVS). These studies demonstrate an immune stimulating function for a bacterial LysM domain-containing polypeptide and suggest that recombinant versions of L1S or other p60 derivatives can be used to promote NK cell activation in therapeutic contexts.