Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Ecological and Evolutionary Science

Untangling Species-Level Composition of Complex Bacterial Communities through a Novel Metagenomic Approach

Christian Milani, Giulia Alessandri, Marta Mangifesta, Leonardo Mancabelli, Gabriele Andrea Lugli, Federico Fontana, Giulia Longhi, Rosaria Anzalone, Alice Viappiani, Sabrina Duranti, Francesca Turroni, Renato Costi, Alfredo Annicchiarico, Andrea Morini, Leopoldo Sarli, Maria Cristina Ossiprandi, Douwe van Sinderen, Marco Ventura
J. Gregory Caporaso, Editor
Christian Milani
aLaboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, Parma, Italy
bMicrobiome Research Hub, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Giulia Alessandri
cDepartment of Veterinary Medical Science, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marta Mangifesta
aLaboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Leonardo Mancabelli
aLaboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gabriele Andrea Lugli
aLaboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Federico Fontana
aLaboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Giulia Longhi
aLaboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, Parma, Italy
dGenProbio srl, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rosaria Anzalone
dGenProbio srl, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alice Viappiani
dGenProbio srl, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sabrina Duranti
aLaboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Francesca Turroni
aLaboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, Parma, Italy
bMicrobiome Research Hub, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Renato Costi
eDepartment of Medicine and Surgery, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alfredo Annicchiarico
eDepartment of Medicine and Surgery, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrea Morini
eDepartment of Medicine and Surgery, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Leopoldo Sarli
eDepartment of Medicine and Surgery, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Maria Cristina Ossiprandi
bMicrobiome Research Hub, University of Parma, Parma, Italy
cDepartment of Veterinary Medical Science, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Douwe van Sinderen
fAPC Microbiome Institute and School of Microbiology, Bioscience Institute, National University of Ireland, Cork, Ireland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Douwe van Sinderen
Marco Ventura
aLaboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, Parma, Italy
bMicrobiome Research Hub, University of Parma, Parma, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Marco Ventura
J. Gregory Caporaso
Northern Arizona University
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00404-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

16S small-subunit (SSU) rRNA gene-based bacterial profiling is the gold standard for cost-effective taxonomic reconstruction of complex bacterial populations down to the genus level. However, it has been proven ineffective in clinical and research settings requiring higher taxonomic resolution. We therefore developed a bacterial profiling method based on the internal transcribed spacer (ITS) region employing optimized primers and a comprehensive ITS database for accurate cataloguing of bacterial communities at (sub)species resolution. Performance of the microbial ITS profiling pipeline was tested through analysis of host-associated, food, and environmental matrices, while its efficacy in clinical settings was assessed through analysis of mucosal biopsy specimens of colorectal cancer, leading to the identification of putative novel biomarkers. The data collected indicate that the proposed pipeline represents a major step forward in cost-effective identification and screening of microbial biomarkers at (sub)species level, with relevant impact in research, industrial, and clinical settings.

IMPORTANCE We developed a novel method for accurate cataloguing of bacterial communities at (sub)species level involving amplification of the internal transcribed spacer (ITS) region through optimized primers, followed by next-generation sequencing and taxonomic classification of amplicons by means of a comprehensive database of bacterial ITS sequences. Host-associated, food, and environmental matrices were employed to test the performance of the microbial ITS profiling pipeline. Moreover, mucosal biopsy samples from colorectal cancer patients were analyzed to demonstrate the scientific relevance of this profiling approach in a clinical setting through identification of putative novel biomarkers. The results indicate that the ITS-based profiling pipeline proposed here represents a key metagenomic tool with major relevance for research, industrial, and clinical settings.

INTRODUCTION

The relevance of microbial communities, i.e., the microbiota, has been demonstrated for a range of environmental and host-associated ecological niches, including soil, water, and air, as well as in animal and human body sites (1–6). Due to their impact on many physiological and immunological aspects of their hosts, microbial populations harbored by humans have been associated with various diseases and disorders. Limitations in the taxonomic resolution of the commonly used microbial profiling approach based on the 16S rRNA gene have hampered identification of microbial biomarkers and cost-effective, population-wide screening for bacterial species of clinical relevance.

Currently, shotgun and long-read sequencing methods provide higher taxonomic resolution, although their higher cost represents an important limitation compared to the use of 16S rRNA-based microbial profiling approach (7–9). Moreover, the absence of a PCR amplification step may prevent retrieval of sufficient microbial DNA for analysis of host-associated, low-abundance bacterial communities due to interfering levels of eukaryotic DNA (7).

To overcome the limitations of 16S rRNA gene microbial profiling and shotgun metagenomics, we present here a pipeline for reliable and cost-effective profiling of microbial communities at the (sub)species level. In detail, we extended the internal transcribed spacer (ITS)-profiling approach described previously for bifidobacteria and lactobacilli (10–12) to all members of the domain Bacteria. ITS sequences are characterized by higher variability than 16S rRNA genes, thus allowing (sub)species taxonomic resolution when employed for metagenomic profiling purposes (10–12). Remarkably, the ITS profiling method combines the lower cost and high sensitivity of a marker gene amplification approach with the resolution of shotgun metagenomics (7, 8).

RESULTS AND DISCUSSION

The bacterial ITS profiling pipeline.Our ITS profiling pipeline encompasses universal ITS primers for bacteria, a database including ITS sequences retrieved from bacterial genomes currently available in the NCBI database, and a ready-to-use bioinformatic script for QIIME v.2.

