Multi-omics analysis to decipher the molecular link between chronic exposure to pollution and human skin dysfunction.

Environmental pollution is composed of several factors, namely particulate matter (PM2.5, PM10), ozone and Ultra Violet (UV) rays among others and first and the most exposed tissue to these substances is the skin epidermis. It has been established that several skin disorders such as eczema, acne, lentigines and wrinkles are aggravated by exposure to atmospheric pollution. While pollutants can interact with skin surface, contamination of deep skin by ultrafine particles or Polycyclic aromatic hydrocarbons (PAH) might be explained by their presence in blood and hair cortex. Molecular mechanisms leading to skin dysfunction due to pollution exposure have been poorly explored in humans. In addition to various host skin components, cutaneous microbiome is another target of these environment aggressors and can actively contribute to visible clinical manifestation such as wrinkles and aging. The present study aimed to investigate the association between pollution exposure, skin microbiota, metabolites and skin clinical signs in women from two cities with different pollution levels. Untargeted metabolomics and targeted proteins were analyzed from D-Squame samples from healthy women (n = 67 per city), aged 25-45 years and living for at least 15 years in the Chinese cities of Baoding (used as a model of polluted area) and Dalian (control area with lower level of pollution). Additional samples by swabs were collected from the cheeks from the same population and microbiome was analysed using bacterial 16S rRNA as well as fungal ITS1 amplicon sequencing and metagenomics analysis. The level of exposure to pollution was assessed individually by the analysis of polycyclic aromatic hydrocarbons (PAH) and their metabolites in hair samples collected from each participant. All the participants of the study were assessed for the skin clinical parameters (acne, wrinkles, pigmented spots etc.). Women from the two cities (polluted and less polluted) showed distinct metabolic profiles and alterations in skin microbiome. Profiling data from 350 identified metabolites, 143 microbes and 39 PAH served to characterize biochemical events that correlate with pollution exposure. Finally, using multiblock data analysis methods, we obtained a potential molecular map consisting of multi-omics signatures that correlated with the presence of skin pigmentation dysfunction in individuals living in a polluted environment. Overall, these signatures point towards macromolecular alterations by pollution that could manifest as clinical sign of early skin pigmentation and/or other imperfections.

Fast computation of genome-metagenome interaction effects.

MOTIVATION: Association studies have been widely used to search for associations between common genetic variants observations and a given phenotype. However, it is now generally accepted that genes and environment must be examined jointly when estimating phenotypic variance. In this work we consider two types of biological markers: genotypic markers, which characterize an observation in terms of inherited genetic information, and metagenomic marker which are related to the environment. Both types of markers are available in their millions and can be used to characterize any observation uniquely. OBJECTIVE: Our focus is on detecting interactions between groups of genetic and metagenomic markers in order to gain a better understanding of the complex relationship between environment and genome in the expression of a given phenotype. CONTRIBUTIONS: We propose a novel approach for efficiently detecting interactions between complementary datasets in a high-dimensional setting with a reduced computational cost. The method, named SICOMORE, reduces the dimension of the search space by selecting a subset of supervariables in the two complementary datasets. These supervariables are given by a weighted group structure defined on sets of variables at different scales. A Lasso selection is then applied on each type of supervariable to obtain a subset of potential interactions that will be explored via linear model testing. RESULTS: We compare SICOMORE with other approaches in simulations, with varying sample sizes, noise, and numbers of true interactions. SICOMORE exhibits convincing results in terms of recall, as well as competitive performances with respect to running time. The method is also used to detect interaction between genomic markers in Medicago truncatula and metagenomic markers in its rhizosphere bacterial community. SOFTWARE AVAILABILITY: An R package is available [4], along with its documentation and associated scripts, allowing the reader to reproduce the results presented in the paper.

Changes of the human skin microbiota upon chronic exposure to polycyclic aromatic hydrocarbon pollutants.

