家畜基因组选择研究进展

近年来,随着基因芯片技术的发展与育种技术的进步,动植物的基因组选择成为研究热点。在家畜育种中,基因组选择凭借其准确性高、世代间隔短和育种成本低等优势被应用于各种经济动物的种畜选择中。本文详细介绍了基因分型技术和基因组育种值估计方法(最小二乘法、RR-BLUP法、GBLUP法、ssGBLUP法、贝叶斯A法、贝叶斯B法等),并对这些育种方法选用的标记范围、准确性以及计算速度进行了比较,总结了我国和其他国家基因组选择在种畜选择中的应用情况及存在的问题,展望了目前国内外在基因组选择上的最新研究动态及进展,以期为其他育种工作者进一步了解基因组选择提供参考。     

Selection for increased body length in Subantarctic fur seals on Amsterdam Island

The traditional distinction between ecological and evolutionary times is eroding, calling for tighter links between ecology and evolution. An example of such a brigde between the two disciplines is the so-called 'animal model', a methodology initially developed by animal breeders, which has become very popular among ecologists studying contemporary microevolution. Using a Bayesian multi-trait 'animal model', we investigated the quantitative genetics of body size, a fitness-related trait, in Subantarctic fur seals (Arctocephalus tropicalis) breeding on Amsterdam Island, Southern Ocean. Our approach jointly modelled the growth and selection processes at work in this population. Body length is heritable for both sexes, and females are under selection for increased body length in this population. We strongly suspect the peculiar ecological context of impoverished, suitable prey availability exacerbated by density-dependence phenomena to be an important selective agent on females breeding on Amsterdam Island.     

Bayesian Multi-Trait Analysis Reveals a Useful Tool to Increase Oil Concentration and to Decrease Toxicity in Jatropha curcas L.

The biggest challenge for jatropha breeding is to identify superior genotypes that present high seed yield and seed oil content with reduced toxicity levels. Therefore, the objective of this study was to estimate genetic parameters for three important traits (weight of 100 seed, oil seed content, and phorbol ester concentration), and to select superior genotypes to be used as progenitors in jatropha breeding. Additionally, the genotypic values and the genetic parameters estimated under the Bayesian multi-trait approach were used to evaluate different selection indices scenarios of 179 half-sib families. Three different scenarios and economic weights were considered. It was possible to simultaneously reduce toxicity and increase seed oil content and weight of 100 seed by using index selection based on genotypic value estimated by the Bayesian multi-trait approach. Indeed, we identified two families that present these characteristics by evaluating genetic diversity using the Ward clustering method, which suggested nine homogenous clusters. Future researches must integrate the Bayesian multi-trait methods with realized relationship matrix, aiming to build accurate selection indices models.     

Multivariate selection shapes environment-dependent variation in the clonal morphology of a red seaweed

Within-individual strategies of variation (e.g., phenotypic plasticity) are particularly relevant to modular organisms, in which ramets of the same genetic individual may encounter diverse environments imposing diverse patterns of selection. Hence, measuring selection in heterogeneous environments is essential to understanding whether environment-dependent phenotypic change enhances the fitness of modular individuals. In sublittoral marine habitats, competition for light and space among modular taxa generates extreme patchiness in resource availability. Little is known, however, of the potential for plasticity within individuals to arise from spatially-variable selection in such systems. We tested whether plasticity enhances genet-level fitness in Asparagopsis armata, a clonal seaweed in which correlated traits mediate morphological responses to variation in light. Using the capacity for rapid, clonal growth to measure fitness, we identified aspects of ramet morphology targeted by selection in two contrasting light environments and compared patterns of selection across environments. We found that directional selection on single traits, coupled with linear and nonlinear selection on multi-trait interactions, shape ramet morphology within environments and favor different phenotypes in each. Evidence of environment-dependent, multivariate selection on correlated traits is novel for any marine modular organism and demonstrates that seaweeds, such as A. armata, may potentially adapt to environmental heterogeneity via plasticity in clonal morphology.     

Multi-trait genomic prediction for nitrogen response indices in tropical maize hybrids

