Bayesian methods, such as BayesR, for predicting the genetic value or risk of individuals from their genotypes, such as Single Nucleotide Polymorphisms (SNP), are often implemented using a Markov Chain Monte Carlo (MCMC) process. However, the generation of Markov chains is computationally slow. We introduce a form of blocked Gibbs sampling for estimating SNP effects from Markov chains that greatly reduces computational time by sampling each SNP effect iteratively n-times from conditional block posteriors. Subsequent iteration over all blocks m-times produces chains of length m × n. We use this strategy to solve large-scale genomic prediction and fine-mapping problems using the Bayesian MCMC mixed-effects genetic model, BayesR3. We validate the method using simulated data, followed by analysis of empirical dairy cattle data using high dimension milk mid infra-red spectra data as an example of “omics” data and show its use to increase the precision of mapping variants affecting milk, fat, and protein yields relative to a univariate analysis of milk, fat, and protein.

The extraction of automated plant phenomics from digital images has advanced in recent years. However, the accuracy of extracted phenomics, especially for individual plants in a field environment, requires improvement. In this paper, a new and efficient method of extracting individual plant areas and their mean normalized difference vegetation index from high-resolution digital images is proposed. The algorithm was applied on perennial ryegrass row field data multispectral images taken from the top view. First, the center points of individual plants from digital images were located to exclude plant positions without plants. Second, the accurate area of each plant was extracted using its center point and radius. Third, the accurate mean normalized difference vegetation index of each plant was extracted and adjusted for overlapping plants. The correlation between the extracted individual plant phenomics and fresh weight ranged between 0.63 and 0.75 across four time points. The methods proposed are applicable to other crops where individual plant phenotypes are of interest.

Background: Lentil (Lens culinaris Medik.) is a globally-significant agricultural crop used to feed millions of people. Lentils have been cultivated in the Australian states of Victoria and South Australia for several decades, but efforts are now being made to expand their cultivation into Western Australia and New South Wales. Plant architecture plays a pivotal role in adaptation, leading to improved and stable yields especially in new expansion regions. Image-based high-throughput phenomics technologies provide opportunities for an improved understanding of plant development, architecture, and trait genetics. This paper describes a novel method for mapping and quantifying individual branch structures on immature glasshouse-grown lentil plants grown using a LemnaTec Scanalyser 3D high-throughput phenomics platform, which collected side-view RGB images at regular intervals under controlled photographic conditions throughout the experiment. A queue and distance-based algorithm that analysed morphological skeletons generated from images of lentil plants was developed in Python. This code was incorporated into an image analysis pipeline using open-source software (PlantCV) to measure the number, angle, and length of individual branches on lentil plants. Results: Branching structures could be accurately identified and quantified in immature plants, which is sufficient for calculating early vigour traits, however the accuracy declined as the plants matured. Absolute accuracy for branch counts was 77.9% for plants at 22 days after sowing (DAS), 57.9% at 29 DAS and 51.9% at 36 DAS. Allowing for an error of ± 1 branch, the associated accuracies for the same time periods were 97.6%, 90.8% and 79.2% respectively. Occlusion in more mature plants made the mapping of branches less accurate, but the information collected could still be useful for trait estimation. For branch length calculations, the amount of variance explained by linear mixed-effects models was 82% for geodesic length and 87% for Euclidean branch lengths. Within these models, both the mean geodesic and Euclidean distance measurements of branches were found to be significantly affected by genotype, DAS and their interaction. Two informative metrices were derived from the calculations of branch angle; ‘splay’ is a measure of how far a branch angle deviates from being fully upright whilst ‘angle-difference’ is the difference between the smallest and largest recorded branch angle on each plant. The amount of variance explained by linear mixed-effects models was 38% for splay and 50% for angle difference. These lower R² values are likely due to the inherent difficulties in measuring these parameters, nevertheless both splay and angle difference were found to be significantly affected by cultivar, DAS and their interaction. When 276 diverse lentil genotypes with varying degrees of salt tolerance were grown in a glasshouse-based experiment where a portion were subjected to a salt treatment, the branching algorithm was able to distinguish between salt-treated and untreated lentil lines based on differences in branch counts. Likewise, the mean geodesic and Euclidean distance measurements of branches were both found to be significantly affected by cultivar, DAS and salt treatment. The amount of variance explained by the linear mixed-effects models was 57.8% for geodesic branch length and 46.5% for Euclidean branch length. Conclusion: The methodology enabled the accurate quantification of the number, angle, and length of individual branches on glasshouse-grown lentil plants. This methodology could be applied to other dicotyledonous species.