Bacterial Profiling and Dynamic Succession Analysis of Phlebopus portentosus Casing Soil Using MiSeq Sequencing

Phlebopusportentosus (Berk. and Broome) Boedijin is a popular edible mushroom found in China and Thailand. To date, P. portentosus is the only species in the order Boletales that can be successfully cultivated worldwide. The use of a casing layer or casing soil overlaying the substrate is a crucial step in the production of this mushroom. In this study, bacterial profiling and dynamic succession analyses of casing soil during the cultivation of P. portentosus were performed. One hundred and fifty samples were collected, and MiSeq sequencing of the V3-V4 region of the 16S rRNA gene was conducted. After performing a decontamination procedure, only 38 samples were retained, including 6 casing soil-originated samples (OS), 6 casing soil samples (FHCS) and 5 upper substrate samples (FHCU) from the period of complete colonization by mycelia; 6 casing soil samples (PCS) and 5 upper substrate samples (PCU) from the primordium period; and 6 casing soil samples (FCS) and 4 upper substrate samples (FCU) from fruit body period. The results revealed that bacterial diversity increased sharply from the hyphal to the primordium stage and then decreased during harvesting. The non-metric multidimensional scaling (NMDS) ordination and analysis of similarities (ANOSIM) analysis suggested that the community composition during different stages was significantly different in casing soil. The most abundant phyla in all of the samples were Proteobacteria, Chloroflexi, Acidobacteria, Actinobacteria, Saccharibacteria, and Bacteroidetes. Burkholderia was the most abundant genus in all the samples except the OS samples. The relative abundance of Burkholderia in the FHCS samples (55.79%) decreased to 35.14% in the PCS samples and then increased to 45.60% in the FCS samples. The abundances of Acidobacterium, Rhizobium, Acidisphaera, Bradyrhizobium, and Bacillus increased from the FHCS to PCS samples. The linear discriminant analysis (LDA) effect size (LEfSe) suggested that Acidobacterium and Acidisphaera are micromarkers for PCS, whereas Bradyrhizobium, Roseiarcus, and Pseudolabrys were associated with fruit body stages. The network analyses resulted in 23 edges, including 4 negative and 19 positive edges. Extensive mutualistic interactions may occur among casing soil bacteria. Furthermore, these bacteria play important roles in mycelial elongation, primordium formations, and the production of increased yields.

Fouling Microbial Communities on Plastics Compared with Wood and Steel: Are They Substrate- or Location-Specific?

Although marine biofouling has been widely studied on different substrates, information on biofouling on plastics in the Arabian Gulf is limited. Substrate- and location-specific effects were investigated by comparing the microbial communities developed on polyethylene terephthalate (PET) and polyethylene (PE) with those on steel and wood, at two locations in the Sea of Oman. Total biomass was lower on PET and PE than on steel and wood. PET had the highest bacterial abundance at both locations, whereas chlorophyll a concentrations did not vary between substrates. MiSeq 16S ribosomal RNA sequencing revealed comparable operational taxonomic unit (OTU) richness on all substrates at one location but lower numbers on PET and PE at the other location. Non-metric multidimensional scaling (NMDS) showed distinct clusters of the bacterial communities based on substrate (analysis of similarity (ANOSIM), R = 0.45-0.97, p < 0.03) and location (ANOSIM, R = 0.56, p < 0.0001). The bacterial genera Microcystis and Hydrogenophaga and the diatoms Licmophora and Mastogloia were specifically detected on plastics. Desulfovibrio and Pseudomonas spp. exhibited their highest abundance on steel and Corynebacterium spp. on wood. Scanning electron microscopy (SEM) revealed fissure formation on PET and PE, indicating physical degradation. The presence of free radicals on PET and carbonyl bonds (C=O) on PE, as revealed by Fourier transform infrared (FTIR) spectroscopy, indicated abiotic degradation while hydroxyl groups and spectral peaks for proteins and polysaccharides on PE indicated biotic degradation. We conclude that fouling microbial communities are not only substrate-specific but also location-specific and microbes developing on plastics could potentially contribute to their degradation in the marine environment.

