<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Aging</journal-id>
<journal-title>Frontiers in Aging</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Aging</abbrev-journal-title>
<issn pub-type="epub">2673-6217</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1644012</article-id>
<article-id pub-id-type="doi">10.3389/fragi.2025.1644012</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Aging</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Aging-dependent skin microbiome alterations across body sites in a United Kingdom cohort</article-title>
<alt-title alt-title-type="left-running-head">Swaney et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fragi.2025.1644012">10.3389/fragi.2025.1644012</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Swaney</surname>
<given-names>Mary Hannah</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3094850/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Newman</surname>
<given-names>Duncan J.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3095785/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Mao</surname>
<given-names>Junhong</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1204710/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Hilton</surname>
<given-names>Anthony C.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Worthington</surname>
<given-names>Tony</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Li</surname>
<given-names>Min</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2200916/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Colgate-Palmolive Company</institution>, <addr-line>Piscataway</addr-line>, <addr-line>NJ</addr-line>, <country>United States</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>College of Health and Life Sciences, Aston University</institution>, <addr-line>Birmingham</addr-line>, <country>United Kingdom</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/37407/overview">John Tower</ext-link>, University of Southern California, United States</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/381929/overview">Danay Saavedra</ext-link>, University of Miami, United States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/959610/overview">Judy Malas</ext-link>, University of Illinois Chicago, United States</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Min Li, <email>min_li@colpal.com</email>
</corresp>
<fn fn-type="equal" id="fn001">
<label>
<sup>&#x2020;</sup>
</label>
<p>These authors have contributed equally to this work and share first authorship</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>19</day>
<month>09</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="collection">
<year>2025</year>
</pub-date>
<volume>6</volume>
<elocation-id>1644012</elocation-id>
<history>
<date date-type="received">
<day>09</day>
<month>06</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>09</day>
<month>09</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2025 Swaney, Newman, Mao, Hilton, Worthington and Li.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Swaney, Newman, Mao, Hilton, Worthington and Li</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>The aging process profoundly influences not only the health and visual appearance of the skin, but also the composition of the microbial communities residing on its surface.</p>
</sec>
<sec>
<title>Methods</title>
<p>To investigate these microbial changes, we employed a comprehensive, multi-scale approach that probes community composition, species interactions, and predicted metabolic function of the skin microbiome of the face and forearm in young and old age individuals from the United Kingdom using 16S rRNA gene sequencing.</p>
</sec>
<sec>
<title>Results</title>
<p>Our findings revealed significant and site-specific age-related shifts in the microbiome involving diversity, interpersonal heterogeneity, network connectivity, and metabolic potential, suggesting loss of microbiome robustness and a shift towards a hyperdiversified, fragile microbial community in old age. Furthermore, we applied Dirichlet Multinomial Mixtures to uncover novel age-driven microbiome profiles unique across each skin site, highlighting <italic>Cutibacterium acnes</italic>, <italic>Staphylococcus hominis</italic>, and microbial community diversity as key differentiating biomarkers of the skin microbiome across the lifespan.</p>
</sec>
<sec>
<title>Discussion</title>
<p>Overall, through examining the aging skin microbiome from a systems perspective, our study reinforces and enhances the findings from previous aging microbiome studies and underscores the importance of site-specific differences in skin microbiome dynamics with age. These insights suggest that microbial interventions could mitigate age-related changes, enhancing skin health and wellbeing throughout life.</p>
</sec>
</abstract>
<kwd-group>
<kwd>skin microbiome</kwd>
<kwd>aging skin</kwd>
<kwd>aging skin microbiome</kwd>
<kwd>microbial diversity</kwd>
<kwd>skin</kwd>
</kwd-group>
<counts>
<page-count count="15"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Genetics, Genomics and Epigenomics of Aging</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>The human skin microbiome is the collection of microorganisms that inhabits the skin and carries out beneficial and necessary processes to support skin health (<xref ref-type="bibr" rid="B7">Byrd et al., 2018</xref>). Starting at birth, the skin microbiota play a critical role in the development and homeostatic maintenance of the immune system and skin physiology over one&#x2019;s lifetime. For example, they influence key skin functions that include acidification, antimicrobial defense, lipid synthesis, and barrier integrity (<xref ref-type="bibr" rid="B7">Byrd et al., 2018</xref>; <xref ref-type="bibr" rid="B3">Almoughrabie et al., 2023</xref>; <xref ref-type="bibr" rid="B85">Zheng et al., 2022</xref>; <xref ref-type="bibr" rid="B42">Liu et al., 2023</xref>). Balance of the skin microbiome is crucial; an imbalance, or dysbiosis, can disrupt these functions and is associated with various conditions including atopic dermatitis and acne vulgaris (<xref ref-type="bibr" rid="B54">Nakatsuji and Gallo, 2018</xref>; <xref ref-type="bibr" rid="B28">Huang et al., 2023</xref>). Therefore, supporting the host-microbiome interactions at the skin surface is critical for maintaining the delicate balance between health and dermatological conditions.</p>
<p>Aging is a time-dependent process characterized by a decline in functional capacity that affects nearly all living organisms (<xref ref-type="bibr" rid="B43">L&#xf3;pez-Ot&#xed;n et al., 2013</xref>; <xref ref-type="bibr" rid="B24">Guo et al., 2022</xref>). It involves not only a decline in tissue and cellular function, but also broader deterioration of physiological processes, which can lead to increased inflammation and visible alterations in appearance, particularly in the skin (<xref ref-type="bibr" rid="B44">L&#xf3;pez-Ot&#xed;n et al., 2023</xref>; <xref ref-type="bibr" rid="B81">Tobin, 2017</xref>). The skin experiences significant transformations during aging due to intrinsic factors like reduced epidermal thickness, slower cell turnover, decreased collagen production, and changes in immune function, as well as extrinsic factors such as sun exposure, pollutants, and smoking (<xref ref-type="bibr" rid="B19">Farage et al., 2008</xref>). These factors collectively alter the physiological properties of the skin, leading to variations in pH, lipid and amino acid composition, and reduced sweat and sebum production (<xref ref-type="bibr" rid="B10">Choi, 2019</xref>; <xref ref-type="bibr" rid="B73">Schreiner et al., 2000</xref>; <xref ref-type="bibr" rid="B83">Waller and Maibach, 2006</xref>). Such differences affect the skin microbiota, which rely on the skin as a microbial niche and function to support a healthy skin barrier.</p>
<p>Aging skin undergoes significant alterations in its microbiome compared to that of younger individuals, with the skin microbiome reported to be an accurate predictor of chronological age, more so than the gut or oral microbiomes (<xref ref-type="bibr" rid="B27">Huang et al., 2020</xref>). With age, there is a notable rise in alpha diversity, which often coincides with a reduction in <italic>Cutibacterium</italic> (<xref ref-type="bibr" rid="B75">Shibagaki et al., 2017</xref>; <xref ref-type="bibr" rid="B30">Jug&#xe9; et al., 2018</xref>; <xref ref-type="bibr" rid="B26">Howard et al., 2022</xref>; <xref ref-type="bibr" rid="B84">Ying et al., 2015</xref>; <xref ref-type="bibr" rid="B34">Kim et al., 2022</xref>; <xref ref-type="bibr" rid="B37">Larson et al., 2022</xref>). This shift is further marked by a rise in <italic>Corynebacterium</italic> and <italic>Streptococcus</italic> populations (<xref ref-type="bibr" rid="B75">Shibagaki et al., 2017</xref>; <xref ref-type="bibr" rid="B30">Jug&#xe9; et al., 2018</xref>; <xref ref-type="bibr" rid="B26">Howard et al., 2022</xref>). The specific bacteria at the species level that differentiate young and aged skin can vary significantly from study to study, reflecting the substantial impact of population demographics, environmental factors, and methodological differences on the identification of age-related species of interest (<xref ref-type="bibr" rid="B34">Kim et al., 2022</xref>; <xref ref-type="bibr" rid="B37">Larson et al., 2022</xref>; <xref ref-type="bibr" rid="B51">Myers et al., 2023</xref>; <xref ref-type="bibr" rid="B22">Garlet et al., 2024</xref>). Despite this variability, numerous studies consistently identify <italic>Cutibacterium acnes</italic> and <italic>Corynebacterium kroppenstedtii</italic> as key markers of age-related microbiome changes (<xref ref-type="bibr" rid="B37">Larson et al., 2022</xref>; <xref ref-type="bibr" rid="B22">Garlet et al., 2024</xref>; <xref ref-type="bibr" rid="B16">Dimitriu et al., 2019</xref>; <xref ref-type="bibr" rid="B87">Zhou et al., 2023</xref>). <italic>C. acnes</italic> tends to decrease with age, while <italic>C. kroppenstedtii</italic> increases, underscoring their roles as key indicators of microbiome shifts that could influence skin health and the aging process. Furthermore, the skin microbiome is generally stable and resilient to environmental exposures during adulthood (<xref ref-type="bibr" rid="B56">Oh et al., 2016</xref>). However, in old age, there is an observed increase in interpersonal variation in the skin microbiome, indicating decline of microbiome stability and robustness (<xref ref-type="bibr" rid="B30">Jug&#xe9; et al., 2018</xref>; <xref ref-type="bibr" rid="B26">Howard et al., 2022</xref>; <xref ref-type="bibr" rid="B84">Ying et al., 2015</xref>; <xref ref-type="bibr" rid="B37">Larson et al., 2022</xref>; <xref ref-type="bibr" rid="B22">Garlet et al., 2024</xref>).</p>
<p>While many studies have thoroughly evaluated the microbiome of the face as it ages (<xref ref-type="bibr" rid="B27">Huang et al., 2020</xref>; <xref ref-type="bibr" rid="B75">Shibagaki et al., 2017</xref>; <xref ref-type="bibr" rid="B30">Jug&#xe9; et al., 2018</xref>; <xref ref-type="bibr" rid="B26">Howard et al., 2022</xref>; <xref ref-type="bibr" rid="B34">Kim et al., 2022</xref>; <xref ref-type="bibr" rid="B37">Larson et al., 2022</xref>; <xref ref-type="bibr" rid="B51">Myers et al., 2023</xref>; <xref ref-type="bibr" rid="B22">Garlet et al., 2024</xref>; <xref ref-type="bibr" rid="B16">Dimitriu et al., 2019</xref>; <xref ref-type="bibr" rid="B87">Zhou et al., 2023</xref>), there is a gap in comprehensive understanding of the aging microbiome of non-facial skin sites. Characterizing the aging skin microbiome across the body is crucial for uncovering insights into susceptibility to infections, skin conditions, and overall skin health. Moreover, as consumer interest in the &#x201c;skinification&#x201d; of beauty and personal care products grows, there is a demand for targeted solutions that address diverse aspects of skin health, including anti-aging treatments and personalized care. As the skin microbiome is a pivotal component of skin health, evaluating the microbiome variations with age across different body sites is essential. This broader understanding presents a unique opportunity for developing innovative interventions and products that can enhance skin vitality and wellbeing across the lifespan. Expanding this research not only contributes to a more holistic understanding of aging and skin health but also underscores the potential for microbiome-focused interventions to mitigate age-related skin changes.</p>
<p>In this study, we characterized the skin microbiome of the face and forearm in young and old aged individuals from a UK-based population using 16S rRNA gene sequencing. The objective of our study was to evaluate the differences in composition, species interactions, and predicted metabolic function of the bacterial community with age across two distinct skin sites. By focusing on a UK-based population, we aimed to provide insights into aging-related skin microbiome differences of this geographical region, which is relatively understudied compared to North America- and Asia-based populations. We employed comprehensive techniques including differential abundance testing, microbial association network analysis, microbiome clustering, and functional prediction to identify age-related patterns in the skin microbiome and to enhance our understanding of the skin microbiome&#x2019;s interplay with the aging process.</p>
</sec>
<sec sec-type="results" id="s2">
<title>Results</title>
<p>A total of 59 participants from Birmingham, United Kingdom were enrolled in the study, including thirty young age (YA) (26.7 &#xb1; 4.45 years) and twenty-nine old age (OA) (72.3 &#xb1; 4.04 years) individuals (<xref ref-type="sec" rid="s12">Supplementary Table S1</xref>). The gender distribution was similar for both age groups (YA: n &#x3d; 12 male and n &#x3d; 18 female; OA: n &#x3d; 12 male and n &#x3d; 17 female). Other participant demographics and clinical characteristics are summarized in <xref ref-type="sec" rid="s12">Supplementary Table S1</xref>. Skin microbiome samples were collected from the antecubital fossa (arm) and cheek (face) of the participants and were sequenced using V1-V3 16S rRNA gene sequencing, yielding a median of 105,646 high-quality sequences per sample across both sites.</p>
<sec id="s2-1">
<title>Microbiome diversity and composition vary by age and skin site</title>
<p>Alpha diversity metrics Shannon diversity and Pielou evenness were significantly higher in the old age group compared to the young age group on the face (<italic>p</italic> &#x3d; 0.043 and <italic>p</italic> &#x3d; 0.023, respectively), with no significant differences in richness (number of observed amplicon sequence variants (ASVs)) (<xref ref-type="fig" rid="F1">Figure 1</xref>). This indicates that the increase in diversity is due to a more even microbial community rather than a greater number of taxa. In contrast, no significant differences in richness, evenness, or Shannon diversity were observed on the arm with age (<xref ref-type="fig" rid="F1">Figure 1</xref>), highlighting that the general pattern of increased alpha diversity of the skin microbiome with age may be skin site-specific. Indeed, the majority of skin microbiome studies reporting significant increases in alpha diversity with age observed these differences on the face (<xref ref-type="bibr" rid="B75">Shibagaki et al., 2017</xref>; <xref ref-type="bibr" rid="B30">Jug&#xe9; et al., 2018</xref>; <xref ref-type="bibr" rid="B34">Kim et al., 2022</xref>; <xref ref-type="bibr" rid="B51">Myers et al., 2023</xref>; <xref ref-type="bibr" rid="B22">Garlet et al., 2024</xref>; <xref ref-type="bibr" rid="B40">Li et al., 2020</xref>; <xref ref-type="bibr" rid="B76">Sun et al., 2024</xref>; <xref ref-type="bibr" rid="B33">Kim et al., 2019</xref>). We also assessed differences in alpha diversity attributed to gender. For the arm, there were no significant age-related differences for either gender, similar to the findings for the overall cohort (<xref ref-type="sec" rid="s12">Supplementary Figures S1A&#x2013;C</xref>). Conversely, for the face, only female participants exhibited a significant increase in Shannon diversity and Pielous evenness, indicating that gender may be a driving factor for the observed facial microbial diversity increase with age.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Increased alpha diversity in the face microbiome with age. Alpha diversity metrics were calculated and compared between the old age (OA) group and young age (YA) group for the face and arm. Calculations were performed on the ASV feature table and include: <bold>(A)</bold> Observed ASVs, <bold>(B)</bold> Pielou evenness, and <bold>(C)</bold> Shannon diversity. YA (blue) represents the young age group, and OA (red) represents the old age group. Points represent data outliers. Statistical significance was calculated using the Wilcoxon rank sum test (&#x2a; <italic>p</italic>-value &#x3c;0.05).</p>
</caption>
<graphic xlink:href="fragi-06-1644012-g001.tif">
<alt-text content-type="machine-generated">Box plots showing microbial diversity comparisons between young adults (YA) and older adults (OA) for face and arm samples. Panel A displays observed ASVs, Panel B shows Pielou Evenness with significant difference marked, and Panel C illustrates the Shannon Diversity Index. Blue represents YA and red represents OA groups.</alt-text>
</graphic>
</fig>
<p>In assessing beta diversity, microbial community composition differed significantly between age groups for both the face and the arm using Bray-Curtis dissimilarity and unweighted UniFrac (p &#x3c; 0.05), with notably higher interpersonal variation in the old age group compared to the young age group (Bray-Curtis, p &#x3c; 0.001) (<xref ref-type="fig" rid="F2">Figure 2</xref>). Weighted UniFrac exhibited significant differences in community composition for the face but not the arm between age groups (<xref ref-type="sec" rid="s12">Supplementary Figures S2A, B</xref>). These findings suggest that differences in microbial communities between age groups for the face may be driven by variations in species abundances and evolutionary relationships among abundant taxa, as well as the presence/absence of low-abundance, phylogenetically distinct taxa. In contrast, for the arm, the compositional differences may be primarily due to shifts in the abundances of closely related taxa and the presence/absence of specific, phylogenetically unique taxa. When examining community composition across all samples, significant differences in Bray-Curtis dissimilarity were observed between the face and arm (p &#x3d; 0.001, R<sup>2</sup> &#x3d; 0.0679) as well as between young age and old age (p &#x3d; 0.001, R<sup>2</sup> &#x3d; 0.0287), with skin site having a larger impact on community composition differences than age group (<xref ref-type="sec" rid="s12">Supplementary Figure S2C</xref>). Furthermore, evaluation of microbial community composition by gender revealed significant differences in Bray-Curtis dissimilarity for female participants across both the face and arm, however no significant differences were observed for male participants (<xref ref-type="sec" rid="s12">Supplementary Figures S1D&#x2013;G</xref>), suggesting that gender plays an important role in driving skin microbial community differences with age.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Significant differences in microbial community composition with age in the face and arm microbiome. Microbial community composition shows significant variation between young age (YA) and old age (OA) groups in both the face and arm. <bold>(A,D)</bold> Bray-Curtis dissimilarities between individuals within the same and different age groups for the face and arm, respectively. Statistical significance was calculated using the Wilcoxon rank sum test (&#x2a;&#x2a;&#x2a;&#x2a; <italic>p</italic>-value &#x3c;0.0001, ns &#x3d; not significant). <bold>(B,E)</bold> NMDS ordination of the Bray-Curtis dissimilarity matrix for the face and arm, respectively. <bold>(C,F)</bold> NMDS ordination of the unweighted UniFrac distance matrix for the face and arm, respectively. Statistical significance for compositional differences was calculated using PERMANOVA, with the <italic>p</italic>-value indicated within each ordination plot. Ellipses represent a 95% confidence interval for each age group. Blue circles represent the young age (YA) group and red triangles represent the old age (OA) group.</p>
</caption>
<graphic xlink:href="fragi-06-1644012-g002.tif">
<alt-text content-type="machine-generated">Box plots and scatter plots comparing skin microbial diversity between young adults (YA) and older adults (OA) on the face and arm. Plots A and D show Bray-Curtis dissimilarity indices, indicating significant differences on the face but not on the arm. Plots B and E depict NMDS plots using Bray-Curtis, showing clustering differences with significance values (p=0.011 for face, p=0.003 for arm). Plots C and F show NMDS plots using Unweighted UniFrac, with face and arm showing distinct clustering (p=0.03 for face, p=0.011 for arm).</alt-text>
</graphic>
</fig>
</sec>
<sec id="s2-2">
<title>Age-related taxonomic shifts across skin sites</title>
<p>At the phylum level, the face demonstrated a significant decrease with age in Actinobacteriota (Actinobacteria) alongside a significant increase in Pseodomonadota (Proteobacteria) and taxa outside of the Actinobacteriota, Pseudomonadota, and Bacillota (Firmicutes) phyla (<xref ref-type="sec" rid="s12">Supplementary Figure S3</xref>). The arm did not exhibit significant age-related phylum-level differences. At the genus level, microbial community composition exhibited high interpersonal variation (<xref ref-type="fig" rid="F3">Figure 3A</xref>). Among the top three most abundant genera on the skin, <italic>Cutibacterium</italic>, <italic>Staphylococcus</italic>, and <italic>Corynebacterium</italic>, both skin sites showed a significant decrease in <italic>Cutibacterium</italic> in the old age group (Face <italic>p</italic> &#x3d; 0.033, Arm <italic>p</italic> &#x3d; 0.007), with minimal differences in <italic>Staphylococcus</italic> and <italic>Corynebacterium</italic> abundances (<xref ref-type="fig" rid="F3">Figure 3B</xref>). Interestingly, this pattern was gender-dependent; only female participants exhibited a significant decrease in <italic>C. acnes</italic> with age for both skin sites, in addition to a significant decrease in <italic>Staphylococcus</italic> on the face (<xref ref-type="sec" rid="s12">Supplementary Figure S4</xref>).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Genus- and species-level variation of abundant skin taxa between age groups. Significant differences are observed between age groups for abundant genera and species on the skin. <bold>(A)</bold> Relative proportions of the top twenty genera are presented as stacked bar plots for each sample within each age group and skin site, with samples ordered by <italic>Cutibacterium</italic> abundance from highest to lowest. Genera not included within the top twenty genera are grouped into &#x201c;Other&#x201d;. <bold>(B)</bold> Relative abundance differences between <italic>Cutibacterium, Staphylococcus</italic>, and <italic>Corynebacterium</italic> within each age group and skin site. Genera are colored according to the legend in <bold>(A)</bold>. <bold>(C)</bold> Relative abundance differences for species of interest for the <bold>(C)</bold> face and <bold>(D)</bold> arm. Statistical significance was calculated using the Wilcoxon rank sum test &#x2a; <italic>p</italic>-value &#x3c;0.05, &#x2a;&#x2a; <italic>p</italic>-value &#x3c;0.01. YA &#x3d; young age group, OA &#x3d; old age group.</p>
</caption>
<graphic xlink:href="fragi-06-1644012-g003.tif">
<alt-text content-type="machine-generated">Panel A shows stacked bar charts of microbial genus distribution on face and arm samples from young adults (YA) and older adults (OA). Panel B features violin plots of relative abundance of Cutibacterium, Staphylococcus, and Corynebacterium across age and body sites. Panels C and D include violin plots for specific species like Cutibacterium acnes and Corynebacterium kroppenstedtii, detailing differences between YA and OA for face and arm. Genus and species are color-coded.</alt-text>
</graphic>
</fig>
<p>Differential abundance testing identified 4 ASVs that were differentially abundant between age groups on the face: <italic>C. kroppenstedtii, Phocaicola vulgatus, Streptococcus</italic> sp., and <italic>Streptococcus thermophilus</italic> (<italic>q</italic> &#x3c; 0.05) (<xref ref-type="sec" rid="s12">Supplementary Figure S5A</xref>)<italic>.</italic> On the arm, 19 ASVs were differentially abundant: <italic>Corynebacterium amycolatum, Corynebacterium dentalis,</italic> five <italic>Cutibacterium acnes</italic> ASVs, one <italic>Lawsonella</italic> ASV, <italic>Moraxella cinereus</italic>, three <italic>Staphylococcus capitis</italic> ASVs, four <italic>Staphylococcus hominis</italic> ASVs, two <italic>Staphylococcus</italic> ASVs, and a <italic>Weeksellaceae</italic> ASV (<italic>q</italic> &#x3c; 0.05) (<xref ref-type="sec" rid="s12">Supplementary Figure S5B</xref>).</p>
<p>All ASVs assigned to <italic>C. acnes, C. kroppenstedtii, S. hominis,</italic> and <italic>S. capitis,</italic> which were selected based on their differential abundance results as well as relevance to skin health, were aggregated to the species level for further analysis. <italic>C. acnes</italic> has been previously reported to decrease with age, attributed to a decrease in sebum present on the skin (<xref ref-type="bibr" rid="B37">Larson et al., 2022</xref>; <xref ref-type="bibr" rid="B22">Garlet et al., 2024</xref>; <xref ref-type="bibr" rid="B87">Zhou et al., 2023</xref>). We observed that <italic>C. acnes</italic> is significantly reduced both on the face and the arm (&#x2212;16.1% and &#x2212;11.4% relative abundance, respectively) in the old age group compared to the young age group (Face <italic>p</italic> &#x3d; 0.024, Arm <italic>p</italic> &#x3d; 0.0078) (<xref ref-type="fig" rid="F3">Figures 3C,D</xref>). Similarly, <italic>C. kroppenstedtii</italic> has also been previously found to increase on the skin with age and has been reported to be associated with rosacea (<xref ref-type="bibr" rid="B22">Garlet et al., 2024</xref>; <xref ref-type="bibr" rid="B16">Dimitriu et al., 2019</xref>; <xref ref-type="bibr" rid="B87">Zhou et al., 2023</xref>; <xref ref-type="bibr" rid="B64">Rainer et al., 2020</xref>). On the face, <italic>C. kroppenstedtii</italic> was significantly increased in the old age compared to the young age group (<italic>p</italic> &#x3d; 0.022, &#x2b;3.51% mean relative abundance) (<xref ref-type="fig" rid="F3">Figure 3C</xref>). Lastly, we examined levels of <italic>S. capitis</italic> and <italic>S. hominis</italic>, as several ASVs assigned to these species were differentially abundant between age groups on the arm (<xref ref-type="sec" rid="s12">Supplementary Figure S5B</xref>)<italic>. S. capitis</italic> showed a significant decrease in old age compared to young age (<italic>p</italic> &#x3d; 0.026, &#x2212;2.12% mean relative abundance), while <italic>S. hominis</italic> trended towards an increase in old age (<italic>p</italic> &#x3d; 0.15, &#x2b;9.03% mean relative abundance) (<xref ref-type="fig" rid="F3">Figure 3D</xref>). Assessing differences by gender revealed that the significant decrease with age observed for <italic>C. acnes</italic> on the face and arm and <italic>S. capitis</italic> on the arm is largely driven by the female participants (<xref ref-type="sec" rid="s12">Supplementary Figures S4B, C</xref>). However, <italic>S. hominis</italic> shows an increase in old age for both genders, suggesting a gender-independent mechanism. Overall, numerous taxa demonstrate age-specific changes, which notably are distinct between skin sites.</p>
</sec>
<sec id="s2-3">
<title>Loss of microbial community network structure in old age</title>
<p>To evaluate the impact of age and skin site on microbial interactions and the overall microbial community network, we employed the statistical method SPIEC-EASI to infer microbial associations between skin microbiome members of the different age groups and skin sites. We selected features found in at least 30% of all face or arm samples, with a feature in this instance referring to ASVs agglomerated to the species level or retained at their original taxonomic level if a species-level assignment was unavailable. This resulted in 194 features for the face and 144 features for the arm being included in the analysis.</p>
<p>For both skin sites, the microbial networks in old age had fewer nodes and edges compared to the young age networks, indicating fewer associating taxa within the network with age (<xref ref-type="fig" rid="F4">Figures 4A&#x2013;C,F&#x2013;H</xref>). In addition, the old age networks showed reduced degree assortativity and transitivity [a measure of the preference for a node to attach to other nodes with similar degree count and the overall probability for the network to have adjacent nodes interconnected, respectively (<xref ref-type="bibr" rid="B68">R&#xf6;ttjers and Faust, 2018</xref>)], increased modularity [a measure of the strength of division of a network into modules (<xref ref-type="bibr" rid="B68">R&#xf6;ttjers and Faust, 2018</xref>)], and a lower frequency of high-degree nodes compared to the young age group (<xref ref-type="fig" rid="F4">Figures 4C,D,H&#x2013;I</xref>). Both old age networks had significantly reduced node degrees compared to the young age networks (<italic>p</italic> &#x3c; 0.001) (<xref ref-type="fig" rid="F4">Figures 4E,J</xref>). The face network in young age showed particularly high connectivity, with all 194 features (nodes) having at least one association with other features. In contrast, the old age network for the face exhibited a less dense structure, with only 171 associating features. Similarly, the arm network in old age showed a marked shift from the young-age network, presenting a highly sparse network composed of numerous modules with limited connectivity (121 vs. 75 associating features in young versus old age, respectively). This reduction in connectivity could indicate a microbial community in old age that exhibits reduced interaction and cooperation among microbial members. Interestingly, differentially-abundant and skin health-related taxa that showed significant differences with age (<xref ref-type="fig" rid="F3">Figure 3</xref>; <xref ref-type="sec" rid="s12">Supplementary Figure S5</xref>) also were more prevalent and connected within the network in young age compared to old age for both skin sites. For example, <italic>C. acnes</italic> associated with numerous other taxa within the young age networks, but only formed a single association in the old age networks. Overall, these findings suggest that in old age, the microbial communities of the face and arm become less robust, less stable, and more fragile.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Microbial association networks show increased sparsity and reduced connectivity with age. The SPIEC-EASI statistical method was used to identify microbial associations between skin microbiome members for each age group and skin site. Networks for the face in the <bold>(A)</bold> young age group and <bold>(B)</bold> old age group are shown, alongside <bold>(C)</bold> network statistics (number of nodes, number of edges, degree assortativity, transitivity, and modularity), <bold>(D)</bold> degree frequency distributions, and <bold>(E)</bold> degrees per node comparison for the face. Similarly, networks for the arm in the <bold>(F)</bold> young age group and <bold>(G)</bold> old age group are presented, with corresponding <bold>(H)</bold> network statistics, <bold>(I)</bold> degree frequency distributions, and <bold>(J)</bold> degrees per node comparison for the arm. Nodes represent microbial features (ASVs agglomerated to the species level or ASVs retained at their original taxonomy level if a species-level assignment was unavailable), which are colored by their phylum assignment and sized by their mean relative abundance within the skin site and age group collectively. Edges represent associations identified between features. Only features with at least one association are included in the networks. Nodes representing taxa of interest or differentially abundant taxa are indicated in black text. Statistical significance was calculated using the Wilcoxon rank sum test (&#x2a;&#x2a;&#x2a;&#x2a; <italic>p</italic>-value &#x3c;0.0001). YA &#x3d; young age group, OA &#x3d; old age group.</p>
</caption>
<graphic xlink:href="fragi-06-1644012-g004.tif">
<alt-text content-type="machine-generated">Microbial community networks on face and arm for young and old age groups. Panels A, B, F, and G display network diagrams showing microbial interactions. Panel C and H present bar graphs comparing nodes, edges, assortativity, transitivity, and modularity between age groups. Panels D and I feature frequency distribution plots for degrees of connection. Panels E and J show box plots comparing degrees between young and old age groups. Bubble colors represent different microbial phyla, and bubble sizes indicate relative abundance.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s2-4">
<title>DMM clustering reveals age-related microbiome profiles</title>
<p>To explore the relationship between microbial community composition and age, we performed mathematical modeling of ASV frequencies using Dirichlet multinomial mixtures (DMM) to uncover potential age-related microbiome structures across participants. DMM clustering separated the face samples into two community types with significantly different microbial community compositions (Bray-Curtis dissimilarity, PERMANOVA <italic>p</italic> &#x3d; 1E-04) (<xref ref-type="fig" rid="F5">Figure 5A</xref>). Community Type 1 exhibited high Shannon diversity, with <italic>C. acnes</italic> and <italic>Staphylococcus epidermidis</italic> being the primary drivers of this community type (<xref ref-type="fig" rid="F5">Figures 5B,C</xref>; <xref ref-type="sec" rid="s12">Supplementary Figure S6A</xref>). In contrast, Community Type 2 was significantly less diverse than Community Type 1 (<italic>p</italic> &#x3c; 0.001) and was dominated by <italic>C. acnes</italic> (<xref ref-type="sec" rid="s12">Supplementary Figure S6B</xref>). Community Type 2 had a predominance of young age group individuals (66.7%) compared to 33.3% from the old age group, whereas Community Type 1 was more balanced with 41.2% from the young age group and 58.8% from the old age group (<xref ref-type="fig" rid="F5">Figure 5D</xref>). When comparing facial community type distribution by gender for both Community Type 1 and Community Type 2, female participants mirrored the overall cohort&#x2019;s age group distribution, while the male participants had a relatively equal age group distribution (<xref ref-type="sec" rid="s12">Supplementary Figure S7A</xref>), suggesting that gender is an important driving factor in face microbiome shifts with age.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Dirichlet multinomial mixtures (DMM) clustering identifies age-specific microbiome profiles for the face and arm. <bold>(A,E)</bold> NMDS ordination of the Bray-Curtis dissimilarity matrix for the face and arm, respectively. Points are colored and shaped by their DMM community type. Ellipses represent a 95% confidence interval for each age group. <bold>(B,F)</bold> Stacked bar plot representing the genus-level microbiome community composition of each community type as determined through Dirichlet multinomial mixtures (DMM) in the face and arm samples, respectively. Taxa not within the top twenty most abundant genera are grouped into &#x201c;Other&#x201d;. <bold>(C,G)</bold> Shannon diversity for the community types in face and arm samples, respectively. Statistical significance was calculated using the Wilcoxon rank sum test (&#x2a;&#x2a;&#x2a;&#x2a; <italic>p</italic>-value &#x3c;0.0001). <bold>(D,H)</bold> Number of samples that are assigned to each community type, colored by age group, in the face and arm samples, respectively. YA (blue) &#x3d; young age, OA (red) &#x3d; old age.</p>
</caption>
<graphic xlink:href="fragi-06-1644012-g005.tif">
<alt-text content-type="machine-generated">This figure presents four panels in two rows. The top row (A-D) depicts data for facial microbial communities: A) NMDS scatter plot shows two community types, B) bar chart illustrates bacterial genus relative proportions for each community type, C) box plot of Shannon diversity index comparing the community types, D) bar chart of sample numbers by age group. The bottom row (E-H) shows arm microbial communities: E) NMDS scatter plot with three community types, F) bar chart of bacterial genus proportions for each community, G) box plot of Shannon diversity index among the community types, H) bar chart of sample numbers by age group.</alt-text>
</graphic>
</fig>
<p>For the arm, DMM identified three distinct community types, which exhibited significant differences in microbial community composition (Bray-Curtis dissimilarity, PERMANOVA <italic>p</italic> &#x3d; 1E-04) (<xref ref-type="fig" rid="F5">Figure 5E</xref>). Community Type 1 exhibited high Shannon diversity and was predominantly driven by <italic>C. acnes</italic>, <italic>S. hominis</italic>, and <italic>S. epidermidis,</italic> appearing in 68% of samples from the old age group (<xref ref-type="fig" rid="F5">Figures 5F&#x2013;H</xref>; <xref ref-type="sec" rid="s12">Supplementary Figure S6C</xref>). Community Type 2 demonstrated similar diversity to Community Type 1, however was largely driven by <italic>C. acnes</italic> alone (<xref ref-type="sec" rid="s12">Supplementary Figure S6D</xref>). Notably, all but one sample assigned to Community Type 2 was from the young age group. Lastly, Community Type 3, dominated by <italic>S. homini</italic>s, showed significantly reduced Shannon diversity compared to the other community types (<italic>p</italic> &#x3c; 0.001), with 83.3% of samples from the old age group (<xref ref-type="sec" rid="s12">Supplementary Figure S6E</xref>). When comparing by gender, the age distribution patterns for both male and female participants mirrored the overall distribution, with Community Types 1 and 3 being old age-dominant and Community Type 2 being young-age dominant (<xref ref-type="sec" rid="s12">Supplementary Figure S7B</xref>). This suggests that for the arm, age-related microbial community shifts are gender-independent.</p>
</sec>
<sec id="s2-5">
<title>The face microbiome exhibits a shift in predicted functional potential with age</title>
<p>PICRUSt2 was utilized to assess the functional potential of the face and arm samples across age groups, predicting MetaCyc metabolic pathways based on ASV-level data. For the face, there was a significant shift in the predicted functional profile between young age and old age (<italic>p</italic> &#x3d; 0.006), and numerous metabolic pathways were found to be more abundant in old age compared to young age (<italic>q</italic> &#x3c; 0.05) (<xref ref-type="fig" rid="F6">Figure 6A</xref>). These include functions related to aromatic hydrocarbon and pollutant degradation (aerobic toluene degradation III via p-cresol, superpathway of salicylate degradation, protocatechuate degradation II, aerobic benzoyl-CoA degradation I, and 4-methylcatechol degradation), amino acid degradation (L-tyrosine and L-histidine), and lipopolysaccharide biosynthesis (ADP-L-glycero-&#x3b2;-D-manno-heptose biosynthesis) (<xref ref-type="fig" rid="F6">Figure 6B</xref>). In contrast, there was not an observed shift in predicted metabolic function of the microbial communities between age groups for the arm (<xref ref-type="fig" rid="F6">Figure 6C</xref>). Two pathways, coenzyme M biosynthesis and chitin derivatives degradation, were significantly more abundant in old age compared to young age (<italic>q</italic> &#x3c; 0.05), although the relative abundances of these pathways were very low (<xref ref-type="fig" rid="F6">Figure 6D</xref>).</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Differences in predicted microbial functional potential are more prominent in the face than the arm in old age. <bold>(A,C)</bold> NMDS ordination of the Bray-Curtis dissimilarity matrix computed on MetaCyc pathway relative abundances for the face and arm samples, respectively. Density plots within the plot margins represent the density of points along the x and y-axes for each age group. <bold>(B,D)</bold> Relative abundances of MetaCyc pathways determined to be differentially abundant between young age (YA) and old age (OA) groups for the face and arm, respectively. All pathways shown are statistically significant with <italic>q</italic>-value &#x3c;0.05. YA (blue) &#x3d; young age, OA (red) &#x3d; old age. Statistical significance for pathway compositional differences was calculated using PERMANOVA, with the <italic>p</italic>-value indicated within each ordination plot.</p>
</caption>
<graphic xlink:href="fragi-06-1644012-g006.tif">
<alt-text content-type="machine-generated">Panel A shows an NMDS plot comparing two age groups, YA (blue circles) and OA (red triangles), with distinct clustering and a p-value of 0.006. Panel B presents a bar chart of Metacyc Pathway Descriptions and relative abundances for pathways like toluene degradation and nitrate reduction, with differences between age groups. Panel C depicts another NMDS plot with no significant clustering (p-value 0.24). Panel D displays a bar chart for the relative abundance of pathways like coenzyme M biosynthesis, showing visible differences between the age groups.</alt-text>
</graphic>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="s3">
<title>Discussion</title>
<p>The aging process directly impacts skin appearance and function, resulting in distinct changes that range from visible signs of aging, such as fine lines and wrinkles, to increased susceptibility to skin conditions and infections (<xref ref-type="bibr" rid="B81">Tobin, 2017</xref>). These physiological changes are accompanied by shifts in the microbial populations that inhabit our skin, necessitating a better understanding of the consequences that these altered microbial communities have on skin health, particularly across different body sites. In our study, we find that the skin microbiome of the face and arm undergo significant yet distinct shifts with age, highlighting the site-specific variation in microbiome dynamics that are impacted by the aging process.</p>
<p>The face is a highly exposed skin region that is chronically exposed to external factors such as UV radiation and air pollution, both of which significantly contribute to signs of skin aging (<xref ref-type="bibr" rid="B35">Krutmann et al., 2017</xref>). We observed strong shifts in the microbiome of the face between young and old age groups, including increased alpha diversity, significant microbiome community composition shifts characterized by a large reduction in <italic>C. acnes</italic>, altered network topology, and predicted functional changes. The arm, while being a less exposed skin site than the face (<xref ref-type="bibr" rid="B5">Bulliard et al., 2007</xref>), also exhibited significant changes in the microbiome between age groups, albeit less pronounced. These included shifts in microbiome community composition characterized by a decrease in <italic>C. acnes</italic> and shifted <italic>Staphylococcus</italic> populations, as well as an altered network structure. Thus, skin sites with differential exposure to environmental aggressors show distinct skin microbiome shifts, suggesting that these differences could be driven in part by differences in intrinsic versus extrinsic skin aging factors. Indeed, differential exposure to polycyclic aromatic hydrocarbons (PAH), a class of organic pollutants in ambient air, has been shown to impact the skin microbiome, metabolome, and clinical skin parameters, with high exposure to PAH associated with biodegrading skin bacteria, dry skin, and hyperpigmentation (<xref ref-type="bibr" rid="B39">Leung et al., 2023</xref>; <xref ref-type="bibr" rid="B38">Leung et al., 2020</xref>; <xref ref-type="bibr" rid="B49">Misra et al., 2021</xref>). Similarly, exposure to UVA and UVB radiation, which have profound impacts on skin aging (<xref ref-type="bibr" rid="B35">Krutmann et al., 2017</xref>), has been observed to immediately alter microbiome composition (<xref ref-type="bibr" rid="B6">Burns et al., 2019</xref>). Therefore, the interplay between the skin, its microbiome, and aging are likely to be multi-faceted and dynamic, influenced by a combination of genetic factors, environmental exposures, and lifestyle choices.</p>
<p>A notable reduction in <italic>C. acnes</italic> was observed in old age for both skin sites (<xref ref-type="fig" rid="F3">Figures 3C,D</xref>), supporting previous studies that similarly found this species to be a biomarker negatively associated with aging (<xref ref-type="bibr" rid="B37">Larson et al., 2022</xref>; <xref ref-type="bibr" rid="B22">Garlet et al., 2024</xref>; <xref ref-type="bibr" rid="B87">Zhou et al., 2023</xref>). Reduced sebum availability likely leads to a decrease in <italic>C. acnes</italic>, which plays a key role in maintaining skin pH through production of free fatty acids from sebum metabolism, contributing to the skin&#x2019;s acidic nature and its protective resistance against unwanted microorganisms (<xref ref-type="bibr" rid="B77">Swaney and Kalan, 2021</xref>). While often associated as an etiological contributor to acne, <italic>C. acnes</italic> is a ubiquitous skin commensal that has key roles in immunomodulation, epithelial barrier regulation, lipid synthesis, and protection against pathogens (<xref ref-type="bibr" rid="B3">Almoughrabie et al., 2023</xref>; <xref ref-type="bibr" rid="B12">Claesen et al., 2020</xref>; <xref ref-type="bibr" rid="B52">Nagy et al., 2006</xref>; <xref ref-type="bibr" rid="B1">Agak et al., 2018</xref>; <xref ref-type="bibr" rid="B53">Nakamura et al., 2020</xref>; <xref ref-type="bibr" rid="B71">Sanford et al., 2016</xref>). As a result, depletion of <italic>C. acnes</italic> is likely to lead to hyperdiversification as a result of increased niche availability and reduced colonization resistance. This can provide rationale for the observed increase in interpersonal variation seen on the face with old age; diminished colonization resistance could allow for an increase in presence and abundance of non-specific taxa, particularly those that are transient and environmental, leading to high interpersonal variability of the microbiome. Indeed, predictive functional analysis of the face microbiome revealed increased abundances of functions involved in aromatic hydrocarbon/pollutant degradation (<xref ref-type="fig" rid="F6">Figure 6</xref>), functions that are distributed among environmental microorganisms due to the widespread presence of these compounds in nature (<xref ref-type="bibr" rid="B21">Fuchs et al., 2011</xref>; <xref ref-type="bibr" rid="B32">Ke et al., 2023</xref>; <xref ref-type="bibr" rid="B61">Porter and Young, 2014</xref>; <xref ref-type="bibr" rid="B58">Parales et al., 2008</xref>). Increased abundance of bacteria with biodegradation capabilities has also been associated with high-level PAH exposure (<xref ref-type="bibr" rid="B39">Leung et al., 2023</xref>), providing a possible link in the present study between pollution exposure, the skin microbiome, and aging.</p>
<p>Our study also revealed gender-specific differences in how the microbiome differs with aging across different body sites. For the face, the microbiome changes observed in female participants largely mirrored the trends seen in the overall cohort results, characterized by a significant shift in community composition and diversity with age (<xref ref-type="sec" rid="s12">Supplementary Figure S1</xref>). This may be attributed to the hormonal shifts experienced by women during menopause, which lead to a notable reduction in sebum production (<xref ref-type="bibr" rid="B88">Zouboulis and Boschnakow, 2001</xref>). Given that sebum availability influences growth of the skin microbiota (<xref ref-type="bibr" rid="B79">Swaney et al., 2023</xref>), these physiological changes are likely important drivers of the pronounced differences observed in the female microbiome with age. Conversely, the typically higher and more stable sebum production in males (<xref ref-type="bibr" rid="B82">Townsend and Kalan, 2023</xref>; <xref ref-type="bibr" rid="B60">Pochi et al., 1979</xref>) may provide a more consistent and robust microbial niche on the face. This could buffer against the age-related microbial changes seen in females, potentially explaining the more limited age-related differences observed in the male facial microbiome. However, on the arm, the aging-related microbial changes were similar across genders, suggesting that the less sebaceous environment of the arm is less influenced by gender-specific physiological factors. These findings underscore the importance of considering both skin site and gender in future studies on the skin microbiome and aging.</p>
<p>Our clustering analysis revealed the existence of age-driven microbiome profiles for both the face and arm. In particular, two facial microbiome types were identified using DMM, with high diversity/low <italic>C. acnes</italic> more frequent in old age and low diversity/high <italic>C. acnes</italic> more prevalent in young age (<xref ref-type="fig" rid="F5">Figures 5A&#x2013;D</xref>). These findings support previous efforts to assign facial skin microbiomes to cutotypes by applying Partitioning Around Medoids (PAM) clustering, which similarly identified a <italic>Cutibacterium-</italic>cutotype that exhibited low Shannon diversity and whose proportion declined with age across individuals (<xref ref-type="bibr" rid="B76">Sun et al., 2024</xref>). Therefore, <italic>C. acnes</italic> may play a significant role in maintaining a core, balanced microbial community in younger individuals, with its decline contributing to increased diversification and susceptibility to perturbation and dysbiosis with age (<xref ref-type="bibr" rid="B37">Larson et al., 2022</xref>). Interestingly, we observed two distinct old age-associated microbiome types for the arm; their compositional profiles were both characterized by a predominance of <italic>S. hominis</italic>, but differed in respect to their overall community diversity (<xref ref-type="fig" rid="F5">Figures 5E&#x2013;H</xref>). This divergence suggests that aging may not only shift the dominant taxa, but may also influence ecological stability and resilience of these communities. Indeed, microbial association analysis revealed a stark shift in structure of the microbial networks between young and old age for both skin sites (<xref ref-type="fig" rid="F4">Figure 4</xref>), highlighting a less robust and more fragile microbiome with age. This has been observed in previous studies, which have similarly observed diminished microbial community stability and resilience in older populations (<xref ref-type="bibr" rid="B34">Kim et al., 2022</xref>; <xref ref-type="bibr" rid="B33">Kim et al., 2019</xref>). Overall, these findings suggest that aging not only disrupts the core, <italic>C. acnes</italic>-dominant microbiome, but also promotes a less predictable and more perturbation-susceptible microbial community that notably is skin-site specific.</p>
<p>Our findings highlight the potential role of <italic>S. capitis</italic> and <italic>S. hominis</italic> as key differentiating species between young and old age on the arm, with <italic>S. capitis</italic> more abundant in young age and <italic>S. hominis</italic> more abundant in old age (<xref ref-type="fig" rid="F3">Figure 3D</xref>; <xref ref-type="sec" rid="s12">Supplementary Figure S6</xref>). <italic>S. capitis</italic>, similar to <italic>C. acnes</italic>, thrives in lipid-rich skin regions (<xref ref-type="bibr" rid="B11">Chong et al., 2022</xref>), which may explain the observed decrease in abundance of this species with age. <italic>S. hominis</italic> is the second most frequently isolated coagulase-negative staphylococci from healthy human skin and is known for its ability to produce antimicrobial peptides as well as thioalcohol from apocrine gland secretions, contributing to body odor (<xref ref-type="bibr" rid="B74">Severn et al., 2022</xref>; <xref ref-type="bibr" rid="B69">Rudden et al., 2020</xref>). Interestingly, old-age individuals have a characteristic body odor (<xref ref-type="bibr" rid="B50">Mitro et al., 2012</xref>), suggesting a link between elevated levels of this bacterial species in old age and dermal body odor. While the effect of a <italic>S. capitis</italic>-to-<italic>S. hominis</italic> shift on skin health is unknown, it suggests that <italic>S. hominis</italic> may opportunistically colonize the skin in response to decreased competition and altered ecological niches resulting from aging, potentially exerting influence on microbial homeostasis and the skin&#x2019;s local immune response. In support of this notion, <italic>S. hominis</italic> was previously found to be inversely correlated with sebum content and positively correlated with staphylococcal alpha diversity (<xref ref-type="bibr" rid="B2">Ahle et al., 2022</xref>), highlighting the potential importance of skin physiology and microbial community context in <italic>S. hominis</italic> expansion. <italic>S. hominis</italic> has also shown to have inhibitory properties against <italic>C. acnes</italic> (<xref ref-type="bibr" rid="B2">Ahle et al., 2022</xref>), suggesting interspecies competition may also play a role. Our study was limited to species-level interrogation of these bacteria, however, high strain heterogeneity and intra-species diversity has been well-documented for <italic>S. capitis</italic>, <italic>S. hominis</italic>, and <italic>C. acnes</italic> (<xref ref-type="bibr" rid="B29">Joglekar et al., 2023</xref>; <xref ref-type="bibr" rid="B13">Conwill et al., 2022</xref>), warranting future investigation into the role of strain diversity in the microbiome and aging.</p>
<p>Aging has been shown to impact not only the skin microbiome, but also the microbiome of other body sites such as the oral cavity and gastrointestinal tract (<xref ref-type="bibr" rid="B37">Larson et al., 2022</xref>; <xref ref-type="bibr" rid="B27">Huang et al., 2020</xref>; <xref ref-type="bibr" rid="B72">Sarafidou et al., 2024</xref>; <xref ref-type="bibr" rid="B23">Ghosh et al., 2022</xref>). As a part of this study, we also concurrently collected saliva samples for microbiome analysis from the same cohort (manuscript in progress). Characterizing the microbiota of these two distinct body sites will allow for even more in-depth insight into the impact of aging on the human microbiome, notably within the same population. Furthermore, this will provide the unique opportunity to integrate microbial data across body sites to investigate relationships between the skin and oral microbiota, particularly given the observation that oral bacteria have previously been associated with skin microbiome diversification in aging (<xref ref-type="bibr" rid="B75">Shibagaki et al., 2017</xref>).</p>
<p>While our study provided a detailed overview of how aging may modulate the bacterial members of the skin microbiome, it was limited by the use of 16S rRNA gene sequencing, which excludes non-bacterial microorganisms and does not provide strain-level resolution. Future efforts should employ metagenomic sequencing and additional techniques such as metatranscriptomics and metabolomics to gain a greater depth of understanding into the mechanisms of microbial function and strain specificity in skin aging. Furthermore, our study did not measure parameters of skin health such as moisture levels, elasticity, and the presence of visible signs of aging. Inclusion of these metrics may provide more defined relationships between microbiome shifts and aging. For example, relationships have been identified between the skin microbiota and clinical signs of aging such as grade of Crow&#x2019;s feet wrinkles (<xref ref-type="bibr" rid="B51">Myers et al., 2023</xref>; <xref ref-type="bibr" rid="B22">Garlet et al., 2024</xref>) and dark spot area (<xref ref-type="bibr" rid="B75">Shibagaki et al., 2017</xref>; <xref ref-type="bibr" rid="B40">Li et al., 2020</xref>), as well as skin characteristics including sebum content, sebaceous gland area, transepidermal water loss (TEWL), natural moisturizing factor (NMF), antimicrobial peptides (AMPs), and lipids (<xref ref-type="bibr" rid="B26">Howard et al., 2022</xref>; <xref ref-type="bibr" rid="B34">Kim et al., 2022</xref>; <xref ref-type="bibr" rid="B76">Sun et al., 2024</xref>; <xref ref-type="bibr" rid="B33">Kim et al., 2019</xref>). An additional limitation of the study is the potential influence of confounding variables that we did not control for, which include lifestyle factors (e.g., cosmetic use, timing of last wash) and health-related factors (e.g., underlying health conditions, medication). Lastly, our study was limited by collection of the microbiome at a single time point. Future efforts should employ longitudinal studies to better understand temporal dynamics of the microbiome during the aging process, which may help to elucidate how changes in the skin microbial communities correlate with the progression of skin aging.</p>
<p>In summary, these findings provide comprehensive insights into the microbiome of aging skin in a UK-based population, highlighting the shift from a core, stable community driven by <italic>C. acnes</italic> in young individuals, to a more variable and less robust community in older individuals. This work also supports and complements previous studies of the aging skin microbiome carried out on North America- and Asia-based populations, highlighting the connection between the skin microbiota and the physiology of skin aging across ethnicities and geography. By comparing the microbiome of distinct skin sites, our study also underscores the importance of considering skin site specificity when studying the impact of aging on skin microbiome dynamics, with both intrinsic and extrinsic factors of aging likely playing significant roles. We also demonstrate the potential of microbiome typing techniques to unravel complex microbial interactions and microbiome profiles associated with age, emphasizing the opportunity for advanced techniques such as machine learning in generating more personalized and precision interventions tailored to the needs of an individual&#x2019;s aging skin. Expanding these advanced methods to include other omics data, skin health parameters, and longitudinal studies has the potential to further our understanding of the intricate interplay between the skin microbiota, skin health, and the aging process. Ultimately, our findings highlight a significant opportunity for developing innovative and targeted next-generation anti-aging solutions, such as prebiotic or postbiotics treatments, to restore a more youthful, stable skin microbiome and thereby enhance overall skin health.</p>
</sec>
<sec sec-type="methods" id="s4">
<title>Methods</title>
<sec id="s4-1">
<title>Study design</title>
<p>Ethical approval for the study was granted by Aston University&#x2019;s College of Health and Life Sciences ethics committee (reference number HLS21008). Participants were recruited via advertisement through the Aston Research Centre for Health in Ageing (ARCHA), Aston University, and online social media channels. All participants gave informed written consent prior to sample collection. Initial participant screening was performed using a survey hosted on Qualtrics, during which basic demographic information was collected including age, gender, ethnicity, and postcode. Participants were required to speak and understand written English and were required to be aged between 20 and 40 or 60&#x2013;80&#xa0;years old. Participants who did not meet these criteria or who were pregnant were excluded. The specified age groups were selected to create distinct contrast between young and old populations and minimize the potential confounding effects of variables that might be present in a more heterogeneous middle-aged cohort, such as significant lifestyle changes and hormonal shifts (e.g., perimenopause). Additional metadata was collected at the time of sampling that included health, lifestyle, and geographical information, however this information was not used for study inclusion or exclusion. In total, 59 participants from Birmingham, United Kingdom were included in the study. Thirty participants were included in the young age (YA) group, and 29 participants were included in the old age (OA) group.</p>
</sec>
<sec id="s4-2">
<title>Sample collection</title>
<p>Sample collection and processing was performed at Aston University (Birmingham, United Kingdom). Participants continued their normal daily health and hygiene routine and did not follow any washout guidelines before sampling, reflecting a non-controlled, more representative microbiome state. A sampling solution was made using phosphate buffered saline modified with Tween<sup>&#xae;</sup> 20 to a final concentration of 0.1% (v/v). Sampling templates with a 5&#xa0;cm &#xd7; 5&#xa0;cm area were made from SILASTIC&#x2122; RTV-4136-M Liquid Silicone Rubber (Dow, Midland, Michigan). The sampling solution and templates were sterilized by autoclaving at 121 &#xb0;C for 15&#xa0;min prior to use.</p>
<p>Microbiome samples were collected from the right cheek and from the right antecubital fossa (interior of the inner elbow). For each area a sterile viscose swab (Scientific Laboratory Supplies, Nottingham, United Kingdom, catalog number: SWA3112) was soaked in 2&#xa0;mL of sampling solution and used to swab the participants skin within the sample area in a constant &#x201c;zigzag&#x201d; pattern. The swab head was then stored inside the 2&#xa0;mL sampling solution for storage at 2&#x2013;8 &#xb0;C prior to processing. Blank control swabs were collected by soaking the swab in sampling solution and processed in a manner identical to the skin swabs.</p>
</sec>
<sec id="s4-3">
<title>Preservation and processing of samples</title>
<p>Initial processing of samples occurred within 24&#xa0;h of sampling. DNA extraction was performed on the samples using Qiagen DNeasy PowerSoil Pro Kit (QIAGEN, Hilden, Germany, catalog number: 47014) following the manufacturer&#x2019;s protocol. DNA concentration was measured using a NanoDrop 1000 spectrophotometer (Thermo Fisher, United Kingdom), following storage of the DNA extracts at &#x2212;80 &#xb0;C. For transit prior to sequencing the DNA extracts were dried using an Eppendorf Concentrator plus at 30 &#xb0;C for 45&#xa0;min.</p>
</sec>
<sec id="s4-4">
<title>16S rRNA gene sequencing</title>
<p>Sequencing of the V1-V3 hypervariable region of the 16S rRNA gene was carried out by the Rutgers University Genomic Center. The V1-V3 region was amplified using NEB HiFidelity DNA polymerase (New England Biolabs) with the 28F (5&#x2032; GAGTTTGATCNTGGCTCAG) and 519R primers (5&#x2032; GTNTTACNGCGGCKGCTG). Barcoding was performed using Nextera XT dual index primers (Illumina) and NEB HiFidelity DNA polymerase (New England Biolabs). Libraries were sequenced on an Illumina MiSeq 2 x 300 flow cell.</p>
</sec>
<sec id="s4-5">
<title>Microbiome data analysis</title>
<sec id="s4-5-1">
<title>Sequence processing</title>
<p>Raw FASTQ sequencing data was imported into QIIME2 (version 2024.2.0) (<xref ref-type="bibr" rid="B4">Bolyen et al., 2019</xref>). Quality assessment identified primers and low-quality regions within the forward and reverse reads, therefore quality trimming was applied as follows: forward reads were trimmed 19&#xa0;bp and truncated at 291&#xa0;bp, and reverse reads were trimmed 18&#xa0;bp and truncated at 261&#xa0;bp, providing 22&#xa0;bp of read overlap between forward and reverse reads. The DADA2 plugin within the QIIME2 environment was used to generate an ASV feature table with 23,004 ASVs (<xref ref-type="bibr" rid="B8">Callahan et al., 2016</xref>). Four samples with read counts below 10,000 were removed following quality processing, resulting in 114 total microbiome samples with an initial median read count of 114,226 per sample for downstream filtering.</p>
</sec>
<sec id="s4-5-2">
<title>Taxonomic classification</title>
<p>To achieve improved species-level resolution with taxonomic assignment, a custom Naive Bayes classifier was trained on the extracted V1-V3 region of the Greengenes2 16S database (<xref ref-type="bibr" rid="B47">McDonald et al., 2024</xref>). Briefly, the V1-V3 region was extracted from the preformatted Greengenes2 backbone sequences (2022.10.backbone.full-length.fna.qza) using the 28F forward (GAGTTTGATCNTGGCTCAG) and 519R reverse (GTNTTACNGCGGCKGCTG) primers, with a 400&#xa0;bp minimum length and 550&#xa0;bp maximum length required for the extracted sequences (<italic>qiime feature-classifier extract-reads</italic>) (<xref ref-type="bibr" rid="B47">McDonald et al., 2024</xref>). Extracted sequences were then dereplicated using the &#x201c;uniq&#x201d; method from the RESCRIPT QIIME2 plugin (<italic>qiime rescript dereplicate</italic>). A Naive Bayes classifier was trained on this set of dereplicated V1-V3 sequences (<italic>qiime rescript evaluate-fit-classifier</italic>) (<xref ref-type="bibr" rid="B66">Robeson et al., 2021</xref>), followed by taxonomic assignment of ASVs using the classifier (<italic>qiime feature-classifier classify-sklearn</italic>).</p>
</sec>
<sec id="s4-5-3">
<title>Decontamination</title>
<p>To ensure a high-quality, robust dataset, decontamination of the samples was performed. Because the Greengenes2 2022.10 backbone taxonomy does not explicitly represent mitochondria and chloroplast, a prefiltering step was required to eliminate these contaminant ASVs. Here, ASVs were classified using a V1-V3 Naive Bayes classifier trained on the Silva 138 SSURef NR99 full-length database (<xref ref-type="bibr" rid="B63">Quast et al., 2013</xref>). This step identified 190 ASVs classified as chloroplast or mitochondria, which were then removed from the Greengenes2-classified ASV table. Subsequently, manual filtering was applied to remove ASVs unassigned at the domain level, ASVs assigned as Archaea or Eukaryota, and ASVs assigned as <italic>Pelomonas</italic> or <italic>Bradyrhizobium</italic>&#x2013;common contaminants, especially in low biomass samples (<xref ref-type="bibr" rid="B70">Salter et al., 2014</xref>; <xref ref-type="bibr" rid="B18">Eisenhofer et al., 2019</xref>) &#x2013; which collectively accounted for 3,964 ASVs removed. A final decontamination step using the combined frequency- and prevalence-based approach within decontam (v1.22.0) eliminated 71 more contaminant ASVs (<xref ref-type="bibr" rid="B15">Davis et al., 2018</xref>). In total, 4,035 unique ASVs were identified as contaminants and removed, resulting in 18,969 unique ASVs and a final median read count of 105,646 per sample for downstream analysis.</p>
</sec>
<sec id="s4-5-4">
<title>Diversity analysis</title>
<p>A phylogenetic tree was constructed from the decontaminated ASV sequences using MAFFT and FastTree within the QIIME2 environment (<italic>qiime phylogeny align-to-tree-mafft-fasttree)</italic> to use for phylogenetic-based diversity analysis (<xref ref-type="bibr" rid="B31">Katoh and Standley, 2013</xref>; <xref ref-type="bibr" rid="B62">Price et al., 2010</xref>). Using the R phyloseq package (v1.46.0), samples were then subjected to repeated rarefaction (n &#x3d; 100) without replacement to a depth of 13,074 counts (the minimum library size among all samples) for alpha diversity comparisons (<italic>rarefy_even_depth</italic>) (<xref ref-type="bibr" rid="B48">McMurdie and Holmes, 2013</xref>). For beta diversity comparisons, raw count data were normalized to relative abundances. Beta diversity indices were then computed using Bray-Curtis and Jaccard dissimilarities (via the vegan v2.6&#x2013;6.1&#xa0;R package (<xref ref-type="bibr" rid="B57">Oksanen et al., 2024</xref>)) and weighted and unweighted UniFrac distances (via phyloseq v1.46.0). The resulting dissimilarity matrices were ordinated using nonmetric multidimensional scaling (NMDS). Statistical differences were assessed using PERMANOVA (<italic>adonis2</italic> function from the vegan package). Finally, the top ASVs with abundances correlated to the NMDS coordinates were identified using Spearman correlation.</p>
</sec>
<sec id="s4-5-5">
<title>Taxonomic and differential abundance testing analysis</title>
<p>Differential abundance (DA) testing was performed on a filtered ASV feature table (raw counts as input), excluding ASVs present in fewer than 10% of samples or with fewer than 100 counts per sample. This filtering was performed to manage the burden of multiple-test correction and reduce tool-specific variation across multiple DA test methods (<xref ref-type="bibr" rid="B55">Nearing et al., 2022</xref>). This resulted in 4,823 ASVs as input. To ensure robust results, multiple methods were employed for DA testing, including Maaslin2 (v1.16.0, normalization &#x3d; &#x201c;TSS&#x201d;, transform &#x3d; &#x201c;AST&#x201d;) (<xref ref-type="bibr" rid="B55">Nearing et al., 2022</xref>; <xref ref-type="bibr" rid="B46">Mallick et al., 2021</xref>), DESeq2 (v1.42.1, default parameters) (<xref ref-type="bibr" rid="B45">Love et al., 2014</xref>), ANCOM-BC2 (v2.4.0, default parameters) (<xref ref-type="bibr" rid="B41">Lin and Peddada, 2024</xref>), and ALDEx2 (v1.34.0, default parameters) (<xref ref-type="bibr" rid="B20">Fernandes et al., 2013</xref>). ASVs identified as statistically significant (q &#x3c; 0.05) by at least two methods were considered differentially abundant between age groups. For additional taxonomic comparisons of taxa of interest, the ASV feature table was converted to relative abundances and collapsed to the specified taxonomic level (species, genus, phylum). Statistical comparisons between age groups were conducted using the non-parametric Wilcoxon rank sum test. Hierarchical clustering (method &#x3d; &#x201c;ward.D2&#x2033;) was performed on the Bray-Curtis dissimilarity matrix to explore the relationship between age and microbiome composition. Microbiome types based on microbiome composition at the ASV level were identified using Dirichlet multinomial mixtures (DMM), as implemented in the <italic>DirichletMultinomial</italic> R package (v1.44.0) (<xref ref-type="bibr" rid="B25">Holmes et al., 2012</xref>). The ASV feature table was first filtered to exclude ASVs present in &#x3c;30% of samples and those with &#x3c;100 counts. This stringent filter was applied to eliminate noise from rare taxa and focus on the core, defining members of the community. The lowest Laplace approximation was used to calculate model fit and determine the optimal number of microbiome types per skin site.</p>
</sec>
<sec id="s4-5-6">
<title>SPIEC-EASI microbial network inference</title>
<p>The statistical method SPIEC-EASI, implemented in the SpiecEasi R package (v1.0.7) (<xref ref-type="bibr" rid="B36">Kurtz et al., 2015</xref>), was used to identify associations between microbial species within each skin site and age group. Briefly, the ASV feature table was collapsed to the species level, retaining ASVs without species-level assignment to minimize data loss. Within each skin site, features present in more than 30% of samples were included in the network analysis, resulting in 194 features for the face (148, 42, and 4 assigned at the species, genus, or family level, respectively) and 144 features for the arm (133, 25, and 6 assigned at the species, genus, or family level, respectively). This filtering step was performed to reduce noise and computational burden, as including too many rare ASVs can yield uninterpretable networks. This prevalence-based approach is consistent with other studies employing SPIEC-EASI (<xref ref-type="bibr" rid="B36">Kurtz et al., 2015</xref>; <xref ref-type="bibr" rid="B80">Tipton et al., 2018</xref>; <xref ref-type="bibr" rid="B78">Swaney et al., 2022</xref>). SPIEC-EASI (method &#x3d; &#x201c;mb&#x201d;, lambda.min.ratio &#x3d; 1e-2, nlambda &#x3d; 20) was run on samples grouped by skin site and age group, generating two networks (younger and older age) for each skin site. Network visualization and calculation of topological properties were conducted using the R package igraph (v2.0.3) (<xref ref-type="bibr" rid="B14">Csardi and Nepusz, 2005</xref>).</p>
</sec>
<sec id="s4-5-7">
<title>Microbial function prediction</title>
<p>To predict microbial community function, the standalone version of PICRUSt2 (v2.5.2) (<xref ref-type="bibr" rid="B17">Douglas et al., 2020</xref>) was used to generate MetaCyc pathway abundances (<xref ref-type="bibr" rid="B9">Caspi et al., 2014</xref>). Predicted pathways were normalized by relative abundance, followed by calculation of Bray-Curtis dissimilarities and ordination using NMDS. Statistically significant differences were computed using PERMANOVA (<italic>vegan adonis2</italic>). Differential abundance testing of the pathways by age group and skin site was performed using multiple DA methods to ensure robustness of the identified pathways. The following methods were used from the ggpicrust2 R package (v1.7.3) using default parameters: DESeq2, edgeR (<xref ref-type="bibr" rid="B67">Robinson et al., 2010</xref>), Maaslin2, limma voom (<xref ref-type="bibr" rid="B65">Ritchie et al., 2015</xref>), metagenomeSeq (<xref ref-type="bibr" rid="B59">Paulson et al., 2013</xref>), and LinDA (<xref ref-type="bibr" rid="B86">Zhou et al., 2022</xref>). Of the methods that produced significant results, only the pathways that were identified by all methods were considered differentially abundant.</p>
</sec>
</sec>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: <ext-link ext-link-type="uri" xlink:href="https://www.ncbi.nlm.nih.gov/">https://www.ncbi.nlm.nih.gov/</ext-link>, PRJNA1273048.</p>
</sec>
<sec sec-type="ethics-statement" id="s6">
<title>Ethics statement</title>
<p>The studies involving humans were approved by Aston University&#x2019;s College of Health and Life Sciences ethics committee (reference number HLS21008). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.</p>
</sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>MS: Writing &#x2013; review and editing, Methodology, Data curation, Formal Analysis, Visualization, Writing &#x2013; original draft. DN: Writing &#x2013; original draft, Data curation, Writing &#x2013; review and editing, Investigation, Methodology. JM: Writing &#x2013; review and editing, Supervision. AH: Conceptualization, Supervision, Writing &#x2013; review and editing. TW: Supervision, Writing &#x2013; review and editing, Conceptualization. ML: Project administration, Supervision, Conceptualization, Writing &#x2013; review and editing.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research and/or publication of this article.</p>
</sec>
<ack>
<p>The authors would like to thank Fatima Saravia and Jiapeng Chen for their assistance with data curation, Lyndsay Schaeffer-Korbylo for contribution to study design and discussion of data analysis, and Kristen Mellenbruch for assistance in sample preparation, processing, and shipping. Additionally, we would like to thank Matthew Marples and John Simms for their work in designing and creating the mold used for the creation of the sampling templates.</p>
</ack>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>Authors MS, JM, and ML were employed by Colgate-Palmolive Company.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="ai-statement" id="s10">
<title>Generative AI statement</title>
<p>The author(s) declare that no Generative AI was used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec sec-type="supplementary-material" id="s12">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fragi.2025.1644012/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fragi.2025.1644012/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material>
<label>SUPPLEMENTARY FIGURE S1</label>
<caption>
<p>Alpha and beta diversity age group and skin site comparisons by gender. Alpha diversity metrics were calculated and compared between the old age (OA) group and young age (YA) group for the face and arm, split by gender. Calculations were performed on the ASV feature table and include: <bold>(A)</bold> Observed ASVs, <bold>(B)</bold> Pielou evenness, and <bold>(C)</bold> Shannon diversity. YA (blue) represents the young age group, and OA (red) represents the old age group. Points represent data outliers. Statistical significance for alpha diversity metrics was calculated using the Wilcoxon rank sum test (&#x2a; <italic>p</italic>-value &#x3c;0.05). NMDS ordination of the Bray-Curtis dissimilarity matrix for the face of <bold>(D)</bold> female and <bold>(E)</bold> male participants and for the arm of <bold>(F)</bold> female and <bold>(G)</bold> male participants. Statistical significance for compositional differences between age groups was calculated using PERMANOVA, with the <italic>p</italic>-value indicated within each ordination plot. Ellipses represent a 95% confidence interval for each age group. Blue circles represent the young age (YA) group and red triangles represent the old age (OA) group.</p>
</caption>
</supplementary-material>
<supplementary-material>
<label>SUPPLEMENTARY FIGURE S2</label>
<caption>
<p>Weighted UniFrac and Bray-Curtis NMDS ordination of face and arm samples. NMDS ordination of the weighted UniFrac distance matrix for the <bold>(A)</bold> face and <bold>(B)</bold> arm. <bold>(C)</bold> NMDS ordination of the Bray-Curtis dissimilarity matrix for all samples. Statistical significance for composition differences was calculated using PERMANOVA, with the <italic>p</italic>-value indicated within each ordination plot. Ellipses represent a 95% confidence interval for each age group. One old age sample from the face was excluded from the weighted UniFrac calculation due to a <italic>Pseudomonas</italic>-dominant microbial community profile that highly skewed the ordination.</p>
</caption>
</supplementary-material>
<supplementary-material>
<label>SUPPLEMENTARY FIGURE S3</label>
<caption>
<p>Relative abundances of the top four phyla present by skin site and age group. Phyla not included within the top four phyla are grouped into &#x201c;Other&#x201d;. Statistical significance was calculated using the Wilcoxon rank sum test (&#x2a; <italic>p</italic>-value &#x3c;0.05). YA &#x3d; young age group, OA &#x3d; old age group.</p>
</caption>
</supplementary-material>
<supplementary-material>
<label>SUPPLEMENTARY FIGURE S4</label>
<caption>
<p>Genus- and species-level variation of abundant skin taxa between gender and age groups. <bold>(A)</bold> Relative abundance differences between <italic>Cutibacterium, Staphylococcus</italic>, and <italic>Corynebacterium</italic> within each age group, skin site, and gender. Relative abundance differences for species of interest for the <bold>(B)</bold> face and <bold>(C)</bold> arm, split by gender and age group. Statistical significance was calculated using the Wilcoxon rank sum test &#x2a; <italic>p</italic>-value &#x3c;0.05, &#x2a;&#x2a; <italic>p</italic>-value &#x3c;0.01. YA &#x3d; young age group, OA &#x3d; old age group.</p>
</caption>
</supplementary-material>
<supplementary-material>
<label>SUPPLEMENTARY FIGURE S5</label>
<caption>
<p>Differentially abundant taxa between skin sites and age groups. Shown are all differentially abundant ASVs identified for the <bold>(A)</bold> face and <bold>(B)</bold> arm across age groups that were identified by at least two differential abundance testing methods. The horizontal line indicates the mean relative abundance value of each ASV. Multiple ASVs assigned to the same taxon are distinguished with sequential numbering starting from 1. YA &#x3d; young age group, OA &#x3d; old age group.</p>
</caption>
</supplementary-material>
<supplementary-material>
<label>SUPPLEMENTARY FIGURE S6</label>
<caption>
<p>Top driving ASVs of DMM community types. ASVs identified as the top drivers of Dirichlet multinomial mixtures (DMM) community types are indicated for <bold>(A)</bold> face community type 1, <bold>(B)</bold> face community type 2, <bold>(C)</bold> arm community type 1, <bold>(D)</bold> arm community type 2, and <bold>(E)</bold> arm community type 3. The top drivers are defined as the top 5% of ASVs that contribute the highest absolute weight within each model of community type.</p>
</caption>
</supplementary-material>
<supplementary-material>
<label>SUPPLEMENTARY FIGURE S7</label>
<caption>
<p>Age group distribution of Dirichlet multinomial mixtures (DMM) cluster profiles by gender. Number of samples that are assigned to each community type in <xref ref-type="fig" rid="F5">Figure 5</xref>, colored by age group and split by gender, for the <bold>(A)</bold> face and <bold>(B)</bold> arm. YA (blue) &#x3d; young age, OA (red) &#x3d; old age.</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Image3.jpeg" id="SM1" mimetype="application/jpeg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table1.docx" id="SM2" mimetype="application/docx" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image1.jpeg" id="SM3" mimetype="application/jpeg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image4.jpeg" id="SM4" mimetype="application/jpeg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image7.jpeg" id="SM5" mimetype="application/jpeg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image2.jpeg" id="SM6" mimetype="application/jpeg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image5.jpeg" id="SM7" mimetype="application/jpeg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image6.jpeg" id="SM8" mimetype="application/jpeg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="DataSheet1.xlsx" id="SM9" mimetype="application/xlsx" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Agak</surname>
<given-names>G. W.</given-names>
</name>
<name>
<surname>Kao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ouyang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Moon</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Butt</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Phenotype and antimicrobial activity of Th17 cells induced by Propionibacterium acnes strains associated with healthy and acne skin</article-title>. <source>J. Invest Dermatol</source> <volume>138</volume> (<issue>2</issue>), <fpage>316</fpage>&#x2013;<lpage>324</lpage>. <pub-id pub-id-type="doi">10.1016/j.jid.2017.07.842</pub-id>
<pub-id pub-id-type="pmid">28864077</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ahle</surname>
<given-names>C. M.</given-names>
</name>
<name>
<surname>St&#xf8;dkilde</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Poehlein</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>B&#xf6;meke</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Streit</surname>
<given-names>W. R.</given-names>
</name>
<name>
<surname>Wenck</surname>
<given-names>H.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Interference and co-existence of staphylococci and Cutibacterium acnes within the healthy human skin microbiome</article-title>. <source>Commun. Biol.</source> <volume>5</volume> (<issue>1</issue>), <fpage>923</fpage>. <pub-id pub-id-type="doi">10.1038/s42003-022-03897-6</pub-id>
<pub-id pub-id-type="pmid">36071129</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Almoughrabie</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Cau</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Cavagnero</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>O&#x2019;Neill</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Roso-Mares</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Commensal Cutibacterium acnes induce epidermal lipid synthesis important for skin barrier function</article-title>. <source>Sci. Adv.</source> <volume>9</volume>, <fpage>eadg6262</fpage>. <pub-id pub-id-type="doi">10.1126/sciadv.adg6262</pub-id>
<pub-id pub-id-type="pmid">37595033</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bolyen</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Rideout</surname>
<given-names>J. R.</given-names>
</name>
<name>
<surname>Dillon</surname>
<given-names>M. R.</given-names>
</name>
<name>
<surname>Bokulich</surname>
<given-names>N. A.</given-names>
</name>
<name>
<surname>Abnet</surname>
<given-names>C. C.</given-names>
</name>
<name>
<surname>Al-Ghalith</surname>
<given-names>G. A.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2</article-title>. <source>Nat. Biotechnol.</source> <volume>37</volume> (<issue>8</issue>), <fpage>852</fpage>&#x2013;<lpage>857</lpage>. <pub-id pub-id-type="doi">10.1038/s41587-019-0209-9</pub-id>
<pub-id pub-id-type="pmid">31341288</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bulliard</surname>
<given-names>J. L.</given-names>
</name>
<name>
<surname>De Weck</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Fisch</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Bordoni</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Levi</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Detailed site distribution of melanoma and sunlight exposure: aetiological patterns from a Swiss series</article-title>. <source>Ann. Oncol.</source> <volume>18</volume> (<issue>4</issue>), <fpage>789</fpage>&#x2013;<lpage>794</lpage>. <pub-id pub-id-type="doi">10.1093/annonc/mdl490</pub-id>
<pub-id pub-id-type="pmid">17237475</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Burns</surname>
<given-names>E. M.</given-names>
</name>
<name>
<surname>Ahmed</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Isedeh</surname>
<given-names>P. N.</given-names>
</name>
<name>
<surname>Kohli</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Van Der Pol</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Shaheen</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Ultraviolet radiation, both UVA and UVB, influences the composition of the skin microbiome</article-title>. <source>Exp. Dermatol.</source> <volume>28</volume> (<issue>2</issue>), <fpage>136</fpage>&#x2013;<lpage>141</lpage>. <pub-id pub-id-type="doi">10.1111/exd.13854</pub-id>
<pub-id pub-id-type="pmid">30506967</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Byrd</surname>
<given-names>A. L.</given-names>
</name>
<name>
<surname>Belkaid</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Segre</surname>
<given-names>J. A.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The human skin microbiome</article-title>. <source>Nat. Rev. Microbiol.</source> <volume>16</volume> (<issue>3</issue>), <fpage>143</fpage>&#x2013;<lpage>155</lpage>. <pub-id pub-id-type="doi">10.1038/nrmicro.2017.157</pub-id>
<pub-id pub-id-type="pmid">29332945</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Callahan</surname>
<given-names>B. J.</given-names>
</name>
<name>
<surname>McMurdie</surname>
<given-names>P. J.</given-names>
</name>
<name>
<surname>Rosen</surname>
<given-names>M. J.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>A. W.</given-names>
</name>
<name>
<surname>Johnson</surname>
<given-names>A. J. A.</given-names>
</name>
<name>
<surname>Holmes</surname>
<given-names>S. P.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>DADA2: high-resolution sample inference from Illumina amplicon data</article-title>. <source>Nat. Methods</source> <volume>13</volume> (<issue>7</issue>), <fpage>581</fpage>&#x2013;<lpage>583</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.3869</pub-id>
<pub-id pub-id-type="pmid">27214047</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Caspi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Altman</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Billington</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Dreher</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Foerster</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Fulcher</surname>
<given-names>C. A.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases</article-title>. <source>Nucleic Acids Res.</source> <volume>42</volume> (<issue>Database issue</issue>), <fpage>D459</fpage>&#x2013;<lpage>D471</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gkt1103</pub-id>
<pub-id pub-id-type="pmid">24225315</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Choi</surname>
<given-names>E. H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Aging of the skin barrier</article-title>. <source>Clin. Dermatology</source> <volume>37</volume> (<issue>4</issue>), <fpage>336</fpage>&#x2013;<lpage>345</lpage>. <pub-id pub-id-type="doi">10.1016/j.clindermatol.2019.04.009</pub-id>
<pub-id pub-id-type="pmid">31345321</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chong</surname>
<given-names>C. E.</given-names>
</name>
<name>
<surname>Bengtsson</surname>
<given-names>R. J.</given-names>
</name>
<name>
<surname>Horsburgh</surname>
<given-names>M. J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Comparative genomics of Staphylococcus capitis reveals species determinants</article-title>. <source>Front. Microbiol.</source> <volume>13</volume>, <fpage>1005949</fpage>. <pub-id pub-id-type="doi">10.3389/fmicb.2022.1005949</pub-id>
<pub-id pub-id-type="pmid">36246238</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Claesen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Spagnolo</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Ramos</surname>
<given-names>S. F.</given-names>
</name>
<name>
<surname>Kurita</surname>
<given-names>K. L.</given-names>
</name>
<name>
<surname>Byrd</surname>
<given-names>A. L.</given-names>
</name>
<name>
<surname>Aksenov</surname>
<given-names>A. A.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>A Cutibacterium acnes antibiotic modulates human skin microbiota composition in hair follicles</article-title>. <source>Sci. Transl. Med.</source> <volume>12</volume>, <fpage>eaay5445</fpage>. <pub-id pub-id-type="doi">10.1126/scitranslmed.aay5445</pub-id>
<pub-id pub-id-type="pmid">33208503</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Conwill</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Kuan</surname>
<given-names>A. C.</given-names>
</name>
<name>
<surname>Damerla</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Poret</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Baker</surname>
<given-names>J. S.</given-names>
</name>
<name>
<surname>Tripp</surname>
<given-names>A. D.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Anatomy promotes neutral coexistence of strains in the human skin microbiome</article-title>. <source>Cell Host Microbe</source> <volume>30</volume> (<issue>2</issue>), <fpage>171</fpage>&#x2013;<lpage>82.e7</lpage>. <pub-id pub-id-type="doi">10.1016/j.chom.2021.12.007</pub-id>
<pub-id pub-id-type="pmid">34995483</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Csardi</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Nepusz</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>The igraph software package for complex network research. InterJournal</article-title>. <source>Complex Syst.</source> 1695, 1&#x2010;9. Available online at: <ext-link ext-link-type="uri" xlink:href="http://igraph.org.">http://igraph.org.</ext-link>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Davis</surname>
<given-names>N. M.</given-names>
</name>
<name>
<surname>Proctor</surname>
<given-names>D. M.</given-names>
</name>
<name>
<surname>Holmes</surname>
<given-names>S. P.</given-names>
</name>
<name>
<surname>Relman</surname>
<given-names>D. A.</given-names>
</name>
<name>
<surname>Callahan</surname>
<given-names>B. J.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data</article-title>. <source>Microbiome</source> <volume>6</volume> (<issue>1</issue>), <fpage>226</fpage>. <pub-id pub-id-type="doi">10.1186/s40168-018-0605-2</pub-id>
<pub-id pub-id-type="pmid">30558668</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dimitriu</surname>
<given-names>P. A.</given-names>
</name>
<name>
<surname>Iker</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Malik</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Leung</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Mohn</surname>
<given-names>W. W.</given-names>
</name>
<name>
<surname>Hillebrand</surname>
<given-names>G. G.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>New insights into the intrinsic and extrinsic factors that shape the human skin microbiome</article-title>. <source>MBio</source> <volume>10</volume> (<issue>4</issue>), <fpage>e00839-19</fpage>. <pub-id pub-id-type="doi">10.1128/mBio.00839-19</pub-id>
<pub-id pub-id-type="pmid">31266865</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Douglas</surname>
<given-names>G. M.</given-names>
</name>
<name>
<surname>Maffei</surname>
<given-names>V. J.</given-names>
</name>
<name>
<surname>Zaneveld</surname>
<given-names>J. R.</given-names>
</name>
<name>
<surname>Yurgel</surname>
<given-names>S. N.</given-names>
</name>
<name>
<surname>Brown</surname>
<given-names>J. R.</given-names>
</name>
<name>
<surname>Taylor</surname>
<given-names>C. M.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>PICRUSt2 for prediction of metagenome functions</article-title>. <source>Nat. Biotechnol.</source> <volume>38</volume> (<issue>6</issue>), <fpage>685</fpage>&#x2013;<lpage>688</lpage>. <pub-id pub-id-type="doi">10.1038/s41587-020-0548-6</pub-id>
<pub-id pub-id-type="pmid">32483366</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Eisenhofer</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Minich</surname>
<given-names>J. J.</given-names>
</name>
<name>
<surname>Marotz</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Cooper</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Knight</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Weyrich</surname>
<given-names>L. S.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Contamination in low microbial biomass microbiome studies: issues and recommendations</article-title>. <source>Trends Microbiol.</source> <volume>27</volume> (<issue>2</issue>), <fpage>105</fpage>&#x2013;<lpage>117</lpage>. <pub-id pub-id-type="doi">10.1016/j.tim.2018.11.003</pub-id>
<pub-id pub-id-type="pmid">30497919</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Farage</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Miller</surname>
<given-names>K. W.</given-names>
</name>
<name>
<surname>Elsner</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Maibach</surname>
<given-names>H. I.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Intrinsic and extrinsic factors in skin ageing: a review</article-title>. <source>Int. J. Cosmet. Sci.</source> <volume>30</volume> (<issue>2</issue>), <fpage>87</fpage>&#x2013;<lpage>95</lpage>. <pub-id pub-id-type="doi">10.1111/j.1468-2494.2007.00415.x</pub-id>
<pub-id pub-id-type="pmid">18377617</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fernandes</surname>
<given-names>A. D.</given-names>
</name>
<name>
<surname>Macklaim</surname>
<given-names>J. M.</given-names>
</name>
<name>
<surname>Linn</surname>
<given-names>T. G.</given-names>
</name>
<name>
<surname>Reid</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Gloor</surname>
<given-names>G. B.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>ANOVA-like differential expression (ALDEx) analysis for mixed population RNA-Seq</article-title>. <source>PLoS One</source> <volume>8</volume> (<issue>7</issue>), <fpage>e67019</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0067019</pub-id>
<pub-id pub-id-type="pmid">23843979</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fuchs</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Boll</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Heider</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Microbial degradation of aromatic compounds - from one strategy to four</article-title>. <source>Nat. Rev. Microbiol.</source> <volume>9</volume> (<issue>11</issue>), <fpage>803</fpage>&#x2013;<lpage>816</lpage>. <pub-id pub-id-type="doi">10.1038/nrmicro2652</pub-id>
<pub-id pub-id-type="pmid">21963803</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garlet</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Andre-Frei</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Del Bene</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Cameron</surname>
<given-names>H. J.</given-names>
</name>
<name>
<surname>Samuga</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Rawat</surname>
<given-names>V.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Facial skin microbiome composition and functional shift with aging</article-title>. <source>Microorganisms</source> <volume>12</volume> (<issue>5</issue>), <fpage>1021</fpage>. <pub-id pub-id-type="doi">10.3390/microorganisms12051021</pub-id>
<pub-id pub-id-type="pmid">38792850</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ghosh</surname>
<given-names>T. S.</given-names>
</name>
<name>
<surname>Shanahan</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>O&#x2019;Toole</surname>
<given-names>P. W.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>The gut microbiome as a modulator of healthy ageing</article-title>. <source>Nat. Rev. Gastroenterol. Hepatol.</source> <volume>19</volume> (<issue>9</issue>), <fpage>565</fpage>&#x2013;<lpage>584</lpage>. <pub-id pub-id-type="doi">10.1038/s41575-022-00605-x</pub-id>
<pub-id pub-id-type="pmid">35468952</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Dou</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Aging and aging-related diseases: from molecular mechanisms to interventions and treatments</article-title>. <source>Signal Transduct. Target. Ther.</source> <volume>7</volume> (<issue>1</issue>), <fpage>391</fpage>&#x2013;<lpage>40</lpage>. <pub-id pub-id-type="doi">10.1038/s41392-022-01251-0</pub-id>
<pub-id pub-id-type="pmid">36522308</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Holmes</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Harris</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Quince</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Dirichlet multinomial mixtures: generative models for microbial metagenomics</article-title>. <source>PLoS One</source> <volume>7</volume> (<issue>2</issue>), <fpage>e30126</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0030126</pub-id>
<pub-id pub-id-type="pmid">22319561</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Howard</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Bascom</surname>
<given-names>C. C.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Binder</surname>
<given-names>R. L.</given-names>
</name>
<name>
<surname>Fadayel</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Huggins</surname>
<given-names>T. G.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Aging-associated changes in the adult human skin microbiome and the host factors that affect skin microbiome composition</article-title>. <source>J. Invest Dermatol.</source> <volume>142</volume> (<issue>7</issue>), <fpage>1934</fpage>&#x2013;<lpage>1946.e21</lpage>. <pub-id pub-id-type="doi">10.1016/j.jid.2021.11.029</pub-id>
<pub-id pub-id-type="pmid">34890626</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Haiminen</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Carrieri</surname>
<given-names>A. P.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Parida</surname>
<given-names>L.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Human skin, oral, and gut microbiomes predict chronological age</article-title>. <source>mSystems</source> <volume>5</volume> (<issue>1</issue>), <fpage>e00630-19</fpage>. <pub-id pub-id-type="doi">10.1128/msystems.00630-19</pub-id>
<pub-id pub-id-type="pmid">32047061</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhuo</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>The updates and implications of cutaneous microbiota in acne</article-title>. <source>Cell and Biosci.</source> <volume>13</volume> (<issue>1</issue>), <fpage>113</fpage>&#x2013;<lpage>118</lpage>. <pub-id pub-id-type="doi">10.1186/s13578-023-01072-w</pub-id>
<pub-id pub-id-type="pmid">37344849</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Joglekar</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Conlan</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lee-Lin</surname>
<given-names>S. Q.</given-names>
</name>
<name>
<surname>Deming</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Kashaf</surname>
<given-names>S. S.</given-names>
</name>
</person-group>
<collab>NISC Comparative Sequencing Program</collab>, et al. (<year>2023</year>). <article-title>Integrated genomic and functional analyses of human skin-associated Staphylococcus reveal extensive inter- and intra-species diversity</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>120</volume> (<issue>47</issue>), <fpage>e2310585120</fpage>. <pub-id pub-id-type="doi">10.1073/pnas.2310585120</pub-id>
<pub-id pub-id-type="pmid">37956283</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jug&#xe9;</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Rouaud-Tinguely</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Breugnot</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Servaes</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Grimaldi</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Roth</surname>
<given-names>M. P.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Shift in skin microbiota of Western European women across aging</article-title>. <source>J. Appl. Microbiol.</source> <volume>125</volume> (<issue>3</issue>), <fpage>907</fpage>&#x2013;<lpage>916</lpage>. <pub-id pub-id-type="doi">10.1111/jam.13929</pub-id>
<pub-id pub-id-type="pmid">29791788</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Katoh</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Standley</surname>
<given-names>D. M.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>MAFFT multiple sequence alignment software version 7: improvements in performance and usability</article-title>. <source>Mol. Biol. Evol.</source> <volume>30</volume> (<issue>4</issue>), <fpage>772</fpage>&#x2013;<lpage>780</lpage>. <pub-id pub-id-type="doi">10.1093/molbev/mst010</pub-id>
<pub-id pub-id-type="pmid">23329690</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ke</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Unveiling the regulatory mechanisms of salicylate degradation gene cluster cehGHIR4 in Rhizobium sp. strain X9</article-title>. <source>Appl. Environ. Microbiol.</source> <volume>89</volume> (<issue>10</issue>), <fpage>e0080223</fpage>. <pub-id pub-id-type="doi">10.1128/aem.00802-23</pub-id>
<pub-id pub-id-type="pmid">37800922</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>H. J.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>J. J.</given-names>
</name>
<name>
<surname>Myeong</surname>
<given-names>N. R.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Segregation of age-related skin microbiome characteristics by functionality</article-title>. <source>Sci. Rep.</source> <volume>9</volume> (<issue>1</issue>), <fpage>16748</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-019-53266-3</pub-id>
<pub-id pub-id-type="pmid">31727980</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>H. J.</given-names>
</name>
<name>
<surname>Oh</surname>
<given-names>H. N.</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>H. G.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Aged related human skin microbiome and mycobiome in Korean women</article-title>. <source>Sci. Rep.</source> <volume>12</volume> (<issue>1</issue>), <fpage>2351</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-022-06189-5</pub-id>
<pub-id pub-id-type="pmid">35149745</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Krutmann</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Bouloc</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Sore</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Bernard</surname>
<given-names>B. A.</given-names>
</name>
<name>
<surname>Passeron</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>The skin aging exposome</article-title>. <source>J. Dermatological Sci.</source> <volume>85</volume> (<issue>3</issue>), <fpage>152</fpage>&#x2013;<lpage>161</lpage>. <pub-id pub-id-type="doi">10.1016/j.jdermsci.2016.09.015</pub-id>
<pub-id pub-id-type="pmid">27720464</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kurtz</surname>
<given-names>Z. D.</given-names>
</name>
<name>
<surname>M&#xfc;ller</surname>
<given-names>C. L.</given-names>
</name>
<name>
<surname>Miraldi</surname>
<given-names>E. R.</given-names>
</name>
<name>
<surname>Littman</surname>
<given-names>D. R.</given-names>
</name>
<name>
<surname>Blaser</surname>
<given-names>M. J.</given-names>
</name>
<name>
<surname>Bonneau</surname>
<given-names>R. A.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Sparse and compositionally robust inference of microbial ecological networks</article-title>. <source>PLoS Comput. Biol.</source> <volume>11</volume> (<issue>5</issue>), <fpage>e1004226</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pcbi.1004226</pub-id>
<pub-id pub-id-type="pmid">25950956</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Larson</surname>
<given-names>P. J.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Santiago</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Driscoll</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Fleming</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Voigt</surname>
<given-names>A. Y.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Associations of the skin, oral and gut microbiome with aging, frailty and infection risk reservoirs in older adults</article-title>. <source>Nat. Aging</source> <volume>2</volume> (<issue>10</issue>), <fpage>941</fpage>&#x2013;<lpage>955</lpage>. <pub-id pub-id-type="doi">10.1038/s43587-022-00287-9</pub-id>
<pub-id pub-id-type="pmid">36398033</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Leung</surname>
<given-names>M. H. Y.</given-names>
</name>
<name>
<surname>Tong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Bastien</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Guinot</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Tenenhaus</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Appenzeller</surname>
<given-names>B. M. R.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Changes of the human skin microbiota upon chronic exposure to polycyclic aromatic hydrocarbon pollutants</article-title>. <source>Microbiome</source> <volume>8</volume> (<issue>1</issue>), <fpage>100</fpage>&#x2013;<lpage>117</lpage>. <pub-id pub-id-type="doi">10.1186/s40168-020-00874-1</pub-id>
<pub-id pub-id-type="pmid">32591010</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Leung</surname>
<given-names>M. H. Y.</given-names>
</name>
<name>
<surname>Tong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Bastien</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Appenzeller</surname>
<given-names>B. M. R.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Skin microbiome differentiates into distinct cutotypes with unique metabolic functions upon exposure to polycyclic aromatic hydrocarbons</article-title>. <source>Microbiome</source> <volume>11</volume> (<issue>1</issue>), <fpage>124</fpage>&#x2013;<lpage>14</lpage>. <pub-id pub-id-type="doi">10.1186/s40168-023-01564-4</pub-id>
<pub-id pub-id-type="pmid">37264459</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Bai</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yi</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>New insights into the skin microbial communities and skin aging</article-title>. <source>Front. Microbiol.</source> <volume>11</volume>, <fpage>565549</fpage>. <pub-id pub-id-type="doi">10.3389/fmicb.2020.565549</pub-id>
<pub-id pub-id-type="pmid">33193154</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Peddada</surname>
<given-names>S. D.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Multigroup analysis of compositions of microbiomes with covariate adjustments and repeated measures</article-title>. <source>Nat. Methods</source> <volume>21</volume> (<issue>1</issue>), <fpage>83</fpage>&#x2013;<lpage>91</lpage>. <pub-id pub-id-type="doi">10.1038/s41592-023-02092-7</pub-id>
<pub-id pub-id-type="pmid">38158428</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Ranallo</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Rios</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Grice</surname>
<given-names>E. A.</given-names>
</name>
<name>
<surname>Moon</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Gallo</surname>
<given-names>R. L.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Crosstalk between skin microbiota and immune system in health and disease</article-title>. <source>Nat. Immunol.</source> <volume>24</volume> (<issue>6</issue>), <fpage>895</fpage>&#x2013;<lpage>898</lpage>. <pub-id pub-id-type="doi">10.1038/s41590-023-01500-6</pub-id>
<pub-id pub-id-type="pmid">37081151</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>L&#xf3;pez-Ot&#xed;n</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Blasco</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Partridge</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Serrano</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kroemer</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>The hallmarks of aging</article-title>. <source>Cell</source> <volume>153</volume> (<issue>6</issue>), <fpage>1194</fpage>&#x2013;<lpage>1217</lpage>. <pub-id pub-id-type="doi">10.1016/j.cell.2013.05.039</pub-id>
<pub-id pub-id-type="pmid">23746838</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>L&#xf3;pez-Ot&#xed;n</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Blasco</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Partridge</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Serrano</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kroemer</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Hallmarks of aging: an expanding universe</article-title>. <source>Cell</source> <volume>186</volume> (<issue>2</issue>), <fpage>243</fpage>&#x2013;<lpage>278</lpage>. <pub-id pub-id-type="doi">10.1016/j.cell.2022.11.001</pub-id>
<pub-id pub-id-type="pmid">36599349</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Love</surname>
<given-names>M. I.</given-names>
</name>
<name>
<surname>Huber</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Anders</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2</article-title>. <source>Genome Biol.</source> <volume>15</volume> (<issue>12</issue>), <fpage>550</fpage>. <pub-id pub-id-type="doi">10.1186/s13059-014-0550-8</pub-id>
<pub-id pub-id-type="pmid">25516281</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mallick</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Rahnavard</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>McIver</surname>
<given-names>L. J.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>L. H.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Multivariable association discovery in population-scale meta-omics studies</article-title>. <source>PLoS Comput. Biol.</source> <volume>17</volume> (<issue>11</issue>), <fpage>e1009442</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pcbi.1009442</pub-id>
<pub-id pub-id-type="pmid">34784344</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McDonald</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Balaban</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Cantrell</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Gonzalez</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Greengenes2 unifies microbial data in a single reference tree</article-title>. <source>Nat. Biotechnol.</source> <volume>42</volume> (<issue>5</issue>), <fpage>715</fpage>&#x2013;<lpage>718</lpage>. <pub-id pub-id-type="doi">10.1038/s41587-023-01845-1</pub-id>
<pub-id pub-id-type="pmid">37500913</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McMurdie</surname>
<given-names>P. J.</given-names>
</name>
<name>
<surname>Holmes</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data</article-title>. <source>PLoS One</source> <volume>8</volume> (<issue>4</issue>), <fpage>e61217</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0061217</pub-id>
<pub-id pub-id-type="pmid">23630581</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Misra</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Clavaud</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Guinot</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Bourokba</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Nouveau</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Mezzache</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Multi-omics analysis to decipher the molecular link between chronic exposure to pollution and human skin dysfunction</article-title>. <source>Sci. Rep.</source> <volume>11</volume> (<issue>1</issue>), <fpage>18302</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-021-97572-1</pub-id>
<pub-id pub-id-type="pmid">34526566</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mitro</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Gordon</surname>
<given-names>A. R.</given-names>
</name>
<name>
<surname>Olsson</surname>
<given-names>M. J.</given-names>
</name>
<name>
<surname>Lundstr&#xf6;m</surname>
<given-names>J. N.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>The smell of age: perception and discrimination of body odors of different ages</article-title>. <source>PLoS One</source> <volume>7</volume> (<issue>5</issue>), <fpage>e38110</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0038110</pub-id>
<pub-id pub-id-type="pmid">22666457</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Myers</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Bouslimani</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Hansen</surname>
<given-names>S. T.</given-names>
</name>
<name>
<surname>Clavaud</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Azouaoui</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>A multi-study analysis enables identification of potential microbial features associated with skin aging signs</article-title>. <source>Front. Aging.</source> <volume>4</volume>, <fpage>1304705</fpage>. <pub-id pub-id-type="doi">10.3389/fragi.2023.1304705</pub-id>
<pub-id pub-id-type="pmid">38362046</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nagy</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Pivarcsi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Kis</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Koreck</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Bodai</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>McDowell</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2006</year>). <article-title>Propionibacterium acnes and lipopolysaccharide induce the expression of antimicrobial peptides and proinflammatory cytokines/chemokines in human sebocytes</article-title>. <source>Microbes Infect.</source> <volume>8</volume> (<issue>8</issue>), <fpage>2195</fpage>&#x2013;<lpage>2205</lpage>. <pub-id pub-id-type="doi">10.1016/j.micinf.2006.04.001</pub-id>
<pub-id pub-id-type="pmid">16797202</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nakamura</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>O&#x2019;Neill</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Williams</surname>
<given-names>M. R.</given-names>
</name>
<name>
<surname>Cau</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Nakatsuji</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Horswill</surname>
<given-names>A. R.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Short chain fatty acids produced by Cutibacterium acnes inhibit biofilm formation by Staphylococcus epidermidis</article-title>. <source>Sci. Rep.</source> <volume>10</volume> (<issue>1</issue>), <fpage>21237</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-020-77790-9</pub-id>
<pub-id pub-id-type="pmid">33277548</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nakatsuji</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Gallo</surname>
<given-names>R. L.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The role of the skin microbiome in atopic dermatitis</article-title>. <source>Ann. allergy, asthma and Immunol. official Publ. Am. Coll. Allergy, Asthma, and Immunol.</source> <volume>122</volume> (<issue>3</issue>), <fpage>263</fpage>&#x2013;<lpage>269</lpage>. <pub-id pub-id-type="doi">10.1016/j.anai.2018.12.003</pub-id>
<pub-id pub-id-type="pmid">30550810</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nearing</surname>
<given-names>J. T.</given-names>
</name>
<name>
<surname>Douglas</surname>
<given-names>G. M.</given-names>
</name>
<name>
<surname>Hayes</surname>
<given-names>M. G.</given-names>
</name>
<name>
<surname>MacDonald</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Desai</surname>
<given-names>D. K.</given-names>
</name>
<name>
<surname>Allward</surname>
<given-names>N.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Microbiome differential abundance methods produce different results across 38 datasets</article-title>. <source>Nat. Commun.</source> <volume>13</volume> (<issue>1</issue>), <fpage>342</fpage>. <pub-id pub-id-type="doi">10.1038/s41467-022-28034-z</pub-id>
<pub-id pub-id-type="pmid">35039521</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oh</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Byrd</surname>
<given-names>A. L.</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kong</surname>
<given-names>H. H.</given-names>
</name>
<name>
<surname>Segre</surname>
<given-names>J. A.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Temporal stability of the human skin microbiome</article-title>. <source>Cell</source> <volume>165</volume> (<issue>4</issue>), <fpage>854</fpage>&#x2013;<lpage>866</lpage>. <pub-id pub-id-type="doi">10.1016/j.cell.2016.04.008</pub-id>
<pub-id pub-id-type="pmid">27153496</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oksanen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Simpson</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Blanchet</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Kindt</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Legendre</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Minchin</surname>
<given-names>P.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Vegan: community ecology package</article-title>.</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Parales</surname>
<given-names>R. E.</given-names>
</name>
<name>
<surname>Parales</surname>
<given-names>J. V.</given-names>
</name>
<name>
<surname>Pelletier</surname>
<given-names>D. A.</given-names>
</name>
<name>
<surname>Ditty</surname>
<given-names>J. L.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Diversity of microbial toluene degradation pathways</article-title>. <source>Adv. Appl. Microbiol.</source> <volume>64</volume> (<issue>1&#x2013;73</issue>), <fpage>1</fpage>&#x2013;<lpage>73</lpage>. <pub-id pub-id-type="doi">10.1016/S0065-2164(08)00401-2</pub-id>
<pub-id pub-id-type="pmid">18485280</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Paulson</surname>
<given-names>J. N.</given-names>
</name>
<name>
<surname>Stine</surname>
<given-names>O. C.</given-names>
</name>
<name>
<surname>Bravo</surname>
<given-names>H. C.</given-names>
</name>
<name>
<surname>Pop</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Differential abundance analysis for microbial marker-gene surveys</article-title>. <source>Nat. Methods</source> <volume>10</volume> (<issue>12</issue>), <fpage>1200</fpage>&#x2013;<lpage>1202</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.2658</pub-id>
<pub-id pub-id-type="pmid">24076764</pub-id>
</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pochi</surname>
<given-names>P. E.</given-names>
</name>
<name>
<surname>Strauss</surname>
<given-names>J. S.</given-names>
</name>
<name>
<surname>Downing</surname>
<given-names>D. T.</given-names>
</name>
</person-group> (<year>1979</year>). <article-title>Age-related changes in sebaceous gland activity</article-title>. <source>J. Invest Dermatol.</source> <volume>73</volume> (<issue>1</issue>), <fpage>108</fpage>&#x2013;<lpage>111</lpage>. <pub-id pub-id-type="doi">10.1111/1523-1747.ep12532792</pub-id>
<pub-id pub-id-type="pmid">448169</pub-id>
</citation>
</ref>
<ref id="B61">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Porter</surname>
<given-names>A. W.</given-names>
</name>
<name>
<surname>Young</surname>
<given-names>L. Y.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Benzoyl-CoA, a universal biomarker for anaerobic degradation of aromatic compounds</article-title>. <source>Adv. Appl. Microbiol.</source> <volume>88</volume>, <fpage>167</fpage>&#x2013;<lpage>203</lpage>. <pub-id pub-id-type="doi">10.1016/B978-0-12-800260-5.00005-X</pub-id>
<pub-id pub-id-type="pmid">24767428</pub-id>
</citation>
</ref>
<ref id="B62">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Price</surname>
<given-names>M. N.</given-names>
</name>
<name>
<surname>Dehal</surname>
<given-names>P. S.</given-names>
</name>
<name>
<surname>Arkin</surname>
<given-names>A. P.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>FastTree 2--approximately maximum-likelihood trees for large alignments</article-title>. <source>PLoS One</source> <volume>5</volume> (<issue>3</issue>), <fpage>e9490</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0009490</pub-id>
<pub-id pub-id-type="pmid">20224823</pub-id>
</citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Quast</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Pruesse</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Yilmaz</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Gerken</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Schweer</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yarza</surname>
<given-names>P.</given-names>
</name>
<etal/>
</person-group> (<year>2013</year>). <article-title>The SILVA ribosomal RNA gene database project: improved data processing and web-based tools</article-title>. <source>Nucleic Acids Res.</source> <volume>41</volume> (<issue>Database issue</issue>), <fpage>D590</fpage>&#x2013;<lpage>D596</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gks1219</pub-id>
<pub-id pub-id-type="pmid">23193283</pub-id>
</citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rainer</surname>
<given-names>B. M.</given-names>
</name>
<name>
<surname>Thompson</surname>
<given-names>K. G.</given-names>
</name>
<name>
<surname>Antonescu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Florea</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Mongodin</surname>
<given-names>E. F.</given-names>
</name>
<name>
<surname>Bui</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Characterization and analysis of the skin microbiota in rosacea: a Case&#x2013;control study</article-title>. <source>Am. J. Clin. dermatology</source> <volume>21</volume> (<issue>1</issue>), <fpage>139</fpage>&#x2013;<lpage>147</lpage>. <pub-id pub-id-type="doi">10.1007/s40257-019-00471-5</pub-id>
<pub-id pub-id-type="pmid">31502207</pub-id>
</citation>
</ref>
<ref id="B65">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ritchie</surname>
<given-names>M. E.</given-names>
</name>
<name>
<surname>Phipson</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Law</surname>
<given-names>C. W.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Limma powers differential expression analyses for RNA-sequencing and microarray studies</article-title>. <source>Nucleic Acids Res.</source> <volume>43</volume> (<issue>7</issue>), <fpage>e47</fpage>. <pub-id pub-id-type="doi">10.1093/nar/gkv007</pub-id>
<pub-id pub-id-type="pmid">25605792</pub-id>
</citation>
</ref>
<ref id="B66">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robeson</surname>
<given-names>M. S.</given-names>
</name>
<name>
<surname>O&#x2019;Rourke</surname>
<given-names>D. R.</given-names>
</name>
<name>
<surname>Kaehler</surname>
<given-names>B. D.</given-names>
</name>
<name>
<surname>Ziemski</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Dillon</surname>
<given-names>M. R.</given-names>
</name>
<name>
<surname>Foster</surname>
<given-names>J. T.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>RESCRIPt: reproducible sequence taxonomy reference database management</article-title>. <source>PLoS Comput. Biol.</source> <volume>17</volume> (<issue>11</issue>), <fpage>e1009581</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pcbi.1009581</pub-id>
<pub-id pub-id-type="pmid">34748542</pub-id>
</citation>
</ref>
<ref id="B67">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robinson</surname>
<given-names>M. D.</given-names>
</name>
<name>
<surname>McCarthy</surname>
<given-names>D. J.</given-names>
</name>
<name>
<surname>Smyth</surname>
<given-names>G. K.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>edgeR: a Bioconductor package for differential expression analysis of digital gene expression data</article-title>. <source>Bioinformatics</source> <volume>26</volume> (<issue>1</issue>), <fpage>139</fpage>&#x2013;<lpage>140</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btp616</pub-id>
<pub-id pub-id-type="pmid">19910308</pub-id>
</citation>
</ref>
<ref id="B68">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>R&#xf6;ttjers</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Faust</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>From hairballs to hypotheses-biological insights from microbial networks</article-title>. <source>FEMS Microbiol. Rev.</source> <volume>42</volume> (<issue>6</issue>), <fpage>761</fpage>&#x2013;<lpage>780</lpage>. <pub-id pub-id-type="doi">10.1093/femsre/fuy030</pub-id>
<pub-id pub-id-type="pmid">30085090</pub-id>
</citation>
</ref>
<ref id="B69">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rudden</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Herman</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Rose</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bawdon</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Cox</surname>
<given-names>D. S.</given-names>
</name>
<name>
<surname>Dodson</surname>
<given-names>E.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>The molecular basis of thioalcohol production in human body odour</article-title>. <source>Sci. Rep.</source> <volume>10</volume> (<issue>1</issue>), <fpage>12500</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-020-68860-z</pub-id>
<pub-id pub-id-type="pmid">32719469</pub-id>
</citation>
</ref>
<ref id="B70">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Salter</surname>
<given-names>S. J.</given-names>
</name>
<name>
<surname>Cox</surname>
<given-names>M. J.</given-names>
</name>
<name>
<surname>Turek</surname>
<given-names>E. M.</given-names>
</name>
<name>
<surname>Calus</surname>
<given-names>S. T.</given-names>
</name>
<name>
<surname>Cookson</surname>
<given-names>W. O.</given-names>
</name>
<name>
<surname>Moffatt</surname>
<given-names>M. F.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>Reagent and laboratory contamination can critically impact sequence-based microbiome analyses</article-title>. <source>BMC Biol.</source> <volume>12</volume>, <fpage>87</fpage>. <pub-id pub-id-type="doi">10.1186/s12915-014-0087-z</pub-id>
<pub-id pub-id-type="pmid">25387460</pub-id>
</citation>
</ref>
<ref id="B71">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sanford</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L. J.</given-names>
</name>
<name>
<surname>Williams</surname>
<given-names>M. R.</given-names>
</name>
<name>
<surname>Gangoiti</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>C. M.</given-names>
</name>
<name>
<surname>Gallo</surname>
<given-names>R. L.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Inhibition of HDAC8 and HDAC9 by microbial short-chain fatty acids breaks immune tolerance of the epidermis to TLR ligands</article-title>. <source>Sci. Immunol.</source> <volume>1</volume> (<issue>4</issue>), <fpage>eaah4609</fpage>. <pub-id pub-id-type="doi">10.1126/sciimmunol.aah4609</pub-id>
<pub-id pub-id-type="pmid">28783689</pub-id>
</citation>
</ref>
<ref id="B72">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sarafidou</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Alexakou</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Talioti</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Bakopoulou</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Anastassiadou</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>The oral microbiome in older adults &#x2013;a state-of-the-art review</article-title>. <source>Archives Gerontology Geriatrics Plus</source> <volume>1</volume> (<issue>4</issue>), <fpage>100061</fpage>. <pub-id pub-id-type="doi">10.1016/j.aggp.2024.100061</pub-id>
</citation>
</ref>
<ref id="B73">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schreiner</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Gooris</surname>
<given-names>G. S.</given-names>
</name>
<name>
<surname>Pfeiffer</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lanzend&#xf6;rfer</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Wenck</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Diembeck</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2000</year>). <article-title>Barrier characteristics of different human skin types investigated with X-ray Diffraction, lipid analysis, and electron microscopy imaging</article-title>. <source>J. Investigative Dermatology</source> <volume>114</volume> (<issue>4</issue>), <fpage>654</fpage>&#x2013;<lpage>660</lpage>. <pub-id pub-id-type="doi">10.1046/j.1523-1747.2000.00941.x</pub-id>
<pub-id pub-id-type="pmid">10733669</pub-id>
</citation>
</ref>
<ref id="B74">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Severn</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Williams</surname>
<given-names>M. R.</given-names>
</name>
<name>
<surname>Shahbandi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Bunch</surname>
<given-names>Z. L.</given-names>
</name>
<name>
<surname>Lyon</surname>
<given-names>L. M.</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>The ubiquitous human skin commensal Staphylococcus hominis protects against opportunistic pathogens</article-title>. <source>MBio</source> <volume>13</volume> (<issue>3</issue>), <fpage>e0093022</fpage>. <pub-id pub-id-type="doi">10.1128/mbio.00930-22</pub-id>
<pub-id pub-id-type="pmid">35608301</pub-id>
</citation>
</ref>
<ref id="B75">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shibagaki</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Suda</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Clavaud</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Bastien</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Takayasu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Iioka</surname>
<given-names>E.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Aging-related changes in the diversity of women&#x2019;s skin microbiomes associated with oral bacteria</article-title>. <source>Sci. Rep.</source> <volume>7</volume> (<issue>1</issue>), <fpage>10567</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-10834-9</pub-id>
<pub-id pub-id-type="pmid">28874721</pub-id>
</citation>
</ref>
<ref id="B76">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Yi</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Ge</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Integrated analysis of facial microbiome and skin physio-optical properties unveils cutotype-dependent aging effects</article-title>. <source>Microbiome</source> <volume>12</volume> (<issue>1</issue>), <fpage>163</fpage>&#x2013;<lpage>19</lpage>. <pub-id pub-id-type="doi">10.1186/s40168-024-01891-0</pub-id>
<pub-id pub-id-type="pmid">39232827</pub-id>
</citation>
</ref>
<ref id="B77">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Swaney</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Kalan</surname>
<given-names>L. R.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Living in your skin: microbes, molecules, and mechanisms</article-title>. <source>Infect. Immun.</source> <volume>89</volume>, <fpage>e00695-20</fpage>. <pub-id pub-id-type="doi">10.1128/iai.00695-20</pub-id>
<pub-id pub-id-type="pmid">33468585</pub-id>
</citation>
</ref>
<ref id="B78">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Swaney</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Sandstrom</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Kalan</surname>
<given-names>L. R.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Cobamide sharing is predicted in the human skin microbiome</article-title>. <source>mSystems</source> <volume>7</volume> (<issue>5</issue>), <fpage>e0067722</fpage>. <pub-id pub-id-type="doi">10.1128/msystems.00677-22</pub-id>
<pub-id pub-id-type="pmid">35968974</pub-id>
</citation>
</ref>
<ref id="B79">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Swaney</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Nelsen</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Sandstrom</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Kalan</surname>
<given-names>L. R.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Sweat and sebum preferences of the human skin Microbiota</article-title>. <source>Microbiol. Spectr.</source> <volume>11</volume> (<issue>1</issue>), <fpage>e0418022</fpage>. <pub-id pub-id-type="doi">10.1128/spectrum.04180-22</pub-id>
<pub-id pub-id-type="pmid">36602383</pub-id>
</citation>
</ref>
<ref id="B80">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tipton</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>M&#xfc;ller</surname>
<given-names>C. L.</given-names>
</name>
<name>
<surname>Kurtz</surname>
<given-names>Z. D.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Kleerup</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Morris</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Fungi stabilize connectivity in the lung and skin microbial ecosystems</article-title>. <source>Microbiome</source> <volume>6</volume> (<issue>1</issue>), <fpage>12</fpage>. <pub-id pub-id-type="doi">10.1186/s40168-017-0393-0</pub-id>
<pub-id pub-id-type="pmid">29335027</pub-id>
</citation>
</ref>
<ref id="B81">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tobin</surname>
<given-names>D. J.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Introduction to skin aging</article-title>. <source>J. Tissue Viability</source> <volume>26</volume> (<issue>1</issue>), <fpage>37</fpage>&#x2013;<lpage>46</lpage>. <pub-id pub-id-type="doi">10.1016/j.jtv.2016.03.002</pub-id>
<pub-id pub-id-type="pmid">27020864</pub-id>
</citation>
</ref>
<ref id="B82">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Townsend</surname>
<given-names>E. C.</given-names>
</name>
<name>
<surname>Kalan</surname>
<given-names>L. R.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>The dynamic balance of the skin microbiome across the lifespan</article-title>. <source>Biochem. Soc. Trans.</source> <volume>51</volume> (<issue>1</issue>), <fpage>71</fpage>&#x2013;<lpage>86</lpage>. <pub-id pub-id-type="doi">10.1042/BST20220216</pub-id>
<pub-id pub-id-type="pmid">36606709</pub-id>
</citation>
</ref>
<ref id="B83">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Waller</surname>
<given-names>J. M.</given-names>
</name>
<name>
<surname>Maibach</surname>
<given-names>H. I.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Age and skin structure and function, a quantitative approach (II): protein, glycosaminoglycan, water, and lipid content and structure</article-title>. <source>Skin Res. Technol.</source> <volume>12</volume> (<issue>3</issue>), <fpage>145</fpage>&#x2013;<lpage>154</lpage>. <pub-id pub-id-type="doi">10.1111/j.0909-752X.2006.00146.x</pub-id>
<pub-id pub-id-type="pmid">16827688</pub-id>
</citation>
</ref>
<ref id="B84">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ying</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zeng</surname>
<given-names>D. N.</given-names>
</name>
<name>
<surname>Chi</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Tan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Galzote</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Cardona</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>The influence of age and gender on skin-associated microbial communities in urban and rural human populations</article-title>. <source>PLOS ONE</source> <volume>10</volume> (<issue>10</issue>), <fpage>e0141842</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0141842</pub-id>
<pub-id pub-id-type="pmid">26510185</pub-id>
</citation>
</ref>
<ref id="B85">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Hunt</surname>
<given-names>R. L.</given-names>
</name>
<name>
<surname>Villaruz</surname>
<given-names>A. E.</given-names>
</name>
<name>
<surname>Fisher</surname>
<given-names>E. L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Q.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Commensal Staphylococcus epidermidis contributes to skin barrier homeostasis by generating protective ceramides</article-title>. <source>Cell host and microbe</source> <volume>30</volume> (<issue>3</issue>), <fpage>301</fpage>&#x2013;<lpage>313.e9</lpage>. <pub-id pub-id-type="doi">10.1016/j.chom.2022.01.004</pub-id>
<pub-id pub-id-type="pmid">35123653</pub-id>
</citation>
</ref>
<ref id="B86">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>LinDA: linear models for differential abundance analysis of microbiome compositional data</article-title>. <source>Genome Biol.</source> <volume>23</volume> (<issue>1</issue>), <fpage>95</fpage>. <pub-id pub-id-type="doi">10.1186/s13059-022-02655-5</pub-id>
<pub-id pub-id-type="pmid">35421994</pub-id>
</citation>
</ref>
<ref id="B87">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Fleming</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Legendre</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Roux</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Latreille</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Gendronneau</surname>
<given-names>G.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Skin microbiome attributes associate with biophysical skin ageing</article-title>. <source>Exp. Dermatol</source> <volume>32</volume> (<issue>9</issue>), <fpage>1546</fpage>&#x2013;<lpage>1556</lpage>. <pub-id pub-id-type="doi">10.1111/exd.14863</pub-id>
<pub-id pub-id-type="pmid">37350224</pub-id>
</citation>
</ref>
<ref id="B88">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zouboulis</surname>
<given-names>C. C.</given-names>
</name>
<name>
<surname>Boschnakow</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Chronological ageing and photoageing of the human sebaceous gland</article-title>. <source>Clin. Exp. Dermatol</source> <volume>26</volume> (<issue>7</issue>), <fpage>600</fpage>&#x2013;<lpage>607</lpage>. <pub-id pub-id-type="doi">10.1046/j.1365-2230.2001.00894.x</pub-id>
<pub-id pub-id-type="pmid">11696064</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>