Estimated breeding values and association mapping for persistency and total milk yield using natural cubic smoothing splines

BackgroundFor dairy producers, a reliable description of lactation curves is a valuable tool for management and selection. From a breeding and production viewpoint, milk yield persistency and total milk yield are important traits. Understanding the genetic drivers for the phenotypic variation of both these traits could provide a means for improving these traits in commercial production.MethodsIt has been shown that Natural Cubic Smoothing Splines (NCSS) can model the features of lactation curves with greater flexibility than the traditional parametric methods. NCSS were used to model the sire effect on the lactation curves of cows. The sire solutions for persistency and total milk yield were derived using NCSS and a whole-genome approach based on a hierarchical model was developed for a large association study using single nucleotide polymorphisms (SNP).ResultsEstimated sire breeding values (EBV) for persistency and milk yield were calculated using NCSS. Persistency EBV were correlated with peak yield but not with total milk yield. Several SNP were found to be associated with both traits and these were used to identify candidate genes for further investigation.ConclusionNCSS can be used to estimate EBV for lactation persistency and total milk yield, which in turn can be used in whole-genome association studies.Klara L. Verbyla and Arunas P. Verbyl

Quantitative trait loci analysis of zinc efficiency and grain zinc concentration in wheat using whole genome average interval mapping

Zinc (Zn) deficiency is a widespread problem which reduces yield and grain nutritive value in many cereal growing regions of the world. While there is considerable genetic variation in tolerance to Zn deficiency (also known as Zn efficiency), phenotypic selection is difficult and would benefit from the development of molecular markers. A doubled haploid population derived from a cross between the Zn inefficient genotype RAC875-2 and the moderately efficient genotype Cascades was screened in three experiments to identify QTL linked to growth under low Zn and with the concentrations of Zn and iron (Fe) in leaf tissue and in the grain. Two experiments were conducted under controlled conditions while the third examined the response to Zn in the field. QTL were identified using an improved method of analysis, whole genome average interval mapping. Shoot biomass and shoot Zn and Fe concentrations showed significant negative correlations, while there were significant genetic correlations between grain Zn and Fe concentrations. Shoot biomass, tissue and grain Zn concentrations were controlled by a number of genes, many with a minor effect. Depending on the traits and the site, the QTL accounted for 12–81% of the genetic variation. Most of the QTL linked to seedling growth under Zn deficiency and to Zn and Fe concentrations were associated with height genes with greater seedling biomass associated with lower Zn and Fe concentrations. Four QTL for grain Zn concentration and a single QTL for grain Fe concentration were also identified. A cluster of adjacent QTL related to the severity of symptoms of Zn deficiency, shoot Zn concentration and kernel weight was found on chromosome 4A and a cluster of QTL associated with shoot and grain Fe concentrations and kernel weight was found on chromosome 3D. These two regions appear promising areas for further work to develop markers for enhanced growth under low Zn and for Zn and Fe uptake. Although there was no significant difference between the parents, the grain Zn concentration ranged from 29 to 43 mg kg−1 within the population and four QTL associated with grain Zn concentration were identified. These were located on chromosomes 3D, 4B, 6B and 7A and they described 92% of the genetic variation. Each QTL had a relatively small effect on grain Zn concentration but combining the four high Zn alleles increased the grain Zn by 23%. While this illustrates the potential for pyramiding genes to improve grain Zn, breeding for increased grain Zn concentration requires identification of individual QTL with large effects, which in turn requires construction and testing of new mapping populations in the future

Advanced backcross-QTL analysis in spring barley (H. vulgare ssp. spontaneum) comparing a REML versus a Bayesian model in multi-environmental field trials

A common difficulty in mapping quantitative trait loci (QTLs) is that QTL effects may show environment specificity and thus differ across environments. Furthermore, quantitative traits are likely to be influenced by multiple QTLs or genes having different effect sizes. There is currently a need for efficient mapping strategies to account for both multiple QTLs and marker-by-environment interactions. Thus, the objective of our study was to develop a Bayesian multi-locus multi-environmental method of QTL analysis. This strategy is compared to (1) Bayesian multi-locus mapping, where each environment is analysed separately, (2) Restricted Maximum Likelihood (REML) single-locus method using a mixed hierarchical model, and (3) REML forward selection applying a mixed hierarchical model. For this study, we used data on multi-environmental field trials of 301 BC2DH lines derived from a cross between the spring barley elite cultivar Scarlett and the wild donor ISR42-8 from Israel. The lines were genotyped by 98 SSR markers and measured for the agronomic traits “ears per m²,” “days until heading,” “plant height,” “thousand grain weight,” and “grain yield”. Additionally, a simulation study was performed to verify the QTL results obtained in the spring barley population. In general, the results of Bayesian QTL mapping are in accordance with REML methods. In this study, Bayesian multi-locus multi-environmental analysis is a valuable method that is particularly suitable if lines are cultivated in multi-environmental field trials

