Flexible modelling of spatial variation in agricultural field trials with the R package INLA

The objective of this paper was to fit different established spatial models for analysing agricultural field trials using the open-source R package INLA. Spatial variation is common in field trials, and accounting for it increases the accuracy of estimated genetic effects. However, this is still hindered by the lack of available software implementations. We compare some established spatial models and show possibilities for flexible modelling with respect to field trial design and joint modelling over multiple years and locations. We use a Bayesian framework and for statistical inference the integrated nested Laplace approximations (INLA) implemented in the R package INLA. The spatial models we use are the well-known independent row and column effects, separable first-order autoregressive ( AR1⊗AR1 ) models and a Gaussian random field (Matérn) model that is approximated via the stochastic partial differential equation approach. The Matérn model can accommodate flexible field trial designs and yields interpretable parameters. We test the models in a simulation study imitating a wheat breeding programme with different levels of spatial variation, with and without genome-wide markers and with combining data over two locations, modelling spatial and genetic effects jointly. The results show comparable predictive performance for both the AR1⊗AR1 and the Matérn models. We also present an example of fitting the models to a real wheat breeding data and simulated tree breeding data with the Nelder wheel design to show the flexibility of the Matérn model and the R package INLA

Photoacclimation in Dunaliella tertiolecta reveals a unique NPQ pattern upon exposure to irradiance

Highly time-resolved photoacclimation patterns of the chlorophyte microalga Dunaliella tertiolecta during exposure to an off–on–off (block) light pattern of saturating photon flux, and to a regime of consecutive increasing light intensities are presented. Non-photochemical quenching (NPQ) mechanisms unexpectedly responded with an initial decrease during dark–light transitions. NPQ values started to rise after light exposure of approximately 4 min. State-transitions, measured as a change of PSII:PSI fluorescence emission at 77 K, did not contribute to early NPQ oscillations. Addition of the uncoupler CCCP, however, caused a rapid increase in fluorescence and showed the significance of qE for NPQ. Partitioning of the quantum efficiencies showed that constitutive NPQ was (a) higher than qE-driven NPQ and (b) responded to light treatment within seconds, suggesting an active role of constitutive NPQ in variable energy dissipation, although it is thought to contribute statically to NPQ. The PSII connectivity parameter p correlated well with F′, Fm′ and NPQ during the early phase of the dark–light transients in sub-saturating light, suggesting a plastic energy distribution pattern within energetically connected PSII centres. In consecutive increasing photon flux experiments, correlations were weaker during the second light increment. Changes in connectivity can present an early photoresponse that are reflected in fluorescence signals and NPQ and might be responsive to the short-term acclimation state, and/or to the actinic photon flux

Intra- and Inter-clade Cross-reactivity by HIV-1 Gag Specific T-Cells Reveals Exclusive and Commonly Targeted Regions: Implications for Current Vaccine Trials

The genetic diversity of HIV-1 across the globe is a major challenge for developing an HIV vaccine. To facilitate immunogen design, it is important to characterize clusters of commonly targeted T-cell epitopes across different HIV clades. To address this, we examined 39 HIV-1 clade C infected individuals for IFN-γ Gag-specific T-cell responses using five sets of overlapping peptides, two sets matching clade C vaccine candidates derived from strains from South Africa and China, and three peptide sets corresponding to consensus clades A, B, and D sequences. The magnitude and breadth of T-cell responses against the two clade C peptide sets did not differ, however clade C peptides were preferentially recognized compared to the other peptide sets. A total of 84 peptides were recognized, of which 19 were exclusively from clade C, 8 exclusively from clade B, one peptide each from A and D and 17 were commonly recognized by clade A, B, C and D. The entropy of the exclusively recognized peptides was significantly higher than that of commonly recognized peptides (p = 0.0128) and the median peptide processing scores were significantly higher for the peptide variants recognized versus those not recognized (p = 0.0001). Consistent with these results, the predicted Major Histocompatibility Complex Class I IC50 values were significantly lower for the recognized peptide variants compared to those not recognized in the ELISPOT assay (p<0.0001), suggesting that peptide variation between clades, resulting in lack of cross-clade recognition, has been shaped by host immune selection pressure. Overall, our study shows that clade C infected individuals recognize clade C peptides with greater frequency and higher magnitude than other clades, and that a selection of highly conserved epitope regions within Gag are commonly recognized and give rise to cross-clade reactivities

Global warming and recurrent mass bleaching of corals

During 2015–2016, record temperatures triggered a pan-tropical episode of coral bleaching, the third global-scale event since mass bleaching was first documented in the 1980s. Here we examine how and why the severity of recurrent major bleaching events has varied at multiple scales, using aerial and underwater surveys of Australian reefs combined with satellite-derived sea surface temperatures. The distinctive geographic footprints of recurrent bleaching on the Great Barrier Reef in 1998, 2002 and 2016 were determined by the spatial pattern of sea temperatures in each year. Water quality and fishing pressure had minimal effect on the unprecedented bleaching in 2016, suggesting that local protection of reefs affords little or no resistance to extreme heat. Similarly, past exposure to bleaching in 1998 and 2002 did not lessen the severity of bleaching in 2016. Consequently, immediate global action to curb future warming is essential to secure a future for coral reefs

Evaluating alternate models to estimate genetic parameters of calving traits in United Kingdom Holstein-Friesian dairy cattle