Genetic diversity and quantification of human mastadenoviruses in wastewater from Sydney and Melbourne, Australia

Human mastadenoviruses (HAdVs) are DNA viruses that can cause a wide range of clinical diseases, including gastroenteritis, respiratory illnesses, conjunctivitis, and in more severe cases hepatitis, pancreatitis and disseminated diseases. HAdV infections are generally asymptomatic or self-limiting, but can cause adverse outcomes within vulnerable populations. Since most HAdV serotypes replicate within the human gastrointestinal tract, high levels of HAdV DNA are excreted into wastewater systems. In this study, we identified the genetic diversity of HAdV at a population level using wastewater samples collected from Sydney and Melbourne from 2016 to 2017, with the use of next generation sequencing (NGS) technologies. In addition, HAdV DNA levels were quantified using quantitative polymerase chain reaction (qPCR) based methods to better understand the health risks involved if wastewater contamination occurs. An average of 1.8 × 10⁷ genome copies of HAdV DNA was detected in one litre of wastewater collected in Sydney and Melbourne, over the two-year study period. A total of six major groups of HAdV were identified in wastewater samples using MiSeq, which included 19 different serotypes. Of those, the most prevalent was F41 (83.5%), followed by F40 (11.0%) and A31 (3.7%). In contrast, five groups of HAdV were identified in clinical samples with F41 as the most dominant serotype, (52.5% of gastroenteritis cases), followed by C1 and C2 (each responsible for 15.0%), and B3 was the fourth most common serotype (7.5%). This study demonstrated the practicability of using amplicon based NGS to identify HAdV diversity and quantify HAdV genome levels in environmental water samples, as well as broadening our current understanding of circulating HAdV in the Australian population.

Short-term succession of marine microbial fouling communities and the identification of primary and secondary colonizers

Microbial succession during the initial stages of marine biofouling has been rarely studied, especially in the Arabian Gulf. This study was undertaken to follow temporal shifts in biofouling communities in order to identify primary and secondary colonizers. Quantitative analysis revealed a significant increase in total biomass, coverage of macrofoulers, chlorophyll a concentrations, and bacterial counts with time. The relative abundance of the adnate diatoms increased with time, whereas it decreased in the case of the plocon diatoms. Non-metric multidimensional scaling (NMDS) ordination based on MiSeq data placed the bacterial communities in three distinct clusters, depending on the time of sampling. While the relative abundance of Alphaproteobacteria and Flavobacteriia decreased with time, suggesting their role as primary colonizers, the relative abundance of Actinobacteria and Planctomycetia increased with time, suggesting their role as secondary colonizers. Biofouling is a dynamic process that involves temporal quantitative and qualitative shifts in the micro- and macrofouling communities.

Presence of Newcastle disease viruses of sub-genotypes Vc and VIn in backyard chickens and in apparently healthy wild birds from Mexico in 2017

Virulent Newcastle disease viruses (NDV) have been present in Mexico since 1946, and recently, multiple outbreaks have been reported in the country. Here, we characterized eleven NDV isolated from apparently healthy wild birds and backyard chickens in three different locations of Jalisco, Mexico in 2017. Total RNA from NDV was reverse-transcribed, and 1285 nucleotides, which includes 3/4 of the fusion gene, was amplified and sequenced using a long-read MinION sequencing method. The sequences were 99.99-100% identical to the corresponding region obtained using the Illumina MiSeq. Phylogenetic analysis using MinION sequences demonstrated that nine virulent NDV from wild birds belonged to sub-genotypes Vc and VIn, and two backyard chicken isolates were of sub-genotype Vc. The sub-genotype Vc viruses had nucleotide sequence identity that ranged from 97.7 to 98% to a virus of the same sub-genotype isolated from a chicken in Mexico in 2010. Three viruses from pigeons had 96.3-98.7% nucleotide identity to sub-genotype VIn pigeon viruses, commonly referred to as pigeon paramyxovirus, isolated in the USA during 2000-2016. This study demonstrates that viruses of sub-genotype Vc are still present in Mexico, and the detection of this sub-genotype in both chickens and wild birds suggests that transmission among these species may represent a biosecurity risk. This is the first detection and complete genome sequencing of genotype VI NDV from Mexico. In addition, the utilization of an optimized long-read sequencing method for rapid virulence and genotype identification using the Oxford nanopore MinION system is demonstrated.

