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.     

Structural prediction and comparative docking studies of psychrophilic ��- Galactosidase with lactose, ONPG and PNPG against its counter parts of mesophilic and thermophilic enzymes

Enzymes from psychrophiles catalyze the reactions at low temperatures with higher specific activity. Among all the psychrophilic enzymes produced, cold active β-galactosidase from marine psychrophiles revalorizes a new arena in numerous areas at industrial level. The hydrolysis of lactose in to glucose and galactose by cold active β-galactosidase offers a new promising approach in removal of lactose from milk to overcome the problem of lactose intolerance. Herein we propose, a 3D structure of cold active β-galactosidase enzyme sourced from Pseudoalteromonas haloplanktis by using Modeler 9v8 and best model was developed having 88% of favourable region in ramachandran plot. Modelling was followed by docking studies with the help of Auto dock 4.0 against the three substrates lactose, ONPG and PNPG. In addition, comparative docking studies were also performed for the 3D model of psychrophilic β-galactosidase with mesophilic and thermophilic enzymes. Docking studies revealed that binding affinity of enzyme towards the three different substrates is more for psychrophilic enzyme when compared with mesophilic and thermophilic enzymes. It indicates that the enzyme has high specific activity at low temperature when compared with mesophilic and thermophilic enzymes.     

Utilizing laboratory and field studies to determine physiological responses of cattle to multiple environmental stressors

Heat stress studies are often conducted using controlled laboratory exposures or field exposures. Each approach has limitations and provides a partial understanding of complex interactions between simultaneous environmental stressors. The question is how similar the responses are in each situation. Several physiological measures of thermal status were used to compare heat stress responses of cattle in controlled chamber stress tests and fluctuating field conditions. Angus steers (N=23; 318±8 kg BW) were first placed on either endophyte-infected or -uninfected tall fescue pastures for the field exposure, followed by a controlled heat challenge, which exacerbates the condition known as fescue toxicosis. During the controlled heat challenge, steers were assigned to diets of either 0 or 40 μg ergovaline/kg/d to maintain the treatment states. Respiration rate (RR) was measured via flank counting and telemetric temperature transmitters in the rumen of each animal monitored core temperature (T_rum). Linear regression fit models for RR, T_rum, and air temperature (T_a) were utilized to compare relationships between field and chamber exposure. Correlation coefficients for RR were similar during both chamber (R=0.69) and field exposures (R=0.72). Respiration rate showed greater responsiveness to change in T_a under field conditions having twice the slope (4.40 versus 1.75 bpm/°C) and a lower Y-intercept (−42.14 versus +30.97 bpm) compared to the chamber run. Ruminal temperature was consistent between exposures showing a similar slope (0.04 versus 0.03 °C T_rum/°C T_a) and Y-intercept (38.40 versus 39.30 °C) for its relationship with T_a. Despite respiration rate being the more sensitive indicator of heat stress, ruminal temperature proved to be the most consistent between environments.     

Patterns of heat response and adaptation on summer pasture: A comparison of heat-sensitive (Angus) and -tolerant (Romosinuano) cattle

Heat stress in Bos taurus cattle is a problem that affects many regions of the world. Numerous studies have focused on heat stress in feedlots or environmental chambers; but few have looked at undisturbed cattle on pasture. The present study followed two Bos taurus cattle breeds throughout a mid-Missouri summer to determine thermoregulatory responses to fluctuating summer air temperature (T_a), as well as differences in adaptation to heat. Heat-sensitive Angus steers (ANG; n=22; 480±7.15 kg BW), and heat-tolerant Romosinuano steers (RO; n=11; 352±6 kg BW) were monitored on 12 day from June through August of 2009 in an endophyte free tall fescue pasture. Data were grouped into two, six-day periods representing peak (Period 1) and late (Period 2) summer for determination of adaptation. Respiration rate (RR) was measured via flank counting and telemetric temperature transmitters in the rumen of each animal monitored core temperature (T_rum). Romosinuano sustained a lower (P<0.05) RR and T_rum compared to ANG during both periods. Linear relationships for RR and T_rum, compared against T_a for both Periods were determined. Slopes of RR to T_a from Period 1 to Period 2 decreased (P<0.05) from 2.63 to 1.08 bpm/°C and 2.25 to 0.49 bpm/°C for ANG and RO, respectively. Slopes of T_rum to T_a also decreased (P<0.05) from Periods 1 to 2 from 0.12 to 0.02 °C T_rum/°C T_a for ANG; however, RO showed no differences between periods. Although Romosinuano have a lower respiration rate and ruminal temperature than Angus, they share a similar pattern of respiration rate adaptation from early to late summer periods.     