In maize breeding, genomic prediction may be an efficient tool for selecting single-crosses evaluated under abiotic stress conditions. In addition, a promising strategy is applying multiple-trait genomic prediction using selection indices (SIs), increasing genetics gains and reducing time per cycles. In this study, we aimed (i) to compare accuracy of single- and multi-trait genomic prediction (STGP; MTGP) in two maize datasets, (ii) to evaluate prediction of four selection indices that could contribute to the selection of tropical maize hybrids under contrasting nitrogen conditions, and (iii) to compare the use of linear (GBLUP) and nonlinear (RKHS/GK) kernels in STGP and MTGP analyses. For either single-trait GBLUP and RKHS analyses, the highest values obtained for accuracy were 0.40 and 0.41 using harmonic mean (HM), respectively. From multi-trait GBLUP and GK, using the combination of selection indices in MTGP seems to be suitable, increasing the accuracy. Adding grain yield and plant height in MTGP showed a slight improvement in accuracy compared to STGP. In general, there was a modest benefit of using single-trait RKHS and GK multi-trait, rather than GBLUP.     

Genetic Group Index on Multiple Traits for Homogeneous Covariance Structure

A multi-trait index for ranking different genetic groups/populations with homogeneous covariance structure is developed and illustrated with live data.     

Efficiency of multi-trait,indirect, and trait-assisted genomic selection for improvement of biomass sorghum

Key message

We compare genomic selection methods that use correlated traits to help predict biomass yield in sorghum, and find that trait-assisted genomic selection performs best.

Abstract

Genomic selection (GS) is usually performed on a single trait, but correlated traits can also help predict a focal trait through indirect or multi-trait GS. In this study, we use a pre-breeding population of biomass sorghum to compare strategies that use correlated traits to improve prediction of biomass yield, the focal trait. Correlated traits include moisture, plant height measured at monthly intervals between planting and harvest, and the area under the growth progress curve. In addition to single- and multi-trait direct and indirect GS, we test a new strategy called trait-assisted GS, in which correlated traits are used along with marker data in the validation population to predict a focal trait. Single-trait GS for biomass yield had a prediction accuracy of 0.40. Indirect GS performed best using area under the growth progress curve to predict biomass yield, with a prediction accuracy of 0.37, and did not differ from indirect multi-trait GS that also used moisture information. Multi-trait GS and single-trait GS yielded similar results, indicating that correlated traits did not improve prediction of biomass yield in a standard GS scenario. However, trait-assisted GS increased prediction accuracy by up to \(50\%\) when using plant height in both the training and validation populations to help predict yield in the validation population. Coincidence between selected genotypes in phenotypic and genomic selection was also highest in trait-assisted GS. Overall, these results suggest that trait-assisted GS can be an efficient strategy when correlated traits are obtained earlier or more inexpensively than a focal trait.

     

An interval estimation of expected response to selection

Summary A method is presented here for obtaining an interval estimate of expected response to selection based on results of a progeny test experiment. The structure of the constructed confidence limits is then examined for the influence of the numbers of lines and replicates on the precision of predicting the expected response to selection.     

Genetic analysis of milking ability in Lacaune dairy ewes

The milking ability of Lacaune ewes was characterised by derived traits of milk flow patterns, in an INRA experimental farm, from a divergent selection experiment in order to estimate the correlated effects of selection for protein and fat yields. The analysis of selected divergent line effects (involving 34 616 data and 1204 ewes) indicated an indirect improvement of milking traits (+17% for maximum milk flow and -10% for latency time) with a 25% increase in milk yield. Genetic parameters were estimated by multi-trait analysis with an animal model, on 751 primiparous ewes. The heritabilities of the traits expressed on an annual basis were high, especially for maximum flow (0.54) and for latency time (0.55). The heritabilities were intermediate for average flow (0.30), time at maximum flow (0.42) and phase of increasing flow (0.43), and low for the phase of decreasing flow (0.16) and the plateau of high flow (0.07). When considering test-day data, the heritabilities of maximum flow and latency time remained intermediate and stable throughout the lactation. Genetic correlations between milk yield and milking traits were all favourable, but latency time was less milk yield dependent (-0.22) than maximum flow (+0.46). It is concluded that the current dairy ewe selection based on milk solid yield is not antagonistic to milking ability.     