Differences in Soil Bacterial Community Compositions in Paddy Fields under Organic and Conventional Farming Conditions

Soil bacterial community compositions and temporal changes in organic paddy fields were elucidated using a 16S rRNA gene amplicon sequencing analysis with a high-throughput next generation sequencer. At transplanting, bacterial community compositions in organic and conventional paddy fields were mostly similar despite differences in field management. The bacterial community composition in organic fields differed from that under conventional management during the rice growth period, possibly as a result of the decomposition process of organic fertilizers. However, differences in the frequency of tillage and photosynthetic bacterial inoculations in the organic plots had less of an impact on bacterial communities.

In-depth comparative analysis of Illumina® MiSeq run metrics: Development of a wet-lab quality assessment tool

Whole genome sequencing of bacterial isolates has become a daily task in many laboratories, generating incredible amounts of data. However, data acquisition is not an end in itself; the goal is to acquire high-quality data useful for understanding genetic relationships. Having a method that could rapidly determine which of the many available run metrics are the most important indicators of overall run quality and having a way to monitor these during a given sequencing run would be extremely helpful to this effect. Therefore, we compared various run metrics across 486 MiSeq runs, from five different machines. By performing a statistical analysis using principal components analysis and a K-means clustering algorithm of the metrics, we were able to validate metric comparisons among instruments, allowing for the development of a predictive algorithm, which permits one to observe whether a given MiSeq run has performed adequately. This algorithm is available in an Excel spreadsheet: that is, MiSeq Instrument & Run (In-Run) Forecast. Our tool can help verify that the quantity/quality of the generated sequencing data consistently meets or exceeds recommended manufacturer expectations. Patterns of deviation from those expectations can be used to assess potential run problems and plan preventative maintenance, which can save valuable time and funding resources.

Microsatellite development from genome skimming and transcriptome sequencing: comparison of strategies and lessons from frog species

Background

Even though microsatellite loci frequently have been isolated using recently developed next-generation sequencing (NGS) techniques, this task is still difficult because of the subsequent polymorphism screening requires a substantial amount of time. Selecting appropriate polymorphic microsatellites is a critical issue for ecological and evolutionary studies. However, the extent to which assembly strategy, read length, sequencing depth, and library layout produce a measurable effect on microsatellite marker development remains unclear. Here, we use six frog species for genome skimming and two frog species for transcriptome sequencing to develop microsatellite markers, and investigate the effect of different isolation strategies on the yield of microsatellites.

Results

The results revealed that the number of isolated microsatellites increases with increased data quantity and read length. Assembly strategy could influence the yield and the polymorphism of microsatellite development. Larger k-mer sizes produced fewer total number of microsatellite loci, but these loci had a longer repeat length, suggesting greater polymorphism. However, the proportion of each type of nucleotide repeats was not affected; dinucleotide repeats were always the dominant type. Finally, the transcriptomic microsatellites displayed lower levels of polymorphisms and were less abundant than genomic microsatellites, but more likely to be functionally linked loci.

Conclusions

These observations provide deep insight into the evolution and distribution of microsatellites and how different isolation strategies affect microsatellite development using NGS.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5329-y) contains supplementary material, which is available to authorized users.

Bacterial community response to a preindustrial-to-future CO2 gradient is limited and soil specific in Texas Prairie grassland

