Kinetic bias in estimates of coastal picoplankton community structure obtained by measurements of small-subunit rRNA gene PCR amplicon length heterogeneity

Factors structuring microbial communities in highly impacted coastal marine sediments (Mar Menor lagoon, SE Spain)

Coastal marine lagoons are environments highly vulnerable to anthropogenic pressures such as agriculture nutrient loading or runoff from metalliferous mining. Sediment microorganisms, which are key components in the biogeochemical cycles, can help attenuate these impacts by accumulating nutrients and pollutants. The Mar Menor, located in the southeast of Spain, is an example of a coastal lagoon strongly altered by anthropic pressures, but the microbial community inhabiting its sediments remains unknown. Here, we describe the sediment prokaryotic communities along a wide range of environmental conditions in the lagoon, revealing that microbial communities were highly heterogeneous among stations, although a core microbiome was detected. The microbiota was dominated by Delta- and Gammaproteobacteria and members of the Bacteroidia class. Additionally, several uncultured groups such as Asgardarchaeota were detected in relatively high proportions. Sediment texture, the presence of Caulerpa or Cymodocea, depth, and geographic location were among the most important factors structuring microbial assemblages. Furthermore, microbial communities in the stations with the highest concentrations of potentially toxic elements (Fe, Pb, As, Zn, and Cd) were less stable than those in the non-contaminated stations. This finding suggests that bacteria colonizing heavily contaminated stations are specialists sensitive to change.

Past, Present, and Future of DNA Typing for Analyzing Human and Non-Human Forensic Samples

Forensic DNA analysis has vastly evolved since the first forensic samples were evaluated by restriction fragment length polymorphism (RFLP). Methodologies advanced from gel electrophoresis techniques to capillary electrophoresis and now to next generation sequencing (NGS). Capillary electrophoresis was and still is the standard method used in forensic analysis. However, dependent upon the information needed, there are several different techniques that can be used to type a DNA fragment. Short tandem repeat (STR) fragment analysis, Sanger sequencing, SNapShot, and capillary electrophoresis-single strand conformation polymorphism (CE-SSCP) are a few of the techniques that have been used for the genetic analysis of DNA samples. NGS is the newest and most revolutionary technology and has the potential to be the next standard for genetic analysis. This review briefly encompasses many of the techniques and applications that have been utilized for the analysis of human and nonhuman DNA samples.

Comparative study of the hemolymph microbiome between live and recently dead American lobsters Homarus americanus

Lobsters and other crustaceans do not have sterile hemolymph. Despite this, little is known about the microbiome in the hemolymph of the lobster Homarus americanus. The purpose of this study was to characterize the hemolymph microbiome in lobsters. The lobsters were part of a larger study on the effect of temperature on epizootic shell disease, and several died during the course of the study, providing an opportunity to examine differences in the microbiomes between live and recently dead (1-24 h) animals. The hemolymph microbiomes of live lobsters was different from those in dead animals and both were different from the tank microbiome in which the animals had been held. The microbiomes of live lobsters were more diverse and had a different suite of bacteria than those from dead animals. The dominant taxa in live lobsters belonged to Flavobacteriaceae and Rhodobacteraceae, whereas Vibrionaceae and Enterobacteriaceae were predominant in the dead lobsters. Although aquarium microbiomes overlapped with the hemolymph microbiomes, there was less overlap and lower abundance of taxa in comparison with hemolymph from live lobsters. Previous studies reporting bacteria in the digestive tract of lobsters suggested that Vibrionaceae and Enterobacteriaceae had invaded the hemolymph via the gut. Our study suggests that hemolymph bacteria abundant in live lobsters do not originate from the tank milieu and comprise a rich, natural, or native background of bacterial constituents.

Botanical microbiomes on the cheap: Inexpensive molecular fingerprinting methods to study plant-associated communities of bacteria and fungi

High-throughput sequencing technologies have revolutionized the study of plant-associated microbial populations, but they are relatively expensive. Molecular fingerprinting techniques are more affordable, yet yield considerably less information about the microbial community. Does this mean they are no longer useful for plant microbiome research? In this paper, we review the past 10 years of studies on plant-associated microbiomes using molecular fingerprinting methodologies, including single-strand conformation polymorphism (SSCP), denaturing gradient gel electrophoresis (DGGE), amplicon length heterogeneity PCR (LH-PCR), ribosomal intergenic spacer analysis (RISA) and automated ribosomal intergenic spacer analysis (ARISA), and terminal restriction fragment length polymorphism (TRFLP). We also present data juxtaposing results from TRFLP methods with those generated using Illumina sequencing in the comparison of rhizobacterial populations of Brazilian maize and fungal surveys in Canadian tomato roots. In both cases, the TRFLP approach yielded the desired results at a level of resolution comparable to that of the MiSeq method, but at a fraction of the cost. Community fingerprinting methods (especially TRFLP) remain relevant for the identification of dominant microbes in a population, the observation of shifts in plant microbiome community diversity, and for screening samples before their use in more sensitive and expensive approaches.

Technical pitfalls that bias comparative microbial community analyses of aquatic disease Ian Hewson

The accessibility of high-throughput DNA sequencing technologies has attracted the application of comparative microbial analyses to study diseases. These studies present a window into host microbiome diversity and composition that can be used to address ecological theory in the context of host biology and behavior. Recently, comparative microbiome studies have been used to study non-vertebrate aquatic diseases to elucidate microorganisms potentially involved in disease processes or in disease prevention. These investigations suffer from many well-described biases, especially prior to sequence analyses, that could lead to misleading conclusions. Microbiome-focused studies of aquatic metazoan diseases provide valuable documentation of microbial ecology, although, they are only a starting point for establishing disease etiology, which demands quantitative validation through targeted approaches. The microbiome approach to understanding disease is most useful after laboratory diagnostics guided by pathology have failed to identify a causative agent. This opinion piece presents several technical pitfalls which may affect wider interpretation of microbe-host interactions through comparative microbial community analyses and provides recommendations, based on studies in non-aquatic systems, for incorporation into future aquatic disease research.

Patterns of PCR Amplification Artifacts of the Fungal Barcode Marker in a Hybrid Mushroom