Economic evaluation of genomic selection in small ruminants: a sheep meat breeding program

Recent genomic evaluation studies using real data and predicting genetic gain by modeling breeding programs have reported moderate expected benefits from the replacement of classic selection schemes by genomic selection (GS) in small ruminants. The objectives of this study were to compare the cost, monetary genetic gain and economic efficiency of classic selection and GS schemes in the meat sheep industry. Deterministic methods were used to model selection based on multi-trait indices from a sheep meat breeding program. Decisional variables related to male selection candidates and progeny testing were optimized to maximize the annual monetary genetic gain (AMGG), that is, a weighted sum of meat and maternal traits annual genetic gains. For GS, a reference population of 2000 individuals was assumed and genomic information was available for evaluation of male candidates only. In the classic selection scheme, males breeding values were estimated from own and offspring phenotypes. In GS, different scenarios were considered, differing by the information used to select males (genomic only, genomic+own performance, genomic+offspring phenotypes). The results showed that all GS scenarios were associated with higher total variable costs than classic selection (if the cost of genotyping was 123 euros/animal). In terms of AMGG and economic returns, GS scenarios were found to be superior to classic selection only if genomic information was combined with their own meat phenotypes (GS-Pheno) or with their progeny test information. The predicted economic efficiency, defined as returns (proportional to number of expressions of AMGG in the nucleus and commercial flocks) minus total variable costs, showed that the best GS scenario (GS-Pheno) was up to 15% more efficient than classic selection. For all selection scenarios, optimization increased the overall AMGG, returns and economic efficiency. As a conclusion, our study shows that some forms of GS strategies are more advantageous than classic selection, provided that GS is already initiated (i.e. the initial reference population is available). Optimizing decisional variables of the classic selection scheme could be of greater benefit than including genomic information in optimized designs.     

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.     

Using a network framework to quantitatively select ecological indicators

Ecological indicators are often constructed as an integrated set to represent key information and characteristics of the ecosystem which are tightly linked to management objectives. As an effective tool, ecological indicators play an increasingly important role in ecosystem monitoring, assessment and management. Reasonable selection of an indicator is a prerequisite for effectively using it. A defined protocol with scientific rigor to select ecological indicators is imperative to solve the challenges in ecological indicator selection. This paper compares the Causal Network (CN) with the Ecological Hierarchy Network (EHN) as a framework to select ecological indicators. These frameworks are not exclusive but interdependent in constructing a network framework. Based on the network framework, a quantitative ecological indicator selection method is demonstrated through a theoretical example. In the selection process, the criteria and requirements considering the balance of science and utility are proposed and translated into quantitative constraints of a selection model. By resolving the model under a mathematical operation, the human arbitrary disturbance will be reduced and random selection minimized.     

Sharp simultaneous confidence intervals for the means of selected populations with application to microarray data analysis

Summary Simultaneous inference for a large number, N, of parameters is a challenge. In some situations, such as microarray experiments, researchers are only interested in making inference for the K parameters corresponding to the K most extreme estimates. Hence it seems important to construct simultaneous confidence intervals for these K parameters. The naïve simultaneous confidence intervals for the K means (applied directly without taking into account the selection) have low coverage probabilities. We take an empirical Bayes approach (or an approach based on the random effect model) to construct simultaneous confidence intervals with good coverage probabilities. For N= 10,000 and K= 100, typical for microarray data, our confidence intervals could be 77% shorter than the naïve K‐dimensional simultaneous intervals.     

Using path analysis to measure natural selection

We expand current methods for calculating selection coefficients using path analysis and demonstrate how to analyse nonlinear selection. While this incorporation is a straightforward extension of current procedures, the rules for combining these traits to calculate selection coefficients can be complex. We demonstrate our method with an analysis of selection in an experimental population of Arabidopsis thaliana consisting of 289 individuals. Multiple regression analyses found positive directional selection and positive nonlinear selection only for inflorescence height. In contrast, the path analyses also revealed positive directional selection for number of rosette leaves and positive nonlinear selection for leaf number and time of inflorescence initiation. These changes in conclusions came about because indirect selection was converted into direct selection with the change in causal structure. Path analysis has great promise for improving our understanding of natural selection but must be used with caution since coefficient estimates depend on the assumed causal structure.     