Long non‐coding RNAs in melanoma

Melanoma is the most lethal cutaneous cancer with a highly aggressive and metastatic phenotype. While recent genetic and epigenetic studies have shed new insights into the mechanism of melanoma development, the involvement of regulatory non‐coding RNAs remain unclear. Long non‐coding RNAs (lncRNAs) are a group of endogenous non‐protein‐coding RNAs with the capacity to regulate gene expression at multiple levels. Recent evidences have shown that lncRNAs can regulate many cellular processes, such as cell proliferation, differentiation, migration and invasion. In the melanoma, deregulation of a number of lncRNAs, such as HOTAIR, MALAT1, BANCR, ANRIL, SPRY‐IT1 and SAMMSON, have been reported. Our review summarizes the functional role of lncRNAs in melanoma and their potential clinical application for diagnosis, prognostication and treatment.     

Changes in odor quality discrimination following recovery from olfactory nerve transection

Following recovery from olfactory nerve transection, animals regain their ability to discriminate between odors. Odor discrimination is restored after new neurons establish connections with the olfactory bulb. However, it is not known if the new connections alter odor quality perception. To address this question, 20 adult hamsters were first trained to discriminate between cinnamon and strawberry odors. After reaching criterion (> or = 90% correct response), half of the animals received a bilateral nerve transection (BTX) and half a surgical sham procedure. Animals were not tested again until day 40, a point in recovery when connections are re-established with the bulb. When BTX animals were tested without food reinforcement, they could not perform the odor discrimination task. Sham animals, however, could discriminate, demonstrating that the behavioral response had not been extinguished during the 40 day period. When reinforcement was resumed, BTX animals were able to discriminate between cinnamon and strawberry after four test sessions. In addition, their ability to discriminate between these two familiar odors was no different than that of BTX and sham animals tested with two novel odors, baby powder and coffee. These findings suggest that, after recovery from nerve transection, there are alterations in sensory perception and that restoration of odor quality discrimination requires that the animal must again learn to associate individual odor sensations with a behavioral response.      

Multiplex PCR and Next Generation Sequencing for the Non-Invasive Detection of Bladder Cancer

Background

Highly sensitive and specific urine-based tests to detect either primary or recurrent bladder cancer have proved elusive to date. Our ever increasing knowledge of the genomic aberrations in bladder cancer should enable the development of such tests based on urinary DNA.

Methods

DNA was extracted from urine cell pellets and PCR used to amplify the regions of the TERT promoter and coding regions of FGFR3, PIK3CA, TP53, HRAS, KDM6A and RXRA which are frequently mutated in bladder cancer. The PCR products were barcoded, pooled and paired-end 2 x 250 bp sequencing performed on an Illumina MiSeq. Urinary DNA was analysed from 20 non-cancer controls, 120 primary bladder cancer patients (41 pTa, 40 pT1, 39 pT2+) and 91 bladder cancer patients post-TURBT (89 cancer-free).

Results

Despite the small quantities of DNA extracted from some urine cell pellets, 96% of the samples yielded mean read depths >500. Analysing only previously reported point mutations, TERT mutations were found in 55% of patients with bladder cancer (independent of stage), FGFR3 mutations in 30% of patients with bladder cancer, PIK3CA in 14% and TP53 mutations in 12% of patients with bladder cancer. Overall, these previously reported bladder cancer mutations were detected in 86 out of 122 bladder cancer patients (70% sensitivity) and in only 3 out of 109 patients with no detectable bladder cancer (97% specificity).

Conclusion

This simple, cost-effective approach could be used for the non-invasive surveillance of patients with non-muscle-invasive bladder cancers harbouring these mutations. The method has a low DNA input requirement and can detect low levels of mutant DNA in a large excess of normal DNA. These genes represent a minimal biomarker panel to which extra markers could be added to develop a highly sensitive diagnostic test for bladder cancer.     

Young patients with colorectal cancer have poor survival in the first twenty months after operation and predictable survival in the medium and long-term: Analysis of survival and prognostic markers

Introduction

Male breast cancer (MBC) is a rare, yet potentially aggressive disease. Although literature regarding female breast cancer (FBC) is extensive, little is known about the etiopathogenesis of male breast cancer. Studies from our laboratory show that MBCs have a distinct immunophenotypic profile, suggesting that the etiopathogenesis of MBC is different from FBCs. The aim of this study was to evaluate and correlate the immunohistochemical expression of cell cycle proteins in male breast carcinoma to significant clinico-biological endpoints.

