A Rank-Based Nonparametric Method for Mapping Quantitative Trait Loci in Outbred Half-Sib Pedigrees: Application to Milk Production in a Granddaughter Design

We describe the development of a multipoint nonparametric quantitative trait loci mapping method based on the Wilcoxon rank-sum test applicable to outbred half-sib pedigrees. The method has been evaluated on a simulated dataset and its efficiency compared with interval mapping by using regression. It was shown that the rank-based approach is slightly inferior to regression when the residual variance is homoscedastic normal; however, in three out of four other scenarios envisaged, i.e., residual variance heteroscedastic normal, homoscedastic skewed, and homoscedastic positively kurtosed, the latter outperforms the former one. Both methods were applied to a real data set analyzing the effect of bovine chromosome 6 on milk yield and composition by using a 125-cM map comprising 15 microsatellites and a granddaughter design counting 1158 Holstein-Friesian sires

Dominance and parent-of-origin effects of coding and non-coding alleles at the acylCoA-diacylglycerol-acyltransferase (DGAT1) gene on milk production traits in German Holstein cows

<p>Abstract</p> <p>Background</p> <p>Substantial gene substitution effects on milk production traits have formerly been reported for alleles at the K232A and the promoter VNTR loci in the bovine acylCoA-diacylglycerol-acyltransferase 1 (<it>DGAT1</it>) gene by using data sets including sires with accumulated phenotypic observations of daughters (breeding values, daughter yield deviations). However, these data sets prevented analyses with respect to dominance or parent-of-origin effects, although an increasing number of reports in the literature outlined the relevance of non-additive gene effects on quantitative traits.</p> <p>Results</p> <p>Based on a data set comprising German Holstein cows with direct trait measurements, we first confirmed the previously reported association of <it>DGAT1 </it>promoter VNTR alleles with milk production traits. We detected a dominant mode of effects for the <it>DGAT1 </it>K232A and promoter VNTR alleles. Namely, the contrasts between the effects of heterozygous individuals at the <it>DGAT1 </it>loci differed significantly from the midpoint between the effects for the two homozygous genotypes for several milk production traits, thus indicating the presence of dominance. Furthermore, we identified differences in the magnitude of effects between paternally and maternally inherited <it>DGAT1 </it>promoter VNTR – K232A haplotypes indicating parent-of-origin effects on milk production traits.</p> <p>Conclusion</p> <p>Non-additive effects like those identified at the bovine <it>DGAT1 </it>locus have to be accounted for in more specific QTL detection models as well as in marker assisted selection schemes. The <it>DGAT1 </it>alleles in cattle will be a useful model for further investigations on the biological background of non-additive effects in mammals due to the magnitude and consistency of their effects on milk production traits.</p

Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle

Background Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals. Results Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle. Conclusions This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought

Deconstructing the DGAT1 enzyme: membrane interactions at substrate binding sites

Diacylglycerol acyltransferase 1 (DGAT1) is a key enzyme in the triacylglyceride synthesis pathway. Bovine DGAT1 is an endoplasmic reticulum membrane-bound protein associated with the regulation of fat content in milk and meat. The aim of this study was to evaluate the interaction of DGAT1 peptides corresponding to putative substrate binding sites with different types of model membranes. Whilst these peptides are predicted to be located in an extramembranous loop of the membrane-bound protein, their hydrophobic substrates are membrane-bound molecules. In this study, peptides corresponding to the binding sites of the two substrates involved in the reaction were examined in the presence of model membranes in order to probe potential interactions between them that might influence the subsequent binding of the substrates. Whilst the conformation of one of the peptides changed upon binding several types of micelles regardless of their surface charge, suggesting binding to hydrophobic domains, the other peptide bound strongly to negatively-charged model membranes. This binding was accompanied by a change in conformation, and produced leakage of the liposome-entrapped dye calcein. The different hydrophobic and electrostatic interactions observed suggest the peptides may be involved in the interactions of the enzyme with membrane surfaces, facilitating access of the catalytic histidine to the triacylglycerol substrates

Scans for signatures of selection in Russian cattle breed genomes reveal new candidate genes for environmental adaptation and acclimation

Domestication and selective breeding has resulted in over 1000 extant cattle breeds. Many of these breeds do not excel in important traits but are adapted to local environments. These adaptations are a valuable source of genetic material for efforts to improve commercial breeds. As a step toward this goal we identified candidate regions to be under selection in genomes of nine Russian native cattle breeds adapted to survive in harsh climates. After comparing our data to other breeds of European and Asian origins we found known and novel candidate genes that could potentially be related to domestication, economically important traits and environmental adaptations in cattle. The Russian cattle breed genomes contained regions under putative selection with genes that may be related to adaptations to harsh environments (e.g., AQP5, RAD50, and RETREG1). We found genomic signatures of selective sweeps near key genes related to economically important traits, such as the milk production (e.g., DGAT1, ABCG2), growth (e.g., XKR4), and reproduction (e.g., CSF2). Our data point to candidate genes which should be included in future studies attempting to identify genes to improve the extant breeds and facilitate generation of commercial breeds that fit better into the environments of Russia and other countries with similar climates