Selection of nest-site habitat by interior least terns in relation to sandbar construction

Federally endangered interior least terns (Sternula antillarum) nest on bare or sparsely vegetated sandbars on midcontinent river systems. Loss of nesting habitat has been implicated as a cause of population declines, and managing these habitats is a major initiative in population recovery. One such initiative involves construction of mid-channel sandbars on the Missouri River, where natural sandbar habitat has declined in quantity and quality since the late 1990s. We evaluated nest-site habitat selection by least terns on constructed and natural sandbars by comparing vegetation, substrate, and debris variables at nest sites (n = 798) and random points (n = 1,113) in bare or sparsely vegetated habitats. Our logistic regression models revealed that a broader suite of habitat features was important in nest-site selection on constructed than on natural sandbars. Odds ratios for habitat variables indicated that avoidance of habitat features was the dominant nest-site selection process on both sandbar types, with nesting terns being attracted to nest-site habitat features (gravel and debris) and avoiding vegetation only on constructed sandbars, and avoiding silt and leaf litter on both sandbar types. Despite the seemingly uniform nature of these habitats, our results suggest that a complex suite of habitat features influences nest-site choice by least terns. However, nest-site selection in this social, colonially nesting species may be influenced by other factors, including spatial arrangement of bare sand habitat, proximity to other least terns, and prior habitat occupancy by piping plovers (Charadrius melodus). We found that nest-site selection was sensitive to subtle variation in habitat features, suggesting that rigor in maintaining habitat condition will be necessary in managing sandbars for the benefit of least terns. Further, management strategies that reduce habitat features that are avoided by least terns may be the most beneficial to nesting least terns. © 2011 The Wildlife Society.     

Genetic evaluation and selection response for growth in meat-type quail through random regression models using B-spline functions and Legendre polynomials

The objective was to estimate (co)variance functions using random regression models (RRM) with Legendre polynomials, B-spline function and multi-trait models aimed at evaluating genetic parameters of growth traits in meat-type quail. A database containing the complete pedigree information of 7000 meat-type quail was utilized. The models included the fixed effects of contemporary group and generation. Direct additive genetic and permanent environmental effects, considered as random, were modeled using B-spline functions considering quadratic and cubic polynomials for each individual segment, and Legendre polynomials for age. Residual variances were grouped in four age classes. Direct additive genetic and permanent environmental effects were modeled using 2 to 4 segments and were modeled by Legendre polynomial with orders of fit ranging from 2 to 4. The model with quadratic B-spline adjustment, using four segments for direct additive genetic and permanent environmental effects, was the most appropriate and parsimonious to describe the covariance structure of the data. The RRM using Legendre polynomials presented an underestimation of the residual variance. Lesser heritability estimates were observed for multi-trait models in comparison with RRM for the evaluated ages. In general, the genetic correlations between measures of BW from hatching to 35 days of age decreased as the range between the evaluated ages increased. Genetic trend for BW was positive and significant along the selection generations. The genetic response to selection for BW in the evaluated ages presented greater values for RRM compared with multi-trait models. In summary, RRM using B-spline functions with four residual variance classes and segments were the best fit for genetic evaluation of growth traits in meat-type quail. In conclusion, RRM should be considered in genetic evaluation of breeding programs.     

Integration of genomic information into sport horse breeding programs for optimization of accuracy of selection