The polymerase chain reaction (PCR) is widely used in modern biology and medicine. However, PCR artifacts can complicate the interpretation of PCR-based results. The internal transcribed spacer (ITS) region of the ribosomal RNA gene cluster is the consensus fungal barcode marker and suspected PCR artifacts have been reported in many studies, especially for the analyses of environmental fungal samples. At present, the patterns of PCR artifacts in the whole fungal ITS region (ITS1+5.8S+ITS2) are not known. In this study, we analyzed the error rates of PCR at three template complexity levels using the divergent copies of ITS from the mushroom Agaricus subrufescens. Our results showed that PCR using the Phusion^® High-Fidelity DNA Polymerase has a per nucleotide error rate of about 4 × 10^–6 per replication. Among the detected mutations, transitions were much more frequent than transversions, insertions, and deletions. When divergent alleles were mixed as templates in the same reaction, a significant proportion (∼30%) of recombinant molecules were detected. The in vitro mixed-template results were comparable to those obtained from using the genomic DNA of the original mushroom specimen as template. Our results indicate that caution should be in place when interpreting ITS sequences from individual fungal specimens, especially those containing divergent ITS copies. Similar results could also happen to PCR-based analyses of other multicopy DNA fragments as well as single-copy DNA sequences with divergent alleles in diploid organisms.

Advancement in LH-PCR methodology for multiple microbial species detections in fermented foods

The length-heterogeneity PCR is a low throughput molecular biology methods explored to monitor bacteria populations in different environments. It could be more used in food microbiology analysis, not only for fingerprinting analysis, but it has been hampered until now by a limiting factor which relates to the high percentage of secondary peaks. With the aim to overcome this problem, different experiments were performed focusing on changing PCR parameters in order to obtain more specific amplicon patterns and also to reduce the complexity of community patterns. With this purpose, different annealing temperatures were tested on complex fermented food matrices taken from both animal and vegetable origin and also on the bacteria isolated from the same food source. In particular, the optimal annealing temperature identified for the fermented food samples is 59 °C and the optimal for bacterial strains varied between 63 °C and 65 °C. The approach allowed the modification of the LH-PCR protocol increasing the amplification efficiency and therefore the bacteria species discrimination. These temperatures also allowed the implementation of the previous LH-PCR published database. The modification in the level of accuracy of the LH-PCR technique could also allow an improvement in the relative species quantification by the peak area evaluation.

Assessment of Fecal Indicator Bacteria and Potential Pathogen Co-Occurrence at a Shellfish Growing Area

Routine monitoring of shellfish growing waters for bacteria indicative of human sewage pollution reveals little about the bacterial communities that co-occur with these indicators. This study investigated the bacterial community, potential pathogens, and fecal indicator bacteria in 40 water samples from a shellfish growing area in the Chesapeake Bay, USA. Bacterial community composition was quantified with deep sequencing of 16S rRNA gene amplicons, and absolute gene abundances were estimated with an internal standard (Thermus thermophilus genomes). Fecal coliforms were quantified by culture, and Vibrio vulnificus and V. parahaemolyticus with quantitative PCR. Fecal coliforms and V. vulnificus were detected in most samples, and a diverse assemblage of potential human pathogens were detected in all samples. These taxa followed two general patterns of abundance. Fecal coliforms and 16S rRNA genes for Enterobacteriaceae, Aeromonas, Arcobacter, Staphylococcus, and Bacteroides increased in abundance after a 1.3-inch rain event in May, and, for some taxa, after smaller rain events later in the season, suggesting that these are allochthonous organisms washed in from land. Clostridiaceae and Mycobacterium 16S rRNA gene abundances increased with day of the year and were not positively related to rainfall, suggesting that these are autochthonous organisms. Other groups followed both patterns, such as Legionella. Fecal coliform abundance did not correlate with most other taxa, but were extremely high following the large rainstorm in May when they co-occurred with a broad range of potential pathogen groups. V. vulnificus were absent during the large rainstorm, and did not correlate with 16S rRNA abundances of Vibrio spp. or most other taxa. These results highlight the complex nature of bacterial communities and the limited utility of using specific bacterial groups as indicators of pathogen presence.

Bioinformatics Approach to Assess the Biogeographical Patterns of Soil Communities: The Utility for Soil Provenance

Soil DNA profiling has potential as a forensic tool to establish a link between soil collected at a crime scene and soil recovered from a suspect. However, a quantitative measure is needed to investigate the spatial/temporal variability across multiple scales prior to their application in forensic science. In this study, soil DNA profiles across Miami-Dade, FL, were generated using length heterogeneity PCR to target four taxa. The objectives of this study were to (i) assess the biogeographical patterns of soils to determine whether soil biota is spatially correlated with geographic location and (ii) evaluate five machine learning algorithms for their predictive ability to recognize biotic patterns which could accurately classify soils at different spatial scales regardless of seasonal collection. Results demonstrate that soil communities have unique patterns and are spatially autocorrelated. Bioinformatic algorithms could accurately classify soils across all scales with Random Forest significantly outperforming all other algorithms regardless of spatial level.

Methylophilaceae and Hyphomicrobium as target taxonomic groups in monitoring the function of methanol-fed denitrification biofilters in municipal wastewater treatment plants

Molecular monitoring of bacterial communities can explain and predict the stability of bioprocesses in varying physicochemical conditions. To study methanol-fed denitrification biofilters of municipal wastewater treatment plants, bacterial communities of two full-scale biofilters were compared through fingerprinting and sequencing of the 16S rRNA genes. Additionally, 16S rRNA gene fingerprinting was used for 10-week temporal monitoring of the bacterial community in one of the biofilters. Combining the data with previous study results, the family Methylophilaceae and genus Hyphomicrobium were determined as suitable target groups for monitoring. An increase in the relative abundance of Hyphomicrobium-related biomarkers occurred simultaneously with increases in water flow, NOx  − load, and methanol addition, as well as a higher denitrification rate, although the dominating biomarkers linked to Methylophilaceae showed an opposite pattern. The results indicate that during increased loading, stability of the bioprocess is maintained by selection of more efficient denitrifier populations, and this progress can be analyzed using simple molecular fingerprinting.

Environmental Sequencing Provides Reasonable Estimates of the Relative Abundance of Specific Picoeukaryotes