Genomic sequences of 120,748 bacterial strains were retrieved from NCBI and processed for prediction of rRNA loci through RNAmmer (13) in order to build a database encompassing all retrieved 16S-ITS-23S genomic regions. The resulting universal microbial ITS database (UMID) comprises 131,795 sequences covering 10,361 bacterial species. UMID is available in both FASTA and QIIME v.2 artifact formats (http://probiogenomics.unipr.it/pbi/) and was used for all analyses performed in this study. Moreover, an alternative database named UMID-RefSeq is also provided in QIIME v.2 artifact format as part of the script package (http://probiogenomics.unipr.it/pbi/). The latter was obtained by extracting the 16S-ITS-23S region from genomes available in the NCBI RefSeq database. Notably, the average nucleotide identity (ANI) matrix precomputed and updated by NCBI is used for correction of taxonomic assignments, leading to 91,309 database entries which represent 4,202 bacterial species. A list of bacterial species covered by UMID and UMID-RefSeq is available in Data Set S1. Furthermore, a script to generate an updated UMID-RefSeq is available for download (http://probiogenomics.unipr.it/pbi/).

DATA SET S1

UNI_ITS_fw/UNI_ITS_rv primers and UMID statistics. Download Data Set S1, XLSX file, 0.6 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Universal primers for bacteria were designed through alignment of 16S-ITS-23S sequences in the UMID and processing through PrimerProspector (14) and the SILVA database v.132 (15) Test Probe web application (https://www.arb-silva.de/search/testprobe/). Primers were designed at the 3′ end of the 16S rRNA sequence (UNI_ITS_fw, 5′-KRGGRYKAAGTCGTAACAAG-3′) and the 5′ end of the 23S rRNA sequence (UNI_ITS_rv, 5′-TTTTCRYCTTTCCCTCACGG-3′), corresponding to positions 1484 to 1504 and 460 to 480 of 16S and 23S rRNA genes of Escherichia coli K-12 strain MG1655, respectively, matching conserved regions of ribosomal genes in order to maximize taxonomic coverage (Fig. 1). In this regard, 100 16S rRNA and 100 23S rRNA gene sequences spanning a wide range of bacterial taxonomy were retrieved from the SILVA database and aligned using ClustalX (16) in order to identify a suitable conserved region at the 3′ end of the desired amplicon. Alignments are available in Clustal format as Data Sets S5 and S6 for 16S rRNA and 23S rRNA sequences, respectively. Primer sequences can be found at position 1680 in the 16S rRNA sequence alignment and at position 718 in the 23S rRNA sequence alignment. Furthermore, PrimerProspector (14) and the SILVA database v.132 (15) Test Probe web application (https://www.arb-silva.de/search/testprobe/) were used to identify mismatches in alignments with specific taxonomic ranks (on a phylum and class basis), leading to the iterative introduction of IUPAC bases in the primer sequence in order to maximize taxonomic coverage.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Positions of primers UNI_ITS_fw and UNI_ITS_rv within the ribosomal locus of Escherichia coli K-12 strain MG1655.

In silico PCR amplification of the UMID with UNI_ITS_fw/UNI_ITS_rv primers by means of FastPCR software (17) revealed that the average amplicon size is 972.6 ± 153.7 nucleotides (nt). The sequencing technology used (Illumina, San Diego, CA, USA) has a recommended fragment or amplicon size of 200 to 500 bp, with larger fragment/amplicon sizes, such as 1,000 bp and higher, causing increasingly lower sequencing efficiency and lower sequence yield (18). For this reason, we amended the standard library preparation protocol by employing 15 μl of Agencourt AMPure XP DNA purification beads (Beckman Coulter Genomics GmbH, Bernried, Germany) for the first purification step in order to remove primer dimers. In addition, we performed the second purification step using 30 μl of the above-mentioned purification beads. With these adaptations, we did not observe any sequencing efficiency issues with libraries generated from samples sequenced as part of the current work, especially with the R1/forward reads, confirming recently published data (19). The designed primers target the 16S rRNA-proximal side of the amplicons for taxonomic assignment (Fig. 1). Thus, in the case of paired-end sequencing of the PCR amplicon, only the read starting with the UNI_ITS_fw primer sequence is used for taxonomic assignment in bioinformatic analyses. Alignment of primers against the bacterial sequences contained in the SILVA REF small-subunit (SSU) and large-subunit (LSU) databases, v.132, with five mismatches allowed through the testProbe web application (https://www.arb-silva.de/search/testprobe/) is reported in the bacterial ITS profiling download page of our website (http://probiogenomics.unipr.it/pbi/). Briefly, 89.7% of the aligned 16S rRNA sequences and 85.0% of aligned 23S rRNA sequences reported a perfect match, while an additional 8.2% and 11.8% showed a single mismatch with the 16S rRNA and 23S rRNA genes. Remarkably, just 1.1% of 16S rRNA and 0.3% of 23S rRNA sequences showed ≥3 mismatches, ensuring broad amplification coverage of all known bacterial taxa (Table S1). Moreover, an additional validation step was applied through the PrimerProspector (14) analyze_primers.py and taxa_coverage.py application using default settings and the SILVA NR99 database v.132 (15) in FASTA format. Only 16S rRNA gene sequences of >1,500 bp were used, in order to ensure coverage of the 3′ end. Notably, results revealed average taxonomic coverage of 95.64% and 94.92% for UNI_ITS_fw alignment to16S sequences and UNI_ITS_rv alignment to 23S sequences, respectively, and an absence of rank-specific biases (Data Set S1).

TABLE S1

Percentages of reads mapped to the SILVA SSU and LSU databases. Download Table S1, DOCX file, 0.02 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

In this context, nonmicrobial low-biomass PCR amplification would require a higher number of mismatches, as revealed by PrimerProspector (14) analyze_primers.py for Homo sapiens, Mus musculus, and Arabidopsis thaliana nuclear and mitochondrial DNA as well as Arabidopsis thaliana for chloroplast DNA (Data Set S1).To allow easy integration of the methodology into the QIIME v.2 analysis pipeline, we developed a ready-to-use bash script for Linux and MacOS operating systems along with a preformatted database in artifact format for taxonomic classification in a QIIME 2 environment (http://probiogenomics.unipr.it/pbi/) (20, 21). Notably, all ITS-profiling analyses reported in this study were performed using QIIME 2 v.2018.11, but the package available on our website has also been tested to be compatible with QIIME 2 v. 2020.2 and QIIME 2 v. 2020.6, i.e., the latest release available at this time. A guide to using this script is available in Text S1. Briefly, fastq files are filtered to remove reads without the forward primer UNI_ITS_fw. DADA2 (22) is then employed to denoise sequences, dereplicate sequence variants, and remove chimeras, and the amplicon sequence variants (ASVs) obtained are classified at the species level using the QIIME 2 feature-classifier classify-consensus-vsearch (23) method. Furthermore, alpha and beta diversity analyses using multiple statistical metrics are performed.

TEXT S1

User manual for the analysis script available at http://probiogenomics.unipr.it/pbi/, discussion of the comparison between ITS profiling and shotgun metagenomics results, and detailed discussion of ITS profiling data obtained from the clinical test case encompassing colorectal cancer biopsy specimens. Download Text S1, DOCX file, 0.1 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

A group of six data sets in fastq format for testing purposes is also available in the “Test” package that can be downloaded from our website (http://probiogenomics.unipr.it/pbi/). These data sets consist of ITS amplicon data obtained from analysis of the bacterial communities harbored by human saliva, human vagina, and cheese samples and are accompanied by an Excel file reporting the expected taxonomic profiles.

In silico evaluation of profiling performances.In order to evaluate the efficacy of the proposed analysis pipeline across all taxonomic ranks covered by the reference database, we developed three artificial data sets constituted by UNI_ITS_fw/UNI_ITS_rv 250-bp forward amplicons. Each of these data sets included sequences from 500 different randomly selected species, and they were named 500_A, 500_B, and 500_C (Data Set S2). Analysis of these artificial communities revealed that 97.4% and 87.87% of the amplicons were correctly classified at the genus and species levels, respectively, and confirmed the absence of biases in accuracy of specific taxonomic ranks (Data Set S2).

DATA SET S2

Results obtained from the analysis of artificial datasets. Download Data Set S2, XLSX file, 0.3 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Furthermore, we generated artificial communities constituted by UNI_ITS_fw/UNI_ITS_rv 250 bp forward amplicons obtained from genomes of species representative of bacterial communities typically found in biological matrices of high scientific relevance, i.e., human feces from adults and infants, vagina, sputum, lung, and skin, as well as cow’s milk and cheese. Moreover, the same genomes were used to predict 450-bp amplicons, representing joined 250-bp paired-end reads with 50-bp overlaps, corresponding to multiple 16S rRNA gene variable regions. In detail, 450 bp starting from primers 27F (3′-AGAGTTTGATCMTGGCTCAG-5′), ProbioUni_fw (3′-CCTACGGGRSGCAGCAG-5′), 785F (3′-GGATTAGATACCCTGGTA-5′), 1100F (3′-YAACGAGCGCAACCC-5′), and 1492R (3′-CGGTTACCTTGTTACGACTT-5′) were used to generate artificial data sets covering the V1-V2, V3-V4, V4-V6, V7-V8, and V9-V8 regions of the 16S rRNA gene. The list of species used for each matrix was derived from recent literature (5, 24–35) and shotgun metagenomic analyses reported in this study (see below). Amplicons corresponding to each species were added to each data set at specific relative abundances, reported as expected profiles in Data Set S2, for a total of approximately 50,000 reads. 16S rRNA-based data sets were processed using the 16S rRNA database SILVA v.132 and the sklearn classifier (36) with default settings in a Qiime2 v. 2020.02 environment (20, 21). Remarkably, the results obtained revealed that the profiles generated by ITS profiling match the expected taxonomic profiles, with a species-level classification accuracy of 93.5% (Data Set S2; Fig. 2). In contrast, 16S rRNA gene profiling showed lower accuracy of species-level taxonomic reconstruction. In detail, V1-V2, V3-V4, V4-V6, V7-V8, and V9-V8 regions of the 16S rRNA gene allowed a classification of 30%, 20%, 22%, 18%, and 20%, respectively, of the species across the analyzed artificial data sets (Data Set S2; Fig. 2). While these data clearly demonstrate the higher taxonomic resolution of ITS profiling approach with respect to 16S rRNA gene microbial profiling, it is worth mentioning that synthetic data sets generated and used for test purposes in this study encompass the same sequences present in the databases, since both were retrieved from publicly available genomes. Thus, the reported performances may vary in real settings and the evaluation of the impact of sequence variability with respect to database might require additional ad hoc testing.

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

In silico evaluation of microbial ITS profiling performance and comparison respect to different hypervariable 16S rRNA gene regions. (a) Percentage of bacterial taxa correctly classified by ITS and 16S rRNA gene profiling methods with respect to expected taxonomic profiles of artificial data sets. (b) Taxonomic profiles observed using ITS and 16S rRNA gene profiling methods and their comparison with respect to expected taxonomic profiles of artificial data sets.

Specificity assessment through profiling of an artificial community.To assess specificity, we performed microbial ITS and 16S rRNA gene profiling of an artificial bacterial community composed of 26 species in equal viable-cell numbers, including multiple species of the same genus (Table S2; Data Set S3). 16S rRNA gene-based profiling was performed using a primer set with broad taxonomic coverage, as shown by comparison with other widely used primers (29, 37–39). Remarkably, microbial ITS profiling successfully identified all species within the artificial community (Data Set S3), whereas 16S rRNA gene profiling was unable to classify at the species level or misclassified 16 of the 26 (61.5%) bacterial species (Data Set S3). While a number of misclassifications observed in the 16S rRNA profiling analysis may be due to issues in the SILVA database v.132, overall, these data confirm higher precision and taxonomic resolution of the ITS-based approach versus those achieved by 16S rRNA-based cataloguing.

TABLE S2

Quality filtering table of samples analyzed in this study. Download Table S2, DOCX file, 0.03 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

DATA SET S3

Results obtained from analysis of the artificial community. Download Data Set S3, XLSX file, 0.01 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Evaluation of profiling performance with biological samples.Performance of the microbial ITS profiling pipeline was furthermore tested through analysis of five samples that cover eight different matrices representing both environmental and host-associated microbial populations (Table S2). Alpha diversity analysis through rarefaction curves built up to 10,000 reads revealed that all matrices reach a plateau (Data Set S4), indicating that the microbial ITS profiling approach requires sequencing depths and costs that are comparable with those of 16S rRNA gene profiling. Remarkably, an overall average of 67.8% of the reads, and up to 100%, were successfully classified at the species level (Fig. 3; Data Set S4). Of note, the discriminatory power of microbial ITS-based profiling is dependent on the availability of reference ITSs extracted from genomic sequences, similar to profiling methods based on shotgun metagenomics data (7, 8, 40–42). Therefore, the percentage of classified species is expected to increase due to ongoing genome decoding of novel bacterial species.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Performance of microbial ITS profiling compared to 16S rRNA gene cataloguing and read-based classification of shotgun metagenomic data sets. (a) Average percentage of reads classified to the species level using ITS and 16S rRNA gene profiling for each matrix analyzed. Data obtained for each sample are reported in Data Set S4. (b) Average percentage of reads classified by 16S rRNA gene profiling as a genus also detected by microbial ITS profiling of the same read, for each matrix analyzed. Data obtained for each sample are reported in Data Set S4. (c) Percentage of reads classified at the species level based on read length.

DATA SET S4

Results obtained from the analysis of real biological samples. Download Data Set S4, XLSX file, 1.3 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Microbial compositional data obtained with microbial ITS profiling were compared to those from 16S rRNA gene microbial profiling using de novo ASV generation (at 100% identity) and classification with the feature-classifier classify-sklearn method (43) and the SILVA database v. 132, as suggested by the QIIME v.2 manual (44) (Table S2). Notably, comparison of ITS and 16S rRNA gene-based profiles at the genus level revealed that the two approaches have comparable performances at this taxonomic rank: an average of 84.0% of the reads that were classified at the genus level by 16S rRNA gene profiling were assigned to the same taxonomic rank by the microbial ITS approach, while an average of 85.9% of reads classified at the genus level by microbial ITS profiling showed a taxonomy also reported by 16S rRNA profiling (Fig. 3; Data Set S4). Conversely, the average percentage of reads classified at the species level by 16S rRNA gene profiling is 6.2%, which is substantially lower than that obtained by ITS profiling (67.8%) (Fig. 3; Data Set S4). In this context, one should keep in mind that the limited number of taxa properly profiled by 16S rRNA sequences, yet not detected by ITS, is mostly due to the current wider taxonomic coverage of 16S rRNA data sets than of the ITS database, which was generated from available microbial genomes with a validated ribosomal locus. Nevertheless, as argued above, the number of species classified by the ITS profiling approach is expected to increase due to ongoing genome decoding of novel bacterial species. A comparison of microbial ITS profiling with taxonomic reconstruction from shotgun metagenomics data is available in Text S1.

Furthermore, data sets consisting of 250-bp reads obtained from ITS amplicons of a given sample and representing each of the eight profiled biological matrices were shortened to generate artificial ITS amplicon data sets with read lengths of 200, 150, 100, or 75 bp in order to determine the minimum read length needed to achieve sufficient coverage of the ITS region for optimal taxonomic classification. Intriguingly, the data collected indicated that reads with a length of ≥150 bp provided comparable resolution and discriminatory power, as indicated by average percentages of reads classified at the species level of 80.7%, 78.4%, and 78.3% for data at 250 bp, 200 bp, and 150 bp, respectively (Data Set S4; Fig. 3). In contrast, 100-bp and 75-bp reads showed a percentage of reads classified at the species level lower than 71.1% (Data Set S4; Fig. 3).

The microbial ITS profiling pipeline was employed to assess colon biopsy specimens to identify microbial biomarkers of colorectal cancer (CRC). In detail, two mucosal samples were collected from the same region of the colon from 15 individuals diagnosed with CRC, corresponding to adenocarcinoma and healthy mucosa (Data Set S4; Fig. S1). Analysis of microbial ITS data at the species level revealed that Fusobacterium nucleatum, which has been proposed as a key microbial marker of CRC (45, 46), is present at higher relative abundance in adenocarcinoma mucosa in all cases (Fig. S1). Moreover, Fusobacterium hwasookii, Fusobacterium periodonticum, Clostridium chauvoei, and Lactobacillus oris were also observed to be present at higher relative abundances in adenocarcinoma than in healthy mucosa in all individuals in which these species were detected; thus, they may represent additional microbial biomarkers of CRC (Fig. S1). Though performing colonoscopies for biopsy sample collection obviates the need for profiling resident microbiota for CRC detection, the identification of above-mentioned biomarkers at the species level is important for the establishment of a CRC diagnostic protocol based on fecal microbiota analysis. In this context, profiling of the fecal microbiota followed by quantitative normalization by total bacterial cell counts obtained from qPCR or cytofluorimetric assays (47) or the use of novel approaches such as biosensors for rapid quantitative detection of biomarker-specific DNA in feces may represent key tools for cost-effective population screenings.

FIG S1

Identification of microbial biomarkers of colorectal cancer based on species-level data. The graphical representation reports the bacterial species that, when detected in the mucosal samples analyzed, showed higher abundance in adenocarcinoma than in healthy mucosa of the same individual. Download FIG S1, PDF file, 0.2 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Conclusions.The currently available methods for complex bacterial community profiling do not allow a cost-effective reconstruction of the bacterial population at the species level. For this reason, we developed a comprehensive ITS-based pipeline for profiling at the (sub)species level that encompasses a bacterial ITS database, optimized universal primers for bacteria, and a ready-to-use bioinformatic script. Taxonomic profiling performance was tested through analysis of an artificial community as well as five different host-associated, food, or environmental matrices, and the data obtained were compared with genus-level profiles obtained by means of 16S rRNA gene microbial profiling. Moreover, the efficacy of this novel approach was also validated in a clinical test case consisting of screening of 15 individuals affected by CRC for bacterial species representing known microbial biomarkers. Intriguingly, accurate species-level data also led to the identification of novel biomarkers of CRC.

Based on collected results, the microbial ITS profiling approach represents an innovative tool for cost-effective and accurate identification and screening of microbial biomarkers at the (sub)species level with relevant impact in research, industrial, and clinical settings.

MATERIALS AND METHODS

Ethics approval and consent to participate.All experimental procedures and protocols involving animals were approved by the Veterinarian Animal Care and Use Committee of Parma University and conducted in accordance with the European Community Council Directives dated 22 September 2010 (2010/63/UE). Human participants gave informed written consent before enrollment. Analyses of human samples were performed in the framework of previous studies approved by the Comitato Etico per Parma. All investigations were carried out following the principles of the Declaration of Helsinki.

Sample collection.For the purpose of this study, a total of 70 samples were collected, encompassing 30 colonoscopic biopsy specimens and five biological samples of each of the following matrices: human and animal feces, human vaginal swabs, human sputum, soil, litter from dairy cattle’s husbandries, Parmesan cheese, and dairy cattle milk (Table S2). In detail, by means of endoscopy biopsy forceps, two biopsy samples, one from the adenocarcinoma area and one from the adjacent healthy mucosa, were collected from 15 patients affected by colorectal cancer. Human and animal feces were collected immediately after defecation. Humans and animals included in this study had not taken antibiotics during the previous 6 months. Cow milk samples were harvested directly by hand during milking, after the cow teat ends were cleaned and disinfected, while Parmesan cheese samples were retrieved by trimming fresh Parmesan cheeses. Furthermore, litters were recovered from the ground of different husbandries, while soil samples were collected from different fields. In all cases, immediately after collection, samples were kept on ice and shipped frozen to the laboratory, where they were preserved at −80°C, until they were processed.

Human and animal fecal samples together with sputum and litter samples were subjected to DNA extraction using the QIAamp DNA stool minikit (Qiagen, Germany). A ZymoBIOMICS DNA miniprep kit (Zymo Research Corporation, USA) was used for DNA extraction from vaginal swabs, while DNA extraction from milk and Parmesan cheese samples was performed using the DNeasy mastitis minikit (Qiagen, Germany). DNA from colonoscopic biopsy specimens were extracted using the AllPrep PowerViral DNA/RNA kit (Qiagen, Germany), while the DNeasy PowerSoil kit (Qiagen, Germany) was used for DNA extraction from soil samples. In all cases, the DNA extractions were performed following the manufacturers’ instructions.

16S rRNA gene sequencing.Partial 16S rRNA gene sequences were amplified from extracted DNA using the primer pair Probio_Uni (5′-CCTACGGGRSGCAGCAG-3′) and Probio_Rev (5′-ATTACCGCGGCTGCT-3′) targeting the V3 region of the 16S rRNA gene sequence (37). Illumina adapter overhang nucleotide sequences were added to the partial 16S rRNA gene-specific amplicons, which were further processed by means of the 16S metagenomic sequencing library preparation protocol (part 15044223, rev. B; Illumina). Amplifications were carried out using a Verity thermocycler (Applied Biosystems). The integrity of the PCR amplicons was analyzed by electrophoresis on a 2200 TapeStation instrument (Agilent Technologies, USA). DNA products obtained following PCR-mediated amplification of the 16S rRNA gene sequences were purified by a magnetic purification step employing Agencourt AMPure XP DNA purification beads (Beckman Coulter Genomics GmbH, Bernried, Germany) in order to remove primer dimers. DNA concentration of the amplified sequence library was determined by a fluorometric Qubit quantification system (Life Technologies, USA). Amplicons were diluted to a concentration of 4 nM, and 5-μl quantities of each diluted DNA amplicon sample were mixed to prepare the pooled final library. Sequencing was performed using an Illumina MiSeq sequencer with MiSeq reagent kit v3 chemicals.

Microbial ITS profiling.ITS sequences were amplified from extracted DNA using the primer pair UNI_ITS_fw (5′-KRGGRYKAAGTCGTAACAAG-3′) and UNI_ITS_rv (5′-TTTTCRYCTTTCCCTCACGG-3′), targeting the entire spacer region between the 16S rRNA and 23 rRNA genes within the rRNA locus. The amplification was carried out using GoTaq G2 Hot Start polymerase (Promega, USA) on a Verity thermocycler (Applied Biosystems, USA) according to the following protocol: 95°C for 10 min, followed by 32 cycles of 95°C for 1 min, 52°C for 1 min, and 72°C for 1 min and a final step of 72°C for 5 min. The integrity of PCR amplicons was analyzed by gel electrophoresis. The library of ITS amplicons was prepared according to the 16S metagenomic sequencing library preparation protocol (part 15044223, rev. B; Illumina) with modifications in the purification steps. Specifically, the first purification involved 15 μl of Agencourt AMPure XP DNA purification beads (Beckman Coulter Genomics GmbH, Bernried, Germany) in order to remove primer dimers. Then, the second purification step was performed using 30 μl of the above-mentioned purification beads. Sequencing was performed using an Illumina MiSeq sequencer with MiSeq reagent kit v3 chemicals, using 300 cycles.

16S rRNA microbial profiling analysis.The fastq files were processed using a custom script based on the QIIME2 software suite (20, 44). Paired-end read pairs were assembled to reconstruct the complete Probio_Uni/Probio_Rev amplicons. Quality control retained sequences with a length between 140 and 400 bp and mean sequence quality score of >20, while sequences with homopolymers of >7 bp and mismatched primers were removed to reduce quality issues (though Illumina sequencing is not known to be affected by homopolymer sequences). In order to calculate downstream diversity measures, 16S rRNA amplicon sequence variants (ASVs) were defined at 100% sequence homology using DADA2 (22); ASVs not encompassing at least two sequences of the same sample were removed. Notably, this approach allows highly distinctive taxonomic classification at single nucleotide accuracy (20). All reads were classified to the lowest possible taxonomic rank using QIIME2 (20, 44) and a reference data set from the SILVA database (v. 132) (15). Biodiversity within a given sample (alpha diversity) was calculated considering the observed ASVs for 10 subsamplings of the total read pool.

Microbial ITS profiling analysis.The fastq files were processed using a custom script based on the QIIME software suite (20, 44). The forward read obtained from paired-end sequencing, corresponding to the UNI_ITS_fw amplicon, was used for taxonomic reconstruction. Quality control retained sequences with a quality score of >20, while sequences with homopolymers of >7 bp and mismatched UNI_ITS_fw primer were omitted. In order to calculate downstream diversity measures, ITS ASVs were defined at 100% sequence homology using DADA2 (22); ASVs not encompassing at least two sequences of the same sample were removed. Notably, this approach allows highly distinctive taxonomic classification at single-nucleotide accuracy (20). All reads were classified to the lowest possible taxonomic rank using QIIME2 (20, 44) and the reference data set UMID (probiogenomics.unipr.it/pbi). Biodiversity within a given sample (alpha diversity) was calculated considering the observed ASVs for 10 subsamplings of the total read pool.

Shotgun metagenomics.The extracted DNA was prepared following the Illumina Nextera XT DNA library preparation kit. Briefly, the DNA samples were enzymatically fragmented, barcoded, and purified by using magnetic beads. Samples were then quantified using fluorometric Qubit quantification system (Life Technologies, USA), loaded on a 2200 Tape Station instrument (Agilent Technologies, USA), and normalized to 4 nM. Sequencing was performed using an Illumina NextSeq 500 sequencer with NextSeq High Output v2 kit chemicals for 150 cycles.

Analysis of metagenomic data sets.The obtained fastq files were filtered for reads with a quality of <25, for reads of >80 bp. Moreover, bases were removed from the end of the reads unless the average quality score was >25, in a window of 5 bp. Quality-filtered data were used to further analysis with METAnnotatorX (8) and MetaPhlAn2 (42) for taxonomic profile reconstruction. Both software packages were used with default settings, and the RefSeq database obtained from NCBI in October 2019 was used for taxonomic classification using METAnnotatorX (8).

Data availability.Raw sequences of Universal ITS and 16S rRNA gene profiling data coupled with shotgun metagenomics data are accessible through BioProject accession number PRJNA562817.

DATA SET S5

Alignment file in Clustal format showing alignment of 100 16S rRNA gene sequences used to identify the optimal conserved region for primer design. Download Data Set S5, TXT file, 0.8 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

DATA SET S6

Alignment file in Clustal format showing alignment of 100 23S rRNA gene sequences used to identify the optimal conserved region for primer design. Download Data Set S6, TXT file, 1.1 MB.
Copyright © 2020 Milani et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

ACKNOWLEDGMENTS

This work was primarily funded by the EU Joint Programming Initiative–A Healthy Diet for a Healthy Life (JPI HDHL; http://www.healthydietforhealthylife.eu/) to D.V.S. (in conjunction with Science Foundation Ireland [SFI], grant number 15/JP-HDHL/3280) and to M.V. (in conjunction with MIUR, Italy). D.V.S. is a member of APC Microbiome Ireland, funded by Science Foundation Ireland through the Irish Government’s National Development Plan (grant no. SFI/12/RC/2273-P1 and SFI/12/RC/2273-P2). The study was supported by Fondazione Cariparma, under TeachInParma Project (DV). G.A. is supported by Fondazione Cariparma, Parma, Italy. We furthermore thank GenProbio srl for financial support of the Laboratory of Probiogenomics.

This research benefited from the HPC (High Performance Computing) facility of the University of Parma, Italy. We thank Caseificio La Lovetta for its support in collecting animal, environmental, and Parmesan cheese samples used in this study.

FOOTNOTES

    • Received May 7, 2020.
    • Accepted July 16, 2020.
  • Copyright © 2020 Milani et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Fierer N
    . 2017. Embracing the unknown: disentangling the complexities of the soil microbiome. Nat Rev Microbiol 15:579–590. doi:10.1038/nrmicro.2017.87.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Ling F,
    2. Whitaker R,
    3. LeChevallier MW,
    4. Liu WT
    . 2018. Drinking water microbiome assembly induced by water stagnation. ISME J 12:1520–1531. doi:10.1038/s41396-018-0101-5.
    OpenUrlCrossRef
  3. 3.↵
    1. Dommergue A,
    2. Amato P,
    3. Tignat-Perrier R,
    4. Magand O,
    5. Thollot A,
    6. Joly M,
    7. Bouvier L,
    8. Sellegri K,
    9. Vogel T,
    10. Sonke JE,
    11. Jaffrezo JL,
    12. Andrade M,
    13. Moreno I,
    14. Labuschagne C,
    15. Martin L,
    16. Zhang Q,
    17. Larose C
    . 2019. Methods to investigate the global atmospheric microbiome. Front Microbiol 10:243. doi:10.3389/fmicb.2019.00243.
    OpenUrlCrossRef
  4. 4.↵
    1. Trinh P,
    2. Zaneveld JR,
    3. Safranek S,
    4. Rabinowitz PM
    . 2018. One Health relationships between human, animal, and environmental microbiomes: a mini-review. Front Public Health 6:235. doi:10.3389/fpubh.2018.00235.
    OpenUrlCrossRef
  5. 5.↵
    1. Milani C,
    2. Duranti S,
    3. Bottacini F,
    4. Casey E,
    5. Turroni F,
    6. Mahony J,
    7. Belzer C,
    8. Delgado Palacio S,
    9. Arboleya Montes S,
    10. Mancabelli L,
    11. Lugli GA,
    12. Rodriguez JM,
    13. Bode L,
    14. de Vos W,
    15. Gueimonde M,
    16. Margolles A,
    17. van Sinderen D,
    18. Ventura M
    . 2017. The first microbial colonizers of the human gut: composition, activities, and health implications of the infant gut microbiota. Microbiol Mol Biol Rev 81:e00036-17. doi:10.1128/MMBR.00036-17.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    NIH Human Microbiome Portfolio Analysis Team. 2019. A review of 10 years of human microbiome research activities at the US National Institutes of Health, fiscal years 2007–2016. Microbiome 7:31. doi:10.1186/s40168-019-0620-y.
    OpenUrlCrossRef
  7. 7.↵
    1. Quince C,
    2. Walker AW,
    3. Simpson JT,
    4. Loman NJ,
    5. Segata N
    . 2017. Shotgun metagenomics, from sampling to analysis. Nat Biotechnol 35:833–844. doi:10.1038/nbt.3935.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Milani C,
    2. Casey E,
    3. Lugli GA,
    4. Moore R,
    5. Kaczorowska J,
    6. Feehily C,
    7. Mangifesta M,
    8. Mancabelli L,
    9. Duranti S,
    10. Turroni F,
    11. Bottacini F,
    12. Mahony J,
    13. Cotter PD,
    14. McAuliffe FM,
    15. van Sinderen D,
    16. Ventura M
    . 2018. Tracing mother-infant transmission of bacteriophages by means of a novel analytical tool for shotgun metagenomic datasets: METAnnotatorX. Microbiome 6:145. doi:10.1186/s40168-018-0527-z.
    OpenUrlCrossRef
  9. 9.↵
    1. Cusco A,
    2. Catozzi C,
    3. Vines J,
    4. Sanchez A,
    5. Francino O
    . 2018. Microbiota profiling with long amplicons using Nanopore sequencing: full-length 16S rRNA gene and the 16S-ITS-23S of the rrn operon. F1000Res 7:1755. doi:10.12688/f1000research.16817.2.
    OpenUrlCrossRef
  10. 10.↵
    1. Milani C,
    2. Lugli GA,
    3. Turroni F,
    4. Mancabelli L,
    5. Duranti S,
    6. Viappiani A,
    7. Mangifesta M,
    8. Segata N,
    9. van Sinderen D,
    10. Ventura M
    . 2014. Evaluation of bifidobacterial community composition in the human gut by means of a targeted amplicon sequencing (ITS) protocol. FEMS Microbiol Ecol 90:493–503. doi:10.1111/1574-6941.12410.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Milani C,
    2. Mangifesta M,
    3. Mancabelli L,
    4. Lugli GA,
    5. James K,
    6. Duranti S,
    7. Turroni F,
    8. Ferrario C,
    9. Ossiprandi MC,
    10. van Sinderen D,
    11. Ventura M
    . 2017. Unveiling bifidobacterial biogeography across the mammalian branch of the tree of life. ISME J 11:2834–2847. doi:10.1038/ismej.2017.138.
    OpenUrlCrossRef
  12. 12.↵
    1. Milani C,
    2. Duranti S,
    3. Mangifesta M,
    4. Lugli GA,
    5. Turroni F,
    6. Mancabelli L,
    7. Viappiani A,
    8. Anzalone R,
    9. Alessandri G,
    10. Ossiprandi MC,
    11. van Sinderen D,
    12. Ventura M
    . 2018. Phylotype-level profiling of lactobacilli in highly complex environments by means of an internal transcribed spacer-based metagenomic approach. Appl Environ Microbiol 84:e00706-18. doi:10.1128/AEM.00706-18.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Lagesen K,
    2. Hallin P,
    3. Rodland EA,
    4. Staerfeldt HH,
    5. Rognes T,
    6. Ussery DW
    . 2007. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res 35:3100–3108. doi:10.1093/nar/gkm160.
    OpenUrlCrossRefPubMedWeb of Science
  14. 14.↵
    1. Walters WA,
    2. Caporaso JG,
    3. Lauber CL,
    4. Berg-Lyons D,
    5. Fierer N,
    6. Knight R
    . 2011. PrimerProspector: de novo design and taxonomic analysis of barcoded polymerase chain reaction primers. Bioinformatics 27:1159–1161. doi:10.1093/bioinformatics/btr087.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Quast C,
    2. Pruesse E,
    3. Yilmaz P,
    4. Gerken J,
    5. Schweer T,
    6. Yarza P,
    7. Peplies J,
    8. Glockner FO
    . 2013. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41:D590–6. doi:10.1093/nar/gks1219.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Larkin MA,
    2. Blackshields G,
    3. Brown NP,
    4. Chenna R,
    5. McGettigan PA,
    6. McWilliam H,
    7. Valentin F,
    8. Wallace IM,
    9. Wilm A,
    10. Lopez R,
    11. Thompson JD,
    12. Gibson TJ,
    13. Higgins DG
    . 2007. Clustal W and Clustal X version 2.0. Bioinformatics 23:2947–2948. doi:10.1093/bioinformatics/btm404.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. Kalendar R,
    2. Khassenov B,
    3. Ramankulov Y,
    4. Samuilova O,
    5. Ivanov KI
    . 2017. FastPCR: an in silico tool for fast primer and probe design and advanced sequence analysis. Genomics 109:312–319. doi:10.1016/j.ygeno.2017.05.005.
    OpenUrlCrossRef
  18. 18.↵
    1. Quail MA,
    2. Swerdlow H,
    3. Turner DJ
    . 2009. Improved protocols for the illumina genome analyzer sequencing system. Curr Protoc Hum Genet Chapter 18:Unit 18.2.
    OpenUrlPubMed
  19. 19.↵
    1. Tan G,
    2. Opitz L,
    3. Schlapbach R,
    4. Rehrauer H
    . 2019. Long fragments achieve lower base quality in Illumina paired-end sequencing. Sci Rep 9:2856. doi:10.1038/s41598-019-39076-7.
    OpenUrlCrossRef
  20. 20.↵
    1. Caporaso JG,
    2. Kuczynski J,
    3. Stombaugh J,
    4. Bittinger K,
    5. Bushman FD,
    6. Costello EK,
    7. Fierer N,
    8. Peña AG,
    9. Goodrich JK,
    10. Gordon JI,
    11. Huttley GA,
    12. Kelley ST,
    13. Knights D,
    14. Koenig JE,
    15. Ley RE,
    16. Lozupone CA,
    17. McDonald D,
    18. Muegge BD,
    19. Pirrung M,
    20. Reeder J,
    21. Sevinsky JR,
    22. Turnbaugh PJ,
    23. Walters WA,
    24. Widmann J,
    25. Yatsunenko T,
    26. Zaneveld J,
    27. Knight R
    . 2010. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336. doi:10.1038/nmeth.f.303.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Bolyen E,
    2. Rideout JR,
    3. Dillon MR,
    4. Bokulich NA,
    5. Abnet CC,
    6. Al-Ghalith GA,
    7. Alexander H,
    8. Alm EJ,
    9. Arumugam M,
    10. Asnicar F,
    11. Bai Y,
    12. Bisanz JE,
    13. Bittinger K,
    14. Brejnrod A,
    15. Brislawn CJ,
    16. Brown CT,
    17. Callahan BJ,
    18. Caraballo-Rodríguez AM,
    19. Chase J,
    20. Cope EK,
    21. Da Silva R,
    22. Diener C,
    23. Dorrestein PC,
    24. Douglas GM,
    25. Durall DM,
    26. Duvallet C,
    27. Edwardson CF,
    28. Ernst M,
    29. Estaki M,
    30. Fouquier J,
    31. Gauglitz JM,
    32. Gibbons SM,
    33. Gibson DL,
    34. Gonzalez A,
    35. Gorlick K,
    36. Guo J,
    37. Hillmann B,
    38. Holmes S,
    39. Holste H,
    40. Huttenhower C,
    41. Huttley GA,
    42. Janssen S,
    43. Jarmusch AK,
    44. Jiang L,
    45. Kaehler BD,
    46. Kang KB,
    47. Keefe CR,
    48. Keim P,
    49. Kelley ST,
    50. Knights D,
    51. Koester I,
    52. Kosciolek T,
    53. Kreps J,
    54. Langille MGI,
    55. Lee J,
    56. Ley R,
    57. Liu Y-X,
    58. Loftfield E,
    59. Lozupone C,
    60. Maher M,
    61. Marotz C,
    62. Martin BD,
    63. McDonald D,
    64. McIver LJ,
    65. Melnik AV,
    66. Metcalf JL,
    67. Morgan SC,
    68. Morton JT,
    69. Naimey AT,
    70. Navas-Molina JA,
    71. Nothias LF,
    72. Orchanian SB,
    73. Pearson T,
    74. Peoples SL,
    75. Petras D,
    76. Preuss ML,
    77. Pruesse E,
    78. Rasmussen LB,
    79. Rivers A,
    80. Robeson MS,
    81. Rosenthal P,
    82. Segata N,
    83. Shaffer M,
    84. Shiffer A,
    85. Sinha R,
    86. Song SJ,
    87. Spear JR,
    88. Swafford AD,
    89. Thompson LR,
    90. Torres PJ,
    91. Trinh P,
    92. Tripathi A,
    93. Turnbaugh PJ,
    94. Ul-Hasan S,
    95. van der Hooft JJJ,
    96. Vargas F,
    97. Vázquez-Baeza Y,
    98. Vogtmann E,
    99. von Hippel M,
    100. Walters W,
    101. Wan Y,
    102. Wang M,
    103. Warren J,
    104. Weber KC,
    105. Williamson CHD,
    106. Willis AD,
    107. Xu ZZ,
    108. Zaneveld JR,
    109. Zhang Y,
    110. Zhu Q,
    111. Knight R,
    112. Caporaso JG
    . 2019. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol 37:852–857. doi:10.1038/s41587-019-0209-9.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Callahan BJ,
    2. McMurdie PJ,
    3. Rosen MJ,
    4. Han AW,
    5. Johnson AJ,
    6. Holmes SP
    . 2016. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods 13:581–583. doi:10.1038/nmeth.3869.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Rognes T,
    2. Flouri T,
    3. Nichols B,
    4. Quince C,
    5. Mahe F
    . 2016. VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi:10.7717/peerj.2584.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Arumugam M,
    2. Raes J,
    3. Pelletier E,
    4. Le Paslier D,
    5. Yamada T,
    6. Mende DR,
    7. Fernandes GR,
    8. Tap J,
    9. Bruls T,
    10. Batto J-M,
    11. Bertalan M,
    12. Borruel N,
    13. Casellas F,
    14. Fernandez L,
    15. Gautier L,
    16. Hansen T,
    17. Hattori M,
    18. Hayashi T,
    19. Kleerebezem M,
    20. Kurokawa K,
    21. Leclerc M,
    22. Levenez F,
    23. Manichanh C,
    24. Nielsen HB,
    25. Nielsen T,
    26. Pons N,
    27. Poulain J,
    28. Qin J,
    29. Sicheritz-Ponten T,
    30. Tims S,
    31. Torrents D,
    32. Ugarte E,
    33. Zoetendal EG,
    34. Wang J,
    35. Guarner F,
    36. Pedersen O,
    37. de Vos WM,
    38. Brunak S,
    39. Doré J,
    40. Antolín M,
    41. Artiguenave F,
    42. Blottiere HM,
    43. Almeida M,
    44. Brechot C,
    45. Cara C,
    46. Chervaux C,
    47. Cultrone A,
    48. Delorme C,
    49. Denariaz G,
    50. Dervyn R,
    51. Foerstner KU,
    52. Friss C,
    53. van de Guchte M,
    54. Guedon E,
    55. Haimet F,
    56. Huber W,
    57. van Hylckama-Vlieg J,
    58. Jamet A,
    59. Juste C,
    60. Kaci G,
    61. Knol J,
    62. Lakhdari O,
    63. Layec S,
    64. Le Roux K,
    65. Maguin E,
    66. Mérieux A,
    67. Melo Minardi R,
    68. M'rini C,
    69. Muller J,
    70. Oozeer R,
    71. Parkhill J,
    72. Renault P,
    73. Rescigno M,
    74. Sanchez N,
    75. Sunagawa S,
    76. Torrejon A,
    77. Turner K,
    78. Vandemeulebrouck G,
    79. Varela E,
    80. Winogradsky Y,
    81. Zeller G,
    82. Weissenbach J,
    83. Ehrlich SD,
    84. Bork P, MetaHIT Consortium
    . 2011. Enterotypes of the human gut microbiome. Nature 473:174–180. doi:10.1038/nature09944.
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.↵
    1. Mancabelli L,
    2. Milani C,
    3. Lugli GA,
    4. Turroni F,
    5. Ferrario C,
    6. van Sinderen D,
    7. Ventura M
    . 2017. Meta-analysis of the human gut microbiome from urbanized and pre-agricultural populations. Environ Microbiol 19:1379–1390. doi:10.1111/1462-2920.13692.
    OpenUrlCrossRef
  26. 26.↵
    1. DiGiulio DB,
    2. Callahan BJ,
    3. McMurdie PJ,
    4. Costello EK,
    5. Lyell DJ,
    6. Robaczewska A,
    7. Sun CL,
    8. Goltsman DS,
    9. Wong RJ,
    10. Shaw G,
    11. Stevenson DK,
    12. Holmes SP,
    13. Relman DA
    . 2015. Temporal and spatial variation of the human microbiota during pregnancy. Proc Natl Acad Sci U S A 112:11060–11065. doi:10.1073/pnas.1502875112.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Milani C,
    2. Alessandri G,
    3. Mancabelli L,
    4. Lugli GA,
    5. Longhi G,
    6. Anzalone R,
    7. Viappiani A,
    8. Duranti S,
    9. Turroni F,
    10. Ossiprandi MC,
    11. van Sinderen D,
    12. Ventura M
    . 2019. Bifidobacterial distribution across Italian cheeses produced from raw milk. Microorganisms 7:599. doi:10.3390/microorganisms7120599.
    OpenUrlCrossRef
  28. 28.↵
    1. Milani C,
    2. Fontana F,
    3. Alessandri G,
    4. Mancabelli L,
    5. Lugli GA,
    6. Longhi G,
    7. Anzalone R,
    8. Viappiani A,
    9. Duranti S,
    10. Turroni F,
    11. Ossiprandi MC,
    12. van Sinderen D,
    13. Ventura M
    . 2020. Ecology of lactobacilli present in Italian cheeses produced from raw milk. Appl Environ Microbiol 86:e00139-20. doi:10.1128/AEM.00139-20.
    OpenUrlAbstract/FREE Full Text
  29. 29.↵
    1. Milani C,
    2. Duranti S,
    3. Napoli S,
    4. Alessandri G,
    5. Mancabelli L,
    6. Anzalone R,
    7. Longhi G,
    8. Viappiani A,
    9. Mangifesta M,
    10. Lugli GA,
    11. Bernasconi S,
    12. Ossiprandi MC,
    13. van Sinderen D,
    14. Ventura M,
    15. Turroni F
    . 2019. Colonization of the human gut by bovine bacteria present in Parmesan cheese. Nat Commun 10:1286. doi:10.1038/s41467-019-09303-w.
    OpenUrlCrossRef
  30. 30.↵
    1. Tangedal S,
    2. Nielsen R,
    3. Aanerud M,
    4. Persson LJ,
    5. Wiker HG,
    6. Bakke PS,
    7. Hiemstra PS,
    8. Eagan TM
    . 2019. Sputum microbiota and inflammation at stable state and during exacerbations in a cohort of chronic obstructive pulmonary disease (COPD) patients. PLoS One 14:e0222449. doi:10.1371/journal.pone.0222449.
    OpenUrlCrossRef
  31. 31.↵
    1. Diao W,
    2. Shen N,
    3. Du Y,
    4. Erb-Downward JR,
    5. Sun X,
    6. Guo C,
    7. Ke Q,
    8. Huffnagle GB,
    9. Gyetko MR,
    10. He B
    . 2018. Symptom-related sputum microbiota in stable chronic obstructive pulmonary disease. Int J Chron Obstruct Pulmon Dis 13:2289–2299. doi:10.2147/COPD.S167618.
    OpenUrlCrossRef
  32. 32.↵
    1. Byrd AL,
    2. Belkaid Y,
    3. Segre JA
    . 2018. The human skin microbiome. Nat Rev Microbiol 16:143–155. doi:10.1038/nrmicro.2017.157.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Erb-Downward JR,
    2. Thompson DL,
    3. Han MK,
    4. Freeman CM,
    5. McCloskey L,
    6. Schmidt LA,
    7. Young VB,
    8. Toews GB,
    9. Curtis JL,
    10. Sundaram B,
    11. Martinez FJ,
    12. Huffnagle GB
    . 2011. Analysis of the lung microbiome in the “healthy” smoker and in COPD. PLoS One 6:e16384. doi:10.1371/journal.pone.0016384.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Hilty M,
    2. Burke C,
    3. Pedro H,
    4. Cardenas P,
    5. Bush A,
    6. Bossley C,
    7. Davies J,
    8. Ervine A,
    9. Poulter L,
    10. Pachter L,
    11. Moffatt MF,
    12. Cookson WO
    . 2010. Disordered microbial communities in asthmatic airways. PLoS One 5:e8578. doi:10.1371/journal.pone.0008578.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Beck JM,
    2. Young VB,
    3. Huffnagle GB
    . 2012. The microbiome of the lung. Transl Res 160:258–266. doi:10.1016/j.trsl.2012.02.005.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Hao JG,
    2. Ho TK
    . 2019. Machine learning made easy: a review of Scikit-learn package in Python programming language. J Educ Behav Stat 44:348–361. doi:10.3102/1076998619832248.
    OpenUrlCrossRef
  37. 37.↵
    1. Milani C,
    2. Hevia A,
    3. Foroni E,
    4. Duranti S,
    5. Turroni F,
    6. Lugli GA,
    7. Sanchez B,
    8. Martin R,
    9. Gueimonde M,
    10. van Sinderen D,
    11. Margolles A,
    12. Ventura M
    . 2013. Assessing the fecal microbiota: an optimized ion torrent 16S rRNA gene-based analysis protocol. PLoS One 8:e68739. doi:10.1371/journal.pone.0068739.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Mancabelli L,
    2. Milani C,
    3. Lugli GA,
    4. Turroni F,
    5. Cocconi D,
    6. van Sinderen D,
    7. Ventura M
    . 2017. Identification of universal gut microbial biomarkers of common human intestinal diseases by meta-analysis. FEMS Microbiol Ecol 93. doi:10.1093/femsec/fix153.
    OpenUrlCrossRef
  39. 39.↵
    1. Ticinesi A,
    2. Milani C,
    3. Guerra A,
    4. Allegri F,
    5. Lauretani F,
    6. Nouvenne A,
    7. Mancabelli L,
    8. Lugli GA,
    9. Turroni F,
    10. Duranti S,
    11. Mangifesta M,
    12. Viappiani A,
    13. Ferrario C,
    14. Dodi R,
    15. Dall'Asta M,
    16. Del Rio D,
    17. Ventura M,
    18. Meschi T
    . 2018. Understanding the gut-kidney axis in nephrolithiasis: an analysis of the gut microbiota composition and functionality of stone formers. Gut 67:2097–2106. doi:10.1136/gutjnl-2017-315734.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Wood DE,
    2. Salzberg SL
    . 2014. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol 15:R46. doi:10.1186/gb-2014-15-3-r46.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Segata N,
    2. Waldron L,
    3. Ballarini A,
    4. Narasimhan V,
    5. Jousson O,
    6. Huttenhower C
    . 2012. Metagenomic microbial community profiling using unique clade-specific marker genes. Nat Methods 9:811–814. doi:10.1038/nmeth.2066.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Truong DT,
    2. Franzosa EA,
    3. Tickle TL,
    4. Scholz M,
    5. Weingart G,
    6. Pasolli E,
    7. Tett A,
    8. Huttenhower C,
    9. Segata N
    . 2015. MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat Methods 12:902–903. doi:10.1038/nmeth.3589.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Pedregosa F,
    2. Varoquaux G,
    3. Gramfort A,
    4. Michel V,
    5. Thirion B,
    6. Grisel O,
    7. Blondel M,
    8. Prettenhofer P,
    9. Weiss R,
    10. Dubourg V,
    11. Vanderplas J,
    12. Passos A,
    13. Cournapeau D,
    14. Brucher M,
    15. Perrot M,
    16. Duchesnay E
    . 2011. Scikit-learn: machine learning in Python. J Mach Learn Res 12:2825–2830.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Bokulich NA,
    2. Kaehler BD,
    3. Rideout JR,
    4. Dillon M,
    5. Bolyen E,
    6. Knight R,
    7. Huttley GA,
    8. Caporaso JG
    . 2018. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2's q2-feature-classifier plugin. Microbiome 6:90. doi:10.1186/s40168-018-0470-z.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Brennan CA,
    2. Garrett WS
    . 2019. Fusobacterium nucleatum—symbiont, opportunist and oncobacterium. Nat Rev Microbiol 17:156–166. doi:10.1038/s41579-018-0129-6.
    OpenUrlCrossRef
  46. 46.↵
    1. Yu T,
    2. Guo F,
    3. Yu Y,
    4. Sun T,
    5. Ma D,
    6. Han J,
    7. Qian Y,
    8. Kryczek I,
    9. Sun D,
    10. Nagarsheth N,
    11. Chen Y,
    12. Chen H,
    13. Hong J,
    14. Zou W,
    15. Fang JY
    . 2017. Fusobacterium nucleatum promotes chemoresistance to colorectal cancer by modulating autophagy. Cell 170:548–563.E16. doi:10.1016/j.cell.2017.07.008.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Vandeputte D,
    2. Kathagen G,
    3. D'Hoe K,
    4. Vieira-Silva S,
    5. Valles-Colomer M,
    6. Sabino J,
    7. Wang J,
    8. Tito RY,
    9. De Commer L,
    10. Darzi Y,
    11. Vermeire S,
    12. Falony G,
    13. Raes J
    . 2017. Quantitative microbiome profiling links gut community variation to microbial load. Nature 551:507–511. doi:10.1038/nature24460.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Download PDF
Citation Tools
Untangling Species-Level Composition of Complex Bacterial Communities through a Novel Metagenomic Approach
Christian Milani, Giulia Alessandri, Marta Mangifesta, Leonardo Mancabelli, Gabriele Andrea Lugli, Federico Fontana, Giulia Longhi, Rosaria Anzalone, Alice Viappiani, Sabrina Duranti, Francesca Turroni, Renato Costi, Alfredo Annicchiarico, Andrea Morini, Leopoldo Sarli, Maria Cristina Ossiprandi, Douwe van Sinderen, Marco Ventura
mSystems Jul 2020, 5 (4) e00404-20; DOI: 10.1128/mSystems.00404-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Untangling Species-Level Composition of Complex Bacterial Communities through a Novel Metagenomic Approach
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Untangling Species-Level Composition of Complex Bacterial Communities through a Novel Metagenomic Approach
Christian Milani, Giulia Alessandri, Marta Mangifesta, Leonardo Mancabelli, Gabriele Andrea Lugli, Federico Fontana, Giulia Longhi, Rosaria Anzalone, Alice Viappiani, Sabrina Duranti, Francesca Turroni, Renato Costi, Alfredo Annicchiarico, Andrea Morini, Leopoldo Sarli, Maria Cristina Ossiprandi, Douwe van Sinderen, Marco Ventura
mSystems Jul 2020, 5 (4) e00404-20; DOI: 10.1128/mSystems.00404-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS AND DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

metagenomics
ITS
profiling

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077