Accuracy of genomic breeding values in multi-breed dairy cattle populations

<p>Abstract</p> <p>Background</p> <p>Two key findings from genomic selection experiments are 1) the reference population used must be very large to subsequently predict accurate genomic estimated breeding values (GEBV), and 2) prediction equations derived in one breed do not predict accurate GEBV when applied to other breeds. Both findings are a problem for breeds where the number of individuals in the reference population is limited. A multi-breed reference population is a potential solution, and here we investigate the accuracies of GEBV in Holstein dairy cattle and Jersey dairy cattle when the reference population is single breed or multi-breed. The accuracies were obtained both as a function of elements of the inverse coefficient matrix and from the realised accuracies of GEBV.</p> <p>Methods</p> <p>Best linear unbiased prediction with a multi-breed genomic relationship matrix (GBLUP) and two Bayesian methods (BAYESA and BAYES_SSVS) which estimate individual SNP effects were used to predict GEBV for 400 and 77 young Holstein and Jersey bulls respectively, from a reference population of 781 and 287 Holstein and Jersey bulls, respectively. Genotypes of 39,048 SNP markers were used. Phenotypes in the reference population were de-regressed breeding values for production traits. For the GBLUP method, expected accuracies calculated from the diagonal of the inverse of coefficient matrix were compared to realised accuracies.</p> <p>Results</p> <p>When GBLUP was used, expected accuracies from a function of elements of the inverse coefficient matrix agreed reasonably well with realised accuracies calculated from the correlation between GEBV and EBV in single breed populations, but not in multi-breed populations. When the Bayesian methods were used, realised accuracies of GEBV were up to 13% higher when the multi-breed reference population was used than when a pure breed reference was used. However no consistent increase in accuracy across traits was obtained.</p> <p>Conclusion</p> <p>Predicting genomic breeding values using a genomic relationship matrix is an attractive approach to implement genomic selection as expected accuracies of GEBV can be readily derived. However in multi-breed populations, Bayesian approaches give higher accuracies for some traits. Finally, multi-breed reference populations will be a valuable resource to fine map QTL.</p

Effects of attention on the control of locomotion in individuals with chronic low back pain

<p>Abstract</p> <p>Background</p> <p>People who suffer from low back pain (LBP) exhibit an abnormal gait pattern, characterized by shorter stride length, greater step width, and an impaired thorax-pelvis coordination which may undermine functional walking. As a result, gait in LBP may require stronger cognitive regulation compared to pain free subjects thereby affecting the degree of automaticity of gait control. Conversely, because chronic pain has a strong attentional component, diverting attention away from the pain might facilitate a more efficient walking pattern.</p> <p>Methods</p> <p>Twelve individuals with LBP and fourteen controls participated. Subjects walked on a treadmill at comfortable speed, under varying conditions of attentional load: (a) no secondary task, (b) naming the colors of squares on a screen, (c) naming the colors of color words ("color Stroop task"), and (d) naming the colors of words depicting motor activities. Markers were attached to the thorax, pelvis and feet. Motion was recorded using a three-camera SIMI system with a sample frequency of 100 Hz. To examine the effects of health status and attention on gait, mean and variability of stride parameters were calculated. The coordination between thoracic and pelvic rotations was quantified through the mean and variability of the relative phase between those oscillations.</p> <p>Results</p> <p>LBP sufferers had a lower walking speed, and consequently a smaller stride length and lower mean thorax-pelvis relative phase. Stride length variability was significantly lower in the LBP group but no significant effect of attention was observed. In both groups gait adaptations were found under performance of an attention demanding task, but significantly more so in individuals with LBP as indicated by an interaction effect on relative phase variability.</p> <p>Conclusion</p> <p>Gait in LBP sufferers was characterized by less variable upper body movements. The diminished flexibility in trunk coordination was aggravated under the influence of an attention demanding task. This provides further evidence that individuals with LBP tighten their gait control, and this suggests a stronger cognitive regulation of gait coordination in LBP. These changes in gait coordination reduce the capability to deal with unexpected perturbations, and are therefore maladaptive.</p

Genome-wide association analysis of thirty one production, health, reproduction and body conformation traits in contemporary U.S. Holstein cows