High-throughput sequencing (HTS) is revolutionizing environmental surveys of microbial diversity in the three domains of life by providing detailed information on which taxa are present in microbial assemblages. However, it is still unclear how the relative abundance of specific taxa gathered by HTS correlates with cell abundances. Here, we quantified the relative cell abundance of 6 picoeukaryotic taxa in 13 planktonic samples from 6 European coastal sites using epifluorescence microscopy on tyramide signal amplification-fluorescence in situ hybridization preparations. These relative abundance values were then compared with HTS data obtained in three separate molecular surveys: 454 sequencing of the V4 region of the 18S ribosomal DNA (rDNA) using DNA and RNA extracts (DNA-V4 and cDNA-V4) and Illumina sequencing of the V9 region (cDNA-V9). The microscopic and molecular signals were generally correlated, indicating that a relative increase in specific 18S rDNA was the result of a large proportion of cells in the given taxa. Despite these positive correlations, the slopes often deviated from 1, precluding a direct translation of sequences to cells. Our data highlighted clear differences depending on the nucleic acid template or the 18S rDNA region targeted. Thus, the molecular signal obtained using cDNA templates was always closer to relative cell abundances, while the V4 and V9 regions gave better results depending on the taxa. Our data support the quantitative use of HTS data but warn about considering it as a direct proxy of cell abundances.

IMPORTANCE

Direct studies on marine picoeukaryotes by epifluorescence microscopy are problematic due to the lack of morphological features and due to the limited number and poor resolution of specific phylogenetic probes used in fluorescence in situ hybridization (FISH) routines. As a consequence, there is an increasing use of molecular methods, including high-throughput sequencing (HTS), to study marine microbial diversity. HTS can provide a detailed picture of the taxa present in a community and can reveal diversity not evident using other methods, but it is still unclear what the meaning of the sequence abundance in a given taxon is. Our aim is to investigate the correspondence between the relative HTS signal and relative cell abundances in selected picoeukaryotic taxa. Environmental sequencing provides reasonable estimates of the relative abundance of specific taxa. Better results are obtained when using RNA extracts as the templates, while the region of 18S ribosomal DNA had different influences depending on the taxa assayed.

Evaluation of the effects of intrapartum antibiotic prophylaxis on newborn intestinal microbiota using a sequencing approach targeted to multi hypervariable 16S rDNA regions

Different factors are known to influence the early gut colonization in newborns, among them the perinatal use of antibiotics. On the other hand, the effect on the baby of the administration of antibiotics to the mother during labor, referred to as intrapartum antibiotic prophylaxis (IAP), has received less attention, although routinely used in group B Streptococcus positive women to prevent the infection in newborns. In this work, the fecal microbiota of neonates born to mothers receiving IAP and of control subjects were compared taking advantage for the first time of high-throughput DNA sequencing technology. Seven different 16S rDNA hypervariable regions (V2, V3, V4, V6 + V7, V8, and V9) were amplified and sequenced using the Ion Torrent Personal Genome Machine. The results obtained showed significant differences in the microbial composition of newborns born to mothers who had received IAP, with a lower abundance of Actinobacteria and Bacteroidetes as well as an overrepresentation of Proteobacteria. Considering that the seven hypervariable regions showed different discriminant ability in the taxonomic identification, further analyses were performed on the V4 region evidencing in IAP infants a reduced microbial richness and biodiversity, as well as a lower number of bacterial families with a predominance of Enterobacteriaceae members. In addition, this analysis pointed out a significant reduction in Bifidobacterium spp. strains. The reduced abundance of these beneficial microorganisms, together with the increased amount of potentially pathogenic bacteria, may suggest that IAP infants are more exposed to gastrointestinal or generally health disorders later in age.

Soil pH effects on the interactions between dissolved zinc, non-nano- and nano-ZnO with soil bacterial communities

Zinc oxide nanoparticles (ZnO NPs) are used in an array of products and processes, ranging from personal care products to antifouling paints, textiles, food additives, antibacterial agents and environmental remediation processes. Soils are an environment likely to be exposed to manmade nanoparticles due to the practice of applying sewage sludge as a fertiliser or as an organic soil improver. However, understanding on the interactions between soil properties, nanoparticles and the organisms that live within soil is lacking, especially with regards to soil bacterial communities. We studied the effects of nanoparticulate, non-nanoparticulate and ionic zinc (in the form of zinc chloride) on the composition of bacterial communities in soil with a modified pH range (from pH 4.5 to pH 7.2). We observed strong pH-dependent effects on the interaction between bacterial communities and all forms of zinc, with the largest changes in bacterial community composition occurring in soils with low and medium pH levels (pH 4.8 and 5.9). The high pH soil (pH 7.2) was less susceptible to the effects of zinc exposure. At the highest doses of zinc (2500 mg/kg dw soil), both nano and non-nano particulate zinc applications elicited a similar response in the soil bacterial community, and this differed significantly to the ionic zinc salt treatment. The results highlight the importance of considering soil pH in nanotoxicology studies, although further work is needed to determine the exact mechanisms controlling the toxicity and fate and interactions of nanoparticles with soil microbial communities.

Characterization of successional changes in bacterial community composition during bioremediation of used motor oil-contaminated soil in a boreal climate

The widespread use of motor oil makes it a notable risk factor to cause scattered contamination in soil. The monitoring of microbial community dynamics can serve as a comprehensive tool to assess the ecological impact of contaminants and their disappearance in the ecosystem. Hence, a field study was conducted to monitor the ecological impact of used motor oil under different perennial cropping systems (fodder galega, brome grass, galega-brome grass mixture and bare fallow) in a boreal climate zone. Length heterogeneity PCR characterized a successional pattern in bacterial community following oil contamination over a four-year bioremediation period. Soil pH and electrical conductivity were associated with the shifts in bacterial community composition. Crops had no detectable effect on bacterial community composition or complexity. However, the legume fodder galega increased soil microbial biomass, expressed as soil total DNA. Oil contamination induced an abrupt change in bacterial community composition at the early stage, yet the effect did not last as long as the oil in soil. The successional variation in bacterial community composition can serve as a sensitive ecological indicator of oil contamination and remediation in situ.

A flexible and economical barcoding approach for highly multiplexed amplicon sequencing of diverse target genes

