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.     

Direct probing of ion pair formation using a symmetric triangulenium dye

The 2,6,10-tris(dialkylamino)trioxatriangulenium dyes (ATOTA(+)) are highly stabilised cationic chromophores with D(3h) symmetry. The symmetry gives rise to a degeneracy of the main electronic transition. In low polarity solvents significant splitting of this degenerate transition is observed and assigned to ion pair formation. Ion pairing of the 2,6,10-tris(dioctylamino)trioxatriangulenium ion with Cl(-), BF(4)(-), PF(6)(-) and TRISPHAT anions was studied using absorption spectroscopy. A clear correlation is found between the size of the anion and the splitting of the ATOTA(+) transitions. In benzene the Cl(-) salt displays a splitting of 1955 cm(-1), while the salt of the much larger TRISPHAT ion has a splitting of 1543 cm(-1). TD-DFT calculations confirm the splitting of the states and provide a detailed insight into the electronic structure of the ion pairs. The different degree of splitting in different ion pairs is found to correlate with the magnitude of the electric field generated in each ion pair, thus leading to the conclusion that the effect seen is an internal Stark effect. By insertion of an amphiphilic derivative of the ATOTA(+) chromophore in an oriented lamellar liquid crystal, it was possible to resolve the two bands of the double peak spectrum and show their perpendicular orientation in the molecular framework, as predicted by the calculations.     

Towards systematic identification of Plasmodium essential genes by transposon shuttle mutagenesis

After the deciphering of the genome sequences of several Plasmodium species, efforts must turn to elucidating gene function and identifying essential gene products. However, random approaches are lacking and gene targeting is inefficient in Plasmodium. Here, we established shuttle transposon mutagenesis in Plasmodium berghei. We constructed a mini-Tn5 derivative that can transpose into parasite genes cloned in Escherichia coli, providing an efficient means of generating knockout fragments. A 10⁴-fold increase in frequencies of double-crossover homologous recombination in the parasite using a new electroporation technology permits to reproducibly generate pools of distinct mutants after transfection with mini-Tn5-interrupted sequences. The procedure opens the way to the systematic identification of essential genes in Plasmodium.     

PC12 and HEL cell communication in a reconstituted synapse analyzed with mathematical modeling

A synapse simulating model comprising of the nerve growth factor (NGF)-differentiated PC12 cells releasing neurotransmitter (NT) and sensor 92.1.7.human erythroleukemia (HEL) cells has been used for simulating the connection between neurons and target cells. A Ca(2+) elevation was observed in both cell types when the PC12 cells were challenged with nicotine. The response patterns of individual cell were subsequently analyzed mathematically. The Ca(2+) signals of the PC12 cells were described by an equation representing a simple bi-exponential function. The NT-noradrenaline discharged by the PC12 cells in response to nicotine caused heterogeneous secondary Ca(2+) elevations in the HEL cells after a certain delay. Model fitting of this response disclosed slow "hidden" oscillations and heterogeneous secondary Ca(2+) signals could be grouped on the basis of the oscillation frequency. As determined in control experiments with noradrenaline (NA), the value of oscillation frequency also revealed a good correlation with the NT concentration.     

Substitution of conserved glycine residue by alanine in natural and synthetic neuropeptide ligands causes partial agonism at the stomoxytachykinin receptor

A few naturally occurring insect tachykinin-related peptides, such as stomoxytachykinin (Stc-TK), contain an Ala-residue instead of the highly conserved Gly-residue that is present in most other members of this peptide family. Stc-TK is a potent, partial agonist of the stable fly (Stomoxys calcitrans) tachykinin receptor, STKR. By means of synthetic analogues, the Gly/Ala exchange, representing the addition of a single methyl group in the active core region of these peptides, was shown to be fully responsible for the generation of this partial agonism, which was also accompanied by an increase in agonistic potency. Surprisingly, this Ala-dependent reduction in maximal response levels was only observed for the agonist-induced cellular calcium rise. Stomoxytachykinin, Stc-TK, did not display partial agonism for the STKR-mediated cyclic AMP response. A possible explanation for this differential partial agonism is that the Gly-containing and Ala-replaced peptides recognize and stabilize active receptor conformations that differ in their functional coupling efficacies towards these response pathways. Drosotachykinins, Drm-TK, tachykinin-like peptides encoded in the fruit fly genome, were shown to be STKR-agonists. Interestingly, one of these peptides, which contains an Ala-residue instead of the conserved Gly-residue, also proved to be a potent, partial agonist for STKR.     

Gene expression profiles and genetic damage in benzo(a)pyrene diol epoxide-exposed TK6 cells

