Combining Partially Overlapping Multi-Omics Data in Databases Using Relationship Matrices

An Expectation-Maximization Algorithm for Combining a Sample of Partially Overlapping Covariance Matrices

The generation of unprecedented amounts of data brings new challenges in data management, but also an opportunity to accelerate the identification of processes of multiple science disciplines. One of these challenges is the harmonization of high-dimensional unbalanced and heterogeneous data. In this manuscript, we propose a statistical approach to combine incomplete and partially-overlapping pieces of covariance matrices that come from independent experiments. We assume that the data are a random sample of partial covariance matrices sampled from Wishart distributions and we derive an expectation-maximization algorithm for parameter estimation. We demonstrate the properties of our method by (i) using simulation studies and (ii) using empirical datasets. In general, being able to make inferences about the covariance of variables not observed in the same experiment is a valuable tool for data analysis since covariance estimation is an important step in many statistical applications, such as multivariate analysis, principal component analysis, factor analysis, and structural equation modeling.

Alternative scoring methods of fusarium head blight resistance for genomic assisted breeding

Fusarium head blight (FHB) is a fungal disease of wheat (Triticum aestivum.L) that causes yield losses and produces mycotoxins which could easily exceed the limits of the EU regulations. Resistance to FHB has a complex genetic architecture and accurate evaluation in breeding programs is key to selecting resistant varieties. The Area Under the Disease Progress Curve (AUDPC) is one of the commonly metric used as a standard methodology to score FHB. Although efficient, AUDPC requires significant costs in phenotyping to cover the entire disease development pattern. Here, we show that there are more efficient alternatives to AUDPC (angle, growing degree days to reach 50% FHB severity, and FHB maximum variance) that reduce the number of field assessments required and allow for fair comparisons between unbalanced evaluations across trials. Furthermore, we found that the evaluation method that captures the maximum variance in FHB severity across plots is the most optimal approach for scoring FHB. In addition, results obtained on experimental data were validated on a simulated experiment where the disease progress curve was modeled as a sigmoid curve with known parameters and assessment protocols were fully controlled. Results show that alternative metrics tested in this study captured key components of quantitative plant resistance. Moreover, the new metrics could be a starting point for more accurate methods for measuring FHB in the field. For example, the optimal interval for FHB evaluation could be predicted using prior knowledge from historical weather data and FHB scores from previous trials. Finally, the evaluation methods presented in this study can reduce the FHB phenotyping burden in plant breeding with minimal losses on signal detection, resulting in a response variable available to use in data-driven analysis such as genome-wide association studies or genomic selection.

Integrating and optimizing genomic, weather, and secondary trait data for multiclass classification

Modern plant breeding programs collect several data types such as weather, images, and secondary or associated traits besides the main trait (e.g., grain yield). Genomic data is high-dimensional and often over-crowds smaller data types when naively combined to explain the response variable. There is a need to develop methods able to effectively combine different data types of differing sizes to improve predictions. Additionally, in the face of changing climate conditions, there is a need to develop methods able to effectively combine weather information with genotype data to predict the performance of lines better. In this work, we develop a novel three-stage classifier to predict multi-class traits by combining three data types-genomic, weather, and secondary trait. The method addressed various challenges in this problem, such as confounding, differing sizes of data types, and threshold optimization. The method was examined in different settings, including binary and multi-class responses, various penalization schemes, and class balances. Then, our method was compared to standard machine learning methods such as random forests and support vector machines using various classification accuracy metrics and using model size to evaluate the sparsity of the model. The results showed that our method performed similarly to or better than machine learning methods across various settings. More importantly, the classifiers obtained were highly sparse, allowing for a straightforward interpretation of relationships between the response and the selected predictors.

Genomic selection of parents and crosses beyond the native gene pool of a breeding program

Genomic selection has become a valuable tool for selecting cultivar candidates in many plant breeding programs. Genomic selection of elite parents and crossing combinations with germplasm developed outside a breeding program has, however, hardly been explored until now. The aim of this study was to assess the potential of this method for commonly ranking and selecting elite germplasm developed within and beyond a given breeding program. A winter wheat (Triticum aestivum L.) population consisting of 611 in-house and 87 externally developed lines was used to compare training population compositions and statistical models for genomically predicting baking quality in this framework. Augmenting training populations with lines from other breeding programs had a larger influence on the prediction ability than adding in-house generated lines when aiming to commonly rank both germplasm sets. Exploiting preexisting information of secondary correlated traits resulted likewise in more accurate predictions both in empirical analyses and simulations. Genotyping germplasm developed beyond a given breeding program is moreover a convenient way to clarify its relationships with a breeder's own germplasm because pedigree information is oftentimes not available for this purpose. Genomic predictions can thus support a more informed diversity management, especially when integrating simply to phenotype correlated traits to partly circumvent resource reallocations for a costly phenotyping of germplasm from other programs.