High throughput sequencing of phylogenetic and functional gene amplicons provides tremendous insight into the structure and functional potential of complex microbial communities. Here, we introduce a highly adaptable and economical PCR approach to barcoding and pooling libraries of numerous target genes. In this approach, we replace gene- and sequencing platform-specific fusion primers with general, interchangeable barcoding primers, enabling nearly limitless customized barcode-primer combinations. Compared to barcoding with long fusion primers, our multiple-target gene approach is more economical because it overall requires lower number of primers and is based on short primers with generally lower synthesis and purification costs. To highlight our approach, we pooled over 900 different small-subunit rRNA and functional gene amplicon libraries obtained from various environmental or host-associated microbial community samples into a single, paired-end Illumina MiSeq run. Although the amplicon regions ranged in size from approximately 290 to 720 bp, we found no significant systematic sequencing bias related to amplicon length or gene target. Our results indicate that this flexible multiplexing approach produces large, diverse, and high quality sets of amplicon sequence data for modern studies in microbial ecology.

Temporal Study of the Microbial Diversity of the North Arm of Great Salt Lake, Utah, U.S

We employed a temporal sampling approach to understand how the microbial diversity may shift in the north arm of Great Salt Lake, Utah, U.S. To determine how variations in seasonal environmental factors affect microbial communities, length heterogeneity PCR fingerprinting was performed using consensus primers for the domain Bacteria, and the haloarchaea. The archaeal fingerprints showed similarities during 2003 and 2004, but this diversity changed during the remaining two years of the study, 2005 and 2006. We also performed molecular phylogenetic analysis of the 16S rRNA genes of the whole microbial community to characterize the taxa in the samples. Our results indicated that in the domain, Bacteria, the Salinibacter group dominated the populations in all samplings. However, in the case of Archaea, as noted by LIBSHUFF for phylogenetic relatedness analysis, many of the temporal communities were distinct from each other, and changes in community composition did not track with environmental parameters. Around 20-23 different phylotypes, as revealed by rarefaction, predominated at different periods of the year. Some phylotypes, such as Haloquadradum, were present year-round although they changed in their abundance in different samplings, which may indicate that these species are affected by biotic factors, such as nutrients or viruses, that are independent of seasonal temperature dynamics.

A Multi Size-Level Assessment of Benthic Marine Communities in a Coastal Environment: Are They Different Sides of the Same Coin?

Organism body size has been demonstrated to be a discriminating element in shaping the response of living beings to environmental factors, thus playing a fundamental role in community structuring. Despite the importance of studies elucidating relations among communities of different size levels in ecosystems, the attempts that have been made in this sense are still very scarce and a reliable approach for these research still has to be defined. We characterized the benthic communities of bacteria, microbial eukaryotes, meiofauna and macrofauna in a coastal environment, encompassing a 10000-fold gradient in body size, testing and discussing a mixed approach of molecular fingerprinting for microbes and morphological observations for meio- and macrofauna. We found no correlation among structures of the different size-level communities: this suggests that community composition at one size-level could have no (or very low) influence on the community composition at other size-levels. Moreover, each community responds in a different way to the environmental parameters and with a degree of sensitivity which seems to increase with organism size. Therefore, our data indicate that the characterization of all the different size levels is clearly a necessity in order to study the dynamics really acting in a system.

The Single-Cell Lab or How to Perform Single-Cell Molecular Analysis

The increasing interest towards cellular heterogeneity within cell populations has pushed the development of new protocols to isolate and analyze single cells. PCR-based amplification techniques are widely used in this field. However, setting up an experiment and analyzing the results can sometimes be challenging. The aim of this chapter is to provide a general overview on single-cell PCR analysis focusing on the potential pitfalls and on the possible solutions to successfully perform the analysis.

Marine bacterial communities are resistant to elevated carbon dioxide levels

It is well established that the release of anthropogenic-derived CO2 into the atmosphere will be mainly absorbed by the oceans, with a concomitant drop in pH, a process termed ocean acidification. As such, there is considerable interest in how changes in increased CO2 and lower pH will affect marine biota, such as bacteria, which play central roles in oceanic biogeochemical processes. Set within an ecological framework, we investigated the direct effects of elevated CO2, contrasted with ambient conditions on the resistance and resilience of marine bacterial communities in a replicated temporal seawater mesocosm experiment. The results of the study strongly indicate that marine bacterial communities are highly resistant to the elevated CO2 and lower pH conditions imposed, as demonstrated from measures of turnover using taxa–time relationships and distance–decay relationships. In addition, no significant differences in community abundance, structure or composition were observed. Our results suggest that there are no direct effects on marine bacterial communities and that the bacterial fraction of microbial plankton holds enough flexibility and evolutionary capacity to withstand predicted future changes from elevated CO2 and subsequent ocean acidification.

Bacterial and archaeal communities in long-term contaminated surface and subsurface soil evaluated through coextracted RNA and DNA

Soil RNA and DNA were coextracted along a contamination gradient at a landfarming field with aged crude oil contamination to investigate pollution-dependent differences in 16S rRNA and rRNA gene pools. Microbial biomass correlated with nucleic acid yields as well as bacterial community change, indicating that the same factors controlled community size and structure. In surface soil, bacterial community evenness, estimated through length heterogeneity PCR (LH-PCR) fingerprinting, appeared higher for RNA-based than for DNA-based communities. The RNA-based community profiles resembled the DNA-based communities of soil with a lower contamination level. Cloning-based identification of bacterial hydrocarbon-degrading taxa in the RNA pool, representing the viable community with high protein synthesis potential, indicated that decontamination processes still continue. Analyses of archaea revealed that only Thaumarchaeota were present in the aerobic samples, whereas more diverse communities were found in the compacted subsurface soil with more crude oil. For subsurface bacteria, hydrocarbon concentration explained neither the community structure nor the difference between RNA-based and DNA-based communities. However, rRNA of bacterial taxa associated with syntrophic and sulphate-reducing alkane degradation was detected. Although the same prokaryotic taxa were identified in DNA and RNA, comparison of the two nucleic acid pools can aid in the assessment of past and future restoration success.

Soil mesocosm studies on atrazine bioremediation

