Microbiota fingerprints lose individually identifying features over time
© The Author(s). 2016
Received: 6 May 2016
Accepted: 22 November 2016
Published: 9 January 2017
Humans host individually unique skin microbiota, suggesting that microbiota traces transferred from skin to surfaces could serve as forensic markers analogous to fingerprints. While it is known that individuals leave identifiable microbiota traces on surfaces, it is not clear for how long these traces persist. Moreover, as skin and surface microbiota change with time, even persistent traces may lose their forensic potential as they would cease to resemble the microbiota of the person who left them. We followed skin and surface microbiota within households for four seasons to determine whether accurate microbiota-based matching of individuals to their households could be achieved across long time delays.
While household surface microbiota traces could be matched to the correct occupant or occupants with 67% accuracy, accuracy decreased substantially when skin and surface samples were collected in different seasons, and particularly when surface samples were collected long after skin samples. Most OTUs persisted on skin or surfaces for less than one season, indicating that OTU loss was the major cause of decreased matching accuracy. OTUs that were more useful for individual identification persisted for less time and were less likely to be deposited from skin to surface, suggesting a trade-off between the longevity and identifying value of microbiota traces.
While microbiota traces have potential forensic value, unlike fingerprints they are not static and may degrade in a way that preferentially erases features useful in identifying individuals.
KeywordsBuilt environment Microbiota Skin microbiota Forensics
The transfer of personal skin microbiota traces to surfaces has been called a microbial ‘fingerprint’ [1, 2], reflecting both the individuality of skin microbiota  and the potential of such traces for use in forensic identification. Identification of individuals from their microbiota traces (hereafter ‘microbiota matching’) has been demonstrated with traces left on computer keyboards and mice , mobile phones , and household surfaces [5, 6]. However, an important forensic property of fingerprints is that they can persist unchanged on surfaces for long periods, allowing people to be reliably matched to fingerprints they deposited on a surface some time ago. Little attention has been given to the time for which identifiable microbiota traces persist on surfaces, and given that an individual’s skin microbiota can change significantly within weeks or months [7, 8], it is unclear whether an individual’s current skin microbiota would be similar enough to an older microbiota trace for a reliable identification to be made. As interest in the forensic potential of microbiota matching continues to grow, it is important that this practical limitation be investigated.
This study attempted microbial matching of individuals to their places of residence based on comparison of skin microbiota to household surface microbiota traces, to determine whether accurate microbial matching could be achieved even with large time delays between skin and surface sampling. While sources including outdoor air and pets [9, 10] can contribute to residential surface microbiota, occupant skin contributes a large proportion or majority of surface microbiota [5, 11, 12]. Occupants’ skin microbiota rapidly colonise a newly occupied residence, and a person leaving a residence can cause a decline in microbiota similarity within days , suggesting that household surface assemblages closely track changes in occupant skin microbiota. This study also examined Operational Taxonomic Unit (OTU) stability on skin and surfaces and applied survival analysis to determine whether the value of OTUs in identifying individuals was related to their temporal stability and chance of deposition from skin to surface.
Results and discussion
Summary of collected samples
Values (number of samples)
Air (144), skin (380), surface (288)
Winter 2014 (203), Spring 2014 (203), Summer 2014 (203), Autumn 2014 (203)
Admiralty A (68), Fortress Hill (68), Ma On Shan (108), Quarry Bay (88), Sai Wan (88), Sha Tin Wai (88), Tai Koo (108), Tuen Mun B (108), Wu Kai Sha (88)
Bed headboard surface (36), bedroom air (36), blanket surface (36), forehead skin (76), fridge door seal surface (36), kitchen air (36), kitchen ventilator surface (36), left forearm skin (76), left palm skin (76), living room air (36), remote control surface (36), right forearm skin (76), right palm skin (76), shower curtain surface (36), toilet air (36), toilet flush button surface (36), TV screen surface (36)
Individual (skin only)
Admiralty A z (20), Fortress Hill z (20), Ma On Shan w (20), Ma On Shan x (20), Ma On Shan y (20), Quarry Bay y (20), Quarry Bay z (20), Sai Wan y (20), Sai Wan z (20), Sha Tin Wai x (20), Sha Tin Wai y (20), Tai Koo x (20), Tai Koo y (20), Tai Koo z (20), Tuen Mun B w (20), Tuen Mun B x (20), Tuen Mun B y (20), Wu Kai Sha w (20), Wu Kai Sha y (20)
When we used surface samples collected a season or more after skin samples, the rate fell, with no accurate matches made after a delay of three seasons. Likewise, when we used skin samples collected after surface samples, the accuracy fell to 33% for delays of two or three seasons. While accuracy decreased with delay in both directions, skin-after-surface delays tended to yield a higher accuracy rate than surface-after-skin delays of the same magnitude. This effect of delay on matching is consistent with previous reports that the human skin microbiota changes with time [7, 8] and that microbiota traces deposited on surfaces begin to degrade within hours even in the absence of cleaning or other mechanical removal . In this study, as the residences were continuously occupied, the most likely mechanisms for OTU loss are cleaning and other household activities as well as displacement of older traces by newly acquired skin OTUs. In other circumstances, such as the deposition of a trace in a public space or in a space vacated soon after the deposition, the degradation of identifying traces may occur faster in the absence of continuous deposition from the original host.
We confirmed these results by comparing UniFrac distances between skin and surface samples, another method that has been successfully employed for microbiota matching . The mean distance between skin and surface samples from the same household was lowest for samples taken contemporaneously, and increased with increasing delay between skin and surface sampling in both directions (Fig. 2b, c), with the exception of a small decrease in the weighted UniFrac distance between skin-after-surface delays of two to three seasons. Weighted but not unweighted UniFrac distances were generally lower for skin-after-surface delays than surface-after-skin delays of the same length of time, similar to the skewed distributions of SourceTracker matching accuracy.
We hypothesised that some OTUs are more useful for microbial matching than others, as they uniquely identify or are strongly associated with an individual or subset of individuals. However, by virtue of their specificity, these OTUs may be particularly susceptible to loss from skin microbiota , failure to be deposited on a surface as a microbiota trace, or loss from a deposited trace. To test whether this is the case, we used two measures of the identifying potential of OTUs, indicator value  and hitting set membership, and examined the relationship between these measures and OTU stability and deposition from skin to surface. Hitting sets , as applied to microbiota matching, are algorithmically determined minimal sets of OTUs that uniquely identify individuals within a cohort, thereby making membership of an OTU in a hitting set a useful marker of that OTU’s value in microbiota matching. Franzosa et al. reported that while hitting sets are relatively stable across time in gut and oral microbiota, allowing up to 80% accurate re-identification of individuals over time, skin microbiota OTU hitting sets are relatively unstable and do not permit accurate identification on repeat sampling . We found a similar result in this dataset (Additional file 1: Figure S2), consistent with the hypothesis that the decreasing accuracy of microbiota matching over time is related to the loss of highly identifying OTUs. While the family Moraxellaceae was the most commonly represented in the hitting sets, as it was among all skin microbiota, overall the taxonomic distribution of hitting set OTUs did not closely resemble that of the total skin microbiota, with non-skin-associated families such as the Sphingomonadaceae, Rhodobacteraceae and Weeksellaceae among the most common hitting set members (Additional file 1: Table S1).
One explanation for the preferential loss of identifying OTUs from skin is simply that they tend to be less abundant and are therefore less stable over time. While the positive correlation between abundance and stability would tend to support this, other factors suggest that this explanation is insufficient: the hazard model included terms for OTU abundance; indicator value was significantly positively correlated with abundance across all OTUs, individuals and seasons (Pearson’s r = 0.20, p < 0.05); and the hitting set algorithm prioritises abundant OTUs. It is more likely that identifying OTUs tend to be transients, acquired by an individual through chance environmental encounters to which other individuals were not exposed. Because these transient OTUs are less likely to be adapted to the skin environment and would not have their populations replenished by repeat exposure, they would be more exposed to eventual loss regardless of their abundance.
Abundant OTUs were also more likely to be deposited from skin to a household surface, with OTUs in the 99th percentile of abundance having 3.9× the baseline probability. More abundant OTUs have more cells available for deposition but also tend to persist for longer and thereby have increased time in which a chance deposition event can occur. Hitting set membership and lower indicator value slightly decreased the probability of a deposition event relative to the baseline, although the probability increased slightly with increasing indicator value with values between 0.9–1 having 1.6× baseline probability. This may reflect a slightly higher probability of a deposition event when an OTU is found on a higher proportion of an individual’s body sites, one of the properties measured by indicator value.
This study suggests that the frequently invoked analogy between microbiota traces and fingerprints is misleading: unlike a fingerprint, skin microbiota changes over time both on the host and in microbial traces left on surfaces. Moreover, these changes are not random but select against low-abundance microorganisms as well as those most useful in identifying the individual who left the trace. Additional work is needed to better quantify the persistence of microbiota traces and how persistence depends on environmental factors. This study was only able to detect OTU loss on seasonal timescales; a similar study with a temporal resolution of hours or days could better quantify the expected lifespan of a trace. Further, this study examined only one type of indoor environment, a continually occupied residence where occupants were in frequent contact with household surfaces and input from other sources was comparatively unimportant. Forensic applications of microbiota matching may require the use of traces from other environments, such as public spaces or spaces with a significant microbial input from environmental sources, and for which data on matching accuracy is currently lacking.
Sample collection and DNA sequencing
Sample collection, genomic DNA (gDNA) extraction and sequencing were conducted as previously described . Briefly, skin, surface and air samples were collected from a range of residences in Hong Kong across four seasons in 2014 (Table 1). Participants were instructed to maintain normal household routines, including cleaning schedules, during the study period. The V4 region of the 16S rRNA gene was amplified with the 515F/806R primer pair . Library construction and paired-end sequencing were performed by the Health GeneTech Corporation (Taoyuan, Taiwan).
Read quality control and OTU formation were performed following the UPARSE pipeline  using USEARCH version 8.1.1861. Raw sequencing reads were trimmed to a uniform length of 138 bp and filtered for a maximum of 1 expected error per read using the USEARCH command fastq_filter. These parameters were selected to maximise read retention while minimising expected errors, based on empirical analysis of the raw reads with the USEARCH command fastq_stats. Reads shorter than 138 bp or with ≥1 expected error after trimming were discarded. Forward and reverse reads were separated based on the identity of the first five nucleotides with the primer sequences, with sequences that did not perfectly match either primer in these positions discarded. As the paired ends did not overlap enough to allow reliable merging by alignment, only forward reads were used for subsequent analysis. Reads were dereplicated with a custom perl script, and dereplicated sequences were clustered at 97% sequence similarity with the USEARCH command cluster_OTUs to form OTU representative sequences, with singleton OTUs excluded. Probable chimeric OTU representative sequences were identified with the USEARCH command uchime_ref with the -strand plus option against the RDP classifier training database (downloaded 2 November 2015 from drive5.com/uchime/rdp_gold.fa). Reads were recruited to OTUs using the USEARCH command usearch_global with parameters -strand plus, -id 0.97, -maxaccepts 8, -maxrejects 64 and -top_hit_only. Consensus taxonomic lineages were assigned to OTUs with the QIIME  version 1.9.1 script assign_taxonomy.py against the Greengenes  version 13_8 97% similar 16S database. OTU representative sequences were aligned against the aligned Greengenes 97% similar 16S database with PyNAST  and the QIIME script align_seqs.py, and a tree constructed with FastTree  and the QIIME script make_phylogeny.py. With each batch of samples submitted to the sequencing facility, a kit control (i.e. DNA extracted from an unused swab or filter in parallel with samples ) was also submitted (total of nine kit control samples). Reads from these controls were not included in OTU formation, but control reads were subjected to the same quality control steps as sample reads and recruited against the OTU representative sequences generated from sample reads. A control OTU table was produced using a custom script that excluded OTUs matching any of the following conditions: representative sequence identified as chimeric; singleton OTU; representative sequence failed PyNAST alignment; OTU assigned to class ‘Chloroplast’ or family ‘mitochondria’. Following Flores et al. , any OTU found at ≥5% relative abundance in any control sample was designated a likely contaminant (Additional file 1: Figure S3). The relative abundance threshold (≥5%) was set higher than that used by Flores et al. (≥1%) due to the larger number of kit controls in this study. The final disposition of all reads sequenced for this project (including those removed for quality control purposes) is given in Additional file 1: Figure S4. A sample OTU table was produced that excluded OTUs matching any of the above conditions as well as likely contaminant OTUs (Additional file 2: Table S2). To account for differences in sequencing depth between samples, all samples were normalised by random subsampling to the number of reads in the most depauperate sample (1381 reads). To confirm that skin was the major source of surface OTUs, the Bayesian SourceTracker method  was used to estimate the proportional contributions of occupant skin and household air to surface samples within each household and season. Kit control samples were included as potential sources as a negative control.
To confirm that microbiota comparison could be used to reliably match occupants to residences, and to investigate the effect of sampling delay on matching accuracy, SourceTracker was used to estimate the proportional contributions of occupant skin microbiota (‘sources’) to household surface microbiota (‘sinks’). Source estimation and scoring of match accuracy were performed as previously described . Briefly, input OTU tables were prepared for all pairwise combinations of residence and season, with surface samples from the target residence in all seasons as sinks, and skin samples from occupants of all residences in the target season as potential sources. This arrangement made available all skin samples within each season as potential sources for each surface sample, without permitting sources to be drawn from multiple seasons. OTUs present in <10% of samples within each input table were excluded. Source contributions were estimated using the SourceTracker script sourcetracker_for_qiime.R (v1.0, downloaded 4 November 2015 from https://github.com/danknights/sourcetracker/archive/v1.0.tar.gz) with default settings. Each sink sample was considered to have an accurate match if the source contribution for occupants of that residence was greater than for occupants of any other residence (following Lax et al. ). The proportion of sink samples with accurate matches for each sampling delay (number of seasons between collection of skin and surface samples, expressed as an integer between −3 and +3) was taken as the accuracy rate for that delay. To investigate the effect of the number of possible matches on matching accuracy, microbiota matching with no season delay was repeated with separate random subsets of two, four, six and eight potential individuals or groups of individuals (including the correct match) selected as potential matches for each residence. Accuracy rates were calculated for each subset size as above. Weighted and unweighted UniFrac  distances were calculated between all samples using the QIIME script beta_diversity.py with default settings.
OTU stability and deposition
OTU hitting sets  for each individual and body site were constructed to determine the stability of uniquely identifying OTUs and to explore the utility of hitting sets for microbiota matching. Hitting sets were generated with Franzosa et al.’s  python script idability.py (downloaded 13 December 2015 from https://bitbucket.org/biobakery/idability/get/default.tar.gz), run with the --meta_mode relab option to optimise for metagenome-like data in relative abundance format. The stability of winter hitting sets as unique identifiers across seasons was determined by running idability.py with the --codes flag, which attempts to re-identify individuals based on hitting sets.
Cox proportional hazard models  were used to determine the effect of a set of OTU properties on the stability of OTUs on occupant skin and household surfaces and on the deposition of OTUs from skin to surfaces. Proportional hazard models fit the instantaneous probability of an event (the ‘hazard’, in this case an OTU loss or deposition event) as a function of selected covariates, which can be either fixed or time-dependent. A useful property of these models is that as long as the baseline hazard function does not vary with time (i.e. the hazard varies only on some proportional combination of the covariates; the ‘proportional hazards assumption’), the proportional effect of each covariate on the hazard (the covariate’s ‘hazard ratio’) can be determined even if the baseline hazard is unknown. A second useful property is that these models allow for right-censored observations, where the period of observation may expire before an event is observed. These properties make proportional hazards models well suited for studying OTU stability and deposition over a fixed period.
Three models were created for this study. The first (the ‘skin stability’ models) examined the loss of skin OTUs as a function of time-dependent covariates representing OTU abundance, indicator value and membership in a hitting set for the modelled occupant body site. Each OTU on each body site served as a separate observation. Relative abundances were binned into four percentile ranges of increasing relative abundance (0–69th percentiles, 70–89th percentiles, 90–98th percentiles and 99th percentile) within each model (i.e. within each occupant body site). The OTU indicator value (IndVal ), a weighted measure of the OTU’s specificity (the proportion of the total abundance of that OTU that is found on the occupant) and fidelity (the proportion of body sites on that occupant in which the OTU is found), was calculated for each occupant in each season with the indval function of the R package labdsv  and binned into four ranges (0–0.2, 0.2–0.6, 0.6–0.9, 0.9–1). The hitting set membership covariate was a logical (true/false) value representing whether the OTU was a member of a hitting set for the modelled occupant and body site, in the season in which the OTU was first observed. The final covariate was the OTU’s indicator value (0–1; non-significant values masked to 0) for the modelled occupant. For each OTU, the observation period began in the season in which the OTU was first observed on the modelled occupant and body site and ended either in the last season in which it was observed (loss event) or in autumn (right-censoring with no loss event). If an OTU was not observed for an interim season(s) but then returned, this interim was not counted as a loss event.
The next set of models examined the stability of OTUs on household surfaces. Only OTUs found in at least one skin sample from an occupant of the modelled residence in any season as well as at least one surface sample from the modelled residence and surface in any season were included. This constrained the models to consider only OTUs that likely originated from occupant skin rather than other sources. Each OTU on each household surface served as a separate observation. Covariates were included for OTU abundance, indicator value (maximum value among occupants of that residence) and hitting set membership (of a hitting set for any occupant of that residence). Determination of loss events and right-censoring were as for the skin stability models.
The final model set investigated OTU deposition from skin to surface as a function of the same set of covariates used in the skin stability model. Each OTU on each body site served as a separate observation. The model covariates were set as with the skin stability models, and response variable was a deposition event, defined as an observation of the skin OTU on at least one household surface in the same or a later season as that of first observation on skin. OTUs present on a modelled surface before they were first observed on any occupants’ skin were excluded from the model.
All models were fit with the coxph function from the R library survival (version 2.38 [29, 30]) using survival objects built with the Surv function with right-censoring. Survival objects were constructed in the right-censored counting format, with each transition between seasons represented by an interval with associated covariate values for the season beginning that interval and a response variable representing whether a loss or deposition event occurred at the end of the interval. Significance for each covariate was determined as the two-tailed p < 0.05 of the observed Wald statistic (z, ratio of the fitted coefficient to its standard error) under the null hypothesis of a hazard ratio of 1 (i.e. no effect), with non-significant covariates rejected.
We thank Ho Lam Hedwig Cheung and William Chan for the sample processing and staff of the Health GeneTech Corporation for their expertise in sequencing. Special thanks to the members of the sampling team (Ellen Li, Fred Kong, Wai Shan Chow, Catherine Chung, Ka Yan Ng, Yuet Ying Wong and Flora Yeh). We are grateful for the participation of the cohort in this study.
This research was supported by the Research Grants Council of Hong Kong through Project 11211815. Study design, sample collection, analysis, interpretation and manuscript writing were performed independently by the authors.
Availability of data and materials
Sequencing reads generated for this project have been deposited at the NCBI Sequence Read Archive (SRA), under the study accession SRP094801. A package containing custom scripts and workflows used in this study is available on FigShare, DOI:10.6084/m9.figshare.3968820.
PKHL, DW and MHYL conceived the experiment; MHYL performed the sample collection and DNA extraction; DW performed the data analysis; DW and PKHL wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Consent was given by subjects to release personal and related data for publication as needed.
Ethics approval and consent to participate
All participants were fully informed of the nature of the study and gave written informed consent, and the study was approved by the City University of Hong Kong Human Subjects Ethics Sub-Committee (ref: 3-2-201312 [H000334]).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Fierer N, et al. Forensic identification using skin bacterial communities. Proc Natl Acad Sci U S A. 2010;107:6477–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Wood M, et al. Athletic equipment microbiota are shaped by interactions with human skin. Microbiome. 2015;3:25.View ArticlePubMedPubMed CentralGoogle Scholar
- Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.View ArticleGoogle Scholar
- Meadow JF, Altrichter AE, Green JL. Mobile phones carry the personal microbiome of their owners. Peer J. 2014;2:e447.View ArticlePubMedPubMed CentralGoogle Scholar
- Lax S, et al. Longitudinal analysis of microbial interaction between humans and the indoor environment. Science. 2014;345:1048–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilkins D, Leung MHY, Lee PKH. Indoor air bacterial communities in Hong Kong households assemble independently of occupant skin microbiomes. Environ Microbiol. 2016;18:1754–63.View ArticlePubMedGoogle Scholar
- Franzosa EA, et al. Identifying personal microbiomes using metagenomic codes. Proc Natl Acad Sci U S A. 2015;112:E2930–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Costello EK, et al. Bacterial community variation in human body habitats across space and time. Science. 2009;326:1694–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Dunn RR, Fierer N, Henley JB, Leff JW, Menninger HL. Home life: factors structuring the bacterial diversity found within and between homes. PLoS ONE. 2013;8:e64133.View ArticlePubMedPubMed CentralGoogle Scholar
- Fujimura KE, et al. Man’s best friend? The effect of pet ownership on house dust microbial communities. J Allergy Clin Immunol. 2010;126:410–412.e3.View ArticlePubMedPubMed CentralGoogle Scholar
- Flores GE, et al. Diversity, distribution and sources of bacteria in residential kitchens. Environ Microbiol. 2012;15:588–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Adams RI, Bateman AC, Bik HM, Meadow JF. Microbiota of the indoor environment: a meta-analysis. Microbiome. 2015;3:1–18.View ArticleGoogle Scholar
- Walters WA, et al. PrimerProspector: de novo design and taxonomic analysis of barcoded polymerase chain reaction primers. Bioinformatics. 2011;27:1159–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Knights D, et al. Bayesian community-wide culture-independent microbial source tracking. Nat Meth. 2011;8:761–3.View ArticleGoogle Scholar
- Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Env Microbiol. 2005;71:8228–35.View ArticleGoogle Scholar
- Leung MHY, Wilkins D, Li EKT, Kong FKF, Lee PKH. Indoor-air microbiome in an urban subway network: diversity and dynamics. Appl Env Microbiol. 2014;80:6760–70.View ArticleGoogle Scholar
- Ruiz-Calderon JF, et al. Walls talk: microbial biogeography of homes spanning urbanization. Sci Adv. 2016;2:e1501061.View ArticlePubMedPubMed CentralGoogle Scholar
- Gibbons SM, et al. Ecological succession and viability of human-associated microbiota on restroom surfaces. Appl Env Microbiol. 2015;81:765–73.View ArticleGoogle Scholar
- Dufrêne M, Legendre P. Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol Monogr. 1997;67:345–66.Google Scholar
- Caporaso JG, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Meth. 2013;10:996–8.View ArticleGoogle Scholar
- Caporaso JG, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Meth. 2010;7:335–6.View ArticleGoogle Scholar
- DeSantis TZ, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Env Microbiol. 2006;72:5069–72.View ArticleGoogle Scholar
- Caporaso JG, et al. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010;26:266–7.View ArticlePubMedGoogle Scholar
- Price MN, Dehal PS, Arkin AP. FastTree 2—approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010;5:e9490.View ArticlePubMedPubMed CentralGoogle Scholar
- Flores GE, Henley JB, Fierer N. A direct PCR approach to accelerate analyses of human-associated microbial communities. PLoS ONE. 2012;7:e44563.View ArticlePubMedPubMed CentralGoogle Scholar
- Cox DR. Regression models and life-tables. J R Stat Soc Ser B Methodol. 1972;34:187–220.Google Scholar
- Roberts DW. labdsv: ordination and multivariate analysis for ecology. R package version 1.7-0. (2015).Google Scholar
- Therneau, T. M. survival: a package for survival analysis in S. R package version 2.38. (2015).Google Scholar
- Therneau TM & Grambsch PM. Modeling survival data: extending the Cox model. (Springer-Verlag, New York, 2000).Google Scholar