Identifying the relation between food groups and biological ageing: a data-driven approach.

BACKGROUND: Heterogeneity in ageing rates drives the need for research into lifestyle secrets of successful agers. Biological age, predicted by epigenetic clocks, has been shown to be a more reliable measure of ageing than chronological age. Dietary habits are known to affect the ageing process. However, much remains to be learnt about specific dietary habits that may directly affect the biological process of ageing. OBJECTIVE: To identify food groups that are directly related to biological ageing, using Copula Graphical Models. METHODS: We performed a preregistered analysis of 3,990 postmenopausal women from the Women's Health Initiative, based in North America. Biological age acceleration was calculated by the epigenetic clock PhenoAge using whole-blood DNA methylation. Copula Graphical Modelling, a powerful data-driven exploratory tool, was used to examine relations between food groups and biological ageing whilst adjusting for an extensive amount of confounders. Two food group-age acceleration networks were established: one based on the MyPyramid food grouping system and another based on item-level food group data. RESULTS: Intake of eggs, organ meat, sausages, cheese, legumes, starchy vegetables, added sugar and lunch meat was associated with biological age acceleration, whereas intake of peaches/nectarines/plums, poultry, nuts, discretionary oil and solid fat was associated with decelerated ageing. CONCLUSION: We identified several associations between specific food groups and biological ageing. These findings pave the way for subsequent studies to ascertain causality and magnitude of these relationships, thereby improving the understanding of biological mechanisms underlying the interplay between food groups and biological ageing.

Natural variation in salt-induced root growth phases and their contribution to root architecture plasticity.

The root system architecture of a plant changes during salt stress exposure. Different accessions of Arabidopsis thaliana have adopted different strategies in remodelling their root architecture during salt stress. Salt induces a multiphase growth response in roots, consisting of a stop phase, quiescent phase, recovery phase and eventually a new level of homoeostasis. We explored natural variation in the length of and growth rate during these phases in both main and lateral roots and find that some accessions lack the quiescent phase. Using mathematical models and correlation-based network, allowed us to correlate dynamic traits to overall root architecture and discover that both the main root growth rate during homoeostasis and lateral root appearance are the strongest determinants of overall root architecture. In addition, this approach revealed a trade-off between investing in main or lateral root length during salt stress. By studying natural variation in high-resolution temporal root growth using mathematical modelling, we gained new insights in the interactions between dynamic root growth traits and we identified key traits that modulate overall root architecture during salt stress.

De novo construction of polyploid linkage maps using discrete graphical models.

MOTIVATION: Linkage maps are used to identify the location of genes responsible for traits and diseases. New sequencing techniques have created opportunities to substantially increase the density of genetic markers. Such revolutionary advances in technology have given rise to new challenges, such as creating high-density linkage maps. Current multiple testing approaches based on pairwise recombination fractions are underpowered in the high-dimensional setting and do not extend easily to polyploid species. To remedy these issues, we propose to construct linkage maps using graphical models either via a sparse Gaussian copula or a non-paranormal skeptic approach. RESULTS: We determine linkage groups, typically chromosomes, and the order of markers in each linkage group by inferring the conditional independence relationships among large numbers of markers in the genome. Through simulations, we illustrate the utility of our map construction method and compare its performance with other available methods, both when the data are clean and contain no missing observations and when data contain genotyping errors. Our comprehensive map construction method makes full use of the dosage SNP data to reconstruct linkage map for any bi-parental diploid and polyploid species. We apply the proposed method to three genotype datasets: barley, peanut and potato from diploid and polyploid populations. AVAILABILITY AND IMPLEMENTATION: The method is implemented in the R package netgwas which is freely available at https://cran.r-project.org/web/packages/netgwas. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dietary Intakes of Vegetable Protein, Folate, and Vitamins B-6 and B-12 Are Partially Correlated with Physical Functioning of Dutch Older Adults Using Copula Graphical Models.

BACKGROUND: In nutritional epidemiology, dealing with confounding and complex internutrient relations are major challenges. An often-used approach is dietary pattern analyses, such as principal component analysis, to deal with internutrient correlations, and to more closely resemble the true way nutrients are consumed. However, despite these improvements, these approaches still require subjective decisions in the preselection of food groups. Moreover, they do not make efficient use of multivariate dietary data, because they detect only marginal associations. We propose the use of copula graphical models (CGMs) to model and make statistical inferences regarding complex associations among variables in multivariate data, where associations between all variables can be learned simultaneously. OBJECTIVE: We aimed to reconstruct nutritional intake and physical functioning networks in Dutch older adults by applying a CGM. METHODS: We addressed this issue by uncovering the pairwise associations between variables while correcting for the effect of remaining variables. More specifically, we used a CGM to infer the precision matrix, which contains all the conditional independence relations between nodes in the graph. The nonzero elements of the precision matrix indicate the presence of a direct association. We applied this method to reconstruct nutrient-physical functioning networks from the combined data of 4 studies (Nu-Age, ProMuscle, ProMO, and V-Fit, total n = 662, mean ± SD age = 75 ± 7 y). The method was implemented in the R package nutriNetwork which is freely available at https://cran.r-project.org/web/packages/nutriNetwork. RESULTS: Greater intakes of vegetable protein and vitamin B-6 were partially correlated with higher scores on the total Short Physical Performance Battery (SPPB) and the chair rise test. Greater intakes of vitamin B-12 and folate were partially correlated with higher scores on the chair rise test and the total SPPB, respectively. CONCLUSIONS: We determined that vegetable protein, vitamin B-6, folate, and vitamin B-12 intakes are partially correlated with improved functional outcome measurements in Dutch older adults.