Accumulation of pesticides in the environment causes serious issues of contamination and toxicity. Bioremediation is an ecologically sound method to manage soil pollution, but the bottleneck here, is the successful scale-up of lab-scale experiments to field applications. This study demonstrates pilot-scale bioremediation in tropical soil using atrazine as model pollutant. Mimicking field conditions, three different bioremediation strategies for atrazine degradation were explored. 100 kg soil mesocosms were set-up, with or without atrazine application history. Natural attenuation and enhanced bioremediation were tested, where augmentation with an atrazine degrading consortium demonstrated best pollutant removal. 90% atrazine degradation was observed in six days in soil previously exposed to atrazine, while soil without history of atrazine use, needed 15 days to remove the same amount of amended atrazine. The bacterial consortium comprised of 3 novel bacterial strains with different genetic atrazine degrading potential. The progress of bioremediation was monitored by measuring the levels of atrazine and its intermediate, cyanuric acid. Genes from the atrazine degradation pathway, namely, atzA, atzB, atzD, trzN and trzD were quantified in all mesocosms for 60 days. The highest abundance of all target genes was observed on the 6th day of treatment. trzD was observed in the bioaugmented mesocosms only. The bacterial community profile in all mesocosms was monitored by LH-PCR over a period of two months. Results indicate that the communities changed rapidly after inoculation, but there was no drastic change in microbial community profile after 1 month. Results indicated that efficient bioremediation of atrazine using a microbial consortium could be successfully up-scaled to pilot scale.

The bias associated with amplicon sequencing does not affect the quantitative assessment of bacterial community dynamics

The performance of two sets of primers targeting variable regions of the 16S rRNA gene V1–V3 and V4 was compared in their ability to describe changes of bacterial diversity and temporal turnover in full-scale activated sludge. Duplicate sets of high-throughput amplicon sequencing data of the two 16S rRNA regions shared a collection of core taxa that were observed across a series of twelve monthly samples, although the relative abundance of each taxon was substantially different between regions. A case in point was the changes in the relative abundance of filamentous bacteria Thiothrix, which caused a large effect on diversity indices, but only in the V1–V3 data set. Yet the relative abundance of Thiothrix in the amplicon sequencing data from both regions correlated with the estimation of its abundance determined using fluorescence in situ hybridization. In nonmetric multidimensional analysis samples were distributed along the first ordination axis according to the sequenced region rather than according to sample identities. The dynamics of microbial communities indicated that V1–V3 and the V4 regions of the 16S rRNA gene yielded comparable patterns of: 1) the changes occurring within the communities along fixed time intervals, 2) the slow turnover of activated sludge communities and 3) the rate of species replacement calculated from the taxa–time relationships. The temperature was the only operational variable that showed significant correlation with the composition of bacterial communities over time for the sets of data obtained with both pairs of primers. In conclusion, we show that despite the bias introduced by amplicon sequencing, the variable regions V1–V3 and V4 can be confidently used for the quantitative assessment of bacterial community dynamics, and provide a proper qualitative account of general taxa in the community, especially when the data are obtained over a convenient time window rather than at a single time point.

Use of length heterogeneity polymerase chain reaction (LH-PCR) as non-invasive approach for dietary analysis of Svalbard reindeer, Rangifer tarandus platyrhynchus

To efficiently investigate the forage preference of Svalbard reindeer (Rangifer tarandus platyrhynchus), we applied length-heterogeneity polymerase chain reaction (LH-PCR) based on length differences of internal transcribed spacer (ITS) regions of ribosomal RNA (rRNA) to fecal samples from R. tarandus platyrhynchus. A length-heterogeneity (LH) database was constructed using both collected potential food sources of Svalbard reindeer and fecal samples, followed by PCR, cloning and sequencing. In total, eighteen fecal samples were collected between 2011 and 2012 from 2 geographic regions and 15 samples were successfully amplified by PCR. The LH-PCR analysis detected abundant peaks, 18.6 peaks on an average per sample, ranging from 100 to 500 bp in size and showing distinct patterns associated with both regions and years of sample collection. Principal component analysis (PCA) resulted in clustering of 15 fecal samples into 3 groups by the year of collection and region with a statistically significant difference at 99.9% level. The first 2 principal components (PCs) explained 71.1% of the total variation among the samples. Through comparison with LH database and identification by cloning and sequencing, lichens (Stereocaulon sp. and Ochrolechia sp.) and plant species (Salix polaris and Saxifraga oppositifolia) were detected as the food sources that contributed most to the Svalbard reindeer diet. Our results suggest that the use of LH-PCR analysis would be a non-invasive and efficient monitoring tool for characterizing the foraging strategy of Svalbard reindeer. Additionally, combining sequence information would increase its resolving power in identification of foraged diet components.

Evaluating detection limits of next-generation sequencing for the surveillance and monitoring of international marine pests

Most surveillance programmes for marine invasive species (MIS) require considerable taxonomic expertise, are laborious, and are unable to identify species at larval or juvenile stages. Therefore, marine pests may go undetected at the initial stages of incursions when population densities are low. In this study, we evaluated the ability of the benchtop GS Junior™ 454 pyrosequencing system to detect the presence of MIS in complex sample matrices. An initial in-silico evaluation of the mitochondrial cytochrome c oxidase subunit I (COI) and the nuclear small subunit ribosomal DNA (SSU) genes, found that multiple primer sets (targeting a ca. 400 base pair region) would be required to obtain species level identification within the COI gene. In contrast a single universal primer set was designed to target the V1-V3 region of SSU, allowing simultaneous PCR amplification of a wide taxonomic range of MIS. To evaluate the limits of detection of this method, artificial contrived communities (10 species from 5 taxonomic groups) were created using varying concentrations of known DNA samples and PCR products. Environmental samples (water and sediment) spiked with one or five 160 hr old Asterias amurensis larvae were also examined. Pyrosequencing was able to recover DNA/PCR products of individual species present at greater than 0.64% abundance from all tested contrived communities. Additionally, single A. amurensis larvae were detected from both water and sediment samples despite the co-occurrence of a large array of environmental eukaryotes, indicating an equivalent sensitivity to quantitative PCR. NGS technology has tremendous potential for the early detection of marine invasive species worldwide.

Using boreholes as windows into groundwater ecosystems

