Marco C. A. M. Bink

Using complex plant pedigrees to map valuable genes

Statistical methods pioneered by human and animal geneticists use marker and pedigree information to detect quantitative trait loci within complex pedigrees. These methods, adapted to plants, promise to expand the range of data useful for identifying the genetic factors influencing plant growth, development and evolutionary responses, and to increase the relevance and cost effectiveness of quantitative trait loci mapping in applied contexts.

Detection and use of QTL for complex traits in multiple environments

QTL mapping methods for complex traits are challenged by new developments in marker technology, phenotyping platforms, and breeding methods. In meeting these challenges, QTL mapping approaches will need to also acknowledge the central roles of QTL by environment interactions (QEI) and QTL by trait interactions in the expression of complex traits like yield. This paper presents an overview of mixed model QTL methodology that is suitable for many types of populations and that allows predictive modeling of QEI, both for environmental and developmental gradients. Attention is also given to multi-trait QTL models which are essential to interpret the genetic basis of trait correlations. Biophysical (crop growth) model simulations are proposed as a complement to statistical QTL mapping for the interpretation of the nature of QEI and to investigate better methods for the dissection of complex traits into component traits and their genetic controls.

Genomic selection for fruit quality traits in apple (Malus×domestica Borkh.).

The genome sequence of apple (Malus×domestica Borkh.) was published more than a year ago, which helped develop an 8K SNP chip to assist in implementing genomic selection (GS). In apple breeding programmes, GS can be used to obtain genomic breeding values (GEBV) for choosing next-generation parents or selections for further testing as potential commercial cultivars at a very early stage. Thus GS has the potential to accelerate breeding efficiency significantly because of decreased generation interval or increased selection intensity. We evaluated the accuracy of GS in a population of 1120 seedlings generated from a factorial mating design of four females and two male parents. All seedlings were genotyped using an Illumina Infinium chip comprising 8,000 single nucleotide polymorphisms (SNPs), and were phenotyped for various fruit quality traits. Random-regression best liner unbiased prediction (RR-BLUP) and the Bayesian LASSO method were used to obtain GEBV, and compared using a cross-validation approach for their accuracy to predict unobserved BLUP-BV. Accuracies were very similar for both methods, varying from 0.70 to 0.90 for various fruit quality traits. The selection response per unit time using GS compared with the traditional BLUP-based selection were very high (>100%) especially for low-heritability traits. Genome-wide average estimated linkage disequilibrium (LD) between adjacent SNPs was 0.32, with a relatively slow decay of LD in the long range (r 2 = 0.33 and 0.19 at 100 kb and 1,000 kb respectively), contributing to the higher accuracy of GS. Distribution of estimated SNP effects revealed involvement of large effect genes with likely pleiotropic effects. These results demonstrated that genomic selection is a credible alternative to conventional selection for fruit quality traits.

Towards genomic selection in apple (Malus ×domestica Borkh.) breeding programmes: Prospects, challenges and strategies

The apple genome sequence and the availability of high-throughput genotyping technologies have initiated a new era where SNP markers are abundant across the whole genome. Genomic selection (GS) is a statistical approach that utilizes all available genome-wide markers simultaneously to estimate breeding values or total genetic values. For breeding programmes, GS is a promising alternative to the traditional marker-assisted selection for manipulating complex polygenic traits often controlled by many small-effect genes. Various factors, such as genetic architecture of selection traits, population size and structure, genetic evaluation systems, density of SNP markers and extent of linkage disequilibrium, have been shown to be the key drivers of the accuracy of GS. In this paper, we provide an overview of the status of these aspects in current apple-breeding programmes. Strategies for GS for fruit quality and disease resistance are discussed, and an update on an empirical genomic selection study in a New Zealand apple-breeding programme is provided, along with a foresight of expected accuracy from such selection.

Bayesian analysis of complex traits in pedigreed plant populations

