Prediction accuracies for growth and wood attributes of interior spruce in space using genotyping-by-sequencing

Background

Genomic selection (GS) in forestry can substantially reduce the length of breeding cycle and increase gain per unit time through early selection and greater selection intensity, particularly for traits of low heritability and late expression. Affordable next-generation sequencing technologies made it possible to genotype large numbers of trees at a reasonable cost.

Results

Genotyping-by-sequencing was used to genotype 1,126 Interior spruce trees representing 25 open-pollinated families planted over three sites in British Columbia, Canada. Four imputation algorithms were compared (mean value (MI), singular value decomposition (SVD), expectation maximization (EM), and a newly derived, family-based k-nearest neighbor (kNN-Fam)). Trees were phenotyped for several yield and wood attributes. Single- and multi-site GS prediction models were developed using the Ridge Regression Best Linear Unbiased Predictor (RR-BLUP) and the Generalized Ridge Regression (GRR) to test different assumption about trait architecture. Finally, using PCA, multi-trait GS prediction models were developed. The EM and kNN-Fam imputation methods were superior for 30 and 60% missing data, respectively. The RR-BLUP GS prediction model produced better accuracies than the GRR indicating that the genetic architecture for these traits is complex. GS prediction accuracies for multi-site were high and better than those of single-sites while multi-site predictability produced the lowest accuracies reflecting type-b genetic correlations and deemed unreliable. The incorporation of genomic information in quantitative genetics analyses produced more realistic heritability estimates as half-sib pedigree tended to inflate the additive genetic variance and subsequently both heritability and gain estimates. Principle component scores as representatives of multi-trait GS prediction models produced surprising results where negatively correlated traits could be concurrently selected for using PCA2 and PCA3.

Conclusions

The application of GS to open-pollinated family testing, the simplest form of tree improvement evaluation methods, was proven to be effective. Prediction accuracies obtained for all traits greatly support the integration of GS in tree breeding. While the within-site GS prediction accuracies were high, the results clearly indicate that single-site GS models ability to predict other sites are unreliable supporting the utilization of multi-site approach. Principle component scores provided an opportunity for the concurrent selection of traits with different phenotypic optima.

Electronic supplementary material

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

Activation and inhibition of cyclin-dependent kinase-2 by phosphorylation; a molecular dynamics study reveals the functional importance of the glycine-rich loop

Nanoseconds long molecular dynamics (MD) trajectories of differently active complexes of human cyclin-dependent kinase 2 (inactive CDK2/ATP, semiactive CDK2/Cyclin A/ATP, fully active pT160-CDK2/Cyclin A/ATP, inhibited pT14-; pY15-; and pT14,pY15,pT160-CDK2/Cyclin A/ATP) were compared. The MD simulations results of CDK2 inhibition by phosphorylation at T14 and/or Y15 sites provide insight into the structural aspects of CDK2 deactivation. The inhibitory sites are localized in the glycine-rich loop (G-loop) positioned opposite the activation T-loop. Phosphorylation of T14 and both inhibitory sites T14 and Y15 together causes ATP misalignment for phosphorylation and G-loop conformational change. This conformational change leads to the opening of the CDK2 substrate binding box. The phosphorylated Y15 residue negatively affects substrate binding or its correct alignment for ATP terminal phospho-group transfer to the CDK2 substrate. The MD simulations of the CDK2 activation process provide results in agreement with previous X-ray data.     

Influence of the acetylcholinesterase active site protonation on omega loop and active site dynamics

Existence of alternative entrances in acetylcholinesterase (AChE) could explain the contrast between the very high AChE catalytic efficiency and the narrow and long access path to the active site revealed by X-ray crystallography. Alternative entrances could facilitate diffusion of the reaction products or at least water and ions from the active site. Previous molecular dynamics simulations identified side door and back door as the most probable alternative entrances. The simulations of non-inhibited AChE suggested that the back door opening events occur only rarely (0.8% of the time in the 10ns trajectory). Here we present a molecular dynamics simulation of non-inhibited AChE, where the back door opening appears much more often (14% of the time in the 12ns trajectory) and where the side door opening was observed quite frequently (78% of trajectory time). We also present molecular dynamics, where the back door does not open at all, or where large conformational changes of the AChE omega loop occur together with alternative passage opening events. All these differences in AChE dynamical behavior are caused by different protonation states of two glutamate residues located on bottom of the active site gorge (Glu202 and G450 in Mus musculus AChE). Our results confirm the results of previous molecular dynamics simulations, expand the view and suggest the probable reasons for the overall conformational behavior of AChE omega loop.     

Manipulation of division symmetry and developmental fate in cultures of potato microspores

The effect of media composition on microspore culture was investigated in one tetraploid and two diploid potatoes. The viability of microspores isolated from 4.5 to 5 mm buds was in the range of 33 to 52%. In media for anther culture, microspores showed no further development and lost viability within 2 days. In M1 medium containing mineral components, sucrose, uridine, cytidine, myo-inositol, glutamine and lactalbumin hydrolysate, 18 to 37% of microspores underwent mitosis within 14 days. Up to 95% of the divisions were symmetric and produced equal nuclei. Some symmetrically divided microspores eventually produced structures with 3 to 10 nuclei. The proportion of the total microspore population producing multinuclear structures reached 9% in diploid clones responsive to anther culture and 1 to 2% in recalcitrant cv. Borka. Symmetric mitoses in M1 medium were induced in the presence of glutamine and lactalbumin hydrolysate. Nucleosides and myo-inositol had no effect on microspore division. In the absence of all organic components except sucrose, most mitoses were asymmetric, formation of multinuclear structures was reduced and most pollen accumulated starch indicative of gametophytic fate. In complete M1 medium, starch accumulation was suppressed. Suppression also occurred in asymmetrically divided microspores, indicating a direct inhibition of pollen development independent of the mode of microspore division. This inhibitory effect of M1 medium might present a stress which triggers the induction of symmetric microspore division and subsequent formation of multinuclear structures. This revised version was published online in June 2006 with corrections to the Cover Date.     