Groundwater ecosystems remain poorly understood yet may provide ecosystem services, make a unique contribution to biodiversity and contain useful bio-indicators of water quality. Little is known about ecosystem variability, the distribution of invertebrates within aquifers, or how representative boreholes are of aquifers. We addressed these issues using borehole imaging and single borehole dilution tests to identify three potential aquifer habitats (fractures, fissures or conduits) intercepted by two Chalk boreholes at different depths beneath the surface (34 to 98 m). These habitats were characterised by sampling the invertebrates, microbiology and hydrochemistry using a packer system to isolate them. Samples were taken with progressively increasing pumped volume to assess differences between borehole and aquifer communities. The study provides a new conceptual framework to infer the origin of water, invertebrates and microbes sampled from boreholes. It demonstrates that pumping 5 m³ at 0.4–1.8 l/sec was sufficient to entrain invertebrates from five to tens of metres into the aquifer during these packer tests. Invertebrates and bacteria were more abundant in the boreholes than in the aquifer, with associated water chemistry variations indicating that boreholes act as sites of enhanced biogeochemical cycling. There was some variability in invertebrate abundance and bacterial community structure between habitats, indicating ecological heterogeneity within the aquifer. However, invertebrates were captured in all aquifer samples, and bacterial abundance, major ion chemistry and dissolved oxygen remained similar. Therefore the study demonstrates that in the Chalk, ecosystems comprising bacteria and invertebrates extend from around the water table to 70 m below it. Hydrogeological techniques provide excellent scope for tackling outstanding questions in groundwater ecology, provided an appropriate conceptual hydrogeological understanding is applied.

Corn silage in dairy cow diets to reduce ruminal methanogenesis: effects on the rumen metabolically active microbial communities

Methane produced by the methanogenic Archaea that inhabit the rumen is a potent greenhouse gas and represents an energy loss for the animal. Although several strategies have been proposed to mitigate enteric CH4 production, little is known about the effects of dietary changes on the microbial consortia involved in ruminal methanogenesis. Thus, the current study aimed to examine how the metabolically active microbes are affected when dairy cows were fed diets with increasing proportions of corn silage (CS). For this purpose, 9 ruminally cannulated lactating dairy cows were used in a replicated 3 × 3 Latin square design and fed a total mixed ration (60:40 forage:concentrate ratio on a dry matter basis) with the forage portion being either alfalfa silage (0% CS), corn silage (100% CS), or a 50:50 mixture (50% CS). Enteric CH4 production was determined using respiration chambers and total rumen content was sampled for the determination of fermentation characteristics and molecular biology analyses (cDNA-based length heterogeneity PCR, quantitative PCR). The cDNA-based length heterogeneity PCR targeting active microbes revealed similar bacterial communities in cows fed 0% CS and 50% CS diets, whereas important differences were observed between 0% CS and 100% CS diets, including a reduction in the bacterial richness and diversity in cows fed 100% CS diet. As revealed by quantitative PCR, feeding the 100% CS diet increased the number of total bacteria, Prevotella spp., Archaea, and methanogenic activity, though it reduced protozoal number. Meanwhile, increasing the CS proportion in the diet increased propionate concentration but decreased ruminal pH, CH4 production (L/kg of dry matter intake), and concentrations of acetate and butyrate. Based on these microbial and fermentation changes, and because CH4 production was reduced by feeding 100% CS diet, this study shows that the use of cDNA-based quantitative PCR to estimate archaeal growth and activity is not reliable enough to reflect changes in ruminal methanogenesis. A more robust technique to characterize changes in archaeal community structures will help to better understand the microbial process involved in ruminal methanogenesis and, hence, enabling the development of more effective dietary CH4 mitigation strategies.

Effect of rice cultivation systems on indigenous arbuscular mycorrhizal fungal community structure