A Bayesian approach to analyze complex traits is presented that can help plant eneticists and breeders in exploiting the marker and phenotypic data on pedigreed populations as available from ongoing breeding programs. The statistical model for the quantitative trait may include non-genetic and genetic components. The latter component can be divided into QTL on known marker linkage groups, major genes and a polygenic component. The full probability model, prior assumptions on model variables are presented and criterion for model selection and posterior inferences are given. Simulated data on a known pedigreed population structure of the EU project HiDRAS was used to illustrate the use of the Bayesian approach to analyze complex traits. It was shown that estimates for QTL parameters were more accurate when non-genetic factors were included in the model and when a polygenic component was included when not all linkage groups were analyzed simultaneously. The Bayesian approach has been implemented into the software package FlexQTL and allows plant breeders explore their pedigreed populations for segregating QTL alleles that are relevant in their breeding program.

Multiple QTL mapping in related plant populations via a pedigree-analysis approach

Abstract.QTL mapping experiments in plant breeding may involve multiple populations or pedigrees that are related through their ancestors. These known relationships have often been ignored for the sake of statistical analysis, despite their potential increase in power of mapping. We describe here a Bayesian method for QTL mapping in complex plant populations and reported the results from its application to a (previously analysed) potato data set. This Bayesian method was originally developed for human genetics data, and we have proved that it is useful for complex plant populations as well, based on a sensitivity analysis that was performed here. The method accommodates robustness to complex structures in pedigree data, full flexibility in the estimation of the number of QTL across multiple chromosomes, thereby accounting for uncertainties in the transmission of QTL and marker alleles due to incomplete marker information, and the simultaneous inclusion of non-genetic factors affecting the quantitative trait.

Bayesian Markov Random Field Analysis for Protein Function Prediction Based on Network Data

Inference of protein functions is one of the most important aims of modern biology. To fully exploit the large volumes of genomic data typically produced in modern-day genomic experiments, automated computational methods for protein function prediction are urgently needed. Established methods use sequence or structure similarity to infer functions but those types of data do not suffice to determine the biological context in which proteins act. Current high-throughput biological experiments produce large amounts of data on the interactions between proteins. Such data can be used to infer interaction networks and to predict the biological process that the protein is involved in. Here, we develop a probabilistic approach for protein function prediction using network data, such as protein-protein interaction measurements. We take a Bayesian approach to an existing Markov Random Field method by performing simultaneous estimation of the model parameters and prediction of protein functions. We use an adaptive Markov Chain Monte Carlo algorithm that leads to more accurate parameter estimates and consequently to improved prediction performance compared to the standard Markov Random Fields method. We tested our method using a high quality S.cereviciae validation network with 1622 proteins against 90 Gene Ontology terms of different levels of abstraction. Compared to three other protein function prediction methods, our approach shows very good prediction performance. Our method can be directly applied to protein-protein interaction or coexpression networks, but also can be extended to use multiple data sources. We apply our method to physical protein interaction data from S. cerevisiae and provide novel predictions, using 340 Gene Ontology terms, for 1170 unannotated proteins and we evaluate the predictions using the available literature.

Novel genomic approaches unravel genetic architecture of complex traits in apple

BackgroundUnderstanding the genetic architecture of quantitative traits is important for developing genome-based crop improvement methods. Genome-wide association study (GWAS) is a powerful technique for mining novel functional variants. Using a family-based design involving 1,200 apple (Malus × domestica Borkh.) seedlings genotyped for an 8K SNP array, we report the first systematic evaluation of the relative contributions of different genomic regions to various traits related to eating quality and susceptibility to some physiological disorders. Single-SNP analyses models that accounted for population structure, or not, were compared with models fitting all markers simultaneously. The patterns of linkage disequilibrium (LD) were also investigated.ResultsA high degree of LD even at longer distances between markers was observed, and the patterns of LD decay were similar across successive generations. Genomic regions were identified, some of which coincided with known candidate genes, with significant effects on various traits. Phenotypic variation explained by the loci identified through a whole-genome scan ranged from 3% to 25% across different traits, while fitting all markers simultaneously generally provided heritability estimates close to those from pedigree-based analysis. Results from ‘Q+K’ and ‘K’ models were very similar, suggesting that the SNP-based kinship matrix captures most of the underlying population structure. Correlations between allele substitution effects obtained from single-marker and all-marker analyses were about 0.90 for all traits. Use of SNP-derived realized relationships in linear mixed models provided a better goodness-of-fit than pedigree-based expected relationships. Genomic regions with probable pleiotropic effects were supported by the corresponding higher linkage group (LG) level estimated genetic correlations.ConclusionsThe accuracy of artificial selection in plants species can be increased by using more precise marker-derived estimates of realized coefficients of relationships. All-marker analyses that indirectly account for population- and pedigree structure will be a credible alternative to single-SNP analyses in GWAS. This study revealed large differences in the genetic architecture of apple fruit traits, and the marker-trait associations identified here will help develop genome-based breeding methods for apple cultivar development.