High Resolution Genomic Scans Reveal Genetic Architecture Controlling Alcohol Preference in Bidirectionally Selected Rat Model

Investigations on the influence of nature vs. nurture on Alcoholism (Alcohol Use Disorder) in human have yet to provide a clear view on potential genomic etiologies. To address this issue, we sequenced a replicated animal model system bidirectionally-selected for alcohol preference (AP). This model is uniquely suited to map genetic effects with high reproducibility, and resolution. The origin of the rat lines (an 8-way cross) resulted in small haplotype blocks (HB) with a corresponding high level of resolution. We sequenced DNAs from 40 samples (10 per line of each replicate) to determine allele frequencies and HB. We achieved ~46X coverage per line and replicate. Excessive differentiation in the genomic architecture between lines, across replicates, termed signatures of selection (SS), were classified according to gene and region. We identified SS in 930 genes associated with AP. The majority (50%) of the SS were confined to single gene regions, the greatest numbers of which were in promoters (284) and intronic regions (169) with the least in exon\u27s (4), suggesting that differences in AP were primarily due to alterations in regulatory regions. We confirmed previously identified genes and found many new genes associated with AP. Of those newly identified genes, several demonstrated neuronal function involved in synaptic memory and reward behavior, e.g. ion channels (Kcnf1, Kcnn3, Scn5a), excitatory receptors (Grin2a, Gria3, Grip1), neurotransmitters (Pomc), and synapses (Snap29). This study not only reveals the polygenic architecture of AP, but also emphasizes the importance of regulatory elements, consistent with other complex traits

Visualization and Curve-Parameter Estimation Strategies for Efficient Exploration of Phenotype Microarray Kinetics

The Phenotype MicroArray (OmniLog® PM) system is able to simultaneously capture a large number of phenotypes by recording an organism's respiration over time on distinct substrates. This technique targets the object of natural selection itself, the phenotype, whereas previously addressed '-omics' techniques merely study components that finally contribute to it. The recording of respiration over time, however, adds a longitudinal dimension to the data. To optimally exploit this information, it must be extracted from the shapes of the recorded curves and displayed in analogy to conventional growth curves.The free software environment R was explored for both visualizing and fitting of PM respiration curves. Approaches using either a model fit (and commonly applied growth models) or a smoothing spline were evaluated. Their reliability in inferring curve parameters and confidence intervals was compared to the native OmniLog® PM analysis software. We consider the post-processing of the estimated parameters, the optimal classification of curve shapes and the detection of significant differences between them, as well as practically relevant questions such as detecting the impact of cultivation times and the minimum required number of experimental repeats.We provide a comprehensive framework for data visualization and parameter estimation according to user choices. A flexible graphical representation strategy for displaying the results is proposed, including 95% confidence intervals for the estimated parameters. The spline approach is less prone to irregular curve shapes than fitting any of the considered models or using the native PM software for calculating both point estimates and confidence intervals. These can serve as a starting point for the automated post-processing of PM data, providing much more information than the strict dichotomization into positive and negative reactions. Our results form the basis for a freely available R package for the analysis of PM data

Use of haplotypes to identify regions harbouring lethal recessive variants in pigs

Additional file 2: Figure S1. Heatmap of the linkage disequilibrium r 2 values between each of the SNPs located within regions 1.1. and 1.2 on SSC1

Joint modelling of location and scale parameters of the t distribution

Joint modelling of location and scale parameters has generally been confined to exponential families. In this paper the location and scale parameters of the t distribution are allowed to depend on covariates. The closed form of the likelihood allows inference to proceed in a similar fashion to the Gaussian location and scale model and provides a framework for a simple scoring algorithm to estimate the parameters. The algorithm includes a procedure to estimate the degrees of freedom parameter of the t distribution. Homogeneity and asymptotic tests are discussed and a methodology is derived to detect heteroscedasticity when the response is t distributed. Simulations reveal considerable bias in the estimates of the degrees of freedom parameter and only minor bias in the estimated fixed effects associated with the scale parameter. In comparison, the estimated location effects are well behaved. To illustrate the joint modelling of location and scale parameters of the t distribution the methodology is applied to two data sets