Reliable selection criteria are required for young riding horses to increase genetic gain by increasing accuracy of selection and decreasing generation intervals. In this study, selection strategies incorporating genomic breeding values (GEBVs) were evaluated. Relevant stages of selection in sport horse breeding programs were analyzed by applying selection index theory. Results in terms of accuracies of indices (r_TI) and relative selection response indicated that information on single nucleotide polymorphism (SNP) genotypes considerably increases the accuracy of breeding values estimated for young horses without own or progeny performance. In a first scenario, the correlation between the breeding value estimated from the SNP genotype and the true breeding value (= accuracy of GEBV) was fixed to a relatively low value of r_mg = 0.5. For a low heritability trait (h² = 0.15), and an index for a young horse based only on information from both parents, additional genomic information doubles r_TI from 0.27 to 0.54. Including the conventional information source ‘own performance’ into the before mentioned index, additional SNP information increases r_TI by 40%. Thus, particularly with regard to traits of low heritability, genomic information can provide a tool for well-founded selection decisions early in life. In a further approach, different sources of breeding values (e.g. GEBV and estimated breeding values (EBVs) from different countries) were combined into an overall index when altering accuracies of EBVs and correlations between traits. In summary, we showed that genomic selection strategies have the potential to contribute to a substantial reduction in generation intervals in horse breeding programs.     

Genomic selection for feed efficiency in dairy cattle

Feed is a major component of variable costs associated with dairy systems and is therefore an important consideration for breeding objectives. As a result, measures of feed efficiency are becoming popular traits for genetic analyses. Already, several countries account for feed efficiency in their breeding objectives by approximating the amount of energy required for milk production, maintenance, etc. However, variation in actual feed intake is currently not captured in dairy selection objectives, although this could be possible by evaluating traits such as residual feed intake (RFI), defined as the difference between actual and predicted feed (or energy) intake. As feed intake is expensive to accurately measure on large numbers of cows, phenotypes derived from it are obvious candidates for genomic selection provided that: (1) the trait is heritable; (2) the reliability of genomic predictions are acceptable to those using the breeding values; and (3) if breeding values are estimated for heifers, rather than cows then the heifer and cow traits need to be correlated. The accuracy of genomic prediction of dry matter intake (DMI) and RFI has been estimated to be around 0.4 in beef and dairy cattle studies. There are opportunities to increase the accuracy of prediction, for example, pooling data from three research herds (in Australia and Europe) has been shown to increase the accuracy of genomic prediction of DMI from 0.33 within country to 0.35 using a three-country reference population. Before including RFI as a selection objective, genetic correlations with other traits need to be estimated. Weak unfavourable genetic correlations between RFI and fertility have been published. This could be because RFI is mathematically similar to the calculation of energy balance and failure to account for mobilisation of body reserves correctly may result in selection for a trait that is similar to selecting for reduced (or negative) energy balance. So, if RFI is to become a selection objective, then including it in an overall multi-trait selection index where the breeding objective is net profit is sensible, as this would allow genetic correlations with other traits to be properly accounted for. If genetic parameters are accurately estimated then RFI is a logical breeding objective. If there is uncertainty in these, then DMI may be preferable.     

Kinship rivalry does not trigger specific allocation strategies in Lupinus angustifolius

Background and Aims

Research on the ability of plants to recognize kin and modify plant development to ameliorate competition with coexisting relatives is an area of very active current exploration. Empirical evidence, however, is insufficient to provide a sound picture of this phenomenon.

Methods

An experiment was designed to assess multi-trait phenotypic expression in response to competition with conspecifics of varied degrees of genealogical relatedness. Groups of siblings, cousins and strangers of Lupinus angustifolius were set in competition in a pots assay. Several whole-plant and organ-level traits, directly related to competition for above- and below-ground resources, were measured. In addition, group-level root proliferation was measured as a key response trait to relatedness to neighbours, as identified in previous work.

Key Results

No major significant phenotypic differences were found between individuals and groups that could be assigned to the gradient of relatedness used here. This occurred in univariate models, and also when multi-trait interactions were evaluated through multi-group comparisons of Structural Equation Models. Root proliferation was higher in phenotypically more heterogeneous groups, but phenotypic heterogeneity was independent of the relatedness treatments of the experiment, and root proliferation was alike in the neighbourhoods of siblings, cousins and strangers.

Conclusions

In contrast to recent findings in other species, genealogical relatedness to competing neighbours has a negligible impact on the phenotypic expression of individuals and groups of L. angustifolius. This suggests that kin recognition needs further exploration to assess its generality, the ecological scenarios where it might have been favoured or penalized by natural selection, and its preponderance in different plant lineages.