Synthesis of hexose-related imidazolidinones: Novel glycation products in the Maillard reaction

Carbohydrate-peptide esters which mimic the reactivity of sugar 6-phosphates in nonenzymatic glycations were used as model compounds for the study of the Maillard reaction in vitro. We found that intramolecular cyclization of the monosaccharide ester in which the sugar moiety (D-glucose or D-galactose) is linked, through the C-6 hydroxy group, to the C-terminal carboxy group of the endogenous opioid pentapeptide leucine-enkephalin, in methanol as the solvent, resulted in the formation of imidazolidinone diastereoisomers having cis or trans relative geometry of the substituents at the imidazolidinone ring moiety. The diastereoisomeric imidazolidinones were separated and each transformed by hydrolysis into the corresponding D-gluco- and D-galacto-related imidazolidinone products of leucine-enkephalin. Along with the previous evidence that, from the same sugar-peptide esters by changing the reaction conditions Amadori rearrangement products could be obtained [Horvat et al. (1998) J Chem Soc Perkin Trans 1:99–13], the presented results point to the possibility that similar carbohydrate-related imidazolidinones may also be generated in the early stage of the Maillard reaction in vivo.     

Binding and orientation of conantokins in PL vesicles and aligned PL multilayers

The association of a ligand with its cognate cell surface receptor can be facilitated by interactions between the ligand and the lipid phase of the cell membrane. With respect to the N-methyl-D-aspartate receptor (NMDAR), we have previously established a low affinity, nonreceptor-mediated interaction of the peptidic conantokins with synaptic membranes in conjunction with a high affinity binding to the NMDARs present therein [Klein, R. C., Prorok, M., and Castellino, F. J. (2003) J. Pept. Res. 61, 307-317]. In the current study, several techniques including size-exclusion chromatography, circular dichroism, fluorescence, and NMR spectroscopies were used to investigate the binding, conformation, and orientation of conantokins and their variants to a variety of phospholipid (PL) vesicles and multilayers. We have found that conantokins bind to PLs and that the effectors Ca(2+) and spermine slightly increase this binding ability. The conantokins preserve a high degree of helical conformation when bound to vesicles in the presence of Ca(2+). In the absence of Ca(2+), only conantokin-G (con-G) manifests an increase in conantokin helicity with increasing vesicle concentration. In solution, the conantokins appear to be localized at the headgroup of vesicles and do not insert into the hydrophobic core of the bilayer. On aligned PL films, the helical axis of the conantokins can either reside normal to the membrane surface or partition in a parallel orientation, depending on the nature of the conantokins and the PLs used. These orientation preferences may be conjoined with the biological activities of the conantokins.     

Genomic occurrence of microsatellites containing integral and non-integral repeat numbers

We calculated occurrences of all dinucleotide and trinucleotide microsatellites in the human, mouse, and yeast genomes. The microsatellites were considered separately not only according to the repeated dinucleotide or trinucleotide and the microsatellite length but also according to the starting/terminal nucleotide. The analysis showed that dramatically non-equal amounts occurred in the human genome of microsatellites that differed only by the terminal nucleotides. For example, the 23-mer (TTG)(7)TT occurs 635 times in the human genome whereas (GTT)(7)GT is present only three times in the human genome though the two 23-mers share a 22 nucleotide sequence. The dramatically non-equal occurrences of microsatellites differing only by the terminal nucleotides are observed for most dinucleotide and trinucleotide microsatellites and in all analyzed genomes. We suppose that the strikingly non-equal genomic occurrences of these closely related microsatellites originate from conformational properties of DNA.     

Improved gas chromatographic-mass spectrometric determination of the N-methylcarbamoyl adduct at the N-terminal valine of globin,a metabolic product of the solvent N,N-dimethylformamide

A sensitive method for determination of the N-methylcarbamoyl adduct at the N-terminal valine of globin, a new metabolic product of the industrial solvent N,N-dimethylformamide (DMF), has been developed and validated. The method includes conversion of the adduct by the Edman degradation to 3-methyl-5-isopropylhydantoin (MVH), which is followed by optimized gas chromatographic analysis with mass spectrometric detection at m/z 114. The recovery of MVH from terminal N-methylcarbamoylvaline was determined using a model dipeptide to be 90%. Calibration of the method is done with MVH, employing 3-methyl-5-isobutylhydantoin as the internal standard. The limit of detection is 0.2 nmol MVH/g globin when a 100-mg sample is used. Within- and between-day precision is 4-10%. The method has been used to determine the background levels of MVH in unexposed subjects. Further, toxicokinetic studies in volunteers laid the grounds for setting the reference value for biological monitoring of occupational exposure to DMF.