BACKGROUND: Polycyclic aromatic hydrocarbons (PAHs) are of environmental and public health concerns and contribute to adverse skin attributes such as premature skin aging and pigmentary disorder. However, little information is available on the potential roles of chronic urban PAH pollutant exposure on the cutaneous microbiota. Given the roles of the skin microbiota have on healthy and undesirable skin phenotypes and the relationships between PAHs and skin properties, we hypothesize that exposure of PAHs may be associated with changes in the cutaneous microbiota. In this study, the skin microbiota of over two hundred Chinese individuals from two cities in China with varying exposure levels of PAHs were characterized by bacterial and fungal amplicon and shotgun metagenomics sequencing. RESULTS: Skin site and city were strong parameters in changing microbial communities and their assembly processes. Reductions of bacterial-fungal microbial network structural integrity and stability were associated with skin conditions (acne and dandruff). Multivariate analysis revealed associations between abundances of Propionibacterium and Malassezia with host properties and pollutant exposure levels. Shannon diversity increase was correlated to exposure levels of PAHs in a dose-dependent manner. Shotgun metagenomics analysis of samples (n = 32) from individuals of the lowest and highest exposure levels of PAHs further highlighted associations between the PAHs quantified and decrease in abundances of skin commensals and increase in oral bacteria. Functional analysis identified associations between levels of PAHs and abundance of microbial genes of metabolic and other pathways with potential importance in host-microbe interactions as well as degradation of aromatic compounds. CONCLUSIONS: The results in this study demonstrated the changes in composition and functional capacities of the cutaneous microbiota associated with chronic exposure levels of PAHs. Findings from this study will aid the development of strategies to harness the microbiota in protecting the skin against pollutants. Video Abstract.

Learning the optimal scale for GWAS through hierarchical SNP aggregation.

BACKGROUND: Genome-Wide Association Studies (GWAS) seek to identify causal genomic variants associated with rare human diseases. The classical statistical approach for detecting these variants is based on univariate hypothesis testing, with healthy individuals being tested against affected individuals at each locus. Given that an individual's genotype is characterized by up to one million SNPs, this approach lacks precision, since it may yield a large number of false positives that can lead to erroneous conclusions about genetic associations with the disease. One way to improve the detection of true genetic associations is to reduce the number of hypotheses to be tested by grouping SNPs. RESULTS: We propose a dimension-reduction approach which can be applied in the context of GWAS by making use of the haplotype structure of the human genome. We compare our method with standard univariate and group-based approaches on both synthetic and real GWAS data. CONCLUSION: We show that reducing the dimension of the predictor matrix by aggregating SNPs gives a greater precision in the detection of associations between the phenotype and genomic regions.

Genetic basis of nitrogen use efficiency and yield stability across environments in winter rapeseed.

BACKGROUND: Nitrogen use efficiency is an important breeding trait that can be modified to improve the sustainability of many crop species used in agriculture. Rapeseed is a major oil crop with low nitrogen use efficiency, making its production highly dependent on nitrogen input. This complex trait is suspected to be sensitive to genotype × environment interactions, especially genotype × nitrogen interactions. Therefore, phenotyping diverse rapeseed populations under a dense network of trials is a powerful approach to study nitrogen use efficiency in this crop. The present study aimed to determine the quantitative trait loci (QTL) associated with yield in winter oilseed rape and to assess the stability of these regions under contrasting nitrogen conditions for the purpose of increasing nitrogen use efficiency. RESULTS: Genome-wide association studies and linkage analyses were performed on two diversity sets and two doubled-haploid populations. These populations were densely genotyped, and yield-related traits were scored in a multi-environment design including seven French locations, six growing seasons (2009 to 2014) and two nitrogen nutrition levels (optimal versus limited). Very few genotype × nitrogen interactions were detected, and a large proportion of the QTL were stable across nitrogen nutrition conditions. In contrast, strong genotype × trial interactions in which most of the QTL were specific to a single trial were found. To obtain further insight into the QTL × environment interactions, genetic analyses of ecovalence were performed to identify the genomic regions contributing to the genotype × nitrogen and genotype × trial interactions. Fifty-one critical genomic regions contributing to the additive genetic control of yield-associated traits were identified, and the structural organization of these regions in the genome was investigated. CONCLUSIONS: Our results demonstrated that the effect of the trial was greater than the effect of nitrogen nutrition levels on seed yield-related traits under our experimental conditions. Nevertheless, critical genomic regions associated with yield that were stable across environments were identified in rapeseed.