Rising atmospheric CO₂ concentration directly stimulates plant productivity and affects nutrient dynamics in the soil. However, the influence of CO₂ enrichment on soil bacterial communities remains elusive, likely due to their complex interactions with a wide range of plant and soil properties. Here, we investigated the bacterial community response to a decade long preindustrial-to-future CO₂ gradient (250-500 ppm) among three contrasting soil types using 16S rRNA gene amplicon sequencing. In addition, we examined the effect of seasonal variation and plant species composition on bacterial communities. We found that Shannon index (H') and Faith's phylogenetic diversity (PD) did not change in response to the CO₂ gradient (R²  = 0.01, p > 0.05). CO₂ gradient and season also had a negligible effect on overall community structure, although silty clay soil communities were better structured on a CO₂ gradient (p < 0.001) among three soils. Similarly, CO₂ gradient had no significant effect on the relative abundance of different phyla. However, we observed soil-specific variation of CO₂ effects in a few individual families. For example, the abundance of Pirellulaceae family decreased linearly with CO₂ gradient, but only in sandy loam soils. Conversely, the abundance of Micromonosporaceae and Gaillaceae families increased with CO₂ gradient in clay soils. Soil water content (SWC) and nutrient properties were the key environmental constraints shaping bacterial community structure, one manifestation of which was a decline in bacterial diversity with increasing SWC. Furthermore, the impact of plant species composition on community structure was secondary to the strong influence of soil properties. Taken together, our findings indicate that bacterial communities may be largely unresponsive to indirect effects of CO₂ enrichment through plants. Instead, bacterial communities are strongly regulated by edaphic conditions, presumably because soil differences create distinct environmental niches for bacteria.

Whole Genome Sequencing of Influenza A and B Viruses With the MinION Sequencer in the Clinical Setting: A Pilot Study

Introduction: Whole genome sequencing (WGS) of influenza viruses is important for preparing vaccines and coping with newly emerging viruses. However, WGS is difficult to perform using conventional next-generation sequencers in developing countries, where facilities are often inadequate. In this study, we developed a high-throughput WGS method for influenza viruses in clinical specimens with the MinION portable sequencer.

Methods: Whole genomes of influenza A and B viruses were amplified by multiplex RT-PCR from 13 clinical specimens collected in Tokyo, Japan. Barcode tags for multiplex MinION sequencing were added with each multiplex RT-PCR amplicon by nested PCR with custom barcoded primers. All barcoded amplicons were mixed and multiplex sequencing using the MinION sequencer with 1D² sequencing kit. In addition, multiplex RT-PCR amplicons generated from each clinical specimen were sequenced using the Illumina MiSeq platform to validate the performance of MinION sequencer. The accuracy, recall, and precision rates of MinION sequencing were calculated by comparing the results of variant calling in the Illumina MiSeq platform and MinION sequencer.

Results: Whole genomes of influenza A and B viruses were successfully amplified by multiplex RT-PCR from 13 clinical samples. We identified 6 samples as influenza type A virus H3N2 subtype and 7 as influenza B virus Yamagata lineage using the Illumina MiSeq platform. The overall accuracy, recall, and precision rates of the MinION sequencer were, respectively 99.95%, 89.41%, and 97.88% from 1D reads and 99.97%, 93.28%, and 99.86% from 1D² reads.

Conclusion: We developed a novel WGS method for influenza A and B viruses. It is necessary to improve read accuracy and analytical tools in order to better utilize the MinION sequencer for real-time monitoring of genetic rearrangements and for evaluation of newly emerging viruses.

A Metagenome-Based Investigation of Gene Relationships for Non-Substrate-Associated Microbial Phosphorus Cycling in the Water Column of Streams and Rivers

Phosphorus (P) is a nutrient of primary importance in all living systems, and it is especially important in streams and rivers which are sensitive to anthropogenic P inputs and eutrophication. Microbes are accepted as the primary mineralizers and solubilizers of P improving bioavailability for organisms at all trophic levels. Here, we use a genomics approach with metagenome sequencing of 24 temperate streams and rivers representing a total P (TP) gradient to identify relationships between functional genes, functional gene groupings, P, and organisms within the P biogeochemical cycle. Combining information from network analyses, functional groupings, and system P levels, we have constructed a System Relational Overview of Gene Groupings (SROGG) which is a cohesive system level representation of P cycle gene and nutrient relationships. Using SROGG analysis in concert with other statistical approaches, we found that the compositional makeup of P cycle genes is strongly correlated to environmental P whereas absolute abundance of P genes shows no significant correlation to environmental P. We also found orthophosphate (PO₄^3-) to be the dominant factor correlating with system P cycle gene composition with little evidence of a strong organic phosphorous correlation present even in more oligotrophic streams.

BisAMP: A web-based pipeline for targeted RNA cytosine-5 methylation analysis

RNA cytosine-5 methylation (m⁵C) has emerged as a key epitranscriptomic mark, which fulfills multiple roles in structural modulation, stress signaling and the regulation of protein translation. Bisulfite sequencing is currently the most accurate and reliable method to detect m⁵C marks at nucleotide resolution. Targeted bisulfite sequencing allows m⁵C detection at single base resolution, by combining the use of tailored primers with bisulfite treatment. A number of computational tools currently exist to analyse m⁵C marks in DNA bisulfite sequencing. However, these methods are not directly applicable to the analysis of RNA m⁵C marks, because DNA analysis focuses on CpG methylation, and because artifactual unconversion and misamplification in RNA can obscure actual methylation signals. We describe a pipeline designed specifically for RNA cytosine-5 methylation analysis in targeted bisulfite sequencing experiments. The pipeline is directly applicable to Illumina MiSeq (or equivalent) sequencing datasets using a web interface (https://bisamp.dkfz.de), and is defined by optimized mapping parameters and the application of tailored filters for the removal of artifacts. We provide examples for the application of this pipeline in the unambiguous detection of m⁵C marks in tRNAs from mouse embryonic stem cells and neuron-differentiated stem cells as well as in 28S rRNA from human fibroblasts. Finally, we also discuss the adaptability of BisAMP to the analysis of DNA methylation. Our pipeline provides an accurate, fast and user-friendly framework for the analysis of cytosine-5 methylation in amplicons from bisulfite-treated RNA.

Museum metabarcoding: A novel method revealing gut helminth communities of small mammals across space and time

Natural history collections spanning multiple decades provide fundamental historical baselines to measure and understand changing biodiversity. New technologies such as next generation DNA sequencing have considerably increased the potential of museum specimens to address significant questions regarding the impact of environmental changes on host and parasite/pathogen dynamics. We developed a new technique to identify intestinal helminth parasites and applied it to shrews (Eulipotyphla: Soricidae) because they are ubiquitous, occupy diverse habitats, and host a diverse and abundant parasite fauna. Notably, we included museum specimens preserved in various ways to explore the efficacy of using metabarcoding analyses that may enable identification of helminth symbiont communities from historical archives. We successfully sequenced the parasite communities (using 12S mtDNA, 16S mtDNA, 28S rDNA) of 23 whole gastrointestinal tracts. All gastrointestinal tracts were obtained from the Museum of Southwestern Biology, USA, and from recent field collections, varying both in time since fixation (ranging from 4 months to 16 years) and preservation method (70% or 95% ethanol stored at room temperature, or flash frozen in liquid nitrogen and stored at -80 °C). Our proof of concept demonstrates the feasibility of applying next generation DNA sequencing techniques to authoritatively identify the parasite/pathogen communities within whole gastrointestinal tracts from museum specimens of varying age and fixation, and the value of future preservation of host-associated whole gastrointestinal tracts in public research archives. This powerful approach facilitates future comparative examinations of the distributions and interactions among multiple associated groups of organisms through time and space.

Recombinant GII.P16/GII.4 Sydney 2012 Was the Dominant Norovirus Identified in Australia and New Zealand in 2017

For the past two decades, norovirus pandemic variants have emerged every 3–5 years, and dominate until they are replaced by alternate strains. However, this scenario changed in 2016 with the co-circulation of six prevalent viruses, three of which possessed the pandemic GII.4 Sydney 2012 capsid. An increased number of institutional gastroenteritis outbreaks were reported within the Oceania region in mid-2017. This study identified emerging noroviruses circulating in Australia and New Zealand in 2017 to assess the changing dynamics of the virus infection. RT-PCR-based methods, next generation sequencing, and phylogenetic analyses were used to genotype noroviruses from both clinical and wastewater samples. Antigenic changes were observed between the capsid of pandemic Sydney 2012 variant and the two new Sydney recombinant viruses. The combination of these antigenic changes and the acquisition of a new ORF1 through recombination could both facilitate their ongoing persistence in the population. Overall, an increased prevalence of GII.P16/GII.4 Sydney 2012 viruses was observed in 2017, replacing the GII.P16/GII.2 recombinant that dominated in the region at the end of 2016. This shift in strain dominance was also observed in wastewater samples, demonstrating the reliability of wastewater as a molecular surveillance tool.

Use of next-generation sequencing in microbial risk assessment

Despite the ever increase in rigorous control and monitoring measures to assure safe food along the entire farm‐to‐fork chain, the past decade has also witnessed an increase in microbial food alerts. Hence, research on food safety and quality remain of utmost importance. Complementary, and at least as important, is the necessity to be able to assess the potential microbial risks along the food chain. Risk assessment relies on sound scientific data. Unfortunately, often, quality data are limited if not lacking. High‐throughput tools such as next‐generation sequencing (NGS) could fill this gap. NGS approaches can be used to generate ample qualitative and quantitative data to be used in the risk assessment process. NGS applications are not new in food microbiology with applications ranging from pathogen detection along the food chain, food epidemiology studies, whole genome analysis of food‐associated microorganisms up to describing complete food microbiomes. Yet, its application in the area of microbial risk assessment is still at an early stage and faces important challenges. The possibilities of NGS for risk assessment are ample, but so are the questions on the subject. One of the major strengths of NGS lies in its capacity to generate a lot of data, but to what extend can this wealth be of use in hazard identification, hazard characterisation and exposure assessment to perform a sound risk characterisation, which in turn will make it possible to take substantiated risk management decisions.

Accurate and Rapid Sequence Analysis of Adeno-Associated Virus Plasmids by Illumina Next-Generation Sequencing

Sequence validation of plasmid DNA is a crucial quality control step that must occur prior to adeno-associated virus (AAV) vector packaging through plasmid transfection. AAV cis-plasmids present unique challenges to sequence analysis, as they contain inverted terminal repeats and are prone to sequence rearrangements. An accurate and rapid next-generation sequencing approach has been established to analyze full-length sequences of AAV cis-plasmids within 3.5 days. Here, a step-by-step protocol is described that can reliably detect and identify the location and frequency of sequence variants commonly observed in AAV cis-plasmids.

The complete mitogenome of Caridina indistincta 'sp. A' (Decapoda: Atyidae)

A complete mitochondrial genome sequence was determined for a member of the Caridina indistincta species complex known as C. indistincta 'sp. A'. The 15,461 bp sequence (GenBank: MH189850) was obtained via genome skimming, and contains 13 protein coding genes, 22 tRNA genes, 2 rRNA genes, and a 646 bp control region arranged in the pancrustacean ground pattern. Caridina indistincta sp. A is a freshwater macroinvertebrate important for ecosystem health monitoring in Australia and this reference will be a useful resource for metabarcoding and eDNA studies.

Performance comparison of deep sequencing platforms for detecting HIV-1 variants in the pol gene

The present study compares the performances of an in-house sequencing protocol developed on MiSeq, the Sanger method, and the 454 GS-FLX for detecting and quantifying drug-resistant mutations (DRMs) in the human immunodeficiency virus polymerase gene (reverse transcriptase [RT] and protease [PR]). MiSeq sequencing identified all the resistance mutations detected by bulk sequencing (n = 84). Both the MiSeq and 454 GS-FLX platforms identified 67 DRMs in the RT and PR regions, but a further 25 DRMs were identified by only one or other of them. Pearson's analysis showed good concordance between the percentage of drug-resistant variants determined by MiSeq and 454 GS-FLX sequencing (ρ = .77, P < .0001). The MiSeq platform is as accurate as the 454 GS-FLX Roche system for determining RT and PR DRMs and could be used for monitoring human immunodeficiency virus type 1 drug resistance.

Fine grained compositional analysis of Port Everglades Inlet microbiome using high throughput DNA sequencing

BACKGROUND

Similar to natural rivers, manmade inlets connect inland runoff to the ocean. Port Everglades Inlet (PEI) is a busy cargo and cruise ship port in South Florida, which can act as a source of pollution to surrounding beaches and offshore coral reefs. Understanding the composition and fluctuations of bacterioplankton communities ("microbiomes") in major port inlets is important due to potential impacts on surrounding environments. We hypothesize seasonal microbial fluctuations, which were profiled by high throughput 16S rRNA amplicon sequencing and analysis.

METHODS & RESULTS

Surface water samples were collected every week for one year. A total of four samples per month, two from each sampling location, were used for statistical analysis creating a high sampling frequency and finer sampling scale than previous inlet microbiome studies. We observed significant differences in community alpha diversity between months and seasons. Analysis of composition of microbiomes (ANCOM) tests were run in QIIME 2 at genus level taxonomic classification to determine which genera were differentially abundant between seasons and months. Beta diversity results yielded significant differences in PEI community composition in regard to month, season, water temperature, and salinity. Analysis of potentially pathogenic genera showed presence of Staphylococcus and Streptococcus. However, statistical analysis indicated that these organisms were not present in significantly high abundances throughout the year or between seasons.

DISCUSSION

Significant differences in alpha diversity were observed when comparing microbial communities with respect to time. This observation stems from the high community evenness and low community richness in August. This indicates that only a few organisms dominated the community during this month. August had lower than average rainfall levels for a wet season, which may have contributed to less runoff, and fewer bacterial groups introduced into the port surface waters. Bacterioplankton beta diversity differed significantly by month, season, water temperature, and salinity. The 2013-2014 dry season (October-April), was warmer and wetter than historical averages. This may have driven significant differences in beta diversity. Increased nitrogen and phosphorous concentrations were observed in these dry season months, possibly creating favorable bacterial growth conditions. Potentially pathogenic genera were present in the PEI. However their relatively low, non-significant abundance levels highlight their relatively low risk for public health concerns. This study represents the first to sample a large port at this sampling scale and sequencing depth. These data can help establish the inlet microbial community baseline and supplement the vital monitoring of local marine and recreational environments, all the more poignant in context of local reef disease outbreaks and worldwide coral reef collapse in wake of a harsh 2014-16 El Niño event.

Robust Sub-nanomolar Library Preparation for High Throughput Next Generation Sequencing

BACKGROUND

Current library preparation protocols for Illumina HiSeq and MiSeq DNA sequencers require ≥2 nM initial library for subsequent loading of denatured cDNA onto flow cells. Such amounts are not always attainable from samples having a relatively low DNA or RNA input; or those for which a limited number of PCR amplification cycles is preferred (less PCR bias and/or more even coverage). A well-tested sub-nanomolar library preparation protocol for Illumina sequencers has however not been reported. The aim of this study is to provide a much needed working protocol for sub-nanomolar libraries to achieve outcomes as informative as those obtained with the higher library input (≥ 2 nM) recommended by Illumina's protocols.

RESULTS

Extensive studies were conducted to validate a robust sub-nanomolar (initial library of 100 pM) protocol using PhiX DNA (as a control), genomic DNA (Bordetella bronchiseptica and microbial mock community B for 16S rRNA gene sequencing), messenger RNA, microRNA, and other small noncoding RNA samples. The utility of our protocol was further explored for PhiX library concentrations as low as 25 pM, which generated only slightly fewer than 50% of the reads achieved under the standard Illumina protocol starting with > 2 nM.

CONCLUSIONS

A sub-nanomolar library preparation protocol (100 pM) could generate next generation sequencing (NGS) results as robust as the standard Illumina protocol. Following the sub-nanomolar protocol, libraries with initial concentrations as low as 25 pM could also be sequenced to yield satisfactory and reproducible sequencing results.

Effects of fermented milk treatment on microbial population and metabolomic outcomes in a three-stage semi-continuous culture system

We investigated the impact of a fermented milk product on gut microbiota and their metabolism in 3 different conditions of the colon with a systemic viewpoint. An in vitro semi-continuous anaerobic cultivation was used to assess the colon compartment-specific influence of fermented milk, followed by a multiomics approach combining 16S rDNA amplicon sequencing and nuclear magnetic resonance (NMR) spectroscopy. The microbiome profiling and metabolomic features were significantly different across three colon compartments and after fermented milk treatment. Integrative correlation analysis indicated that the alteration of butyrate-producing microbiota (Veillonella, Roseburia, Lachnospira, and Coprococcus) and some primary metabolites (butyrate, ethanol, lactate, and isobutyrate) in the treatment group had a strong association with the fermented milk microorganisms. Our findings suggested that fermented milk treatment significantly affected microbial population in an in vitro cultivation system as well as the colonic metabolome in different ways in each of colon compartment.

A comparison of computationally predicted functional metagenomes and microarray analysis for microbial P cycle genes in a unique basalt-soil forest

Here we compared microbial results for the same Phosphorus (P) biogeochemical cycle genes from a GeoChip microarray and PICRUSt functional predictions from 16S rRNA data for 20 samples in the four spatially separated Gotjawal forests on Jeju Island in South Korea. The high homogeneity of microbial communities detected at each site allows sites to act as environmental replicates for comparing the two different functional analysis methods. We found that while both methods capture the homogeneity of the system, both differed greatly in the total abundance of genes detected, as well as the diversity of taxa detected. Additionally, we introduce a more comprehensive functional assay that again captures the homogeneity of the system but also captures more extensive community gene and taxonomic information and depth. While both methods have their advantages and limitations, PICRUSt appears better suited to asking questions specifically related to microbial community P as we did here. This comparison of methods makes important distinctions between both the results and the capabilities of each method and can help select the best tool for answering different scientific questions.

Development and characterization of 43 microsatellite markers for the critically endangered primrose Primula reinii using MiSeq sequencing

Primula reinii (Primulaceae), a perennial herb belonging to the Primula section Reinii, occurs on wet, shaded rocky cliffs in the mountains of Japan. This threatened species comprises four varieties; these plants are very localized and rare in the wild. In this study, 43 microsatellite markers were developed using MiSeq sequencing to facilitate conservation genetics of these critically endangered primroses. We developed novel microsatellite markers for three varieties of P. reinii, and tested its polymorphism and genetic diversity using natural populations. These novel markers displayed relatively high polymorphism; the number of alleles and expected heterozygosities ranged from 2 to 6 (mean = 3.2) and 0.13 to 0.82 (mean = 0.45), respectively. All loci were in Hardy-Weinberg equilibrium. These microsatellite markers will be powerful tools to assess P. reinii genetic diversity and develop effective conservation and management strategies.

Comparative Analysis of Gut Microbiota of Culex restuans (Diptera: Culicidae) Females From Different Parents

The potential for gut microbiota to impede or enhance pathogen transmission is well-documented but the factors that shape this microbiota in mosquito vectors are poorly understood. We characterized and compared the gut microbiota of adult females of Culex restuans (Theobald; Diptera: Culicidae) from different parents. Cx. restuans larvae from nine field-collected egg rafts were reared on a common diet and gut microbiota of newly emerged adult females characterized by MiSeq sequencing of the V4 hypervariable region of the 16S rRNA gene. Bacterial diversity and evenness in individuals from one egg raft were significantly lower compared to those of three of the other eight egg rafts. The gut microbiota of adult females reared from seven of the nine egg rafts clustered together suggesting that individuals from most egg rafts had similar profiles of gut microbiota. These findings suggest that the microbiota of adult females from the same parents do not differ appreciably from the microbiota of adult females from different parents. However, additional studies using mosquitoes separated by geographic distances greater than those studied here and estimating the genetic distances between populations from different egg rafts are needed to provide further insights into the influence of host genetics on gut microbiota. Also worthwhile are studies evaluating how individuals from different egg rafts and harboring different gut microbiota compare in relation to vector competence for different pathogens.