Anisotropic matern correlation and spatial prediction using REML

The Mat´ern correlation function provides great flexibility for modeling spatially correlated random processes in two dimensions, in particular via a smoothness parameter, whose estimation allows data to determine the degree of smoothness of a spatial process. The extension to include anisotropy provides a very general and flexible class of spatial covariance functions that can be used in a model-based approach to geostatistics, in which parameter estimation is achieved via REML and prediction is within the E-BLUP framework. In this article we develop a general class of linear mixed models using an anisotropic Mat´ern class with an extended metric. The approach is illustrated by application to soil salinity data in a rice-growing field in Australia, and to fine-scale soil pH data. It is found that anisotropy is an important aspect of both datasets, emphasizing the value of a straightforward and accessible approach to modeling anisotropy

Transcriptomic analysis of wheat near-isogenic lines identifies PM19-A1 and A2 as candidates for a major dormancy QTL

BACKGROUND: Next-generation sequencing technologies provide new opportunities to identify the genetic components responsible for trait variation. However, in species with large polyploid genomes, such as bread wheat, the ability to rapidly identify genes underlying quantitative trait loci (QTL) remains non-trivial. To overcome this, we introduce a novel pipeline that analyses, by RNA-sequencing, multiple near-isogenic lines segregating for a targeted QTL. RESULTS: We use this approach to characterize a major and widely utilized seed dormancy QTL located on chromosome 4AL. It exploits the power and mapping resolution afforded by large multi-parent mapping populations, whilst reducing complexity by using multi-allelic contrasts at the targeted QTL region. Our approach identifies two adjacent candidate genes within the QTL region belonging to the ABA-induced Wheat Plasma Membrane 19 family. One of them, PM19-A1, is highly expressed during grain maturation in dormant genotypes. The second, PM19-A2, shows changes in sequence causing several amino acid alterations between dormant and non-dormant genotypes. We confirm that PM19 genes are positive regulators of seed dormancy. CONCLUSIONS: The efficient identification of these strong candidates demonstrates the utility of our transcriptomic pipeline for rapid QTL to gene mapping. By using this approach we are able to provide a comprehensive genetic analysis of the major source of grain dormancy in wheat. Further analysis across a diverse panel of bread and durum wheats indicates that this important dormancy QTL predates hexaploid wheat. The use of these genes by wheat breeders could assist in the elimination of pre-harvest sprouting in wheat.Jose M. Barrero, Colin Cavanagh, Klara L. Verbyla, Josquin F.G. Tibbits, Arunas P. Verbyla, B. Emma Huang, Garry M. Rosewarne, Stuart Stephen, Penghao Wang, Alex Whan, Philippe Rigault, Matthew J. Hayden, and Frank Guble

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

Safely managed hygiene : a risk-based assessment of handwashing water quality

Sustainable Development Goal (SDG) Indicator 6.2.1 requires household handwashing facilities to have soap and water, but there are no guidelines for handwashing water quality. In contrast, drinking water quality guidelines are defined: water must be "free from contamination" to be defined as "safely managed" (SDG Indicator 6.1.1). We modeled the hypothesized mechanism of infection due to contaminated handwashing water to inform risk-based guidelines for microbial quality of handwashing water. We defined two scenarios that should not occur: (1) if handwashing caused fecal contamination, indicated using Escherichia coli, on a person's hands to increase rather than decrease and (2) if hand-to-mouth contacts following handwashing caused an infection risk greater than an acceptable threshold. We found water containing <1000 E. coli colony-forming units (CFU) per 100 mL removes E. coli from hands with>99.9% probability. However, for the annual probability of infection to be <1:1000, handwashing water must contain <2 × 10; -6; focus-forming units of rotavirus, <1 × 10; -4; CFU of Vibrio cholerae, and <9 × 10; -6; Cryptosporidium oocysts per 100 mL. Our model suggests that handwashing with nonpotable water will generally reduce fecal contamination on hands but may be unable to lower the annual probability of infection risks from hand-to-mouth contacts below 1:1000

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

The Embedding Problem for Markov Models of Nucleotide Substitution

10.1371/journal.pone.0069187PLoS ONE87-POLN