<p>Abstract</p> <p>Background</p> <p>Genome-wide association analysis is a powerful tool for annotating phenotypic effects on the genome and knowledge of genes and chromosomal regions associated with dairy phenotypes is useful for genome and gene-based selection. Here, we report results of a genome-wide analysis of predicted transmitting ability (PTA) of 31 production, health, reproduction and body conformation traits in contemporary Holstein cows.</p> <p>Results</p> <p>Genome-wide association analysis identified a number of candidate genes and chromosome regions associated with 31 dairy traits in contemporary U.S. Holstein cows. Highly significant genes and chromosome regions include: BTA13's <it>GNAS </it>region for milk, fat and protein yields; BTA7's <it>INSR </it>region and BTAX's <it>LOC520057 </it>and <it>GRIA3 </it>for daughter pregnancy rate, somatic cell score and productive life; BTA2's <it>LRP1B </it>for somatic cell score; BTA14's <it>DGAT1-NIBP </it>region for fat percentage; <it>BTA1</it>'s <it>FKBP2 </it>for protein yields and percentage, BTA26's <it>MGMT </it>and BTA6's <it>PDGFRA </it>for protein percentage; BTA18's 53.9-58.7 Mb region for service-sire and daughter calving ease and service-sire stillbirth; BTA18's <it>PGLYRP1</it>-<it>IGFL1 </it>region for a large number of traits; BTA18's <it>LOC787057 </it>for service-sire stillbirth and daughter calving ease; BTA15's <it>CD82</it>, BTA23's <it>DST </it>and the <it>MOCS1</it>-<it>LRFN2 </it>region for daughter stillbirth; and BTAX's <it>LOC520057 </it>and <it>GRIA3 </it>for daughter pregnancy rate. For body conformation traits, BTA11, BTAX, BTA10, BTA5, and BTA26 had the largest concentrations of SNP effects, and <it>PHKA2 </it>of BTAX and <it>REN </it>of BTA16 had the most significant effects for body size traits. For body shape traits, BTAX, BTA19 and BTA3 were most significant. Udder traits were affected by BTA16, BTA22, BTAX, BTA2, BTA10, BTA11, BTA20, BTA22 and BTA25, teat traits were affected by BTA6, BTA7, BTA9, BTA16, BTA11, BTA26 and BTA17, and feet/legs traits were affected by BTA11, BTA13, BTA18, BTA20, and BTA26.</p> <p>Conclusions</p> <p>Genome-wide association analysis identified a number of genes and chromosome regions associated with 31 production, health, reproduction and body conformation traits in contemporary Holstein cows. The results provide useful information for annotating phenotypic effects on the dairy genome and for building consensus of dairy QTL effects.</p

Sensitivity of methods for estimating breeding values using genetic markers to the number of QTL and distribution of QTL variance

The objective of this simulation study was to compare the effect of the number of QTL and distribution of QTL variance on the accuracy of breeding values estimated with genomewide markers (MEBV). Three distinct methods were used to calculate MEBV: a Bayesian Method (BM), Least Angle Regression (LARS) and Partial Least Square Regression (PLSR). The accuracy of MEBV calculated with BM and LARS decreased when the number of simulated QTL increased. The accuracy decreased more when QTL had different variance values than when all QTL had an equal variance. The accuracy of MEBV calculated with PLSR was affected neither by the number of QTL nor by the distribution of QTL variance. Additional simulations and analyses showed that these conclusions were not affected by the number of individuals in the training population, by the number of markers and by the heritability of the trait. Results of this study show that the effect of the number of QTL and distribution of QTL variance on the accuracy of MEBV depends on the method that is used to calculate MEBV

Structure-Function Analysis of Diacylglycerol Acyltransferase Sequences from 70 Organisms

<p>Abstract</p> <p>Background</p> <p>Diacylglycerol acyltransferase families (DGATs) catalyze the final and rate-limiting step of triacylglycerol (TAG) biosynthesis in eukaryotic organisms. Understanding the roles of DGATs will help to create transgenic plants with value-added properties and provide clues for therapeutic intervention for obesity and related diseases. The objective of this analysis was to identify conserved sequence motifs and amino acid residues for better understanding of the structure-function relationship of these important enzymes.</p> <p>Results</p> <p>117 DGAT sequences from 70 organisms including plants, animals, fungi and human are obtained from database search using tung tree DGATs. Phylogenetic analysis separates these proteins into DGAT1 and DGAT2 subfamilies. These DGATs are integral membrane proteins with more than 40% of the total amino acid residues being hydrophobic. They have similar properties and amino acid composition except that DGAT1s are approximately 20 kDa larger than DGAT2s. DGAT1s and DGAT2s have 41 and 16 completely conserved amino acid residues, respectively, although only two of them are shared by all DGATs. These residues are distributed in 7 and 6 sequence blocks for DGAT1s and DGAT2s, respectively, and located at the carboxyl termini, suggesting the location of the catalytic domains. These conserved sequence blocks do not contain the putative neutral lipid-binding domain, mitochondrial targeting signal, or ER retrieval motif. The importance of conserved residues has been demonstrated by site-directed and natural mutants.</p> <p>Conclusions</p> <p>This study has identified conserved sequence motifs and amino acid residues in all 117 DGATs and the two subfamilies. None of the completely conserved residues in DGAT1s and DGAT2s is present in recently reported isoforms in the multiple sequences alignment, raising an important question how proteins with completely different amino acid sequences could perform the same biochemical reaction. The sequence analysis should facilitate studying the structure-function relationship of DGATs with the ultimate goal to identify critical amino acid residues for engineering superb enzymes in metabolic engineering and selecting enzyme inhibitors in therapeutic application for obesity and related diseases.</p