Microarray analysis is a powerful tool to identify the biological effects of drugs or chemicals on cellular gene expression. In this study, we compare the relationships between traditional measures of genetic toxicology and mutagen-induced alterations in gene expression profiles. TK6 cells were incubated with 0.01, 0.1, or 1.0 microM +/-anti-benzo(a)pyrene-trans-7,8-dihydrodiol-9,10-epoxide (BPDE) for 4 h and then cultured for an additional 20 h. Aliquots of the exposed cells were removed at 4 and 24 h in order to quantify DNA adduct levels by 32P post-labeling and measure cell viability by cloning efficiency and flow cytometry. Gene expression profiles were developed by extracting total RNA from the control and exposed cells at 4 and 24 h, labeling with Cy3 or Cy5 and hybridizing to a human 350 gene array. Mutant frequencies in the Thymidine Kinase and Hypoxanthine Phosphoribosyl Transferase genes were also determined. The 10alpha-(deoxyguanosin-N(2)-yl)-7alpha,8beta,9beta-trihydroxy-7,8,9,10-tetrahydrobenzo(a)pyrene (dG-N(2)-BPDE) adduct increased as a function of dose and was the only adduct identified. A dose-related decrease in cell viability was evident at 24 h, but not at 4 h. Cell death occurred by apoptosis. At 4 h, analysis of the gene expression profiles revealed that Glutathione Peroxidase and Gadd45 were consistently upregulated (greater than 1.5-fold and significantly (P < 0.001) greater than the control in two experiments) in response to 1.0 microM BPDE exposure. Fifteen genes were consistently down-regulated (less than 0.67-fold and significantly (P < 0.001) lower than the control in two experiments) at 4 h in cultures exposed to 1.0 microM BPDE. Genes with altered expression at 4 h included genes important in the progression of the cell-cycle and those that inhibit apoptosis. At 24 h post-exposure, 16 genes, involved in cell-cycle control, detoxification, and apoptosis were consistently upregulated; 10 genes were repressed in cultures exposed to the high dose of BPDE. Real-time quantitative PCR confirmed the differential expression of selected genes. These data suggest that changes in gene expression will help to identify effects of drugs and chemicals on molecular pathways in cells, and will provide useful information about the molecular responses associated with DNA damage. Of the endpoints evaluated, DNA adduct formation was the most sensitive indicator of DNA damage. DNA adduct formation was clearly evident at low doses, but the number of genes with significantly altered expression (P < 0.001) was minimal. Alterations in gene expression were more robust at doses associated with cellular toxicity and induction of mutations.     

The essential fatty acid status in phenylketonuria patients under treatment

Phenylketonuric patients are on a special diet that lacks certain essential fatty acids. This study evaluates the essential fatty acid status of a group of phenylketonuric patients in the Netherlands undergoing dietary treatment. To this end, the essential fatty acid status of nine phenylketonuria patients was studied. On the basis of age and gender, two control subjects were selected for each patient. The essential fatty acid composition of duplicate food portions and the essential fatty acid status of plasma and erythrocytes were analyzed. Phenylketonuria subjects had a different essential fatty acid profile from their peers, especially concerning the n-3 fatty acids. N-6 and n-3 fatty long-chain polyenes were hardly consumed by phenylketonuria subjects, in contrast to the control subjects. Linoleic acid, on the other hand, was consumed in significantly higher amounts by phenylketonuria subjects and made up about 40% of their daily fat consumption. The essential fatty acid consumption pattern of the phenylketonuria subjects is mirrored by the essential fatty acid concentrations in blood. The essential fatty acid status of the phenylketonuric diet should be improved in order to prevent deficiency in n-3 fatty acids.     

Lipoteichoic acid is an important microbe-associated molecular pattern of Lactobacillus rhamnosus GG

Background

Receptors with a single transmembrane (TM) domain are essential for the signal transduction across the cell membrane. NMR spectroscopy is a powerful tool to study structure of the single TM domain. The expression and purification of a TM domain in Escherichia coli (E.coli) is challenging due to its small molecular weight. Although ketosteroid isomerase (KSI) is a commonly used affinity tag for expression and purification of short peptides, KSI tag needs to be removed with the toxic reagent cyanogen bromide (CNBr).

Result

The purification of the TM domain of p75 neurotrophin receptor using a KSI tag with the introduction of a thrombin cleavage site is described herein. The recombinant fusion protein was refolded into micelles and was cleaved with thrombin. Studies showed that purified protein could be used for structural study using NMR spectroscopy.

Conclusions

These results provide another strategy for obtaining a single TM domain for structural studies without using toxic chemical digestion or acid to remove the fusion tag. The purified TM domain of p75 neurotrophin receptor will be useful for structural studies.