Genetic architecture of complex traits and accuracy of genomic prediction: coat colour, milk-fat percentage, and type in Holstein cattle as contrasting model traits

Prediction of genetic merit using dense SNP genotypes can be used for estimation of breeding values for selection of livestock, crops, and forage species; for prediction of disease risk; and for forensics. The accuracy of these genomic predictions depends in part on the genetic architecture of the trait, in particular number of loci affecting the trait and distribution of their effects. Here we investigate the difference among three traits in distribution of effects and the consequences for the accuracy of genomic predictions. Proportion of black coat colour in Holstein cattle was used as one model complex trait. Three loci, KIT, MITF, and a locus on chromosome 8, together explain 24% of the variation of proportion of black. However, a surprisingly large number of loci of small effect are necessary to capture the remaining variation. A second trait, fat concentration in milk, had one locus of large effect and a host of loci with very small effects. Both these distributions of effects were in contrast to that for a third trait, an index of scores for a number of aspects of cow confirmation ("overall type"), which had only loci of small effect. The differences in distribution of effects among the three traits were quantified by estimating the distribution of variance explained by chromosome segments containing 50 SNPs. This approach was taken to account for the imperfect linkage disequilibrium between the SNPs and the QTL affecting the traits. We also show that the accuracy of predicting genetic values is higher for traits with a proportion of large effects (proportion black and fat percentage) than for a trait with no loci of large effect (overall type), provided the method of analysis takes advantage of the distribution of loci effects.     

Chemical composition of the wound plug and entire plants for species of the coenocytic green alga, Caulerpa

Levels of ash, NaOH-soluble protein, lipid, and TCA-soluble carbohydrate were determined tor seven species of the green alga Caulerpa, for both the entire plant and wound plugs produced by cutting the rhizomes or blades. Insoluble carbohydrate was estimated by subtraction. Wound plugs had higher ash levels and a distinctly higher level of total carbohydrate. Both entire plants and wound plugs had low levels of lipid and protein. Carbohydrate is the major component of the wound plug in Caulerpa in contrast to protein reported for Bryopsis.     

Differential roles of p55 and p75 tumor necrosis factor receptors on stretch-induced pulmonary edema in mice

Ventilator-induced lung injury plays a crucial role in the outcome of patients with acute lung injury. Previous studies have shown a role for the cytokine tumor necrosis factor-alpha (TNF) in stretch-induced alveolar neutrophil recruitment, but the involvement of TNF in stretch-induced pulmonary edema is unclear. We investigated the effects of TNF through its individual p55 and p75 receptors on early pulmonary edema formation during high stretch ventilation, before neutrophil infiltration. Anesthetized wild-type or TNF receptor single/double knockout mice were ventilated with high tidal volume ( approximately 38 ml/kg) for 2 h or until they developed arterial hypotension. Pulmonary edema was assessed by physiological parameters including respiratory mechanics and blood gases, and by lavage fluid protein, lung wet:dry weight ratio, and lung permeability measurements using fluorescence-labeled albumin. High stretch ventilation in wild-type and TNF receptor double knockout animals induced similar pulmonary edema, and only 25-30% of mice completed the protocol. In contrast, the p55 receptor knockout mice were strongly protected from edema formation, with all animals completing the protocol. Myeloperoxidase assay indicated that this protective effect was not associated with decreased pulmonary neutrophil sequestration. The p75 receptor knockout mice, however, displayed increased susceptibility to edema formation, and no animals survived the full 2 h. These results demonstrate a novel role for TNF signaling (independent from its effects on neutrophil recruitment) specifically through the p55 receptor, in promoting high stretch-induced pulmonary edema, whereas p75 signaling may play an opposing role.     

Impact of QTL properties on the accuracy of multi-breed genomic prediction

Background

Although simulation studies show that combining multiple breeds in one reference population increases accuracy of genomic prediction, this is not always confirmed in empirical studies. This discrepancy might be due to the assumptions on quantitative trait loci (QTL) properties applied in simulation studies, including number of QTL, spectrum of QTL allele frequencies across breeds, and distribution of allele substitution effects. We investigated the effects of QTL properties and of including a random across- and within-breed animal effect in a genomic best linear unbiased prediction (GBLUP) model on accuracy of multi-breed genomic prediction using genotypes of Holstein-Friesian and Jersey cows.

