Best Linear Unbiased Prediction (BLUP) for regional yield trials: a comparison to additive main effects and multiplicative interaction (AMMI) analysis

Multilocation trials are often used to analyse the adaptability of genotypes in different environments and to find for each environment the genotype that is best adapted; i.e. that is highest yielding in that environment. For this purpose, it is of interest to obtain a reliable estimate of the mean yield of a cultivar in a given environment. This article compares two different statistical estimation procedures for this task: the Additive Main Effects and Multiplicative Interaction (AMMI) analysis and Best Linear Unbiased Prediction (BLUP). A modification of a cross validation procedure commonly used with AMMI is suggested for trials that are laid out as a randomized complete block design. The use of these procedure is exemplified using five faba bean datasets from German registration trails. BLUP was found to outperform AMMI in four of five faba bean datasets.     

Decoupling of impact factors reveals the response of German winter wheat yields to climatic changes

Yield development of agricultural crops over time is not merely the result of genetic and agronomic factors, but also the outcome of a complex interaction between climatic and site‐specific soil conditions. However, the influence of past climatic changes on yield trends remains unclear, particularly under consideration of different soil conditions. In this study, we determine the effects of single agrometeorological factors on the evolution of German winter wheat yields between 1958 and 2015 from 298 published nitrogen (N)‐fertilization experiments. For this purpose, we separate climatic from genetic and agronomic yield effects using linear mixed effect models and estimate the climatic influence based on a coefficient of determination for these models. We found earlier occurrence of wheat growth stages, and shortened development phases except for the phase of stem elongation. Agrometeorological factors are defined as climate covariates related to the growth of winter wheat. Our results indicate a general and strong effect of agroclimatic changes on yield development, in particular due to increasing mean temperatures and heat stress events during the grain‐filling period. Except for heat stress days with more than 31°C, yields at sites with higher yield potential were less prone to adverse weather effects than at sites with lower yield potential. Our data furthermore reveal that a potential yield levelling, as found for many West‐European countries, predominantly occurred at sites with relatively low yield potential and about one decade earlier (mid‐1980s) compared to averaged yield data for the whole of Germany. Interestingly, effects related to high precipitation events were less relevant than temperature‐related effects and became relevant particularly during the vegetative growth phase. Overall, this study emphasizes the sensitivity of yield productivity to past climatic conditions, under consideration of regional differences, and underlines the necessity of finding adaptation strategies for food production under ongoing and expected climate change.     

Estimates of dietary threshold levels for crude protein and amino acids to obtain plateau values of apparent ileal crude protein and amino acid digestibilities in newly weaned pigs

This study was conducted to estimate dietary threshold levels for crude protein (CP) and amino acids (AA) at which apparent ileal digestibilities (AID) of CP and AA in assay diets for newly weaned pigs reach plateaus. A total of 14 (12 + 2 for replacement) three-week old barrows were fitted with simple T cannulas at the distal ileum. Corn starch-based diets containing six graded levels of CP from casein, 90, 155, 220, 285, 350, or 415 g CP/kg assay diet (as-fed), were formulated. At 28 days of age, the pigs were randomly allocated to the six dietary treatments with two pigs per CP level in four weekly repeated measurement periods. They were fed twice daily a total of 30 g (as-fed) per kg of individual body weight at 8:00 and 20:00 h. The dietary CP and AA levels affected AID of CP and most AA (p = 0.005 to p = 0.040) in the assay diets. The AID of CP and AA were higher at 155 and 220 compared to 90 g CP/kg assay diet (p < 0.001 to p = 0.047). Initially, the AID of CP and AA increased sharply then gradually reached, at individual break points (p < 0.001 to p = 0.047), their plateaus (p < 0.001), which did not change up to dietary CP levels of 415 g/kg assay diet and the AID values became independent of the dietary AA levels. The piglets' capacity to digest CP and absorb AA was not limiting under these experimental conditions. There was no effect of age on AID of CP and AA (p = 0.056 to p = 0.899) except for a linear increase (p = 0.045) in AID of glycine from Period 1 to 4. Segmented quadratic with plateau relationships between the AID of CP and AA and their dietary contents were fitted for CP and each AA. The lower end points of 95% confidence intervals of the plateau AID values were defined to represent the initial plateau AID. The dietary CP and AA contents, corresponding to the initial plateau AID values, represent the dietary threshold levels. For CP and the indispensable amino acids, the plateau AID [%] and the dietary threshold levels [g/kg DM], respectively, in casein were: CP, 94.2 and 176; arginine, 95.1 and 7; histidine, 96.0 and 5; isoleucine, 96.4 and 8; leucine, 96.8 and 16; lysine, 96.8 and 12; methionine, 97.9 and 5; phenylalanine, 96.2 and 10; threonine, 93.4 and 9; tryptophan, 94.3 and 2; valine, 95.9 and 11. For the determination of plateau AID in piglets, the crude protein and amino acid contents in the assay diets should meet or exceed the corresponding dietary threshold levels.     

Background correction of two-colour cDNA microarray data using spatial smoothing methods

The analysis of two-colour cDNA microarray data usually involves subtracting background values from foreground values prior to normalization and further analysis. This approach has the advantage of reducing bias and the disadvantage of blowing up the variance of lower abundant spots. Whenever background subtraction is considered, it implicitly assumes locally constant background values. In practice, this assumption is often not met, which casts doubts on the usefulness of simple background subtraction. In order to improve background correction, we propose local background smoothing within the pre-processing pipeline of cDNA microarray data prior to background correction. For this purpose, we employ a geostatistical framework with ordinary kriging using both isotropic and anisotropic models of spatial correlation and 2-D locally weighted regression. We show that application of local background smoothing prior to background correction is beneficial in comparison to using raw background estimates. This is done using data of a self-versus-self experiment in Arabidopsis where subsets of differentially expressed genes were simulated. Using locally smoothed background values in conjunction with existing background correction methods increases the power, increases the accuracy and decreases the number of false positive results.     

