- Open Access
Methods for exploring the faecal microbiome of premature infants: a review
Maternal Health, Neonatology and Perinatology volume 7, Article number: 11 (2021)
The premature infant gut microbiome plays an important part in infant health and development, and recognition of the implications of microbial dysbiosis in premature infants has prompted significant research into these issues. The approaches to designing investigations into microbial populations are many and varied, each with its own benefits and limitations. The technique used can influence results, contributing to heterogeneity across studies. This review aimed to describe the most common techniques used in researching the preterm infant microbiome, detailing their various limitations. The objective was to provide those entering the field with a broad understanding of available methodologies, so that the likely effects of their use can be factored into literature interpretation and future study design. We found that although many techniques are used for characterising the premature infant microbiome, 16S rRNA short amplicon sequencing is the most common. 16S rRNA short amplicon sequencing has several benefits, including high accuracy, discoverability and high throughput capacity. However, this technique has limitations. Each stage of the protocol offers opportunities for the injection of bias. Bias can contribute to variability between studies using 16S rRNA high throughout sequencing. Thus, we recommend that the interpretation of previous results and future study design be given careful consideration.
The premature infant gut microbiome has become an important, modifiable factor in the field of neonatal intensive care. Compared with infants born full-term, the characteristic microbiome of premature infants (born <37 weeks gestation) is dysbiotic: highly variable [1,2,3], low in diversity [4,5,6], low in common commensals [1, 6, 7], and harbouring more potential pathogens [8, 9]. This dysbiotic microbiome composition puts immune-compromised premature infants at an increased risk of acute and chronic disease, and developmental abnormalities [10,11,12]. Premature infants are also more likely to be born via caesarean section, be formula fed, receive antibiotics, and spend much of their early life in a clinical environment, all of which have the potential to exacerbate the microbial dysbiosis [1, 2, 13, 14].
Unfortunately, understanding this microbial composition is confounded by the diversity of investigative methods used. Methodologies for examining the microbiome are complex, technically challenging and vary between laboratories. Therefore, it is impossible to rule out protocol bias as a factor contributing to the variability seen between studies. In fact, a number of studies have demonstrated the role of methodological bias in influencing the outcomes of microbiome analysis [15, 16], thus contributing to significant heterogeneity in results between studies.
This review aimed to describe the most common techniques used in researching the preterm infant microbiome. We chose to focus on those studies investigating preterm infants, due to the explosion of interest into this area, although these techniques can also be applied to microbiome study of full-term infants. This magnified interest likely stems from the disproportionate health burden placed on preterm infants and its link to the microbiome. The objective of this review was to provide those entering the field, particularly those in neonatal clinical care, with a broad understanding of the different methods used, so that literature interpretation and future study design can be enhanced. This review identifies and describes the most commonly used methods for examining the premature infant’s microbiome, and maps this information against studies comparing efficacy of techniques. This process is designed to illuminate which techniques will be most appropriate for the examination of the microbiome of premature infants
Search and Eligibility Criteria
The PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) approach (Fig. 1), was taken to search for relevant literature up until August 2020. Studies investigating the premature microbiome were identified via searches in the SCOPUS and PubMed databases using the search terms; “Microbiota” AND “Infant” AND “Premature” AND “Faeces”. Journal articles describing a wide variety of study designs, sample sizes, interventions, comparators and outcomes were included in this review. Articles were excluded if they were not original studies, were case studies, were not in English, were unable to be accessed or did not specifically investigate the premature infant microbiome. Reviews found in the initial search were also used to locate other papers that addressed the review question.
Data Collection Process
A standardised data collection protocol was established to extract all relevant information for qualitative analysis. Author, date of publication, aims/hypotheses, a summary of the methods, a summary of the findings and limitations were recorded. Methodology-specific information was also collected for primary techniques, secondary techniques, storage and DNA extraction. Emphasis was placed on 16S rRNA short amplicon sequencing, as this methodology was the most common primary technique, and further information was collected for target variable regions, platform, pipeline and reference databases.
The review of the literature explored the methodological diversity used in the study of the premature infant microbiome, and a summary of the major techniques used across the studies is presented in Fig. 2. The outcome of the systematic review is summarised in Fig. 1. Two hundred and seventeen articles were identified. Of these 137 articles remained after duplicates were removed, with a further 45 articles removed after assessing the full article for eligibility. A total of 92 articles were reviewed. There was a surprising lack of detail displayed in the methods section of many studies, despite there being several technical choices at each stage of the workflow with the potential to contribute to bias. The summary information is based only on the data that was made available.
Thirteen techniques for characterising the preterm infant gut microbiome were identified (Fig. 2). A wide range of techniques have been used as the primary tool for microbial compositional analysis including:
Shotgun metagenomics ,
Long-read nanopore sequencing , and
Random amplified polymorphic DNA/pulsed-field gel electrophoresis (RAPD/PFGE) .
Molecular techniques dominated, specifically 16S rRNA short amplicon sequencing, which made up 53.2% of primary techniques. qPCR was the second most commonly used technique, used in 13% of studies, with the remaining techniques being used by ≤10% of studies each. Traditional non-molecular techniques using selective and differential agar medium represented a very small fraction of the primary techniques used (5.4%), and a further four studies used culture techniques as a secondary method.
Storage Conditions and Extraction Protocol on DNA for Molecular Techniques
Sample storage protocols were consistent across the molecular techniques. However, the DNA extraction techniques used were highly variable. Freezing at -80°C dominated storage methods for both faeces (71.4%) and DNA (50%), with non-freezing protocols only used in 2.9% of studies. DNA extraction was the area demonstrating the greatest variability, with 15 different methods utilised. The QIAamp DNA Stool Kit, a kit that combines heat, chemical and enzymatic lysis was the most commonly used (40.3%), with using no kit at all (14.9%) the second most common option. The PowerLyzer PowerSoil kit was third (10.4%), with the other twelve kits making up the remaining 34.4%.
16S rRNA Amplicon Sequencing – Specific Techniques
The most common molecular technique was 16S rRNA short amplicon sequencing. However, the methods used were highly variable in sequencing platforms, variable target regions, and pipelines. Roche 454 (57.8%) sequencing platform was the most commonly utilised of the four, with Illumina second (31.1%). The use of this platform has increased in recent years. V4 was the most common variable target region used (22.7%). However, there were sixteen unique combinations used across the 92 studies. Of the eight pipelines used, QIIME/QIIME2 (Quantitative Insights Into Microbial Ecology)  made up half of the pipelines used, with Mothur  a distant second (20%). The Ribosomal Database Project (46.3%) was the predominant reference database used, with SILVA (19.5%) and Greengenes (17.1%) being used in most of the remaining studies. Techniques specific to 16S rRNA short amplicon sequencing varied greatly, despite it being the most common method.
Trends Over Time
Many outdated techniques and tools are being abandoned for newer, more robust methods. High throughput molecular techniques have become more commonly used over time, especially in 16S rRNA amplicon sequencing, which was first used in 2004. This upward trend in 16S rRNA amplicon sequencing is coupled with a decline in both fingerprinting- and culture-based techniques, with all culture-based studies occurring prior to 2015. There is also a trend towards the use of Illumina platforms and pipelines that use error modelling within 16S rRNA amplicon sequencing.
Techniques for examining the microbiome can be categorised into two main groups, molecular and non-molecular. Molecular techniques have become dominant due to their depth of analysis, speed and cost reduction. Nevertheless, there are several techniques to choose from even within molecular methodologies, and within a given molecular technique there is variation possible in protocols. This lack of consistency can contribute to the inconsistencies in results between studies into the premature microbiome.
Culture Based Approaches
Few studies still rely solely on traditional non-molecular methods for microbiome characterisation, with the most recent study occurring in 2014. Non-molecular techniques are based on traditional microbiological methods that involve growing microbial communities on predetermined growth media under strict laboratory conditions designed to optimise growth. Methods vary depending on the type of micro-organisms present and downstream applications. Techniques include broth culture, enrichment and microbial identification. Examples of growth media include Luria Broth, also known as Lysogeny Broth , which is common for the cultivation of Escherichia coli, and selective agars such as blood, MacConkey  or Xylose Lysine Deoxycholate agar , which are specific for other taxa prevalent in the gastrointestinal tract. Microorganisms are placed in growth medium and left to grow under strict conditions, giving them time to grow into individual colonies. Colony morphology can then be used to determine specific taxa and colony counts, and are used for the calculation of concentrations and serial dilutions. These techniques are primarily used to identify specific microorganisms of interest due to their specificity, and employed as diagnostic tools for the detection of pathogenic species.
Non-molecular techniques can be useful despite their limitations in sensitivity and specificity, particularly for anaerobic species, as well as for discovery and scaling. They can improve the robustness of results via identification of specific species of interest or identification of unidentified sequences that may belong to a known organism , when used in combination with molecular techniques, such as 16S rRNA high throughputs sequencing. Other major benefits of non-molecular techniques include that the materials are inexpensive and that the protocol requires limited equipment. However, specific culturing conditions that select for specific microbes, of which there must be prior knowledge, mean that many species can go undetected [102, 103]. Moreover, they are time consuming and labour intensive when performed at large scale. Thus, traditional culture techniques have largely been displaced by molecular techniques due to these time and labour issues, as well as these older techniques’ restricted insight into microbial communities.
Molecular Based Approaches
Molecular techniques, including 16S rRNA high throughput sequencing, fingerprinting, microarrays and quantitative PCRs, are rapid, sensitive and highly specific, particularly for commensal organisms. Molecular techniques have rapidly replaced non-molecular techniques for use in identifying microbiome composition since their advent, due to these benefits (Fig. 2). The utilisation of genetic information to differentiate between taxa has made a more detailed exploration possible, and may provide information on the abundance and composition of these microbial communities beyond those routinely grown in the laboratory. The most described microbiota include bacterial communities, which can be identified through utilisation of the variable regions of the 16S ribosomal RNA (16S rRNA) gene, which is flanked by highly conserved regions. DNA is extracted from faeces in this method, commonly using a commercially derived extraction kit. It is amplified by PCR and then differentiated into groups based on similarity to identify the taxa present, allowing deep community sampling. Samples must first be collected and stored and the DNA extracted for all 16S and other molecular techniques.
Sample collection protocols will vary depending on study design. However, the timing of collection is an important factor to consider when comparing results across studies or during study design. Most studies provide specific time points, based on the gestational or post-gestational age of the infants. However, there are studies that group samples together more broadly, for example, binning samples together as “early-infancy” or meconium. Meconium is the earliest stool of a mammalian infant, comprised of a thick tar-like substance that lines the intestine of the unborn infant. Typically, meconium is not released until after birth. However, sometimes it will be released into the amniotic fluid prior to birth. It can also be released at different time points post-delivery, typically within the first three to five days. It may be that accurate comparisons cannot be made across studies when using such broad definitions, as the infant microbiome is dynamic, with choregraphed abrupt changes in composition [34, 104]. Therefore, the timing of collection is an important factor to take into consideration when interpretating the literature and planning future study design.
Impact of Storage Conditions and Extraction Protocol on DNA for Molecular Techniques
Storage conditions influence the stability and constitution of faecal microbial communities , which could prejudice study conclusions. Inadequate storage can lead to continued growth of specific organisms, altering the proportions of taxa/genera in a sample, and can lead to DNA/RNA fragmentation. Studies show DNA/RNA fragmentation can occur after 24 hours when samples are stored at room temperature , and that significant changes in bacterial communities can occur in samples after this time [105,106,107]. Potentially, changes can occur in as little as 30 minutes, as demonstrated by Gorzalek, et al. . Currently, available storage methods include freezing or refrigeration at different temperatures, and the use of anaerobic incubation systems, aqueous storage/transport mediums and faecal occult blood tests.
Freezing and refrigeration
Optimal sample storage conditions depend upon the duration of storage. If samples are to be processed immediately, storage on ice for up to 48 hours , or 4°C for 24 hours  appear to be sufficient for sample preservation. However, immediate freezing of faecal samples to inhibit bacterial growth is the optimal procedure for longer term storage. Long term storage of faecal samples at -80°C has been shown to yield microbiota similar to that of fresh samples [106, 111, 112]. Storage of faecal samples at -20°C also has shown similar efficacy in sample preservation across several studies [105, 113, 114]. However, this appears to be time limited, with some studies reporting changes in taxa over longer-term storage at -20°C for storage times greater than a week, resulting in significant changes to Bacteroides spp.  and for up 53 days in the Firmicutes to Bacteroidetes ratio . Storage at -80°C produces the most consistent results and appears to be the most common preservation method (Fig. 2). It was common for samples to be stored at -20°C or 4°C temporarily, or to use non-freeze methods until storage at lower temperatures was possible in situations where immediate freezing was not possible, such as with at home collection.
Other storage techniques
Immediate freezing of faecal samples can be logistically difficult, especially for large-scale population-based studies, and freeze-thawing effects may significantly diminish sample integrity. Therefore, other preservation methods may be better suited for some protocols. These preservation techniques include chemical and drying preservation such as DNA/RNA Shield and anaerobic incubation systems. Preservation buffers, aqueous reagents that stabilise and protect cellular DNA/RNA, like DNA/RNA Shield (Zymo Research) and RNAlater (Thermofisher), may also preserve genetic integrity for weeks without refrigeration or freezing [117,118,119,120,121,122,123]. Using these buffers ensures that samples are also protected from the potential stress caused by freeze-thawing effects. These buffers and other non-freeze preservation methods are a good alternative when freezing is not feasible.
Some potential issues have been highlighted with non-freezing preservation methods, despite these methods being a more practical alternative. Preservation buffers may result in lower diversity [124, 125] relative to immediate freezing. Moreover, some older preservation buffers may impede downstream DNA extraction and amplification of target genes . Anaerobic incubation systems, like Anaerocult®, are only effective for storage of anaerobic strains and thus have obvious limitations. Most current research still supports the efficacy of preservation buffers, despite several studies highlighting their limitations. Thus, both freezing at -80°C and suspension in a stabilisation buffer are acceptable practices when considering all the available storage options.
The first step for molecular analysis is DNA extraction, which can be carried out using commercially available kits. Extraction is an important step in molecular techniques, which involves separating DNA from the other cellular material contained within samples of interest. The process involves cell lysis, or the disruption of cell walls, separation of the DNA from the other cell components and its subsequent isolation. DNA extraction can be laborious and carries a high risk of sample contamination, as a significant amount of handling of the biological material is involved. Fortunately, there are several commercially available extraction kits that make the process less laborious, more streamlined and more reproducible due to the widespread interest in the human microbiome.
The amount of tissue needed for DNA extraction and sequencing is dependent on the extraction methods and downstream application respectively. Generally, 100 to 1000 nanograms of DNA is required for whole genome sequencing, and as little as 1 to 10 nanograms for amplicon sequencing . The amount of stool required will be dependent on the efficacy of the protocol in extracting the DNA, which is dependent on the methods used for the given protocol. The QIAamp DNA Stool Kit, the most commonly used in the preterm infant microbiome field, is optimised for 190-220 milligrams of stool, but as mentioned above, the amount of stool required is kit-dependent.
Unfortunately, different extraction protocols and kits can contribute greatly to variation in microbial community structure [127, 128], introducing bias to outcomes [115, 129, 130]. In one study the extraction method was demonstrated to be the second-greatest contributing factor to variation . This variation arises in large part due to different methods of homogenisation and lysis. These steps are critical, as different stool fractions can contain different microbial compositions, and different microbes are lysed better by different techniques.
Microbial cell wall structure differs between Gram-negative and Gram-positive bacteria and require different lysis methods for DNA extraction. In Gram-negative bacteria, the cell wall is thin and made up of both a peptidoglycan and phospholipid bilayer containing lipopolysaccharides, whereas Gram-positive bacteria have a thick peptidoglycan cell wall. As a result, Gram-positive microbes require more vigorous lysis methods and Gram-negative microbes are more easily lysed .
Different forms of homogenisation or lysis can, therefore, contribute to bias by not effectively disrupting the cell wall of all microbes present in a sample or, conversely, by destroying the DNA of easily lysed cells. Mechanical, chemical, and enzymatic lysis methods can also produce different proportions of taxa, with mechanical methods producing higher bacterial numbers and greater diversity . Two comprehensive studies that explored several kits, including the QIAamp DNA Stool Kit, found that the International Human Microbiota Standards (IHMS) Protocol Q, that includes mechanical lysis, performs best across several parameters [128, 131]. Despite this, our review found that the most common extraction method utilised in premature infant microbiome studies was the Qiagen QIAamp DNA Stool Kit. This method uses a combination of heat, chemical and enzymatic lysis, with some studies adapting the protocol to add mechanical lysis through bead beating. Unfortunately, different bead-beating instruments have been shown to produce bias [128, 132, 133]. Despite this, mechanical disruption is essential for comprehensive profiling of the human gut microbiome [133, 134], and until a standardised protocol is established, researchers must be careful to consider the bias generated through different kits.
The increasing usage of 16S rRNA amplicon sequencing has been matched with a reduction in the use of other techniques, like different fingerprinting methods. Fingerprinting methods are more cost effective and faster to perform , although high throughput sequencing techniques provide a broad detailed analysis of microbial communities. These techniques are favoured in comparison to traditional culture methods, as they provide greater sensitivity and specificity for individual organisms, and can be used to analyse large numbers of samples . Broadly speaking, fingerprinting methods provide a profile of microbial communities that uses amplification of a target gene (commonly the 16S rRNA gene) and the utilisation of gel electrophoresis to observe physical separation of amplicons, allowing exploration of highly abundant taxa. Fingerprinting methods have been used in studies exploring the preterm microbiome (18.5%), although they are currently less common. These techniques include denaturing/temperature gradient gel electrophoresis (D/TGGE) and terminal restriction fragment length polymorphism (TRFLP), as well as denaturing high performance liquid chromatography (dHPLC) in a single study. dHPLC uses liquid chromatography to identify polymorphisms , while the others rely on electrophoresis to differentiate between sequences, although all are considered fingerprinting methods.
Denaturing/temperature gradient gel electrophoresis
Gradient electrophoresis is the size dependent movement and separation of dispersed nucleic acids through an acrylamide gel. As DNA has a negative charge, it moves through the acrylamide gel or molecular mesh towards the positive electrodes at a rate that is inversely proportional to the size of the nucleic acid sequences, thus allowing the differentiation of different sized sequences . More detailed exploration is achieved by applying either a temperature (TGGE) or chemical gradient (DGGE) to denature the samples as they move across acrylamide gel, based on the chemical make-up of the sequences [138, 139].
Both DGGE and TGGE differ in the mechanism of DNA denaturation. In DGGE the nucleic acids are exposed to increasingly extreme chemical conditions, leading to the denaturation of the DNA in a stepwise process. This allows for visualisation of the sequence differences by their position on the gel. The method relies on differences in the ability to denature the bases, which is determined by base pair sequences to separate genes by size. In contrast, TGGE uses a temperature gradient in combination with the electrophoresis. Strands separate across the gel depending on base-pair content, with smaller molecules travelling faster  as the temperature increases.
Terminal restriction fragment length polymorphism (TRFLP)
TRFLP also uses electrophoresis to differentiate between sequences based on terminal restriction fragment size, like TGGE and DGGE. This allows sequence identification for microbial community profiling . The method involves PCR amplification of a target gene with fluorescently labelled primers and subsequent digestion with restriction enzymes. The sizes of the different terminal fragments are then determined by separating the fluorescently tagged terminal fragments via capillary or polyacrylamide electrophoresis in a sequencing gel, creating unique banding patterns allowing identification of microorganisms . T-RFLP has high throughput capability and can be highly sensitive, but it also has limited accuracy as incomplete or non-specific digestion can lead to overestimation of diversity. Homology of sequences can contribute to an underestimation of taxa present . Furthermore, libraries must be built prior to analysis.
dHPLC uses liquid chromatography to identify DNA polymorphisms , unlike TGGE, DGGE and T-RFLP. DNA strands are separated into hetero- and homoduplexes using an ion-pair, reverse-phase liquid chromatography on a poly alkyl column matrix ., following partial heat denaturation. The presence of polymorphisms is revealed by the differential retention of these homo- and heteroduplex DNA fragments . Heteroduplexes are double stranded DNA that have formed during PCR amplification that are mismatched at the site of mutation. Mismatched double stranded DNA fragments have reduced retention on the column matrix, and subsequently in a reduced retention time, thus allowing for identification of polymorphisms. As a result, dHPLC can be useful as a screening test for mutations that may be involved in diseases or associated with antibiotic resistance [136, 145]. The method has also been used to differentiate between taxa at species depth by applying the same underlying principle of scanning for mutations to the detection of sequence variations between PCR-amplified bacterial 16S rRNA genes , and also as a tool for re-sequencing of genomes .
Limitations (Pros and Cons)
Some consider dHPLC to be the optimal fingerprinting method , potentially allowing identification of bacteria at the species and/or biotype levels . However, these techniques require extensive downstream processing, can produce PCR bias [102, 147] and have limited detection depth, as it is difficult to relate banding patterns created in gels to species or lineages created by fingerprinting methods . Thus, fingerprinting methods are usually limited to identification at the order/family level , and to only the most abundant organisms. This methodology also makes it difficult to combine data from multiple studies into a single analysis . Fingerprinting techniques can be useful for exploring dominant members of microbial communities, including clustering of communities based on dominant members . However, their application is limited in describing entire microbial communities.
Microarrays were originally developed to monitor gene expression, but their application has been expanded to include comparative genomics, DNA sequencing analysis, single-nucleotide polymorphism (SNP) analysis and microbial detection , including studies on the premature infant microbiome [77, 78]. Microarrays are microscopic slides printed with probes made of predefined oligonucleotide sequences complementary to the small subunit (SSU) rRNA. The oligonucleotide probes detect gene expression or mRNA transcripts expressed by specific genes and extracted from target organisms. Reverse transcriptase converts mRNA into complementary DNA (cDNA), and this cDNA is fragmented and fluorescently labelled and added to the microarray . cDNA then binds complementary oligonucleotide probes via hybridisation, and measurement of the observed fluorescent intensity at a given probe is an indication of the abundance of predetermined sequences that are chosen prior to analysis and are of interest . This makes phylogenetic oligonucleotide arrays (phyloarrays), including HITChip, suited to the analysis of microbial communities.
HITChip is an ecosystem-specific phylogenetic microarray developed for microbial detection in the human gastrointestinal tract [152, 153], and is the only microarray to be used in studies on the premature infant gut microbiota. HITChip is an oligonucleotide microarray that uses 4,800 oligonucleotide probes based on two hypervariable regions (V1 and V6) of the 16S rRNA gene, identifying 1,140 phylotypes. Phylotypes were designed following the analysis of 16,000 human gastrointestinal tract 16S rRNA gene sequences . As a result, HITChip is highly specific to the human gastro-intestinal tract microbiome and provides a high level of diversity.
Limitations (Pros and Cons)
Benefits of microarrays include ease of use, speed and cost , and potential for investigating microbial gene functionality . These intermediate methodologies allow processing of large sample sizes, while providing more taxonomic depth, like fingerprinting methods. Microarrays target the ribosomal RNA gene, allowing comparisons of diversity and taxonomy, and thus display similar robustness [156, 157], like 16S rRNA sequencing. However, when compared to high throughput techniques, phylogenetic arrays are limited when assessing new lineages, as they can only detect predefined taxa [135, 158]. Other methodologies are better suited when there is potential for taxonomic or gene discovery, as microarrays are limited to predefined taxa. Microarrays are not commonly used in studies on the premature infant gut microbiome, as it is a relatively new area of study, and other methodologies may be better suited for characterising this niche.
qPCR and Fluorescent In-situ hybridization
Polymerase chain reaction (PCR) is a highly sensitive molecular technique that was originally developed for detection of DNA/RNA sequences, but has since progressed beyond purely nucleic acid detection. Quantitative polymerase chain reaction (qPCR) builds on standard PCR by providing the quantity of amplified genes. qPCR also differs from standard methods as it monitors the amplification of targeted DNA molecules in real time or during PCR instead of at the end. This process allows not only detection, but also quantification and characterisation of nucleic acids . Fluorescent dye is added to the PCR reaction in dye-based qPCR, and the fluorescent signal increases proportionately to the quantity of DNA being replicated. This allows quantification of DNA after each cycle. However, qPCR only allows one target to be examined at a time, thus throughput is limited. The more accurate probe-based qPCR provides one way around this drawback, by simultaneously examining multiple targets via recognition of sequence-specific probes. The fluorescent signal from the probe in probe-based qPCR is proportional to the target sequence that is present in the reaction , as it is in dye-based qPCR.
Fluorescent in situ hybridisation (FISH) is another probe-based technique. FISH is a molecular technique that uses complimentary binding to identify or quantify cDNA that can be used for microbial identification, like microarrays and qPCR . FISH uses fluorescently labelled DNA probes that match specific DNA sequences that can be observed under a microscope, allowing direct quantification of specific taxa. Fluorescent oligonucleotide probes are created for targets, either 16S or 23S rRNA sequences. The target and probe sequences are denatured with heat or chemicals, and mixed together prior to hybridisation. Hybridisation then occurs between complementary target and probe sequences, with fluorescence microscopy facilitating detection of hybridisation via observations of fluorescently labelled cDNA. This target-specific methodology facilitates high accuracy when targeting specific microbes.
Limitations (Pros and Cons)
Both qPCR and in-situ hybridisation can provide highly accurate quantification , can be highly sensitive , and can produce similar results to metagenomic methods when considering the main intestinal microbial groups . However, they are limited in their application, as prior knowledge of sequences is required, like fingerprinting and microarrays. Thus, these methods have no discovery power and no capacity for assessing diversity. FISH has been designed to examine the major microbial groups present in premature infants, but is based on groups present in full-term infants [79, 80]. However, predefining taxa in this way is a significant limitation, as premature infants are known to have significantly different microbial populations to infants born full-term [24, 26]. Moreover, both qPCR and FISH are not scalable, and therefore are only effective for low target numbers. qPCR or in-situ hybridisation methods may be beneficial when specific populations are being targeted, as they have limited bias and are cheaper compared to sequencing methods, but they are not suitable for projects mapping entire microbial ecosystems, like that of the premature infant gut microbiome.
DNA sequencing is the process of nucleic acid sequence determination, and covers a broad range of techniques across three generations of sequencing. The first generation of sequencing began with a low throughput technique, Sanger sequencing, which only sequenced a single DNA fragment at a time . Sanger uses a labour-intensive cell-based amplification step, involving cloned sequences being placed into plasmids for amplification, prior to extraction and purification . The second generation of sequencing techniques, often referred to as next-generation sequencing (NGS) or high-throughput sequencing (HTS), involved 16S rRNA Metabarcoding and Metagenomics (shotgun sequencing). NGS refers to any sequencing method using the concept of parallel processing. This parallel processing increased the volume of reads per run to millions, vastly improving efficiency, as did the development of a cell-free system. NGS also runs elongation and detection steps in parallel, again improving efficiency . However, NGS technologies are limited in that they use short reads (bp), which create a computational challenge when assembling or mapping to genomes. A third generation of sequencing was developed to overcome this challenge, long read or single molecule direct sequencing. The capacity of all sequencing technologies to produce large volumes of relatively accurate data, coupled with the continual reduction in cost, has led to their adoption across most modern studies investigating microbial populations. 16S rRNA high-throughput amplicon sequencing (metabarcoding) is now the most common method used for studies specifically characterising the preterm infant’s gut microbiome (Fig. 2). However, other methods, including shotgun metagenomics, and third generation single molecule direct (long-read sequencing), have also been applied. All techniques described have their strengths and limitations.
Next Generation Sequencing
16S rRNA amplicon sequencing
16S rRNA amplicon sequencing, or metabarcoding, has become the most common technique for characterisation of the preterm infant microbiome since it was first used in 2004. 16S rRNA metabarcoding uses high throughput sequencing to target variable regions of the 16S rRNA gene, allowing accurate identification of microbial community composition [166,167,168]. The 16S rRNA gene codes for 16S ribosomal RNA, a component of the 30S small subunit of prokaryotic ribosomes. The 16S rRNA gene is highly conserved across taxa, but also has several variable regions allowing differentiation between taxa, due to a slow rate of evolution. The variable regions are conserved enough that most taxa can be characterised, but variable enough that taxa can be differentiated. There are nine of these hypervariable regions that range in base pair length and are involved in the secondary structure of the small ribosomal subunit. The regions vary in conservation, and thus different regions correlate with different levels of taxonomic resolution. The protocol for 16S rRNA gene amplicon sequencing involves DNA extraction, PCR amplification of the variable target region(s), grouping of sequences into OTUs, ASVs or an equivalent, and then mapping these sequence variants to a reference database for taxonomic identification.
16S rRNA metabarcoding is the most common technique for characterising the preterm infant. However, despite this there is no predominant protocol. There are a myriad of options at all stages of the workflow, all of which can introduce bias that alters outputs, which is supported by the observation that samples cluster by study . Technical differences in how samples are collected and stored, how DNA is extracted, the primers that are selected and variable regions targeted, the sequencing platform, bioinformatic pipelines and reference databases could all produce systemic bias that obscures biological differences [15, 16]. As the 16S rRNA protocol is the most common technique for characterising the preterm infant gut microbiome, a more detailed explanation of its varied protocols is discussed below. Caution should be used during both interpretation of the literature and study design until a standardised protocol is agreed upon.
Selection of Variable Regions and Primer Bias
Once DNA is extracted, and prior to sequencing, the target DNA from the variable region of interest must be amplified via PCR. However, there is much debate on which variable sub-region to target and matching primers to use. Indexing primers are complementary base pair sequences that are required to ‘select’ and amplify variable sub-regions. These 16S rRNA variable sub-regions can vary by up to 40% in taxa between samples analysed with the same pipeline , and it has been argued that the most critical step for accurate amplicon analysis is the choice of primers , as primer selection can alter coverage .
Samples with the same extraction and storage protocols have been demonstrated to cluster by primer selection . This is because poor primer selection can influence quantitative abundances , and contribute to under-representation or over-representation of taxa [173, 174] or selection against particular taxa [135, 175, 176]. For example, most primers may inadequately detect Bifidobacterium , possibly over-exaggerating the low levels already observed in preterm infants. For identifying species, targeting of these hypervariable sub-regions is limiting, as different sub-regions show bias in the taxa that they can identify due to limited variability within the sub-region itself. Thus, while V1-V3 may be good for Escherichia and Shigella species, Klebsiella will require V3-V5, and Clostridium and Staphylococcus require V6-V9 sequencing . As a result, studies that target hypervariable sub-regions must settle for taxonomic resolution at the genus level. Thus, the only way to ensure good taxonomic identification would be to sequence the entire 16S gene, given these limitations and the bias that can be introduced through variable region and primer selection. However, high error rates and cost are still major deterrents.
Arguments have been made for targeting the V4-V6 , V4 [172, 179, 180], and V3-V4 [169, 171, 181] sub-regions, with V4 being the most common for characterising the microbiome of preterm infants, when targeting hypervariable sub-regions for high throughput sequencing (Fig. 2). The Earth Microbiome Project  recommends the V4 sub-region, and it has been demonstrated to have low PCR and sequencing errors due to complete overlap of paired end sequences . Other work has also shown that the phylogenetic relationships based on V4 were closest to that entire 16S rRNA gene . However, some evidence suggests that targeting the V4 region may not be as accurate as previously thought . The debate about which region is best is ongoing, but research conducted by Almeida et al. makes a convincing argument for the use of the V3-V4 region above all else. It compared variable regions across different combinations of pipelines and reference databases for both mock communities and simulations , and found the V3-V4 region consistently produced the most reliable taxonomic inferences. Taken together, the frequent use of V4 and the findings across studies, targeting either the V4 or V3-V4 sub regions may be best practice until standardisation occurs. Additionally, consideration should be given when making comparisons across studies that use different variable regions.
Several sequencing platforms are available for 16S rRNA short read amplicon sequencing, with Illumina MiSeq, Roche 454 (originally 454 Life Sciences) and Thermo Fisher’s Ion Torrent Personal Genome Machine (PGM) all being used in the context of the premature infant microbiome. Roche 454 has historically been the dominant platform, as NGS technologies began with it. However, Illumina now dominates the market, with its consistent growth and the eventual abandonment of the Roche 454 sequencing platform in 2016. It is important to understand the differing methods across platforms, and their limitations and biases for accurate interpretation of the literature, although most modern sequencing technologies will opt for Illumina sequencing.
Illumina sequencing technology facilitates massively parallel sequencing by using optical signals to detect base pairs in real time. DNA libraries, containing fragments that vary between 100-150bp, are loaded onto a flow cell and placed in the sequencer for this process. The sequences bind to the flow cell via complementary adaptors. A process called clonal bridge amplification or cluster generation then amplifies each read, creating a spot (cluster) on the flow cell (slide) with thousands of copies of the same DNA strand. Then, through a process coined sequencing by synthesis, fluorescently tagged nucleotides bind to the complementary bases on the DNA strand via repeated cycles of single-base extension. A fluorescent signal (the colour of which is dependent on the base) is emitted upon incorporation of each nucleotide, and a picture taken, indicating what nucleotide was added. Once the forward DNA strand is read, the reads are washed away, and the process is repeated for the reverse strand. Computers then construct the sequence by detecting the base at each site in each image.
Roche 454 relies on the production of sequence clusters, like Illumina, but through a process called clonal emulsion PCR (emPCR). In emPCR, single stranded DNA fragments (up to 1kb) from a DNA library are attached to the surface of a bead, rather than a slide, with one bead for each DNA fragment. The reads bind to the bead via complementary adaptors. The beads are then compartmentalised into single wells containing emulsified oil, and are subjected to thermal cycling to achieve clonal amplification. This process produces many copies of the original template, as in Illumina’s clonal bridge amplification. The slide containing the wells is then flooded with one of the four nucleoside triphosphates (NTP) that bind to their complements, releasing a light signal upon addition. The original NTP mix is washed away and the next NTP is added and the cycle is repeated. The light intensities are then plotted on a graph for each sequence read, with graphs then used to determine the sequence computationally.
Ion Torrents PGM also uses clonal-emPCR, but differs from Roche 454 both in how it determines the nucleotide sequences and the size of DNA fragments. Ion Torrents PGM or proton sequencing uses DNA fragments of ~200bp, which are again bound to beads via adaptors. These then undergo PCR and are washed with different NTPs. It then exploits the release of hydrogen ions, which occurs through the addition of an NTP to a DNA polymer for nucleotide sequence determination. The release of hydrogen ions causes changes in pH that are used to determine the DNA sequences.
Limitations (Pros and Cons)
No platform is without its limitations: limitations that can contribute to platform-associated biases and study-based clustering , despite significant developments in the sequencing field. For example, Roche can have high sequencing error rates associated with A and T bases , high error rates in homopolymer regions resulting from accumulated variance in light intensity [184,185,186], and can have up to 15% of sequences resulting from artificial amplification . Ion Torrent is also subject to high homopolymer error rates [186, 188, 189], as well as organism-specific read truncation, due to the similar methods of Roche and Ion Torrent, in which multiple nucleotides can be incorporated during a single cycle . Illumina still have their own systematic base-calling biases , even though platforms produce comparatively lower error rates . These include production of homopolymer-associated sequencing errors , different quality reads across different sequencing tiles , increased single-base errors associated with GGC motifs  and different sequencing error rates at the different read ends . In the two dominant platforms in preterm infant studies Roche 454 and Illumina, the differences caused by platform are minor [172, 183]. However, the lower error rates, higher throughput  and higher read quality  achieved by Illumina, results in higher quality data. This allows stringent quality control parameters, resulting in more reliable outputs for downstream analyses .
Bioinformatics and Reference Databases
Bioinformatics is an interdisciplinary field of science that combines biology, computer science and statistics, in order to process large amounts of biological data, such as that produced by 16S rRNA gene amplicon sequencing. Bioinformatic tools like QIIME  and Mothur  are required to clean up and make inferences on microbial composition from data that is not human-readable post sequencing and prior to downstream analysis. Bioinformatic tools or pipelines need to be both precise and reliable in order to produce accurate biological conclusions using the vast amounts of genetic data that is being produced with sequencing. These tools convert raw data into interpretable taxonomic abundances by comparing sequencing reads in the form of OTUs, ASVs or an equivalent (sequence variants that represent a true sequence)  to a defined reference database, identifying the taxa present in samples by assigning the most likely taxonomic lineages. The accuracy of the taxonomy classifications produced is then reliant on both the diversity and breadth of annotated sequences in the reference databases , as well as the accuracy of the ever improving algorithms used by bioinformatic pipelines.
There is no agreement on optimal practices, although the bioinformatic pipelines are rapidly changing and improving, and many researchers are unaware of the biases associated with using different tools. Combinations of different software packages, databases and targeted regions can produce vastly different levels of accuracy when examining mock communities and running simulations . When comparing several bioinformatic tools: QIIME, QIIME 2, Mothur and MAPseq , Almeida et al. found that QIIME 2 was the optimal tool in regards to detection sensitivity and composition prediction. QIIME 2 had the largest proportion of classified sequences at the most accurate relative abundances . However, MAPseq was more precise, with fewer genera being miss-assigned. A more recent study which compared the most popular current bioinformatic pipelines for 16S rRNA gene amplicon sequencing, found that DADA2 was the best choice for studies requiring the highest possible biological resolution, but that USEARCH-UNOISE3  had the best overall performance . The common theme running through USEARCH-UNOISE3, DADA2 and QIIME 2 (uses DADA2/deblur plugins) is their denoising or clustering algorithm.
Denoising and clustering are methods for correcting sequencing errors through grouping of similar sequence variants into a bin. This was originally done through OTU clustering, in which sequences are clustered based on a 97% similarity threshold. However, there are several methods for implementing this threshold: closed-reference, open-reference and de novo clustering . The de-novo method clusters reads against one another, based on the threshold, without a reference database, unlike the reference-based approaches. In contrast, the closed-reference method clusters reads against a database and excludes those sequences that do not align. Open-reference clustering also clusters against a database, but then clusters reads that do not align de novo. The most successful method is debatable [200, 201], but may be dependent on the study design. Nonetheless, the quest for more reliable data has seen a shift away from OTU-clustering towards error modelling, which takes into account both abundance and error.
Denoisers, as seen in DADA2  and deblur , generate error models learnt from the reads and use these models for sequence variant assignment with either ASVs (DADA2) or subOTUs (deblur). The error modelling approach allows for clustering down to the level of single-nucleotide differences in the sequence region, improving resolution, and allows consistently reproducible labels with intrinsic biological meaning . These improvements allow researchers to distinguish between true sequences and those generated during PCR amplification and sequences, and result in comparable general community structure across different tools. However, some variability still exists, with differences in the number sequence variant produced and resulting alpha diversity , despite these improvements. These differences should be considered when making cross-study comparisons and during study design.
The most common bioinformatic pipeline for studies exploring the preterm infant gut microbiome is QIIME. QIIME’s use of OTU clustering and its production of a large number of spurious OTUs and inflated alpha diversity  should be taken into account when considering older literature. However, QIIME was succeeded by QIIME2 in 2018 (first published in 2019) , which uses an updated error modelling approach, with either DADA2 or deblur plugins. Most studies on the preterm infant microbiome predate the release of QIIME2, and therefore will have used the older version. It is unclear if newer studies will make the transition, due to the limited number of papers released since the pipeline’s publication. However, in order to produce more robust data, research in this field needs to move towards these new and improved methods.
Along with choosing the best bioinformatics tool, the reference database used is also an important consideration. Pipelines can use a homology based or Bayesian approach to match sequence variants to sequences from the reference databases. This was originally achieved with a similarity threshold of >95% sequence match being considered to represent the same genus and >97% match for species level identification . However, recent work suggests that these thresholds are too low for accurate assignment . The reference databases contain FASTA files with reference sequences assigned to nodes of taxonomy. However, discrepancies exist in both nomenclature and lineages of taxonomy between databases [206, 207], which has obvious implications for the taxa identified in a sample. The Ribosomal Database Project (RDP) is the database used most often for studies on the preterm infant microbiome. SILVA may have better recall and be more precise than the more commonly used RDP (Fig. 2), when examining the human microbiome, based on a benchmarking paper by Almeida et al. . However, more research on best practices, along with standardisation across databases is needed.
Bioinformatics – Contamination
Contamination is another way variability between studies can be introduced. DNA contamination from collection, extraction and sequencing protocols (and kits) can impact upon the interpretation of results [208, 209]. Many of these contaminants could be considered normal inhabitants of the human gastro-intestinal tract, and so it is important that measures are taken to mitigate the risk of contamination, and that steps are taken to account for or remove this contamination. Negative controls, spike-in controls and microbial standards, when coupled with appropriate bioinformatic tools, are effective ways to account for both accuracy of techniques and contamination. Bioinformatic tools like microdecon  can be used to remove homogenous contamination. This is important for both study design and interpretation of the literature, as contamination may produce unusual or novel findings if the appropriate mitigation strategies are not in place .
Analysing microbial data can provide a significant challenge due to the volume and complexity of the data. Additionally, it is difficult to provide a best practice approach for statistical analysis of microbiome data, as it is highly dependent on the hypotheses and objectives of the study. However, generally, there are three main metrics that are considered in microbial analysis: alpha diversity, beta diversity and differential abundance. Unfortunately, finding a meaningful way to conduct these analyses can be a convoluted process.
Alpha diversity refers to the diversity within a sample, summarising the ecological structure with respect to either richness (number of taxonomic groups), evenness (distribution of abundances of the groups), or a combination of the two . Alpha diversity can be represented by richness, the Chao 1 index, the Shannon-Weiner index, the Simpson index, Pielou’s evenness or Faith’s phylogenetic diversity. All these indices differ in what they represent, with richness and the Shannon-Weiner index common in the context of microbiome research. Richness is simply referring to the count of sequence variants, whereas the Shannon-Weiner index takes both richness and evenness into account. The Chao 1 index is a bias-corrected richness estimate  that has become less common with the advent of newer pipelines (e.g. DADA2), due to their handling of singletons, which is based on the assumption that many are spurious sequencing variants.
Beta diversity refers to microbial differences between samples or groups. There are several metrics for beta diversity that represent distances between samples, as there are in alpha diversity. Beta diversity metrics include the Jaccard distance, Bray-Curtis distances and UniFrac. Jaccard distance is the number of sequence variants shared by samples divided by the number shared. Bray-Curtis distances, one of the most widely used in microbiome research, builds on this by taking abundances into account as well, whereas UniFrac distances, either weighted or unweighted, represent the differences between samples based on phylogenetic differences. These distance matrices can all be represented with direct comparisons of distances between samples or groups, hierarchical clustering or ordination techniques. Ordination is the most common as it reduces the complex distance data to a 2D or 3D plot, making for easy interpretation.
Alpha and beta diversity, along with differential abundance, all require normalisation prior to analysis. Differential abundance is simply making comparisons for taxonomic abundance between samples or metadata. However, the process is more complicated than simply counting the number of reads per sample. This is because using read counts as a measure of abundance is flawed, as the number of reads is actually an artefact of sequencing, and is therefore not a good representation of abundance . Alpha and beta diversity also take read counts into account, and so both are also sensitive to sequencing depth. Both diversity metrics require an equal number of reads per sample for valid analysis. This means that if appropriate measures are not taken, library sizes can determine diversity results . Thus, reads must be normalized to account for the differing number of reads per a sample prior to analysis. However, although normalisation is the solution, different methods are required for different analyses.
Older methods of normalisation include Total Sum Scaling/Normalisation (TSS) and rarefying. In TSS, data is transformed to proportions by dividing the reads for each sequence variant by the total number of reads, whereas rarefying adjusts for differences in library sizes by assigning a sequencing depth threshold, and subsequently subsampling samples with a depth above the threshold and discarding those below. However, both methods are poor options for differential abundance testing and can have high type 1 errors [215,216,217,218]. Additionally, TSS doesn’t account for heteroskedasticity  and rarefying discards potentially useful data. As a result, other modern methods have begun to replace the old methods.
Newer methods include variance stabilizing transformation with DESeq2 , upper quantile normalisation , CSS normalisation  and Trimmed Means of M-values (TMM) with EdgeR . Methods like DESeq2 and EdgeR are generally favoured due to their performance across several comparative papers [215,216,217]. However, these comparisons are specific to standardisation of within-sample variance, the ability of data to cluster in ordinations, and their performance in differential abundance testing, but there are several important limitations . These methods tend to focus on standardising within-sample variance across samples, as they were created for differential abundance testing. As a result, the newer methods do not guarantee equal number of reads across samples. They supress species evenness and overestimate the importance of low abundance taxa (through log transformations) . The overestimation of low abundance taxa and suppression of evenness can contribute to an inaccurate representation of the community, and, along with non-equal read counts, can lead to inaccurate comparisons between samples.
So, although proportions, specifically TSS and rarefying, are not suitable for differential abundance testing, they are more suitable for diversity analysis, as they give a more accurate representation of the microbial communities while accounting for differences in read depths [218, 222]. Additionally, methods like the variance stabilising transformation with DESeq2 are favoured for differential abundance testing. It is critical that researchers are aware of the strengths and limitations of normalisation methods for accurate interpretation and robust study design and differences between the subsequent analyses.
Metagenomic Shotgun Sequencing (whole genome sequencing)
Shotgun metagenomics is another NGS approach that has been used a handful of times for characterising the gut microbiome of premature infants. Shotgun metagenomic sequencing targets all DNA in a sample, in contrast to 16S rRNA sequencing, which targets a specific region/gene. The protocols differ from 16S rRNA sequencing slightly, although they use the same sequencing technology. Shotgun metagenomics does not require amplification, as there is no target region/gene, but it does require the removal of host DNA prior to mapping, as all extracted DNA in a sample is sequenced. This alternative NGS method provides greater taxonomic resolution and gene annotation, allowing more comprehensive analyses. As a result, studies using this technology are typically looking to link functional gene profiles or pathogen strains to disease.
Limitations (Pros and Cons)
Shotgun sequencing has its limitations, despite its obvious benefits. There are numerous experimental and computational approaches that can be carried out at each step, as in 16S rRNA sequencing . DNA extraction methods have been shown to affect composition , due to kits and reagents containing microbes  and differences in lysis techniques . Library preparation and sequencing can introduce errors through PCR amplification  and selection of platforms [227, 228]. Furthermore, specifically metagenome profiling can cause protocol-associated variability, having several options for bioinformatics, as all metagenomic profiling techniques have their own limitations .
There are two approaches for metagenome-profiling: assembly-free methods and assembly-based methods. Assembly-free methods, also known as read-based profiling or ‘mapping’, make comparisons to reference databases that contain whole genomes, such as Kraken  or Centrifuge , or to selected marker genes, such as mOTU . Alternatively, assembly-based analysis uses assemblers like Meta-IDBA  and SOAPdenonovo2  to reconstruct genomes de novo. Assembly-based methods can construct multiple whole genomes and resolve novel organisms, but can be a significant computational burden and are limited in assessing complex communities. Alternatively, read-based analysis is computationally efficient and can deal with more complex communities, assuming there are enough sequencing depth and genomes in the reference database. However, identification is limited to those microbes previously defined, and so community structure/function is limited. Both approaches have their strengths and weaknesses, and which is best may depend on the question being asked.
There are also pros and cons when comparing NGS shotgun approaches to other sequencing methods. Metagenomics has more reliable species identification and broader analyses potential relative to metabarcoding, but the bioinformatics is more involved, requiring more time, skill and computational power, and the sequencing is more expensive, as entire genomes are being sequenced instead of a single gene. As a result, older studies using shotgun approaches tend to have lower sample sizes  or only use the technique on a subset of the cohort . Additionally, where fragments of bacterial genomes are mixed in with contamination from host species and other organisms, 16S rRNA amplicon sequencing has specificity for bacteria, does not require full reference genomes, and does not require large quantities of, nor high quality, DNA [234, 235], as opposed to shotgun sequencing. However, the adoption of this technique will likely become more widespread as the price of shotgun metagenomic sequencing continues to drop, in combination with improved computational methods. However, 16S rRNA amplicon sequencing targeting variable regions continues to dominate studies in this field at the present time.
Third Generation Sequencing
Long read sequencing is another sequencing approach. Full length sequencing of the 16S gene was made possible with the advent of third generation sequencing technology, also known as long-read sequencing. This approach is possible with platforms like Oxford Nanopore Technologies Minion  and techniques like Pac Bio’s Circular Consensus and Continuous Long Read Sequencing . These technologies allow discrimination between millions of reads that may only differ by a single nucleotide , and have the capacity to produce reads in excess of 10,000 base pairs (bp) [238, 239]. This allows the sequencing of the entire 1,500 bp 16S gene and increases the resolution in taxonomic profiling to species and strain level.
Third generation sequencing can produce these long reads because their design is distinct from previous sequencing methods. Nanopore technology produces long sequences by passing a single DNA molecule through a DNA pore, measuring changes in current across a membrane. The current passing through the membrane is dictated by the size of the base pairs in the sequence that is passed through the pore. Alternatively, PacBio’s SMRT (Single Molecule, Real-Time) sequencing repeatedly passes a DNA molecule through a DNA polymerase attached to a well, with short sequences being read until there are enough overlapping reads to identify the entire sequence.
Long-read technologies can also be applied to whole genome sequencing (WGS) and shotgun metagenomics, as well as full length sequencing of the 16S rRNA gene. WGS or shotgun metagenomic approaches allow greater sequencing depth, meaning species level detection of the preterm microbiome can be achieved, like sequencing the entire 16S gene. This capacity was demonstrated by Legget et al., who took advantage of Oxford Nanopore’s ability to produce near-real time data in developing a metagenomic screening platform for preterm infant microbiome samples . Moreover, as shotgun metagenomic sequencing targets all genomic DNA in a sample, the data can be used for other analyses, like functional profiling and antibiotic resistance gene profiling. This provides a comprehensive investigation of microbial ecology.
Limitations (Pros and Cons)
The long reads produced from the two technologies are their major advantages. Nanopore can produce reads generally ranging from 10Kbp to 1Mbp, with the longest sequence produced being >2Mbp . These longer reads, along with advances in the associated computational methods, allow for greater sequencing depth than short-read technologies, with potential for greater accuracy , as they can distinguish between sequencing artifacts and actual biological sequences . However, high error rates [242,243,244] are still a problem in TGS, despite the claim of high accuracy. For Oxford Nanopore, this high error rate comes from using changes in current to identify base pairs . In PacBio’s Single Molecule, Real-Time sequencing no current technology can precisely capture the rate of information produced (DNA polymerase adds 100bp/s), , which is one reason why the DNA must be passed through the enzyme multiple times to overcome this issue. So, although these technologies show promise, the high error rates and cost are still deterrents, which is probably why they have been used so seldom in studies on the preterm infant microbiome.
Variability of results will continue to be a limitation when investigating microbial populations in premature infants until there is standardisation of protocols. This review aimed to describe the most common techniques used in researching the preterm infant microbiome, and their limitations. The objective was to provide those entering the field with a broad understanding, so that considerations can be taken for both literature interpretation and future study design. 16S rRNA amplicon sequencing is the most commonly used method, as it is cheaper than both long-read and shotgun metagenomic sequencing, more detailed than non-molecular techniques and allows the characterisation of taxa present across a wide range of samples. This approach, however, has several limitations that can introduce bias. Full length sequencing of the 16S gene or a shotgun metagenomics approach may provide better options, especially as accuracy continues to increase, along with a reduction in cost. However, until these options become more viable, 16S high throughput sequencing targeting a select number of hyper variable sub-regions will continue to dominate.
There are a number of options at different stages within 16S sequencing methods that can contribute to bias, and with the large number of tools and databases available, it can be a difficult task deciding on an optimal approach. In this work we briefly described the bias across methodologies, with emphasis on 16S techniques. The most commonly used techniques within 16S rRNA high throughput sequencing are sample storage at -80°C, QIAamp DNA Stool Kit for extraction, sequencing on the Roche 454 platform, targeting the V4 region, and using the QIIME or QIIME2 pipeline in combination with the Ribosomal Database Project reference database. However, the optimal combination for 16SrRNA sequencing would likely be storage at -80°C, an extraction kit that includes mechanical lysis, such as (IHMS) Protocol Q or the Power Faecal Pro (Qiagen), use of the Illumina platform, targeting of the V3/V4 regions, using the QIIME2 pipeline (or at least error modelling) in combination with the SILVA database. However, the research question, as well as reproducibility and consistency across studies should also be considered. To conclude, until standardisation of microbiome research occurs, significant consideration needs to be given to ensure correct interpretation of the literature and robust study design.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
16S terminal restriction fragment polymorphism
Amplicon sequence variant
Denaturing gradient gel electrophoresis
Denaturing high-performance liquid chromatography
Emulsion polymerase chain reaction
Fluorescence in situ hybridization
- InviMag Kit:
InviMag Stool DNA Kit
- MoBio Kit:
MoBio Powersoil Bacterial DNA Kit
Operational taxonomic unit
- PacBio SMRT:
Pac Bio Single Molecule, Real-Time Sequencing
Polymerase chain reaction
Pulse-field gel electrophoresis
- PowerLyzer Kit:
PowerLyzer PowerSoil Kit
- QIAamp Kit:
QIAamp DNA Stool Kit
Quantitative Insights Into Microbial Ecology
Quantitative polymerase chain reaction
Random amplified polymorphic DNA
Ribosomal Database Project
Temperature gradient gel electrophoresis
Third generations sequencing
Whole genome sequencing
Itani T, Ayoub Moubareck C, Melki I, Rousseau C, Mangin I, Butel MJ, et al. Establishment and development of the intestinal microbiota of preterm infants in a Lebanese tertiary hospital. Anaerobe. 2017;43:4–14.
Raveh-Sadka T, Firek B, Sharon I, Baker R, Brown CT, Thomas BC, et al. Evidence for persistent and shared bacterial strains against a background of largely unique gut colonization in hospitalized premature infants. ISME J. 2016;10(12):2817–30.
Barrett E, Kerr C, Murphy K, O'Sullivan O, Ryan CA, Dempsey EM, et al. The individual-specific and diverse nature of the preterm infant microbiota. Arch Dis Child Fetal Neonatal Ed. 2013;98(4):F334–40.
Ward DV, Scholz M, Zolfo M, Taft DH, Schibler KR, Tett A, et al. Metagenomic Sequencing with Strain-Level Resolution Implicates Uropathogenic E. coli in Necrotizing Enterocolitis and Mortality in Preterm Infants. Cell Rep. 2016;14(12):2912–24.
Chernikova DA, Koestler DC, Hoen AG, Housman ML, Hibberd PL, Moore JH, et al. Fetal exposures and perinatal influences on the stool microbiota of premature infants. J Matern Fetal Neonatal Med. 2016;29(1):99–105.
Hill CJ, Lynch DB, Murphy K, Ulaszewska M, Jeffery IB, O'Shea CA, et al. Evolution of gut microbiota composition from birth to 24 weeks in the INFANTMET Cohort. Microbiome. 2017;5(1):4.
Patel AL, Mutlu EA, Sun Y, Koenig L, Green S, Jakubowicz A, et al. Longitudinal Survey of Microbiota in Hospitalized Preterm Very-Low-Birth-Weight Infants. J Pediatr Gastroenterol Nutr. 2016;62(2):292–303.
Schwiertz A, Gruhl B, Lobnitz M, Michel P, Radke M, Blaut M. Development of the intestinal bacterial composition in hospitalized preterm infants in comparison with breast-fed, full-term infants. Pediatr Res. 2003;54(3):393–9.
Arboleya S, Binetti A, Salazar N, Fernandez N, Solis G, Hernandez-Barranco A, et al. Establishment and development of intestinal microbiota in preterm neonates. FEMS Microbiol Ecol. 2012;79(3):763–72.
Itani T, Ayoub Moubareck C, Melki I, Rousseau C, Mangin I, Butel MJ, et al. Preterm infants with necrotising enterocolitis demonstrate an unbalanced gut microbiota. Acta Paediatr. 2018;107(1):40–7.
Dietert RR. The microbiome-immune-host defense barrier complex (microimmunosome) and developmental programming of noncommunicable diseases. Reprod Toxicol. 2017;68:49–58.
Dinan TG, Cryan JF. Gut instincts: microbiota as a key regulator of brain development, ageing and neurodegeneration. J Physiol. 2017;595(2):489–503.
Taft DH, Ambalavanan N, Schibler KR, Yu Z, Newburg DS, Ward DV, et al. Intestinal microbiota of preterm infants differ over time and between hospitals. Microbiome. 2014;2(1):2-36.
Mshvildadze M, Neu J, Shuster J, Theriaque D, Li N, Mai V. Intestinal microbial ecology in premature infants assessed with non-culture-based techniques. J Pediatr. 2010;156(1):20–5.
Lozupone CA, Stombaugh J, Gonzalez A, Ackermann G, Wendel D, Vazquez-Baeza Y, et al. Meta-analyses of studies of the human microbiota. Genome Res. 2013;23(10):1704–14.
Mao DP, Zhou Q, Chen CY, Quan ZX. Coverage evaluation of universal bacterial primers using the metagenomic datasets. BMC Microbiol. 2012;12:66.
Arboleya S, Salazar N, Solis G, Fernandez N, Hernandez-Barranco AM, Cuesta I, et al. Assessment of intestinal microbiota modulation ability of Bifidobacterium strains in in vitro fecal batch cultures from preterm neonates. Anaerobe. 2013;19:9–16.
Bjorkstrom MV, Hall L, Soderlund S, Hakansson EG, Hakansson S, Domellof M. Intestinal flora in very low-birth weight infants. Acta Paediatr. 2009;98(11):1762–7.
Carl MA, Ndao IM, Springman AC, Manning SD, Johnson JR, Johnston BD, et al. Sepsis from the gut: the enteric habitat of bacteria that cause late-onset neonatal bloodstream infections. Clin Infect Dis. 2014;58(9):1211–8.
Chrzanowska-Liszewska D, Seliga-Siwecka J, Kornacka MK. The effect of Lactobacillus rhamnosus GG supplemented enteral feeding on the microbiotic flora of preterm infants-double blinded randomized control trial. Early Hum Dev. 2012;88(1):57–60.
Ishizeki S, Sugita M, Takata M, Yaeshima T. Effect of administration of bifidobacteria on intestinal microbiota in low-birth-weight infants and transition of administered bifidobacteria: A comparison between one-species and three-species administration. Anaerobe. 2013;23:38–44.
Abdulkadir B, Nelson A, Skeath T, Marrs EC, Perry JD, Cummings SP, et al. Routine Use of Probiotics in Preterm Infants: Longitudinal Impact on the Microbiome and Metabolome. Neonatology. 2016;109(4):239–47.
Andriantsoanirina V, Teolis AC, Xin LX, Butel MJ, Aires J. Bifidobacterium longum and Bifidobacterium breve isolates from preterm and full term neonates: Comparison of cell surface properties. Anaerobe. 2014;28:212–5.
fArboleya S, Ang L, Margolles A, Li YY, Zhang DY, Liang X, et al. Deep 16S rRNA metagenomics and quantitative PCR analyses of the premature infant fecal microbiota. Anaerobe. 2012;18(3):378–80.
Arboleya S, Sánchez B, Milani C, Duranti S, Solís G, Fernández N, et al. Intestinal microbiota development in preterm neonates and effect of perinatal antibiotics. J Pediatr. 2015;166(3):538–44.
Arboleya S, Sanchez B, Solis G, Fernandez N, Suarez M, Hernandez-Barranco AM, et al. Impact of Prematurity and Perinatal Antibiotics on the Developing Intestinal Microbiota: A Functional Inference Study. Int J Mol Sci. 2016;17(5):649.
Ardissone AN, de la Cruz DM, Davis-Richardson AG, Rechcigl KT, Li N, Drew JC, et al. Meconium microbiome analysis identifies bacteria correlated with premature birth. PLoS One. 2014;9(3):e90784.
Barron LK, Warner BB, Tarr PI, Shannon WD, Deych E, Warner BW. Independence of gut bacterial content and neonatal necrotizing enterocolitis severity. J Pediatr Surg. 2017;52(6):993–8.
Cassir N, Benamar S, Khalil JB, Croce O, Saint-Faust M, Jacquot A, et al. Clostridium butyricum Strains and Dysbiosis Linked to Necrotizing Enterocolitis in Preterm Neonates. Clin Infect Dis. 2015;61(7):1107–15.
Chang JY, Shin SM, Chun J, Lee JH, Seo JK. Pyrosequencing-based Molecular Monitoring of the Intestinal Bacterial Colonization in Preterm Infants. J Pediatr Gastroenterol Nutr. 2011;53(5):512–9.
Costello EK, Carlisle EM, Bik EM, Morowitz MJ, Relman DA. Microbiome Assembly across Multiple Body Sites in Low-Birthweight Infants. Mbio. 2013;4(6).
Dardas M, Gill SR, Grier A, Pryhuber GS, Gill AL, Lee YH, et al. The impact of postnatal antibiotics on the preterm intestinal microbiome. Pediatr Res. 2014;76(2):150–8.
Gupta RW, Tran L, Norori J, Ferris MJ, Eren AM, Taylor CM, et al. Histamine-2 receptor blockers alter the fecal microbiota in premature infants. J Pediatr Gastroenterol Nutr. 2013;56(4):397–400.
La Rosa PS, Warner BB, Zhou Y, Weinstock GM, Sodergren E, Hall-Moore CM, et al. Patterned progression of bacterial populations in the premature infant gut. Proc Natl Acad Sci U S A. 2014;111(34):12522–7.
Leach ST, Lui K, Naing Z, Dowd SE, Mitchell HM, Day AS. Multiple Opportunistic Pathogens, but Not Pre-existing Inflammation, May Be Associated with Necrotizing Enterocolitis. Dig Dis Sci. 2015;60(12):3728–34.
Madan JC, Salari RC, Saxena D, Davidson L, O'Toole GA, Moore JH, et al. Gut microbial colonisation in premature neonates predicts neonatal sepsis. Arch Dis Child Fetal Neonatal Ed. 2012;97(6):F456–62.
Magne F, Abely M, Boyer F, Morville P, Pochart P, Suau A. Low species diversity and high interindividual variability in faeces of preterm infants as revealed by sequences of 16S rRNA genes and PCR-temporal temperature gradient gel electrophoresis profiles. FEMS Microbiol Ecol. 2006;57(1):128–38.
Mai V, Torrazza RM, Ukhanova M, Wang X, Sun Y, Li N, et al. Distortions in development of intestinal microbiota associated with late onset sepsis in preterm infants. PLoS One. 2013;8(1):e52876.
Mai V, Young CM, Ukhanova M, Wang XY, Sun YJ, Casella G, et al. Fecal Microbiota in Premature Infants Prior to Necrotizing Enterocolitis. PLoS One. 2011;6(6):e20647.
McMurtry VE, Gupta RW, Tran L, Blanchard EE, Penn D, Taylor CM, et al. Bacterial diversity and Clostridia abundance decrease with increasing severity of necrotizing enterocolitis. Microbiome. 2015;3(1)3-11.
Millar M, Seale J, Greenland M, Hardy P, Juszczak E, Wilks M, et al. The Microbiome of Infants Recruited to a Randomised Placebo-controlled Probiotic Trial (PiPS Trial). EBioMedicine. 2017;20:255–62.
Moles L, Gomez M, Jimenez E, Fernandez L, Bustos G, Chaves F, et al. Preterm infant gut colonization in the neonatal ICU and complete restoration 2 years later. Clin Microbiol Infect. 2015;21(10):936.e1–10.
Morowitz MJ, Denef VJ, Costello EK, Thomas BC, Poroyko V, Relman DA, et al. Strain-resolved community genomic analysis of gut microbial colonization in a premature infant. Proc Natl Acad Sci U S A. 2011;108(3):1128–33.
Arboleya S, Martinez-Camblor P, Solis G, Suarez M, Fernandez N, de los Reyes-Gavilan CG, et al. Intestinal Microbiota and Weight-Gain in Preterm Neonates. Front Microbiol. 2017;8:183.
Normann E, Fahlen A, Engstrand L, Lilja HE. Intestinal microbial profiles in extremely preterm infants with and without necrotizing enterocolitis. Acta Paediatr. 2013;102(2):129–36.
Poroyko V, Morowitz M, Bell T, Ulanov A, Wang M, Donovan S, et al. Diet creates metabolic niches in the "immature gut" that shape microbial communities. Nutr Hosp. 2011;26(6):1283–95.
Rosberg-Cody E, Ross RP, Hussey S, Ryan CA, Murphy BP, Fitzgerald GF, et al. Mining the microbiota of the neonatal gastrointestinal tract for conjugated linoleic acid-producing bifidobacteria. Appl Environ Microbiol. 2004;70(8):4635–41.
Said MB, Hays S, Maucort-Boulch D, Oulmaati A, Hantova S, Loys CM, et al. Gut microbiota in preterm infants with gross blood in stools: A prospective, controlled study. Early Hum Dev. 2014;90(10):579–85.
Shaw AG, Sim K, Randell P, Cox MJ, McClure ZE, Li MS, et al. Late-Onset Bloodstream Infection and Perturbed Maturation of the Gastrointestinal Microbiota in Premature Infants. PLoS One. 2015;10(7):e0132923.
Sherman MP, Sherman J, Arcinue R, Niklas V. Randomized Control Trial of Human Recombinant Lactoferrin: A Substudy Reveals Effects on the Fecal Microbiome of Very Low Birth Weight Infants. J Pediatr. 2016;173:S37–42.
Stewart CJ, Embleton ND, Marrs EC, Smith DP, Nelson A, Abdulkadir B, et al. Temporal bacterial and metabolic development of the preterm gut reveals specific signatures in health and disease. Microbiome. 2016;4(1):67.
Stewart CJ, Embleton ND, Marrs ECL, Smith DP, Fofanova T, Nelson A, et al. Longitudinal development of the gut microbiome and metabolome in preterm neonates with late onset sepsis and healthy controls. Microbiome. 2017;5(1):75.
Stewart CJ, Skeath T, Nelson A, Fernstad SJ, Marrs EC, Perry JD, et al. Preterm gut microbiota and metabolome following discharge from intensive care. Sci Rep. 2015;5:17141.
Taft DH, Ambalavanan N, Schibler KR, Yu Z, Newburg DS, Deshmukh H, et al. Center Variation in Intestinal Microbiota Prior to Late-Onset Sepsis in Preterm Infants. PLoS One. 2015;10(6):e0130604.
Underwood MA, Davis JCC, Kalanetra KM, Gehlot S, Patole S, Tancredi DJ, et al. Digestion of Human Milk Oligosaccharides by Bifidobacterium breve in the Premature Infant. J Pediatr Gastroenterol Nutr. 2017;65(4):449–55.
Underwood MA, Gaerlan S, De Leoz ML, Dimapasoc L, Kalanetra KM, Lemay DG, et al. Human milk oligosaccharides in premature infants: absorption, excretion, and influence on the intestinal microbiota. Pediatr Res. 2015;78(6):670–7.
Young GR, Smith DL, Embleton ND, Berrington JE, Schwalbe EC, Cummings SP, et al. Reducing Viability Bias in Analysis of Gut Microbiota in Preterm Infants at Risk of NEC and Sepsis. Front Cell Infect Microbiol. 2017;7:237.
Zeber-Lubecka N, Kulecka M, Ambrozkiewicz F, Paziewska A, Lechowicz M, Konopka E, et al. Effect of Saccharomyces boulardii and Mode of Delivery on the Early Development of the Gut Microbial Community in Preterm Infants. PLoS One. 2016;11(2):e0150306.
Zhou YJ, Shan G, Sodergren E, Weinstock G, Walker WA, Gregory KE. Longitudinal Analysis of the Premature Infant Intestinal Microbiome Prior to Necrotizing Enterocolitis: A Case-Control Study. PLoS One. 2015;10(3): e0118632.
Costeloe K, Hardy P, Juszczak E, Wilks M, Millar MR, Study PPI. Bifidobacterium breve BBG-001 in very preterm infants: a randomised controlled phase 3 trial. Lancet. 2016;387(10019):649–60.
Morrow AL, Lagomarcino AJ, Schibler KR, Taft DH, Yu Z, Wang B, et al. Early microbial and metabolomic signatures predict later onset of necrotizing enterocolitis in preterm infants. Microbiome. 2013;1(1):13.
Warner BB, Deych E, Zhou Y, Hall-Moore C, Weinstock GM, Sodergren E, et al. Gut bacteria dysbiosis and necrotising enterocolitis in very low birthweight infants: a prospective case-control study. Lancet. 2016;387(10031):1928–36.
Sim K, Shaw AG, Randell P, Cox MJ, McClure ZE, Li MS, et al. Dysbiosis anticipating necrotizing enterocolitis in very premature infants. Clin Infect Dis. 2015;60(3):389–97.
TTB iH, Groer MW, Kane B, Yee AL, Torres BA, Gilbert JA, et al. Dichotomous development of the gut microbiome in preterm infants. Microbiome. 2018;6(1):157.
Liu J, Li Y, Feng Y, Pan L, Xie Z, Yan Z, et al. Patterned progression of gut microbiota associated with necrotizing enterocolitis and late onset sepsis in preterm infants: a prospective study in a Chinese neonatal intensive care unit. PeerJ. 2019;7:e7310.
Aujoulat F, Roudiere L, Picaud JC, Jacquot A, Filleron A, Neveu D, et al. Temporal dynamics of the very premature infant gut dominant microbiota. BMC Microbiol. 2014;14:325.
Butel MJ, Suau A, Campeotto F, Magne F, Aires J, Ferraris L, et al. Conditions of bifidobacterial colonization in preterm infants: A prospective analysis. J Pediatr Gastroenterol Nutr. 2007;44(5):577–82.
Campeotto F, Suau A, Kapel N, Magne F, Viallon V, Ferraris L, et al. A fermented formula in pre-term infants: clinical tolerance, gut microbiota, down-regulation of faecal calprotectin and up-regulation of faecal secretory IgA. Br J Nutr. 2011;105(12):1843–51.
Ferraris L, Butel MJ, Campeotto F, Vodovar M, Roze JC, Aires J. Clostridia in premature neonates' gut: incidence, antibiotic susceptibility, and perinatal determinants influencing colonization. PLoS One. 2012;7(1):e30594.
Rouge C, Goldenberg O, Ferraris L, Berger B, Rochat F, Legrand A, et al. Investigation of the intestinal microbiota in preterm infants using different methods. Anaerobe. 2010;16(4):362–70.
Rouge C, Piloquet H, Butel MJ, Berger B, Rochat F, Ferraris L, et al. Oral supplementation with probiotics in very-low-birth-weight preterm infants: a randomized, double-blind, placebo-controlled trial. Am J Clin Nutr. 2009;89(6):1828–35.
de la Cochetière M-F, Piloquet H. des Robert C, Darmaun D, Galmiche J-P, Rozé J-C. Early Intestinal Bacterial Colonization and Necrotizing Enterocolitis in Premature Infants: The Putative Role of Clostridium. Pediatr Res. 2004;56(3):366–70.
Stewart CJ, Marrs EC, Magorrian S, Nelson A, Lanyon C, Perry JD, et al. The preterm gut microbiota: changes associated with necrotizing enterocolitis and infection. Acta Paediatr. 2012;101(11):1121–7.
Stewart CJ, Marrs ECL, Nelson A, Lanyon C, Perry JD, Embleton ND, et al. Development of the Preterm Gut Microbiome in Twins at Risk of Necrotising Enterocolitis and Sepsis. PLoS One. 2013;8(8):e73465.
Smith B, Bodé S, Skov TH, Mirsepasi H, Greisen G, Krogfelt KA. Investigation of the early intestinal microflora in premature infants with/without necrotizing enterocolitis using two different methods. Pediatr Res. 2012;71(1):115–20.
Millar MR, Linton CJ, Cade A, Glancy D, Hall M, Jalal H. Application of 16S rRNA gene PCR to study bowel flora of preterm infants with and without necrotizing enterocolitis. J Clin Microbiol. 1996;34(10):2506–10.
Moles L, Gomez M, Heilig H, Bustos G, Fuentes S, de Vos W, et al. Bacterial Diversity in Meconium of Preterm Neonates and Evolution of Their Fecal Microbiota during the First Month of Life. PLoS One. 2013;8(6):e66986.
Gomez M, Moles L, Espinosa-Martos I, Bustos G, de Vos WM, Fernandez L, et al. Bacteriological and Immunological Profiling of Meconium and Fecal Samples from Preterm Infants: A Two-Year Follow-Up Study. Nutrients. 2017;9(12).
Van Den Berg JP, Westerbeek EAM, Bröring-Starre T, Garssen J, Van Elburg RM. Neurodevelopment of preterm infants at 24 months after neonatal supplementation of a prebiotic mix: A randomized trial. J Pediatr Gastroenterol Nutr. 2016;63(2):270–6.
Westerbeek EAM, Slump RA, Lafeber HN, Knol J, Georgi G, Fetter WPF, et al. The effect of enteral supplementation of specific neutral and acidic oligosaccharides on the faecal microbiota and intestinal microenvironment in preterm infants. Eur J Clin Microbiol. 2013;32(2):269–76.
Underwood MA, Kalanetra KM, Bokulich NA, Lewis ZT, Mirmiran M, Tancredi DJ, et al. A comparison of two probiotic strains of bifidobacteria in premature infants. J Pediatr. 2013;163(6):1585–91.e9.
Underwood MA, Kalanetra KM, Bokulich NA, Mirmiran M, Barile D, Tancredi DJ, et al. Prebiotic oligosaccharides in premature infants. J Pediatr Gastroenterol Nutr. 2014;58(3):352–60.
Wang Y, Hoenig JD, Malin KJ, Qamar S, Petrof EO, Sun J, et al. 16S rRNA gene-based analysis of fecal microbiota from preterm infants with and without necrotizing enterocolitis. ISME J. 2009;3(8):944–54.
Abdulkadir B, Nelson A, Skeath T, Marrs EC, Perry JD, Cummings SP, et al. Stool bacterial load in preterm infants with necrotising enterocolitis. Early Hum Dev. 2016;95:1–2.
Aly H, Said RN, Wali IE, Elwakkad A, Soliman Y, Awad AR, et al. Medically Graded Honey Supplementation Formula to Preterm Infants as a Prebiotic: A Randomized Controlled Trial. J Pediatr Gastroenterol Nutr. 2017;64(6):966–70.
Arboleya S, Salazar N, Solis G, Fernandez N, Gueimonde M. de los Reyes-Gavilan CG. In vitro evaluation of the impact of human background microbiota on the response to Bifidobacterium strains and fructo-oligosaccharides. Br J Nutr. 2013;110(11):2030–6.
Arboleya S, Solis G, Fernandez N. de los Reyes-Gavilan CG, Gueimonde M. Facultative to strict anaerobes ratio in the preterm infant microbiota: a target for intervention? Gut Microbes. 2012;3(6):583–8.
Forsgren M, Isolauri E, Salminen S, Rautava S. Late preterm birth has direct and indirect effects on infant gut microbiota development during the first six months of life. Acta Paediatr. 2017;106(7):1103–9.
Gregory KE, LaPlante RD, Shan G, Kumar DV, Gregas M. Mode of Birth Influences Preterm Infant Intestinal Colonization With Bacteroides Over the Early Neonatal Period. Adv Neonatal Care. 2015;15(6):386–93.
Grzeskowiak L, Sales Teixeira TF, Bigonha SM, Lobo G, Salminen S, Ferreira CL. Gut Bifidobacterium microbiota in one-month-old Brazilian newborns. Anaerobe. 2015;35(Pt B):54–8.
Hickey L, Garland SM, Jacobs SE, O'Donnell CP, Tabrizi SN. Cross-colonization of infants with probiotic organisms in a neonatal unit. J Hosp Infect. 2014;88(4):226–9.
Luoto R, Rodriguez-Diaz J, Collado MC, Salminen S, Isolauri E, Lehtonen L. Gross blood in stools of premature neonates, a clinical and microbiological follow-up study. Acta Paediatr. 2013;102(5):486–91.
Underwood MA, Salzman NH, Bennett SH, Barman M, Mills DA, Marcobal A, et al. A randomized placebo-controlled comparison of 2 prebiotic/probiotic combinations in preterm infants: impact on weight gain, intestinal microbiota, and fecal short-chain fatty acids. J Pediatr Gastroenterol Nutr. 2009;48(2):216–25.
Leggett RM, Alcon-Giner C, Heavens D, Caim S, Brook TC, Kujawska M, et al. Rapid MinION profiling of preterm microbiota and antimicrobial-resistant pathogens. Nat Microbiol. 2019;5:430-442.
Gomez M, Moles L, Melgar A, Ureta N, Bustos G, Fernandez L, et al. Early Gut Colonization of Preterm Infants: Effect of Enteral Feeding Tubes. J Pediatr Gastroenterol Nutr. 2016;62(6):893–900.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: Open-Source, Platform-Independent, Community-Supported Software for Describing and Comparing Microbial Communities. Appl Environ Microbiol. 2009;75(23):7537–41.
Bertani G. Studies on lysogenesis. I. The mode of phage liberation by lysogenic Escherichia coli. J Bacteriol. 1951;62(3):293–300.
Macconkey A. Lactose-Fermenting Bacteria in Faeces. J Hyg (Lond). 1905;5(3):333–79.
Zajc-Satler J, Gragas AZ. Xylose lysine deoxycholate agar for the isolation of Salmonella and Shigella from clinical specimens. Zentralbl Bakteriol Orig A. 1977;237(2-3):196–200.
Lagier JC, Khelaifia S, Alou MT, Ndongo S, Dione N, Hugon P, et al. Culture of previously uncultured members of the human gut microbiota by culturomics. Nat Microbiol. 2016;1:16203.
Berrington JE, Stewart CJ, Embleton ND, Cummings SP. Gut microbiota in preterm infants: Assessment and relevance to health and disease. Arch Dis Child. 2013;98(4):F286–F90.
Amann RI, Ludwig W, Schleifer KH. Phylogenetic Identification and in-Situ Detection of Individual Microbial-Cells without Cultivation. Microbiol Rev. 1995;59(1):143–69.
Xiong W, Brown CT, Morowitz MJ, Banfield JF, Hettich RL. Genome-resolved metaproteomic characterization of preterm infant gut microbiota development reveals species-specific metabolic shifts and variabilities during early life. Microbiome. 2017;5(1):72.
Cardona S, Eck A, Cassellas M, Gallart M, Alastrue C, Dore J, et al. Storage conditions of intestinal microbiota matter in metagenomic analysis. BMC Microbiol. 2012;12:158.
Carroll IM, Ringel-Kulka T, Siddle JP, Klaenhammer TR, Ringel Y. Characterization of the Fecal Microbiota Using High-Throughput Sequencing Reveals a Stable Microbial Community during Storage. PLoS One. 2012;7(10):e46953.
Roesch LF, Casella G, Simell O, Krischer J, Wasserfall CH, Schatz D, et al. Influence of fecal sample storage on bacterial community diversity. Open Microbiol J. 2009;3:40–6.
Gorzelak MA, Gill SK, Tasnim N, Ahmadi-Vand Z, Jay M, Gibson DL. Methods for Improving Human Gut Microbiome Data by Reducing Variability through Sample Processing and Storage of Stool. PLoS One. 2015;10(8):e0139529.
Wu GD, Lewis JD, Hoffmann C, Chen YY, Knight R, Bittinger K, et al. Sampling and pyrosequencing methods for characterizing bacterial communities in the human gut using 16S sequence tags. BMC Microbiol. 2010;10:206.
Tedjo DI, Jonkers DMAE, Savelkoul PH, Masclee AA, van Best N, Pierik MJ, et al. The Effect of Sampling and Storage on the Fecal Microbiota Composition in Healthy and Diseased Subjects. PLoS One. 2015;10(5):e0126685.
Hale VL, Tan CL, Knight R, Amato KR. Effect of preservation method on spider monkey (Ateles geoffroyi) fecal microbiota over 8 weeks. J Microbiol Methods. 2015;113:16–26.
Shaw AG, Sim K, Powell E, Cornwell E, Cramer T, McClure ZE, et al. Latitude in sample handling and storage for infant faecal microbiota studies: the elephant in the room? Microbiome. 2016;4(1):40.
Walker AW, Martin JC, Scott P, Parkhill J, Flint HJ, Scott KP. 16S rRNA gene-based profiling of the human infant gut microbiota is strongly influenced by sample processing and PCR primer choice. Microbiome. 2015;3:26.
Hang J, Desai V, Zavaljevski N, Yang Y, Lin XX, Satya RV, et al. 16S rRNA gene pyrosequencing of reference and clinical samples and investigation of the temperature stability of microbiome profiles. Microbiome. 2014;2:31.
Maukonen J, Simoes C, Saarela M. The currently used commercial DNA-extraction methods give different results of clostridial and actinobacterial populations derived from human fecal samples. FEMS Microbiol Ecol. 2012;79(3):697–708.
Bahl MI, Bergstrom A, Licht TR. Freezing fecal samples prior to DNA extraction affects the Firmicutes to Bacteroidetes ratio determined by downstream quantitative PCR analysis. FEMS Microbiol Lett. 2012;329(2):193–7.
Nechvatal JM, Ram JL, Basson MD, Namprachan P, Niec SR, Badsha KZ, et al. Fecal collection, ambient preservation, and DNA extraction for PCR amplification of bacterial and human markers from human feces. J Microbiol Methods. 2008;72(2):124–32.
Han M, Hao L, Lin Y, Li F, Wang J, Yang H, et al. A novel affordable reagent for room temperature storage and transport of fecal samples for metagenomic analyses. Microbiome. 2018;6(1):43.
Anderson EL, Li W, Klitgord N, Highlander SK, Dayrit M, Seguritan V, et al. A robust ambient temperature collection and stabilization strategy: Enabling worldwide functional studies of the human microbiome. Sci Rep. 2016;6:31731.
Song SJ, Amir A, Metcalf JL, Amato KR, Xu ZZ, Humphrey G, et al. Preservation Methods Differ in Fecal Microbiome Stability. Affecting Suitability for Field Studies mSystems. 2016;1(3):e00021–16.
Choo JM, Leong LEX, Rogers GB. Sample storage conditions significantly influence faecal microbiome profiles. Sci Rep-Uk. 2015;5:16350.
Flores R, Shi J, Yu G, Ma B, Ravel J, Goedert JJ, et al. Collection media and delayed freezing effects on microbial composition of human stool. Microbiome. 2015;3(1):33.
Dominianni C, Wu J, Hayes RB, Ahn J. Comparison of methods for fecal microbiome biospecimen collection. BMC Microbiol. 2014;14:103.
Menke S, Gillingham MAF, Wilhelm K, Sommer S. Home-Made Cost Effective Preservation Buffer Is a Better Alternative to Commercial Preservation Methods for Microbiome Research. Front Microbiol. 2017;8:102.
Hill CJ, Brown JR, Lynch DB, Jeffery IB, Ryan CA, Ross RP, et al. Effect of room temperature transport vials on DNA quality and phylogenetic composition of faecal microbiota of elderly adults and infants. Microbiome. 2016;4(1):19.
genohub. Whole Genome Sequencing and Re-sequencing Guide [06/01/2021]. Available from: https://genohub.com/dna-seq-library-preparation/#:~:text=DNA%20Sample%20Submission%2D%20Typically%20100,10%20ng%20of%20input%20material.
Mackenzie BW, Waite DW, Taylor MW. Evaluating variation in human gut microbiota profiles due to DNA extraction method and inter-subject differences. Front Microbiol. 2015;6:130.
Costea PI, Zeller G, Sunagawa S, Pelletier E, Alberti A, Levenez F, et al. Towards standards for human fecal sample processing in metagenomic studies. Nat Biotechnol. 2017;35(11):1069–76.
Salonen A, Nikkila J, Jalanka-Tuovinen J, Immonen O, Rajilic-Stojanovic M, Kekkonen RA, et al. Comparative analysis of fecal DNA extraction methods with phylogenetic microarray: Effective recovery of bacterial and archaeal DNA using mechanical cell lysis. J Microbiol Methods. 2010;81(2):127–34.
Claassen S, du Toit E, Kaba M, Moodley C, Zar HJ, Nicol MP. A comparison of the efficiency of five different commercial DNA extraction kits for extraction of DNA from faecal samples. J Microbiol Methods. 2013;94(2):103–10.
Fiedorova K, Radvansky M, Nemcova E, Grombirikova H, Bosak J, Cernochova M, et al. The Impact of DNA Extraction Methods on Stool Bacterial and Fungal Microbiota Community Recovery. Front Microbiol. 2019;10:821.
Vebo HC, Karlsson MK, Avershina E, Finnby L, Rudi K. Bead-beating artefacts in the Bacteroidetes to Firmicutes ratio of the human stool metagenome. J Microbiol Methods. 2016;129:78–80.
Thomas V, Clark J, Doré J. Fecal microbiota analysis: an overview of sample collection methods and sequencing strategies. Future Microbiol. 2015;10(9):1485–504.
Lim MY, Song E-J, Kim SH, Lee J, Nam Y-D. Comparison of DNA extraction methods for human gut microbial community profiling. Syst Appl Microbiol. 2018;41(2):151–7.
Hamady M, Knight R. Microbial community profiling for human microbiome projects: Tools, techniques, and challenges. Genome Res. 2009;19(7):1141–52.
Frueh FW, Noyer-Weidner M. The use of denaturing high-performance liquid chromatography (DHPLC) for the analysis of genetic variations: impact for diagnostics and pharmacogenetics. Clin Chem Lab Med. 2003;41(4):452–61.
Lee PY, Costumbrado J, Hsu C-Y, Kim YH. Agarose gel electrophoresis for the separation of DNA fragments. J Vis Exp. 2012;62:3923.
Viglasky V. Polyacrylamide temperature gradient gel electrophoresis. Methods Mol Biol (Clifton, NJ). 2013;1054:159–71.
Muyzer G, Smalla K. Application of denaturing gradient gel electrophoresis (DGGE) and temperature gradient gel electrophoresis (TGGE) in microbial ecology. Antonie Van Leeuwenhoek. 1998;73(1):127–41.
Liesack W, Dunfield PF. T-RFLP Analysis. In: Walker JM, JFT S, Ragout de Spencer AL, editors. Environmental Microbiology: Methods and Protocols. Totowa: Humana Press; 2004. p. 23–37.
Osborn AM, Moore ERB, Timmis KN. An evaluation of terminal-restriction fragment length polymorphism (T-RFLP) analysis for the study of microbial community structure and dynamics. Environ Microbiol. 2000;2(1):39–50.
Bokulich NA, Mills DA. Differentiation of mixed lactic acid bacteria communities in beverage fermentations using targeted terminal restriction fragment length polymorphism. Food Microbiol. 2012;31(1):126–32.
O'Donovan MC, Oefner PJ, Roberts SC, Austin J, Hoogendoorn B, Guy C, et al. Blind analysis of denaturing high-performance liquid chromatography as a tool for mutation detection. Genomics. 1998;52(1):44–9.
Xiao W, Oefner PJ. Denaturing high-performance liquid chromatography: A review. Hum Mutat. 2001;17(6):439–74.
Hannachi-M’Zali F, Ambler JE, Taylor CF, Hawkey PM. Examination of single and multiple mutations involved in resistance to quinolones in Staphylococcus aureus by a combination of PCR and denaturing high-performance liquid chromatography (DHPLC). J Antimicrob Chemother. 2002;50(5):649–55.
Domann E, Hong G, Imirzalioglu C, Turschner S, Kuhle J, Watzel C, et al. Culture-independent identification of pathogenic bacteria and polymicrobial infections in the genitourinary tract of renal transplant recipients. J Clin Microbiol. 2003;41(12):5500–10.
Lueders T, Friedrich MW. Evaluation of PCR Amplification Bias by Terminal Restriction Fragment Length Polymorphism Analysis of Small-Subunit rRNA and <em>mcrA</em> Genes by Using Defined Template Mixtures of Methanogenic Pure Cultures and Soil DNA Extracts. Appl Environ Microbiol. 2003;69(1):320.
Engelbrektson A, Kunin V, Wrighton KC, Zvenigorodsky N, Chen F, Ochman H, et al. Experimental factors affecting PCR-based estimates of microbial species richness and evenness. ISME J. 2010;4(5):642–7.
Fierer N, Jackson RB. The diversity and biogeography of soil bacterial communities. Proc Natl Acad Sci U S A. 2006;103(3):626–31.
Paliy O, Agans R. Application of phylogenetic microarrays to interrogation of human microbiota. FEMS Microbiol Ecol. 2012;79(1):2–11.
Govindarajan R, Duraiyan J, Kaliyappan K, Palanisamy M. Microarray and its applications. J Pharm Bioallied Sci. 2012;4(Suppl 2):S310–S2.
Brodie, E.L. and Lynch, S.V. (2013). Phyloarrays. In The Human Microbiota, D.N. Fredricks (Ed.). p. 75-84.
Rajilic-Stojanovic M, Heilig HGHJ, Molenaar D, Kajander K, Surakka A, Smidt H, et al. Development and application of the human intestinal tract chip, a phylogenetic microarray: analysis of universally conserved phylotypes in the abundant microbiota of young and elderly adults. Environ Microbiol. 2009;11(7):1736–51.
Wilson KH, Wilson WJ, Radosevich JL, DeSantis TZ, Viswanathan VS, Kuczmarski TA, et al. High-density microarray of small-subunit ribosomal DNA probes. Appl Environ Microbiol. 2002;68(5):2535–41.
He Z, Gentry TJ, Schadt CW, Wu L, Liebich J, Chong SC, et al. GeoChip: a comprehensive microarray for investigating biogeochemical, ecological and environmental processes. ISME J. 2007;1(1):67–77.
Claesson MJ, O'Sullivan O, Wang Q, Nikkila J, Marchesi JR, Smidt H, et al. Comparative analysis of pyrosequencing and a phylogenetic microarray for exploring microbial community structures in the human distal intestine. PLoS One. 2009;4(8):e6669.
van den Bogert B, de Vos WM, Zoetendal EG, Kleerebezem M. Microarray analysis and barcoded pyrosequencing provide consistent microbial profiles depending on the source of human intestinal samples. Appl Environ Microbiol. 2011;77(6):2071–80.
DeSantis TZ, Brodie EL, Moberg JP, Zubieta IX, Piceno YM, Andersen GL. High-density universal 16S rRNA microarray analysis reveals broader diversity than typical clone library when sampling the environment. Microb Ecol. 2007;53(3):371–83.
Deepak S, Kottapalli K, Rakwal R, Oros G, Rangappa K, Iwahashi H, et al. Real-Time PCR: Revolutionizing Detection and Expression Analysis of Genes. Curr Genomics. 2007;8(4):234–51.
Wong W, Farr R, Joglekar M, Januszewski A, Hardikar A. Probe-based Real-time PCR Approaches for Quantitative Measurement of microRNAs. J Vis Exp. 2015;98:52586.
Wallner G, Amann R, Beisker W. Optimizing fluorescent in situ hybridization with rRNA-targeted oligonucleotide probes for flow cytometric identification of microorganisms. Cytometry. 1993;14(2):136–43.
Jian C, Luukkonen P, Yki-Järvinen H, Salonen A, Korpela K. Quantitative PCR provides a simple and accessible method for quantitative microbiota profiling. PLoS One. 2020;15(1):e0227285-e.
Vázquez L, Guadamuro L, Giganto F, Mayo B, Flórez AB. Development and Use of a Real-Time Quantitative PCR Method for Detecting and Quantifying Equol-Producing Bacteria in Human Faecal Samples and Slurry Cultures. Front Microbiol. 2017;(8):1155.
Sanger F, Coulson AR. A rapid method for determining sequences in DNA by primed synthesis with DNA polymerase. J Mol Biol. 1975;94(3):441–8.
Churko JM, Mantalas GL, Snyder MP, Wu JC. Overview of high throughput sequencing technologies to elucidate molecular pathways in cardiovascular diseases. Circ Res. 2013;112(12):1613–23.
Liu Z, DeSantis TZ, Andersen GL, Knight R. Accurate taxonomy assignments from 16S rRNA sequences produced by highly parallel pyrosequencers. Nucleic Acids Res. 2008;36(18):e120.
Liu Z, Lozupone C, Hamady M, Bushman FD, Knight R. Short pyrosequencing reads suffice for accurate microbial community analysis. Nucleic Acids Res. 2007;35(18):e120.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.
Almeida A, Mitchell AL, Tarkowska A, Finn RD. Benchmarking taxonomic assignments based on 16S rRNA gene profiling of the microbiota from commonly sampled environments. Gigascience. 2018;7(5):1-10.
Schloss PD, Gevers D, Westcott SL. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS One. 2011;6(12):e27310.
Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013;41(1):e1.
Tremblay J, Singh K, Fern A, Kirton ES, He SM, Woyke T, et al. Primer and platform effects on 16S rRNA tag sequencing. Front Microbiol. 2015;6:771.
Baker GC, Smith JJ, Cowan DA. Review and re-analysis of domain-specific 16S primers. J Microbiol Methods. 2003;55(3):541–55.
Pammi M, Cope J, Tarr PI, Warner BB, Morrow AL, Mai V, et al. Intestinal dysbiosis in preterm infants preceding necrotizing enterocolitis: a systematic review and meta-analysis. Microbiome. 2017;5(1):31.
Wang Y, Qian PY. Conservative fragments in bacterial 16S rRNA genes and primer design for 16S ribosomal DNA amplicons in metagenomic studies. PLoS One. 2009;4(10):e7401.
Tringe SG, Hugenholtz P. A renaissance for the pioneering 16S rRNA gene. Curr Opin Microbiol. 2008;11(5):442–6.
Sim K, Cox MJ, Wopereis H, Martin R, Knol J, Li MS, et al. Improved detection of bifidobacteria with optimised 16S rRNA-gene based pyrosequencing. PLoS One. 2012;7(3):e32543.
Johnson JS, Spakowicz DJ, Hong BY, Petersen LM, Demkowicz P, Chen L, et al. Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis. Nat Commun. 2019;10:5029.
Yang B, Wang Y, Qian PY. Sensitivity and correlation of hypervariable regions in 16S rRNA genes in phylogenetic analysis. BMC Bioinformatics. 2016;17:135.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79(17):5112–20.
Rausch P, Rühlemann M, Hermes BM, Doms S, Dagan T, Dierking K, et al. Comparative analysis of amplicon and metagenomic sequencing methods reveals key features in the evolution of animal metaorganisms. Microbiome. 2019;7(1):133.
Thompson LR, Sanders JG, McDonald D, Amir A, Ladau J, Locey KJ, et al. A communal catalogue reveals Earth’s multiscale microbial diversity. Nature. 2017;551(7681):457–63.
Luo C, Tsementzi D, Kyrpides N, Read T, Konstantinidis KT. Direct Comparisons of Illumina vs. Roche 454 Sequencing Technologies on the Same Microbial Community DNA Sample. PLOS ONE. 2012;7(2):e30087.
Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, et al. Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005;437(7057):376–80.
Quince C, Lanzén A, Curtis TP, Davenport RJ, Hall N, Head IM, et al. Accurate determination of microbial diversity from 454 pyrosequencing data. Nat Methods. 2009;6(9):639–41.
Loman NJ, Misra RV, Dallman TJ, Constantinidou C, Gharbia SE, Wain J, et al. Performance comparison of benchtop high-throughput sequencing platforms. Nat Biotechnol. 2012;30(5):434-+.
Gomez-Alvarez V, Teal TK, Schmidt TM. Systematic artifacts in metagenomes from complex microbial communities. ISME J. 2009;3(11):1314–7.
Lam HYK, Clark MJ, Chen R, Chen R, Natsoulis G, O'Huallachain M, et al. Performance comparison of whole-genome sequencing platforms. Nat Biotechnol. 2012;30(1):78–82.
Quail MA, Smith M, Coupland P, Otto TD, Harris SR, Connor TR, et al. A tale of three next generation sequencing platforms: comparison of Ion Torrent, Pacific Biosciences and Illumina MiSeq sequencers. BMC Genomics. 2012;13:341.
Salipante SJ, Kawashima T, Rosenthal C, Hoogestraat DR, Cummings LA, Sengupta DJ, et al. Performance comparison of Illumina and ion torrent next-generation sequencing platforms for 16S rRNA-based bacterial community profiling. Appl Environ Microbiol. 2014;80(24):7583–91.
Erlich Y, Mitra PP. delaBastide M, McCombie WR, Hannon GJ. Alta-Cyclic: a self-optimizing base caller for next-generation sequencing. Nat Methods. 2008;5(8):679–82.
Dolan PC, Denver DR. TileQC: A system for tile-based quality control of Solexa data. BMC Bioinformatics. 2008;9(1):250.
Nakamura K, Oshima T, Morimoto T, Ikeda S, Yoshikawa H, Shiwa Y, et al. Sequence-specific error profile of Illumina sequencers. Nucleic Acids Res. 2011;39(13):e90.
Schröder J, Bailey J, Conway T, Zobel J. Reference-Free Validation of Short Read Data. PLoS One. 2010;5(9):e12681.
Caporaso JGKJ, Stombaugh J, Bittinger K, Bushman FD, Costello EK. QIIME allows analysis of high-throughput community sequencing data: Nat Methods; 2010.
Callahan BJ, McMurdie PJ, Holmes SP. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 2017;11(12):2639–43.
Matias Rodrigues JF, Schmidt TSB, Tackmann J, von Mering C. MAPseq: highly efficient k-mer search with confidence estimates, for rRNA sequence analysis. Bioinformatics. 2017;33(23):3808–10.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.
Prodan A, Tremaroli V, Brolin H, Zwinderman AH, Nieuwdorp M, Levin E. Comparing bioinformatic pipelines for microbial 16S rRNA amplicon sequencing. PLoS One. 2020;15(1):e0227434-e.
Westcott SL, Schloss PD. De novo clustering methods outperform reference-based methods for assigning 16S rRNA gene sequences to operational taxonomic units. PeerJ. 2015;3:e1487-e.
Rideout JR, He Y, Navas-Molina JA, Walters WA, Ursell LK, Gibbons SM, et al. Subsampled open-reference clustering creates consistent, comprehensive OTU definitions and scales to billions of sequences. PeerJ. 2014;2:e545.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: High resolution sample inference from Illumina amplicon data. Nature Methods. 2016;13:581-583.
Amir A, McDonald D, Navas-Molina JA, Kopylova E, Morton JT, Zech Xu Z, et al. Deblur Rapidly Resolves Single-Nucleotide Community Sequence Patterns. mSystems. 2017;2(2):e00191–16.
Nearing JT, Douglas GM, Comeau AM, Langille MGI. Denoising the Denoisers: an independent evaluation of microbiome sequence error-correction approaches. PeerJ. 2018;6:e5364.
Schloss PD, Handelsman J. Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. Appl Environ Microbiol. 2005;71(3):1501–6.
Edgar RC. Updating the 97% identity threshold for 16S ribosomal RNA OTUs. Bioinformatics. 2018;34(14):2371–5.
Balvočiūtė M, Huson DH. SILVA, RDP, Greengenes, NCBI and OTT - how do these taxonomies compare? BMC Genomics. 2017;18(Suppl 2):114.
Glassing A, Dowd SE, Galandiuk S, Davis B, Chiodini RJ. Inherent bacterial DNA contamination of extraction and sequencing reagents may affect interpretation of microbiota in low bacterial biomass samples. Gut Pathogens. 2016;8(1):24.
Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF, et al. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 2014;12:87.
McKnight D, Huerlimann R, Bower DS, Schwarzkopf L, Alford RA, Zenger KR. microDecon: A highly accurate read-subtraction tool for the post-sequencing removal of contamination in metabarcoding studies. Environmental DNA. 2019;14(1):25.
Willis AD. Rarefaction, Alpha Diversity, and Statistics. Front Microbiol. 2019;10:2407.
Chao A. Nonparametric Estimation of the Number of Classes in a Population. Scand J Stat. 1984;11(4):265–70.
Vandeputte D, Kathagen G, D’hoe K, Vieira-Silva S, Valles-Colomer M, Sabino J, et al. Quantitative microbiome profiling links gut community variation to microbial load. Nature. 2017;551(7681):507–11.
Lande R, DeVries PJ, Walla TR. When species accumulation curves intersect: implications for ranking diversity using small samples. Oikos. 2000;89(3):601–5.
Bullard JH, Purdom E, Hansen KD, Dudoit S. Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics. 2010;11(1):94.
Dillies MA, Rau A, Aubert J, Hennequet-Antier C, Jeanmougin M, Servant N, et al. A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis. Brief Bioinform. 2013;14(6):671–83.
McMurdie PJ, Holmes S. Waste Not, Want Not: Why Rarefying Microbiome Data Is Inadmissible. PLoS Comput Biol. 2014;10(4):e1003531.
Weiss S, Xu ZZ, Peddada S, Amir A, Bittinger K, Gonzalez A, et al. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome. 2017;5(1):27.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Paulson JN, Stine OC, Bravo HC, Pop M. Differential abundance analysis for microbial marker-gene surveys. Nat Methods. 2013;10(12):1200–2.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
McKnight DT, Huerlimann R, Bower DS, Schwarzkopf L, Alford RA, Zenger KR. Methods for normalizing microbiome data: An ecological perspective. Methods Ecol Evol. 2019;10(3):389–400.
Quince C, Walker AW, Simpson JT, Loman NJ, Segata N. Shotgun metagenomics, from sampling to analysis. Nat Biotechnol. 2017;35(9):833–44.
Wesolowska-Andersen A, Bahl MI, Carvalho V, Kristiansen K, Sicheritz-Ponten T, Gupta R, et al. Choice of bacterial DNA extraction method from fecal material influences community structure as evaluated by metagenomic analysis. Microbiome. 2014;2:19.
Tanner MA, Goebel BM, Dojka MA, Pace NR. Specific ribosomal DNA sequences from diverse environmental settings correlate with experimental contaminants. Appl Environ Microbiol. 1998;64(8):3110–3.
Yuan SQ, Cohen DB, Ravel J, Abdo Z, Forney LJ. Evaluation of Methods for the Extraction and Purification of DNA from the Human Microbiome. PLoS One. 2012;7(3):e33865.
Nelson MC, Morrison HG, Benjamino J, Grim SL, Graf J. Analysis, Optimization and Verification of Illumina-Generated 16S rRNA Gene Amplicon Surveys. PLoS One. 2014;9(4):e94249.
Baym M, Kryazhimskiy S, Lieberman TD, Chung H, Desai MM, Kishony R. Inexpensive Multiplexed Library Preparation for Megabase-Sized Genomes. PLoS One. 2015;10(6):e0128036.
Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019;20(1):257.
Kim D, Song L, Breitwieser FP, Salzberg SL. Centrifuge: rapid and sensitive classification of metagenomic sequences. Genome Res. 2016;26(12):1721–9.
Milanese A, Mende DR, Paoli L, Salazar G, Ruscheweyh H-J, Cuenca M, et al. Microbial abundance, activity and population genomic profiling with mOTUs2. Nat Commun. 2019;10(1):1014.
Peng Y, Leung HCM, Yiu SM, Chin FYL. Meta-IDBA: a de Novo assembler for metagenomic data. Bioinformatics (Oxford, England). 2011;27(13):i94–i101.
Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, et al. De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010;20(2):265–72.
Salipante SJ, Sengupta DJ, Rosenthal C, Costa G, Spangler J, Sims EH, et al. Rapid 16S rRNA next-generation sequencing of polymicrobial clinical samples for diagnosis of complex bacterial infections. PLoS One. 2013;8(5):e65226.
Sunagawa S, Mende DR, Zeller G, Izquierdo-Carrasco F, Berger SA, Kultima JR, et al. Metagenomic species profiling using universal phylogenetic marker genes. Nat Methods. 2013;10(12):1196–9.
Jain M, Olsen HE, Paten B, Akeson M. The Oxford Nanopore MinION: delivery of nanopore sequencing to the genomics community. Genome Biol. 2016;17(1):239.
Hebert PDN, Braukmann TWA, Prosser SWJ, Ratnasingham S. dewaard JR, Ivanova NV, et al. A Sequel to Sanger: amplicon sequencing that scales. BMC Genomics. 2018;19(1):219.
Goodwin S, McPherson JD, McCombie WR. Coming of age: ten years of next-generation sequencing technologies. Nat Rev Genet. 2016;17(6):333–51.
Levy SE, Myers RM. Advancements in Next-Generation Sequencing. Annu Rev Genomics Hum Genet. 2016;17:95–115.
Payne A, Holmes N, Rakyan V, Loose M. BulkVis: a graphical viewer for Oxford nanopore bulk FAST5 files. Bioinformatics. 2019;35(13):2193-2198.
Callahan BJ, Wong J, Heiner C, Oh S, Theriot CM, Gulati AS, et al. High-throughput amplicon sequencing of the full-length 16S rRNA gene with single-nucleotide resolution. Nucleic Acids Res. 2019;47(18):e103.
Jain M, Fiddes IT, Miga KH, Olsen HE, Paten B, Akeson M. Improved data analysis for the MinION nanopore sequencer. Nat Methods. 2015;12(4):351–6.
Jain M, Koren S, Miga KH, Quick J, Rand AC, Sasani TA, et al. Nanopore sequencing and assembly of a human genome with ultra-long reads. Nat Biotechnol. 2018;36(4):338–45.
Carneiro MO, Russ C, Ross MG, Gabriel SB, Nusbaum C, DePristo MA. Pacific biosciences sequencing technology for genotyping and variation discovery in human data. BMC Genomics. 2012;13:375.
Rang FJ, Kloosterman WP, de Ridder J. From squiggle to basepair: computational approaches for improving nanopore sequencing read accuracy. Genome Biol. 2018;19(1):90.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Westaway, J.A.F., Huerlimann, R., Miller, C.M. et al. Methods for exploring the faecal microbiome of premature infants: a review. matern health, neonatol and perinatol 7, 11 (2021). https://doi.org/10.1186/s40748-021-00131-9