Accuracy of imputation to whole-genome sequence data in Holstein Friesian cattle

BackgroundThe use of whole-genome sequence data can lead to higher accuracy in genome-wide association studies and genomic predictions. However, to benefit from whole-genome sequence data, a large dataset of sequenced individuals is needed. Imputation from SNP panels, such as the Illumina BovineSNP50 BeadChip and Illumina BovineHD BeadChip, to whole-genome sequence data is an attractive and less expensive approach to obtain whole-genome sequence genotypes for a large number of individuals than sequencing all individuals. Our objective was to investigate accuracy of imputation from lower density SNP panels to whole-genome sequence data in a typical dataset for cattle.MethodsWhole-genome sequence data of chromosome 1 (1737 471 SNPs) for 114 Holstein Friesian bulls were used. Beagle software was used for imputation from the BovineSNP50 (3132 SNPs) and BovineHD (40 492 SNPs) beadchips. Accuracy was calculated as the correlation between observed and imputed genotypes and assessed by five-fold cross-validation. Three scenarios S40, S60 and S80 with respectively 40%, 60%, and 80% of the individuals as reference individuals were investigated.ResultsMean accuracies of imputation per SNP from the BovineHD panel to sequence data and from the BovineSNP50 panel to sequence data for scenarios S40 and S80 ranged from 0.77 to 0.83 and from 0.37 to 0.46, respectively. Stepwise imputation from the BovineSNP50 to BovineHD panel and then to sequence data for scenario S40 improved accuracy per SNP to 0.65 but it varied considerably between SNPs.ConclusionsAccuracy of imputation to whole-genome sequence data was generally high for imputation from the BovineHD beadchip, but was low from the BovineSNP50 beadchip. Stepwise imputation from the BovineSNP50 to the BovineHD beadchip and then to sequence data substantially improved accuracy of imputation. SNPs with a low minor allele frequency were more difficult to impute correctly and the reliability of imputation varied more. Linkage disequilibrium between an imputed SNP and the SNP on the lower density panel, minor allele frequency of the imputed SNP and size of the reference group affected imputation reliability.

SPICY: towards automated phenotyping of large pepper plants in the greenhouse

Most high-throughput systems for automated plant phenotyping involve a fixed recording cabinet to which plants are transported. However, important greenhouse plants like pepper are too tall to be transported. In this research we developed a system to automatically measure plant characteristics of tall pepper plants in the greenhouse. With a device equipped with multiple cameras, images of plants are recorded at a 5cm interval over a height of 3m. Two types of features are extracted: (1) features from a 3D reconstruction of the plant canopy; and (2) statistical features derived directly from RGB images. The experiment comprised 151 genotypes of a recombinant inbred population of pepper, to examine the heritability and quantitative trait loci (QTL) of the features. Features extracted from the 3D reconstruction of the canopy were leaf size and leaf angle, with heritabilities of 0.70 and 0.56 respectively. Three QTL were found for leaf size, and one for leaf angle. From the statistical features, plant height showed a good correlation (0.93) with manual measurements, and QTL were in accordance with QTL of manual measurements. For total leaf area, the heritability was 0.55, and two of the three QTL found by manual measurement were found by image analysis.