Molecular marker assisted broadening of the Central European heterotic groups in rye with Eastern European germplasm

Broadening the genetic base of heterotic pools is a key to ensure continued genetic gains in hybrid breeding and extend hybrid cultivation to new areas. In the present study, two Central European heterotic pools (Carsten and Petkus) and five Eastern European open-pollinated varieties (OPVs, Pop-1 to Pop-5) were studied with the objectives to (1) investigate the genetic diversity in OPVs and the heterotic pools using molecular and field data, (2) evaluate the molecular diversity among OPVs, (3) examine the combining ability for grain yield of the OPVs when crossed with testers in field trials, and (4) develop a strategy for targeted introgression of OPV germplasm into the heterotic pools. In total, 610 S₀ plants, 347 from OPVs and 263 from heterotic pools, were developed. Clones of the S₀ plants of OPVs were crossed with two testers belonging to each heterotic pool, while clones of heterotic pools were crossed with only the opposite tester. Testcrosses were evaluated for grain yield in multi-location trials. In addition, 589 S₀ plants were fingerprinted with 30 SSR markers. The data revealed that the Carsten pool has a narrow genetic base and should be the primary target for broadening the established heterotic pattern. Mean and genetic variance suggested that Pop-2 and Pop-4 are good candidates for introgression in Petkus pool and Pop-5 in Carsten pool. Nevertheless, introgression of Pop-5 in Carsten could reduce the genetic diversity between heterotic pools. Therefore, we suggest that either selected plants of Pop-5 should be introgressed or more Eastern European germplasm should be fingerprinted and field evaluated to identify promising germplasm for broadening the established heterotic pattern.     

Nonadditive protein accumulation patterns in Maize (Zea mays L.) hybrids during embryo development

Heterosis describes the superior performance of heterozygous F(1)-hybrid plants compared to their homozygous parental inbred lines. In the present study, heterosis was detected for length, weight, and the time point of seminal root primordia initiation in maize (Zea mays L.) embryos of the reciprocal F(1)-hybrids UH005xUH250 and UH250xUH005. A two-dimensional gel electrophoresis (2-DE) proteome survey of the most abundant proteins of the reciprocal hybrids and their parental inbred lines 25 and 35 days after pollination revealed that 141 of 597 detected proteins (24%) exhibited nonadditive accumulation in at least one hybrid. Approximately 44% of all nonadditively accumulated proteins displayed an expression pattern that was not distinguishable from the low parent value. Electrospray ionization-tandem mass spectrometry (ESI-MS/MS) analyses and subsequent functional classification of the 141 proteins revealed that development, protein metabolism, redox-regulation, glycolysis, and amino acid metabolism were the most prominent functional classes among nonadditively accumulated proteins. In 35-day-old embryos of the hybrid UH250xUH005, a significant up-regulation of enzymes related to glucose metabolism which often exceeded the best parent values was observed. A comparison of nonadditive protein accumulation between rice and maize embryo data sets revealed a significant overlap of nonadditively accumulated proteins suggesting conserved organ- or tissue-specific regulatory mechanisms in monocots related to heterosis.     

Temporal changes in allele frequencies in two European F2 flint maize populations under modified recurrent full-sib selection

Selection and random genetic drift are the two main forces affecting the selection response of recurrent selection (RS) programs by changes in allele frequencies. Therefore, detailed knowledge on allele frequency changes attributable to these forces is of fundamental importance for assessing RS programs. The objectives of our study were to (1) estimate the number, position, and genetic effect of quantitative trait loci (QTL) for selection index and its components in the base populations, (2) determine changes in allele frequencies of QTL regions due to the effects of random genetic drift and selection, and (3) predict allele frequency changes by using QTL results and compare these predictions with observed values. We performed QTL analyses, based on restriction fragment length polymorphisms (RFLPs) and simple sequence repeats (SSRs), in 274 F_2:3 lines of cross KW1265 × D146 (A × B) and 133 F_3:4 lines of cross D145 × KW1292 (C × D) originating from two European flint maize populations. Four (A × B) and seven (C × D) cycles of RS were analyzed with SSRs for significant allele frequency changes due to selection. Several QTL regions for selection index were detected with simple and composite interval mapping. In some of them, flanking markers showed a significant allele frequency change after the first and the final selection cycles. The correlation between observed and predicted allele frequencies was significant only in A × B. We attribute these observations mainly to (1) the high dependence of the power of QTL detection on the population size and (2) the occurrence of undetectable QTL in repulsion phase. Assessment of allele frequency changes in RS programs can be used to detect marker alleles linked to QTL regions under selection pressure. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Potential causes of linkage disequilibrium in a European maize breeding program investigated with computer simulations

Knowledge about the forces generating and conserving linkage disequilibrium (LD) is important for drawing conclusions about the prospects and limitations of association mapping. The objectives of our research were to examine the importance of (1) selection, (2) mutation, and (3) genetic drift for generating LD in a typical maize breeding program. We conducted computer simulations based on genotypic data of Central European maize open-pollinated varieties which have played an important role as founders of the European flint heterotic group. The breeding scheme and the dimensioning underlying our simulations reflect essentially the maize breeding program of the University of Hohenheim. Results suggested that in a plant breeding program of the examined dimension and breeding scheme, genetic drift and selection are major forces generating LD. The currently used population-based association mapping tests do not explicitly correct for LD caused by these two forces. Therefore, increased type I error rates are expected if these tests are applied to plant breeding populations. As a consequence, we recommend to use family-based association tests for association mapping approaches in plant breeding populations.