Methods

Genotypes of three classes of variants obtained from whole-genome sequence data, with moderately low, very low or extremely low average minor allele frequencies (MAF), were imputed in 3000 Holstein-Friesian and 3000 Jersey cows that had real high-density genotypes. Phenotypes of traits controlled by QTL with different properties were simulated by sampling 100 or 1000 QTL from one class of variants and their allele substitution effects either randomly from a gamma distribution, or computed such that each QTL explained the same variance, i.e. rare alleles had a large effect. Genomic breeding values for 1000 selection candidates per breed were estimated using GBLUP modelsincluding a random across- and a within-breed animal effect.

Results

For all three classes of QTL allele frequency spectra, accuracies of genomic prediction were not affected by the addition of 2000 individuals of the other breed to a reference population of the same breed as the selection candidates. Accuracies of both single- and multi-breed genomic prediction decreased as MAF of QTL decreased, especially when rare alleles had a large effect. Accuracies of genomic prediction were similar for the models with and without a random within-breed animal effect, probably because of insufficient power to separate across- and within-breed animal effects.

Conclusions

Accuracy of both single- and multi-breed genomic prediction depends on the properties of the QTL that underlie the trait. As QTL MAF decreased, accuracy decreased, especially when rare alleles had a large effect. This demonstrates that QTL properties are key parameters that determine the accuracy of genomic prediction.

Electronic supplementary material

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

Predicted structures of agonist and antagonist bound complexes of adenosine A3 receptor

We used the GEnSeMBLE Monte Carlo method to predict ensemble of the 20 best packings (helix rotations and tilts) based on the neutral total energy (E) from a vast number (10 trillion) of potential packings for each of the four subtypes of the adenosine G protein-coupled receptors (GPCRs), which are involved in many cytoprotective functions. We then used the DarwinDock Monte Carlo methods to predict the binding pose for the human A(3) adenosine receptor (hAA(3)R) for subtype selective agonists and antagonists. We found that all four A(3) agonists stabilize the 15th lowest conformation of apo-hAA(3)R while also binding strongly to the 1st and 3rd. In contrast the four A(3) antagonists stabilize the 2nd or 3rd lowest conformation. These results show that different ligands can stabilize different GPCR conformations, which will likely affect function, complicating the design of functionally unique ligands. Interestingly all agonists lead to a trans χ1 angle for W6.48 that experiments on other GPCRs associate with G-protein activation while all 20 apo-AA(3)R conformations have a W6.48 gauche+ χ1 angle associated experimentally with inactive GPCRs for other systems. Thus docking calculations have identified critical ligand-GPCR structures involved with activation. We found that the predicted binding site for selective agonist Cl-IB-MECA to the predicted structure of hAA(3)R shows favorable interactions to three subtype variable residues, I253(6.58), V169(EL2), and Q167(EL2), while the predicted structure for hAA(2A)R shows weakened to the corresponding amino acids: T256(6.58), E169(EL2), and L167(EL2), explaining the observed subtype selectivity.     

The use of catheter-tipped pressure transducers for chronic measurement of genital tract pressures in the ewe. II. Effects of oxytocin and prostaglandin F(2)alpha

Chronic catheterisation of the uterus, ampulla, and abdomen was performed in five ewes using solid-state, catheter-tipped pressure transducers. The catheters remained in place for up to 129 d, allowing in vivo studies of the effects of oxytocin and prostaglandin F(2)alpha (PGF(2)alpha). These agents did not produce any measurable increase in abdominal pressure. Intravenous (i.v.) oxytocin elicited a rapid increase in work done by both the uterus and ampulla. Intramuscular (i.m.) PGF(2)alpha produced a delayed uterine response but little change in the ampulla; i.v. PGF(2)alpha produced a rapid response at both sites. Low plasma progesterone concentrations (< 0.5 ng/ml) were associated with a greater uterine and ampullary response to oxytocin and with an enhanced uterine response to PGF(2)alpha. However, the uterine tube response to intravenous PGF(2)alpha was greatest when plasma progesterone concentrations were high.