Arbuscular mycorrhizal fungi (AMF) in an agricultural ecosystem are necessary for proper management of beneficial symbiosis. Here we explored how the patterns of the AMF community in rice roots were affected by rice cultivation systems (the system of rice intensification [SRI] and the conventional rice cultivation system [CS]), and by compost application during growth stages. Rice plants harvested from SRI-managed plots exhibited considerably higher total biomass, root dry weight, and seed fill than those obtained from conventionally managed plots. Our findings revealed that all AMF sequences observed from CS plots belonged (only) to the genus Glomus, colonizing in rice roots grown under this type of cultivation, while rice roots sown in SRI showed sequences belonging to both Glomus and Acaulospora. The AMF community was compared between the different cultivation types (CS and SRI) and compost applications by principle component analysis. In all rice growth stages, AMF assemblages of CS management were not separated from those of SRI management. The distribution of AMF community composition based on T-RFLP data showed that the AMF community structure was different among four cultivation systems, and there was a gradual increase of Shannon-Weaver indices of diversity (H') of the AMF community under SRI during growth stages. The results of this research indicated that rice grown in SRI-managed plots had more diverse AMF communities than those grown in CS plots.

Quality control in multi-tag pyrosequencing of microbial communities

Multi-tag pyrosequencing has become a key method in the analysis of microbial community composition. However, it is well known that kinetic bias during the initial PCR amplification of such microbial communities can dramatically distort amplicon abundance prior to downstream emulsion PCR and pyrosequencing. Here we present a simple protocol combining length-heterogeneity PCR fingerprinting with pyrosequencing to ensure the linearity of microbial community amplification. The method employs a fluorescently labeled reverse primer along with multi-tagged forward primers to initially amplify the microbial community. The resulting labeled amplicons are then fingerprinted, purified, and quantitated prior to emulsion PCR and pyrosequencing. Our data demonstrates: (i) use of this protocol results in a distribution of sequences showing linear amplification following emulsion PCR when compared with the initial length-heterogeneity PCR fingerprints, and (ii) that the added tags and labels do not have a negative effect on overall microbial community profiles.

Columnaris disease in fish: a review with emphasis on bacterium-host interactions

Flavobacterium columnare (F. columnare) is the causative agent of columnaris disease. This bacterium affects both cultured and wild freshwater fish including many susceptible commercially important fish species. F. columnare infections may result in skin lesions, fin erosion and gill necrosis, with a high degree of mortality, leading to severe economic losses. Especially in the last decade, various research groups have performed studies aimed at elucidating the pathogenesis of columnaris disease, leading to significant progress in defining the complex interactions between the organism and its host. Despite these efforts, the pathogenesis of columnaris disease hitherto largely remains unclear, compromising the further development of efficient curative and preventive measures to combat this disease. Besides elaborating on the agent and the disease it causes, this review aims to summarize these pathogenesis data emphasizing the areas meriting further investigation.

Diet-induced alterations in total and metabolically active microbes within the rumen of dairy cows

DNA-based techniques are widely used to study microbial populations; however, this approach is not specific to active microbes, because DNA may originate from inactive and/or dead cells. Using cDNA and DNA, respectively, we aimed to discriminate the active microbes from the total microbial community within the rumen of dairy cows fed diets with increasing proportions of corn silage (CS). Nine multiparous lactating Holstein cows fitted with ruminal cannulas were used in a replicated 3×3 Latin square (32-d period; 21-d adaptation) design to investigate diet-induced shifts in microbial populations by targeting the rDNA gene. Cows were fed a total mixed ration with the forage portion being either barley silage (0% CS), a 50∶50 mixture of barley silage and corn silage (50% CS), or corn silage (100% CS). No differences were found for total microbes analyzed by quantitative PCR, but changes were observed within the active ones. Feeding more CS to dairy cows was accompanied by an increase in Prevotella rRNA transcripts (P = 0.10) and a decrease in the protozoal rRNA transcripts (P<0.05). Although they were distributed differently among diets, 78% of the amplicons detected in DNA- and cDNA-based fingerprints were common to total and active bacterial communities. These may represent a bacterial core of abundant and active cells that drive the fermentation processes. In contrast, 10% of amplicons were specific to total bacteria and may represent inactive or dead cells, whereas 12% were only found within the active bacterial community and may constitute slow-growing bacteria with high metabolic activity. It appears that cDNA-based analysis is more discriminative to identify diet-induced shifts within the microbial community. This approach allows the detection of diet-induced changes in the microbial populations as well as particular bacterial amplicons that remained undetected using DNA-based methods.

New developments in the study of the microbiota of raw-milk, long-ripened cheeses by molecular methods: the case of Grana Padano and Parmigiano Reggiano

Microorganisms are an essential component of cheeses and play important roles during both cheese manufacture and ripening. Both starter and secondary flora modify the physical and chemical properties of cheese, contributing and reacting to changes that occur during the manufacture and ripening of cheese. As the composition of microbial population changes under the influence of continuous shifts in environmental conditions and microorganisms interactions during manufacturing and ripening, the characteristics of a given cheese depend also on microflora dynamics. The microbiota present in cheese is complex and its growth and activity represent the most important, but the least controllable steps. In the past, research in this area was dependent on classical microbiological techniques. However, culture-dependent methods are time-consuming and approaches that include a culturing step can lead to inaccuracies due to species present in low numbers or simply uncultivable. Therefore, they cannot be used as a unique tool to monitor community dynamics. For these reasons approaches to cheese microbiology had to change dramatically. To address this, in recent years the focus on the use of culture-independent methods based on the direct analysis of DNA (or RNA) has rapidly increased. Application of such techniques to the study of cheese microbiology represents a rapid, sound, reliable, and effective way for the detection and identification of the microorganisms present in dairy products, leading to major advances in understanding this complex microbial ecosystem and its impact on cheese ripening and quality. In this article, an overview on the recent advances in the use of molecular methods for thorough analysis of microbial communities in cheeses is given. Furthermore, applications of culture-independent approaches to study the microbiology of two important raw-milk, long-ripened cheeses such as Grana Padano and Parmigiano Reggiano, are presented.

Development of a multiplex real time PCR to detect thermophilic lactic acid bacteria in natural whey starters

A multiplex real time PCR (mRealT-PCR) useful to rapidly screen microbial composition of thermophilic starter cultures for hard cooked cheeses and to compare samples with potentially different technological properties was developed. Novel primers directed toward pheS gene were designed and optimized for multiple detection of Lactobacillus helveticus, Lactobacillus delbrueckii, Streptococcus thermophilus and Lactobacillus fermentum. The assay was based on SYBR Green chemistry followed by melting curves analysis. The method was then evaluated for applications in the specific detection of the 4 lactic acid bacteria (LAB) in 29 different natural whey starters for Parmigiano Reggiano cheese production. The results obtained by mRealT-PCR were also compared with those obtained on the same samples by Fluorescence in Situ Hybridization (FISH) and Length-Heterogeneity PCR (LH-PCR). The mRealT-PCR developed in this study, was found to be effective for analyzing species present in the samples with an average sensitivity down to less than 600 copies of DNA and therefore sensitive enough to detect even minor LAB community members of thermophilic starter cultures. The assay was able to describe the microbial population of all the different natural whey starter samples analyzed, despite their natural variability. A higher number of whey starter samples with S. thermophilus and L. fermentum present in their microbial community were revealed, suggesting that these species could be more frequent in Parmigiano Reggiano natural whey starter samples than previously shown. The method was more effective than LH-PCR and FISH and, considering that these two techniques have to be used in combination to detect the less abundant species, the mRealT-PCR was also faster. Providing a single step sensitive detection of L. helveticus, L. delbrueckii, S. thermophilus and L. fermentum, the developed mRealT-PCR could be used for screening thermophilic starter cultures and to follow the presence of those species during ripening of derived dairy products. A major increase in understanding the starter culture contribution to cheese ecosystem could be harnessed to control cheese ripening and flavor formation.

Archaeal Diversity of Upland Rice Field Soils Assessed by the Terminal Restriction Fragment Length Polymorphism Method Combined with Real Time Quantitative-PCR and a Clone Library Analysis

The PCR amplification-based analysis of microbial diversity is subject to potential problems. In this study, to minimize the bias toward a 1:1 ratio in multitemplate PCR, a real-time PCR assay was carried out using a quenching fluorescence dye primer and amplification efficiency was monitored. Then terminal-restriction fragment length polymorphism (T-RFLP) profiling was performed using the PCR product with minimized PCR bias. This method was applied to an analysis of the diversity of the archaeal community in an upland rice field under different tillage systems and winter cover cropping. Terminal restriction fragments (T-RFs) of PCR-amplified archaeal 16S rRNA genes were assigned to the gene sequences recovered from the same soil by using an archaeal 16S rRNA gene clone library. Our results indicated that soil archaeal members were not influenced but the relative abundance of archaeal species particularly those belonging to Crenarchaeota which changed between the tillage and non-tillage treatments.

Groundtruthing next-gen sequencing for microbial ecology-biases and errors in community structure estimates from PCR amplicon pyrosequencing

Analysis of microbial communities by high-throughput pyrosequencing of SSU rRNA gene PCR amplicons has transformed microbial ecology research and led to the observation that many communities contain a diverse assortment of rare taxa-a phenomenon termed the Rare Biosphere. Multiple studies have investigated the effect of pyrosequencing read quality on operational taxonomic unit (OTU) richness for contrived communities, yet there is limited information on the fidelity of community structure estimates obtained through this approach. Given that PCR biases are widely recognized, and further unknown biases may arise from the sequencing process itself, a priori assumptions about the neutrality of the data generation process are at best unvalidated. Furthermore, post-sequencing quality control algorithms have not been explicitly evaluated for the accuracy of recovered representative sequences and its impact on downstream analyses, reducing useful discussion on pyrosequencing reads to their diversity and abundances. Here we report on community structures and sequences recovered for in vitro-simulated communities consisting of twenty 16S rRNA gene clones tiered at known proportions. PCR amplicon libraries of the V3-V4 and V6 hypervariable regions from the in vitro-simulated communities were sequenced using the Roche 454 GS FLX Titanium platform. Commonly used quality control protocols resulted in the formation of OTUs with >1% abundance composed entirely of erroneous sequences, while over-aggressive clustering approaches obfuscated real, expected OTUs. The pyrosequencing process itself did not appear to impose significant biases on overall community structure estimates, although the detection limit for rare taxa may be affected by PCR amplicon size and quality control approach employed. Meanwhile, PCR biases associated with the initial amplicon generation may impose greater distortions in the observed community structure.

Bacterial and phytoplankton responses to nutrient amendments in a boreal lake differ according to season and to taxonomic resolution

Nutrient limitation and resource competition in bacterial and phytoplankton communities may appear different when considering different levels of taxonomic resolution. Nutrient amendment experiments conducted in a boreal lake on three occasions during one open water season revealed complex responses in overall bacterioplankton and phytoplankton abundance and biovolume. In general, bacteria were dominant in spring, while phytoplankton was clearly the predominant group in autumn. Seasonal differences in the community composition of bacteria and phytoplankton were mainly related to changes in observed taxa, while the differences across nutrient treatments within an experiment were due to changes in relative contributions of certain higher- and lower-level phylogenetic groups. Of the main bacterioplankton phyla, only Actinobacteria had a treatment response that was visible even at the phylum level throughout the season. With increasing resolution (from 75 to 99% sequence similarity) major responses to nutrient amendments appeared using 454 pyrosequencing data of 16S rRNA amplicons. This further revealed that OTUs (defined by 97% sequence similarity) annotated to the same highly resolved freshwater groups appeared to occur during different seasons and were showing treatment-dependent differentiation, indicating that OTUs within these groups were not ecologically coherent. Similarly, phytoplankton species from the same genera responded differently to nutrient amendments even though biovolumes of the majority of taxa increased when both nitrogen and phosphorus were added simultaneously. The bacterioplankton and phytoplankton community compositions showed concurrent trajectories that could be seen in synchronous succession patterns over the season. Overall, our data revealed that the response of both communities to nutrient changes was highly dependent on season and that contradictory results may be obtained when using different taxonomic resolutions.

Persistence and degrading activity of free and immobilised allochthonous bacteria during bioremediation of hydrocarbon-contaminated soils

Rhodococcus sp. and Pseudomonas sp. bioremediation experiments were carried out using free and immobilized cells on natural carrier material (corncob powder) in order to evaluate the feasibility of its use in the bioremediation of hydrocarbon-contaminated soils. Terminal restriction fragment length polymorphism analysis was performed on the 16S rRNA gene as molecular fingerprinting method in order to assess the persistence of inoculated strains in the soil over time. Immobilized Pseudomonas cells degraded hydrocarbons more efficiently in the short term compared to the free ones. Immobilization seemed also to increase cell growth and stability in the soil. Free and immobilized Rhodococcus cells showed comparable degradation percentages, probably due to the peculiarity of Rhodococcus cells to aggregate into irregular clusters in the presence of hydrocarbons as sole carbon source. It is likely that the cells were not properly adsorbed on the porous matrix as a result of the small size of its pores. When Rhodococcus and Pseudomonas cells were co-immobilized on the matrix, a competition established between the two strains, that probably ended in the exclusion of Pseudomonas cells from the pores. The organic matrix might act as protective agent, but it also possibly limited cell density. Nevertheless, when the cells were properly adsorbed on the porous matrix, the immobilization became a suitable bioremediation strategy.

Rapid microbial response to the presence of an ancient relic in the Antarctic Dry Valleys

The extreme cold and aridity of the Antarctic McMurdo Dry Valleys have led to the longstanding belief that metabolic rates of soil microbiota are negligible, and that ecosystem changes take place over millennia. Here we report the first direct experimental evidence that soil microbial communities undergo rapid and lasting changes in response to contemporary environmental conditions. Mummified seals, curious natural artifacts found scattered throughout Dry Valleys, alter their underlying soil environment by stabilizing temperatures, elevating relative humidity and reducing ultraviolet exposure. In a unique, multi-year mummified seal transplantation experiment, we found that endemic Dry Valley microbial communities responded to these changes within 3 years, resulting in a sevenfold increase in CO(2) flux and a significant reduction in biodiversity. These findings challenge prevailing ideas about Antarctic Dry Valley ecosystems and indicate that current and future environmental conditions may strongly influence the ecology of the dominant biota in the Dry Valleys.

Analysis of the salivary microbiome using culture-independent techniques

Background

The salivary microbiota is a potential diagnostic indicator of several diseases. Culture-independent techniques are required to study the salivary microbial community since many of its members have not been cultivated.

Methods

We explored the bacterial community composition in the saliva sample using metagenomic whole genome shotgun (WGS) sequencing, the extraction of 16S rRNA gene fragments from metagenomic sequences (16S-WGS) and high-throughput sequencing of PCR-amplified bacterial 16S rDNA gene (16S-HTS) regions V1 and V3.

Results

The hierarchical clustering of data based on the relative abundance of bacterial genera revealed that distances between 16S-HTS datasets for V1 and V3 regions were greater than those obtained for the same V region with different numbers of PCR cycles. Datasets generated by 16S-HTS and 16S-WGS were even more distant. Finally, comparison of WGS and 16S-based datasets revealed the highest dissimilarity.

The analysis of the 16S-HTS, WGS and 16S-WGS datasets revealed 206, 56 and 39 bacterial genera, respectively, 124 of which have not been previously identified in salivary microbiomes. A large fraction of DNA extracted from saliva corresponded to human DNA. Based on sequence similarity search against completely sequenced genomes, bacterial and viral sequences represented 0.73% and 0.0036% of the salivary metagenome, respectively. Several sequence reads were identified as parts of the human herpesvirus 7.

Conclusions

Analysis of the salivary metagenome may have implications in diagnostics e.g. in detection of microorganisms and viruses without designing specific tests for each pathogen.