Discovering the direct relations between nutrients and epigenetic ageing.

BACKGROUND: Along with the ageing of society, the absolute prevalence of age-related diseases is expected to rise, leading to a substantial burden on healthcare systems and society. Thus, there is an urgent need to promote healthy ageing. As opposed to chronological age, biological age was introduced to accurately represent the ageing process, as it considers physiological deterioration that is linked to morbidity and mortality risk. Furthermore, biological age responds to various factors, including nutritional factors, which have the potential to mitigate the risk of age-related diseases. As a result, a promising biomarker of biological age known as the epigenetic clock has emerged as a suitable measure to investigate the direct relations between nutritional factors and ageing, thereby identifying potential intervention targets to improve healthy ageing. METHODS: In this study, we analysed data from 3,969 postmenopausal women from the Women's Health Initiative to identify nutrients that are associated with the rate of ageing by using an accurate measure of biological age called the PhenoAge epigenetic clock. We used Copula Graphical Models, a data-driven exploratory analysis tool, to identify direct relationships between nutrient intake and age-acceleration, while correcting for every variable in the dataset. RESULTS: We revealed that increased dietary intakes of coumestrol, beta-carotene and arachidic acid were associated with decelerated epigenetic ageing. In contrast, increased intakes of added sugar, gondoic acid, behenic acid, arachidonic acid, vitamin A and ash were associated with accelerated epigenetic ageing in postmenopausal women. CONCLUSION: Our study discovered direct relations between nutrients and epigenetic ageing, revealing promising areas for follow-up studies to determine the magnitude and causality of our estimated diet-epigenetic relationships.

Complementary approaches to dissect late leaf rust resistance in an interspecific raspberry population.

Over the last 10 years, global raspberry production has increased by 47.89%, based mainly on the red raspberry species (Rubus idaeus). However, the black raspberry (Rubus occidentalis), although less consumed, is resistant to one of the most important diseases for the crop, the late leaf rust caused by Acculeastrum americanum fungus. In this context, genetic resistance is the most sustainable way to control the disease, mainly because there are no registered fungicides for late leaf rust in Brazil. Therefore, the aim was to understand the genetic architecture that controls resistance to late leaf rust in raspberries. For that, we used an interspecific multiparental population using the species mentioned above as parents, 2 different statistical approaches to associate the phenotypes with markers [GWAS (genome-wide association studies) and copula graphical models], and 2 phenotyping methodologies from the first to the 17th day after inoculation (high-throughput phenotyping with a multispectral camera and traditional phenotyping by disease severity scores). Our findings indicate that a locus of higher effect, at position 13.3 Mb on chromosome 5, possibly controls late leaf rust resistance, as both GWAS and the network suggested the same marker. Of the 12 genes flanking its region, 4 were possible receptors, 3 were likely defense executors, 1 gene was likely part of signaling cascades, and 4 were classified as nondefense related. Although the network and GWAS indicated the same higher effect genomic region, the network identified other different candidate regions, potentially complementing the genetic control comprehension.

Reconstruction of Networks with Direct and Indirect Genetic Effects.

Genetic variance of a phenotypic trait can originate from direct genetic effects, or from indirect effects, i.e., through genetic effects on other traits, affecting the trait of interest. This distinction is often of great importance, for example, when trying to improve crop yield and simultaneously control plant height. As suggested by Sewall Wright, assessing contributions of direct and indirect effects requires knowledge of (1) the presence or absence of direct genetic effects on each trait, and (2) the functional relationships between the traits. Because experimental validation of such relationships is often unfeasible, it is increasingly common to reconstruct them using causal inference methods. However, most current methods require all genetic variance to be explained by a small number of quantitative trait loci (QTL) with fixed effects. Only a few authors have considered the "missing heritability" case, where contributions of many undetectable QTL are modeled with random effects. Usually, these are treated as nuisance terms that need to be eliminated by taking residuals from a multi-trait mixed model (MTM). But fitting such an MTM is challenging, and it is impossible to infer the presence of direct genetic effects. Here, we propose an alternative strategy, where genetic effects are formally included in the graph. This has important advantages: (1) genetic effects can be directly incorporated in causal inference, implemented via our PCgen algorithm, which can analyze many more traits; and (2) we can test the existence of direct genetic effects, and improve the orientation of edges between traits. Finally, we show that reconstruction is much more accurate if individual plant or plot data are used, instead of genotypic means. We have implemented the PCgen-algorithm in the R-package pcgen.