Methods

75 cases of MBC were identified using the records of the Saskatchewan Cancer Agency over 26 years (1970-1996). Cases were reviewed and analyzed for the immunohistochemical expression of PCNA, Ki67, p27, p16, p57, p21, cyclin-D1 and c-myc and correlated to clinico-biological endpoints of tumor size, node status, stage of the disease, and disease free survival (DFS).

Results

Decreased DFS was observed in the majority of tumors that overexpressed PCNA (98%, p = 0.004). The overexpression of PCNA was inversely correlated to the expression of Ki67 which was predominantly negative (78.3%). Cyclin D1 was overexpressed in 83.7% of cases. Cyclin D1 positive tumors were smaller than 2 cm (55.6%, p = 0.005), had a low incidence of lymph node metastasis (38.2%, p = 0.04) and were associated with increased DFS of >150 months (p = 0.04). Overexpression of c-myc (90%) was linked with a higher incidence of node negativity (58.3%, p = 0.006) and increased DFS (p = 0.04). p27 over expression was associated with decreased lymph node metastasis (p = 0.04). P21 and p57 positive tumors were related to decreased DFS (p = 0.04). Though p16 was overexpressed in 76.6%, this did not reach statistical significance with DFS (p = 0.06) or nodal status (p = 0.07).

Conclusion

Aberrant cell cycle protein expression supports our view that these are important pathways involved in the etiopathogenesis of MBC. Tumors with overexpression of Cyclin D1 and c-myc had better outcomes, in contrast to tumors with overexpression of p21, p57, and PCNA with significantly worse outcomes. P27 appears to be a predictive marker for lymph nodal status. Such observation strongly suggests that dysregulation of cell cycle proteins may play a unique role in the initiation and progression of disease in male breast cancer. Such findings open up new avenues for the treatment of MBC as a suitable candidate for novel CDK-based anticancer therapies in the future.     

Evaluation of healing following frenectomy

Aberrant frenum attachment would cause plaque accumulation and malalignment of teeth. It can be managed by frenotomy or frenectomy methods, through a conventional surgical technique or laser technique. Therefore, it is of interest to compare frenectomy healing surgical and laser techniques. Data from 51 outpatients and post-operative healing of frenectomy was assessed by Landry’s healing score index using 3 weeks postoperative photographs followed by statistical analysis. Based on the healing score index, the laser technique showed better outcomes than the surgical technique. Moreover, the association between the management of high frenal attachment and the healing score index was found to be statistically significant.     

QTLs associated with dry matter intake,metabolic mid-test weight,growth and feed efficiency have little overlap across 4 beef cattle studies

Background

The identification of genetic markers associated with complex traits that are expensive to record such as feed intake or feed efficiency would allow these traits to be included in selection programs. To identify large-effect QTL, we performed a series of genome-wide association studies and functional analyses using 50 K and 770 K SNP genotypes scored in 5,133 animals from 4 independent beef cattle populations (Cycle VII, Angus, Hereford and Simmental × Angus) with phenotypes for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake.

Results

A total of 5, 6, 11 and 10 significant QTL (defined as 1-Mb genome windows with Bonferroni-corrected P-value <0.05) were identified for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake, respectively. The identified QTL were population-specific and had little overlap across the 4 populations. The pleiotropic or closely linked QTL on BTA 7 at 23 Mb identified in the Angus population harbours a promising candidate gene ACSL6 (acyl-CoA synthetase long-chain family member 6), and was the largest effect QTL associated with dry matter intake and mid-test body weight explaining 10.39% and 14.25% of the additive genetic variance, respectively. Pleiotropic or closely linked QTL associated with average daily gain and mid-test body weight were detected on BTA 6 at 38 Mb and BTA 7 at 93 Mb confirming previous reports. No QTL for residual feed intake explained more than 2.5% of the additive genetic variance in any population. Marker-based estimates of heritability ranged from 0.21 to 0.49 for residual feed intake across the 4 populations.

Conclusions

This GWAS study, which is the largest performed for feed efficiency and its component traits in beef cattle to date, identified several large-effect QTL that cumulatively explained a significant percentage of additive genetic variance within each population. Differences in the QTL identified among the different populations may be due to differences in power to detect QTL, environmental variation, or differences in the genetic architecture of trait variation among breeds. These results enhance our understanding of the biology of growth, feed intake and utilisation in beef cattle.