<p>Abstract</p> <p>Background</p> <p>The focus in dairy cattle breeding is gradually shifting from production to functional traits and genetic parameters of calving traits are estimated more frequently. However, across countries, various statistical models are used to estimate these parameters. This study evaluates different models for calving ease and stillbirth in United Kingdom Holstein-Friesian cattle.</p> <p>Methods</p> <p>Data from first and later parity records were used. Genetic parameters for calving ease, stillbirth and gestation length were estimated using the restricted maximum likelihood method, considering different models i.e. sire (−maternal grandsire), animal, univariate and bivariate models. Gestation length was fitted as a correlated indicator trait and, for all three traits, genetic correlations between first and later parities were estimated. Potential bias in estimates was avoided by acknowledging a possible environmental direct-maternal covariance. The total heritable variance was estimated for each trait to discuss its theoretical importance and practical value. Prediction error variances and accuracies were calculated to compare the models.</p> <p>Results and discussion</p> <p>On average, direct and maternal heritabilities for calving traits were low, except for direct gestation length. Calving ease in first parity had a significant and negative direct-maternal genetic correlation. Gestation length was maternally correlated to stillbirth in first parity and directly correlated to calving ease in later parities. Multi-trait models had a slightly greater predictive ability than univariate models, especially for the lowly heritable traits. The computation time needed for sire (−maternal grandsire) models was much smaller than for animal models with only small differences in accuracy. The sire (−maternal grandsire) model was robust when additional genetic components were estimated, while the equivalent animal model had difficulties reaching convergence.</p> <p>Conclusions</p> <p>For the evaluation of calving traits, multi-trait models show a slight advantage over univariate models. Extended sire models (−maternal grandsire) are more practical and robust than animal models. Estimated genetic parameters for calving traits of UK Holstein cattle are consistent with literature. Calculating an aggregate estimated breeding value including direct and maternal values should encourage breeders to consider both direct and maternal effects in selection decisions.</p

Mitogen- and Stress-Activated Kinase 1 (MSK1) Regulates Cigarette Smoke-Induced Histone Modifications on NF-κB-dependent Genes

Cigarette smoke (CS) causes sustained lung inflammation, which is an important event in the pathogenesis of chronic obstructive pulmonary disease (COPD). We have previously reported that IKKα (I kappaB kinase alpha) plays a key role in CS-induced pro-inflammatory gene transcription by chromatin modifications; however, the underlying role of downstream signaling kinase is not known. Mitogen- and stress-activated kinase 1 (MSK1) serves as a specific downstream NF-κB RelA/p65 kinase, mediating transcriptional activation of NF-κB-dependent pro-inflammatory genes. The role of MSK1 in nuclear signaling and chromatin modifications is not known, particularly in response to environmental stimuli. We hypothesized that MSK1 regulates chromatin modifications of pro-inflammatory gene promoters in response to CS. Here, we report that CS extract activates MSK1 in human lung epithelial (H292 and BEAS-2B) cell lines, human primary small airway epithelial cells (SAEC), and in mouse lung, resulting in phosphorylation of nuclear MSK1 (Thr581), phospho-acetylation of RelA/p65 at Ser276 and Lys310 respectively. This event was associated with phospho-acetylation of histone H3 (Ser10/Lys9) and acetylation of histone H4 (Lys12). MSK1 N- and C-terminal kinase-dead mutants, MSK1 siRNA-mediated knock-down in transiently transfected H292 cells, and MSK1 stable knock-down mouse embryonic fibroblasts significantly reduced CS extract-induced MSK1, NF-κB RelA/p65 activation, and posttranslational modifications of histones. CS extract/CS promotes the direct interaction of MSK1 with RelA/p65 and p300 in epithelial cells and in mouse lung. Furthermore, CS-mediated recruitment of MSK1 and its substrates to the promoters of NF-κB-dependent pro-inflammatory genes leads to transcriptional activation, as determined by chromatin immunoprecipitation. Thus, MSK1 is an important downstream kinase involved in CS-induced NF-κB activation and chromatin modifications, which have implications in pathogenesis of COPD

High methylmercury in Arctic and subarctic ponds is related to nutrient levels in the warming eastern Canadian Arctic

Permafrost thaw ponds are ubiquitous in the eastern Canadian Arctic, yet little information exists on their potential as sources of methylmercury (MeHg) to freshwaters. They are microbially active and conducive to methylation of inorganic mercury, and are also affected by Arctic warming. This multiyear study investigated thaw ponds in a discontinuous permafrost region in the Subarctic taiga (Kuujjuarapik-Whapmagoostui, QC) and a continuous permafrost region in the Arctic tundra (Bylot Island, NU). MeHg concentrations in thaw ponds were well above levels measured in most freshwater ecosystems in the Canadian Arctic (>0.1 ng L−1). On Bylot, ice-wedge trough ponds showed significantly higher MeHg (0.3−2.2 ng L−1) than polygonal ponds (0.1−0.3 ng L−1) or lakes (<0.1 ng L−1). High MeHg was measured in the bottom waters of Subarctic thaw ponds near Kuujjuarapik (0.1−3.1 ng L−1). High water MeHg concentrations in thaw ponds were strongly correlated with variables associated with high inputs of organic matter (DOC, a320, Fe), nutrients (TP, TN), and microbial activity (dissolved CO2 and CH4). Thawing permafrost due to Arctic warming will continue to release nutrients and organic carbon into these systems and increase ponding in some regions, likely stimulating higher water concentrations of MeHg. Greater hydrological connectivity from permafrost thawing may potentially increase transport of MeHg from thaw ponds to neighboring aquatic ecosystems