A Comparison of Forensic Age Prediction Models Using Data From Four DNA Methylation Technologies

Individual age estimation can be applied to criminal, legal, and anthropological investigations. DNA methylation has been established as the biomarker of choice for age prediction, since it was observed that specific CpG positions in the genome show systematic changes during an individual’s lifetime, with progressive increases or decreases in methylation levels. Subsequently, several forensic age prediction models have been reported, providing average age prediction error ranges of ±3–4 years, using a broad spectrum of technologies and underlying statistical analyses. DNA methylation assessment is not categorical but quantitative. Therefore, the detection platform used plays a pivotal role, since quantitative and semi-quantitative technologies could potentially result in differences in detected DNA methylation levels. In the present study, we analyzed as a shared sample pool, 84 blood-based DNA controls ranging from 18 to 99 years old using four different technologies: EpiTYPER^®, pyrosequencing, MiSeq, and SNaPshot^TM. The DNA methylation levels detected for CpG sites from ELOVL2, FHL2, and MIR29B2 with each system were compared. A restricted three CpG-site age prediction model was rebuilt for each system, as well as for a combination of technologies, based on previous training datasets, and age predictions were calculated accordingly for all the samples detected with the previous technologies. While the DNA methylation patterns and subsequent age predictions from EpiTYPER^®, pyrosequencing, and MiSeq systems are largely comparable for the CpG sites studied, SNaPshot^TM gives bigger differences reflected in higher predictive errors. However, these differences can be reduced by applying a z-score data transformation.

Development of an NGS-Based Workflow for Improved Monitoring of Circulating Plasmids in Support of Risk Assessment of Antimicrobial Resistance Gene Dissemination

Antimicrobial resistance (AMR) is one of the most prominent public health threats. AMR genes localized on plasmids can be easily transferred between bacterial isolates by horizontal gene transfer, thereby contributing to the spread of AMR. Next-generation sequencing (NGS) technologies are ideal for the detection of AMR genes; however, reliable reconstruction of plasmids is still a challenge due to large repetitive regions. This study proposes a workflow to reconstruct plasmids with NGS data in view of AMR gene localization, i.e., chromosomal or on a plasmid. Whole-genome and plasmid DNA extraction methods were compared, as were assemblies consisting of short reads (Illumina MiSeq), long reads (Oxford Nanopore Technologies) and a combination of both (hybrid). Furthermore, the added value of conjugation of a plasmid to a known host was evaluated. As a case study, an isolate harboring a large, low-copy mcr-1-carrying plasmid (>200 kb) was used. Hybrid assemblies of NGS data obtained from whole-genome DNA extractions of the original isolates resulted in the most complete reconstruction of plasmids. The optimal workflow was successfully applied to multidrug-resistant Salmonella Kentucky isolates, where the transfer of an ESBL-gene-containing fragment from a plasmid to the chromosome was detected. This study highlights a strategy including wet and dry lab parameters that allows accurate plasmid reconstruction, which will contribute to an improved monitoring of circulating plasmids and the assessment of their risk of transfer.

Characterization of two MHC II genes (DOB, DRB) in white-tailed deer (Odocoileus virginianus)

BACKGROUND

The major histocompatibility complex (MHC) is responsible for detecting and addressing foreign pathogens inside the body. While the general structure of MHC genes is relatively well conserved among mammalian species, it is notably different among ruminants due to a chromosomal inversion that splits MHC type II genes into two subregions (IIa, IIb). Recombination rates are reportedly high between these subregions, and a lack of linkage has been documented in domestic ruminants. However, no study has yet examined the degree of linkage between these subregions in a wild ruminant. The white-tailed deer (Odocoileus virginianus), a popular ruminant of the Cervidae family, is habitually plagued by pathogens in its natural environment (e.g. Haemonchus contortus, Elaeophora). Due to the association between MHC haplotypes and disease susceptibility, a deeper understanding of MHC polymorphism and linkage between MHC genes can further aid in this species' successful management. We sequenced MHC-DRB exon 2 (IIa) and MHC-DOB exon 2 (IIb) on the MiSeq platform from an enclosed white-tailed deer population located in Alabama.

RESULTS

We identified 12 new MHC-DRB alleles, and resampled 7 alleles, which along with other published alleles brings the total number of documented alleles in white-tailed deer to 30 for MHC-DRB exon 2. The first examination of MHC-DOB in white-tailed deer found significantly less polymorphism (11 alleles), as was expected of a non-classical MHC gene. While MHC-DRB was found to be under positive, diversifying selection, MHC-DOB was found to be under purifying selection for white-tailed deer. We found no significant linkage disequilibrium between MHC-DRB and MHC-DOB, suggesting that these loci are unlikely to be closely linked.

CONCLUSIONS

Overall, this study identified 12 new MHC-DRB exon 2 alleles and characterized a new, non-classical, MHC II gene (MHC-DOB) for white-tailed deer. We also found a lack of significant linkage between these two loci, which supports previous findings of a chromosomal inversion within the MHC type II gene region in ruminants, and suggests that white-tailed deer may have a recombination hotspot between these MHC regions similar to that found for Bos taurus.

Preliminary Assessment of Microbial Community Structure of Wind-Tidal Flats in the Laguna Madre, Texas, USA

Aside from two samples collected nearly 50 years ago, little is known about the microbial composition of wind tidal flats in the hypersaline Laguna Madre, Texas. These mats account for ~42% of the lagoon's area. These microbial communities were sampled at four locations that historically had mats in the Laguna Madre, including Laguna Madre Field Station (LMFS), Nighthawk Bay (NH), and two locations in Kenedy Ranch (KRN and KRS). Amplicon sequencing of 16S genes determined the presence of 51 prokaryotic phyla dominated by Bacteroidota, Chloroflexi, Cyanobacteria, Desulfobacteria, Firmicutes, Halobacteria, and Proteobacteria. The microbial community structure of NH and KR is significantly different to LMFS, in which Bacteroidota and Proteobacteria were most abundant. Twenty-three cyanobacterial taxa were identified via genomic analysis, whereas 45 cyanobacterial taxa were identified using morphological analysis, containing large filamentous forms on the surface, and smaller, motile filamentous and coccoid forms in subsurface mat layers. Sample sites were dominated by species in Oscillatoriaceae (i.e., Lyngbya) and Coleofasciculaceae (i.e., Coleofasciculus). Most cyanobacterial sequences (~35%) could not be assigned to any established taxa at the family/genus level, given the limited knowledge of hypersaline cyanobacteria. A total of 73 cyanobacterial bioactive metabolites were identified using ultra performance liquid chromatography-Orbitrap MS analysis from these commu nities. Laguna Madre seems unique compared to other sabkhas in terms of its microbiology.

Assessment of biomass potentials of microalgal communities in open pond raceways using mass cultivation

Metagenome studies have provided us with insights into the complex interactions of microorganisms with their environments and hosts. Few studies have focused on microalgae-associated metagenomes, and no study has addressed aquatic microalgae and their bacterial communities in open pond raceways (OPRs). This study explored the possibility of using microalgal biomasses from OPRs for biodiesel and biofertilizer production. The fatty acid profiles of the biomasses and the physical and chemical properties of derived fuels were evaluated. In addition, the phenotype-based environmental adaptation ability of soybean plants was assessed. The growth rate, biomass, and lipid productivity of microalgae were also examined during mass cultivation from April to November 2017. Metagenomics analysis using MiSeq identified ∼127 eukaryotic phylotypes following mass cultivation with (OPR 1) or without (OPR 3) a semitransparent film. Of these, ∼80 phylotypes were found in both OPRs, while 23 and 24 phylotypes were identified in OPRs 1 and 3, respectively. The phylotypes belonged to various genera, such as Desmodesmus, Pseudopediastrum, Tetradesmus, and Chlorella, of which, the dominant microalgal species was Desmodesmus sp. On average, OPRs 1 and 3 produced ∼8.6 and 9.9 g m⁻² d⁻¹ (0.307 and 0.309 DW L⁻¹) of total biomass, respectively, of which 14.0 and 13.3 wt% respectively, was lipid content. Fatty acid profiling revealed that total saturated fatty acids (mainly C16:0) of biodiesel obtained from the microalgal biomasses in OPRs 1 and 3 were 34.93% and 32.85%, respectively; total monounsaturated fatty acids (C16:1 and C18:1) were 32.40% and 31.64%, respectively; and polyunsaturated fatty acids (including C18:3) were 32.68% and 35.50%, respectively. Fuel properties determined by empirical equations were within the limits of biodiesel standards ASTM D6751 and EN 14214. Culture solutions with or without microalgal biomasses enhanced the environmental adaptation ability of soybean plants, increasing their seed production. Therefore, microalgal biomass produced through mass cultivation is excellent feedstock for producing high-quality biodiesel and biofertilizer.

A Sporolactobacillus-, Clostridium-, and Paenibacillus- Dominant Microbial Consortium Improved Anaerobic RDX Detoxification by Starch Addition

In the present study, an anaerobic microbial consortium for the degradation of hexahydro-1,3,5- trinitro-1,3,5-triazine (RDX) was selectively enriched with the co-addition of RDX and starch under nitrogen-deficient conditions. Microbial growth and anaerobic RDX biodegradation were effectively enhanced by the co-addition of RDX and starch, which resulted in increased RDX biotransformation to nitroso derivatives at a greater specific degradation rate than those for previously reported anaerobic RDX-degrading bacteria (isolates). The accumulation of the most toxic RDX degradation intermediate (MNX [hexahydro-1-nitroso-3,5-dinitro-1,3,5-triazine]) was significantly reduced by starch addition, suggesting improved RDX detoxification by the co-addition of RDX and starch. The subsequent MiSeq sequencing that targeted the bacterial 16S rRNA gene revealed that the Sporolactobacillus, Clostridium, and Paenibacillus populations were involved in the enhanced anaerobic RDX degradation. These results suggest that these three bacterial populations are important for anaerobic RDX degradation and detoxification. The findings from this work imply that the Sporolactobacillus, Clostridium, and Paenibacillus dominant microbial consortium may be valuable for the development of bioremediation resources for RDX-contaminated environments.

Comparative Molecular Characterization of Three Gallid alphaherpesvirus Type 3 Strains 301B/1, HPRS24, and SB-1

Marek's disease (MD) is a highly contagious lymphoproliferative disease of chickens caused by Gallid alphaherpesvirus type 2. Gallid alphaherpesvirus type 3 (GaHV-3) strain 301B/1 was previously shown to be an effective MD vaccine with synergistic efficacy when used as a bivalent vaccine with turkey herpesvirus. Since the nucleotide sequences of only two GaHV-3 strains have been determined, we sought to sequence the 301B/1 genome using Illumina MiSeq technology. Phylogenomic analysis indicated that 301B/1 is more closely related to other GaHV-3 strains (SB-1 and HPRS24) than to virulent or attenuated strains of GaHV-2. One hundred and twenty-six open reading frames (ORFs) have been identified within the 301B/1 genome with 108 ORFs showing a high degree of similarity to homologs found in the genomes of SB-1 and HPRS24; 14 ORFs are highly homologous (> 90% identity) with the corresponding ORFs within the SB-1 genome. The R-LORF8 and R-LORF9 genes are the most dissimilar to the collinear genes found in the SB-1 genome but are highly homologous (99%-100% identity) with those within the HPRS24 genome. Overall the 301B/1 genome is most similar to the SB-1 virus genome (99.1%) and to a lesser degree with the HPRS24 virus genome (97.7%). However, six 301B/1 ORFs (UL47, UL48, UL52, pp38, ICP4, and US10) have been identified that contain nonsynonymous substitutions relative to homologs found in the SB-1 genome. Notably, unlike the avian retrovirus long terminal repeat sequences found within the SB-1 genome, none were identified within the 301B/1 genome.

Diversity of navel microbiome in young adults

Introduction. Human skin microbial communities represent a tremendous source of genetic diversity that evolves as a function of human age. Microbiota differs between regions of oily and moist skin, and appears to stabilize with age.Aim. We have a minimal understanding of the time frame required for the stabilization of skin microbiota, and the role played by gender. In the current study, we examined the microbiota present in the navel region of college-attending young adults in the age group of 18-25 years and investigated if diversity is associated with gender (male and female).Method. The study involved 16 female and six male subjects. Isolated DNA samples from navel swabs were processed using the Nextera XT library preparation kit and sequenced using the MiSeq platform. Data were analysed using QIIME and statistical analysis performed in R.Results. Microbiota of navel skin is dominated by Corynebacterium and Staphylococcus and includes opportunistic pathogens like Clostridium and Pseudomonas. Also present as the major component of the flora were the organisms normally associated with the gastrointestinal tract such as Acinetobacter, Campylobacter, Klebsiella and organisms from the Enterobacteriaceae and Moraxellaceae families. Comparison of alpha and beta diversity of the microbiota in the male and female navel regions suggests that the flora is not statistically different (P>0.05). However, pairwise comparison suggests that the abundance of 12 specific genera varied with gender, including higher abundance of Klebsiella and Enterobacter in females.Conclusion. Our findings indicate that the navel skin microbiota of young adults has a core microbiota of Corynebacterium and Staphylococcus. We also noted the presence of a significant number of opportunistic pathogens. A minor gender difference in the abundance of individual organisms was also observed.

Supplementation with Chlorella vulgaris, Chlorella protothecoides, and Schizochytrium sp. increases propionate-producing bacteria in in vitro human gut fermentation

BACKGROUND

Gut microbiota are major contributors to host metabolism and are considered as potential targets of novel therapeutics. Microalgae have a strong potential for use as prebiotics because they are a rich source of proteins, fatty acids, fiber, and minerals for nutritional supplementation in humans. Nevertheless, there has been insufficient research into the effect of microalgae on gut microbiota. To investigate the effects of three edible microalgae (Chlorella vulgaris, Chlorella protothecoides, and Schizochytrium sp.) on gut microbiota, simulated digestion and colonic fermentation were examined.

RESULTS

Following in vitro digestion, the microalgae displayed different levels of bioaccessibility and the nutrient analysis revealed that unabsorbed nutrients during the digestion process could be used for colonic fermentation. Following colonic fermentation, the control, inulin, and microalgae groups displayed different metabolite tendencies when investigated with nuclear magnetic resonance (NMR) spectroscopic analysis. In particular, microalgae supplementation increased the proportion of propionate in the colonic culture (control: 19.14%, Inulin: 18.38%, C. vulgaris: 25.80%, C. protothecoides: 25.46%, and Schizochytrium sp.: 25.56%). Microbial profiling analysis using 16S rRNA gene sequencing also disclosed that the relative abundance of Bacteroides (control: 1.91%, inulin: 2.61%, C. vulgaris: 14.77%, C. protothecoides: 11.17%, and Schizochytrium sp.: 5.51%) and Dialister (control: 0.08%, inulin: 2.06%, C. vulgaris: 6.79%, C. protothecoides: 4.45%, and Schizochytrium sp.: 4.48%), involved in propionate metabolism increased more than in the inulin group.

CONCLUSION

Our findings suggest the potential use of microalgae as a functional food to increase propionate generation because propionate has been reported to be effective in weight loss and the inhibition of pathogen infection. © 2020 Society of Chemical Industry.

First detection of a plasmid located carbapenem resistant blaVIM-1 gene in E. coli isolated from meat products at retail in Belgium in 2015

Carbapenemase-producing Enterobacteriaceae (CPE) confer resistance to antibiotics that are of critical importance to human medicine. There have only been a few reported cases of CPEs in the European food chain. We report the first detection of a carbapenemase-producing Escherichia coli (ST 5869) in the Belgian food chain. Our aim was to characterize the origin of the carbapenem resistance in the E. coli isolate. The isolate was detected during the screening of 178 minced pork samples and was shown to contain the carbapenemase gene bla_VIM-1 by PCR and Sanger sequencing. Whole genome short and long read sequencing (MiSeq and MinION) was performed to characterize the isolate. With a hybrid assembly we reconstructed a 190,205 bp IncA/C2 plasmid containing bla_VIM-1 (S15FP06257_p), in addition to other critically important resistance genes. This plasmid showed only low similarity to plasmids containing bla_VIM-1 previously reported in Germany. Moreover, no sequences existed in the NCBI nucleotide database that completely covered S15FP06257_p. Analysis of the bla_VIM-1 gene cassette demonstrated that it likely originated from an integron of a Klebsiella plasmid reported previously in a clinical isolate in Europe, suggesting that the meat could have been contaminated by human handling in one of the steps of the food chain. This study shows the relevance of fully reconstructing plasmids to characterize their genetic content and to allow source attribution. This is especially important in view of the potential risk of antimicrobial resistance gene transmission through mobile elements as was reported here for the of public health concern bla_VIM-1.

Microbiota Analysis Using Sequencing by Synthesis: From Library Preparation to Sequencing

The highly parallel nature of sequencing by synthesis (SBS) allows millions of amplicons to be sequenced simultaneously, which has led to enormous interest in the investigation of bacterial communities (often referred to as the microbiota). In this protocol, we describe a method for the 'universal' amplification of the v4 region of the bacterial 16S rRNA gene from genomic DNA and prepare these amplicons so that they can be sequenced using the MiSeq system (Illumina). The protocol provides instruction on sequencing of 188 genomic DNA samples plus PCR positive and negative controls, which can be applied to any sample type where bacterial DNA may be of interest.

A method to generate multilocus barcodes of pinned insect specimens using MiSeq

For molecular insect identification, amplicon sequencing methods are recommended because they offer a cost-effective approach for targeting small sets of informative genes from multiple samples. In this context, high-throughput multilocus amplicon sequencing has been achieved using the MiSeq Illumina sequencing platform. However, this approach generates short gene fragments of <500 bp, which then have to be overlapped using bioinformatics to achieve longer sequence lengths. This increases the risk of generating chimeric sequences or leads to the formation of incomplete loci. Here, we propose a modified nested amplicon sequencing method for targeting multiple loci from pinned insect specimens using the MiSeq Illumina platform. The modification exists in using a three-step nested PCR approach targeting near full-length loci in the initial PCR and subsequently amplifying short fragments of between 300 and 350 bp for high-throughput sequencing using Illumina chemistry. Using this method, we generated 407 sequences of three loci from 86% of all the specimens sequenced. Out of 103 pinned bee specimens of replicated species, 71% passed the 95% sequence similarity threshold between species replicates. This method worked best for pinned specimens aged between 0 and 5 years, with a limit of 10 years for pinned and 14 years for ethanol-preserved specimens. Hence, our method overcomes some of the challenges of amplicon sequencing using short read next generation sequencing and improves the possibility of creating high-quality multilocus barcodes from insect collections.

The Role of Microbial Mats in the Removal of Hexavalent Chromium and Associated Shifts in Their Bacterial Community Composition

Microbial mats are rarely reported for chromium-polluted ecosystems, hence information on their bacterial diversity and role in chromium removal are very scarce. We investigated the role of nine microbial mats, collected from three quarry sumps of chromium mining sites, in the removal of hexavalent chromium [Cr(VI)]. Bacterial diversity in these mats and community shifts after incubation with Cr(VI) have been investigated using MiSeq sequencing. In nature, a chromium content of 1,911 ± 100 mg kg^–1 was measured in the microbial mats, constituting the third highest source of environmentally available chromium. The mats were able to remove 1 mg l^–1 of Cr(VI) in 7 days under aerobic conditions. MiSeq sequencing of the original mats yielded 46–99% of the sequences affiliated to Proteobacteria, Firmicutes and Actinobacteria. When the mats were incubated with Cr(VI), the bacterial community shifted in the favor of Alphaproteobacteria and Verrucomicrobiae. We conclude that microbial mats in the quarry sumps harbor diverse microorganisms with the ability to remove toxic Cr(VI), hence these mats can be potentially used to remove chromium from polluted waters.

High-Resolution Mapping of 3' Extremities of RNA Exosome Substrates by 3' RACE-Seq

The main 3'-5' exoribonucleolytic activity of eukaryotic cells is provided by the RNA exosome. The exosome is constituted by a core complex of nine subunits (Exo9), which coordinates the recruitment and the activities of distinct types of cofactors. The RNA exosome cofactors confer distributive and processive 3'-5' exoribonucleolytic, endoribonucleolytic, and RNA helicase activities. In addition, several RNA binding proteins and terminal nucleotidyltransferases also participate in the recognition of exosome RNA substrates.To fully understand the biological roles of the exosome, the respective functions of its cofactors must be deciphered. This entails the high-resolution analysis of 3' extremities of degradation or processing intermediates in different mutant backgrounds or growth conditions. Here, we describe a detailed 3' RACE-seq procedure for targeted mapping of exosome substrate 3' ends. This procedure combines a 3' RACE protocol with Illumina sequencing to enable the high-resolution mapping of 3' extremities and the identification of untemplated nucleotides for selected RNA targets.

Degradability and biofouling of oxo-biodegradable polyethylene in the planktonic and benthic zones of the Arabian Gulf

Little is known about the degradability of oxo-biodegradable polyethylene (OXO-PE) and its surface fouling bacterial communities in the marine environment. The degradation of OXO-PE, PE and polyethylene terephthalate (PET) was compared at two depths (2 m and 6 m) in the Arabian Gulf. Scanning electron microcopy (SEM) revealed more fissure formation on OXO-PE and PE than on PET, indicating physical degradation. The formation of hydroxyl groups and carbonyl bonds, by Fourier-transform infrared spectroscopy (FTIR), suggests chemical degradation of OXO-PE. Plastisphere bacterial communities on OXO-PE and PE were different than on PET. Proteobacteria, Bacteriodetes and Planctomycetes were detected on all plastics, however, sequences of Alteromonas and Zoogloea were OXO-PE-specific suggesting a possible involvement of these bacterial genera in OXO-PE degradation. We conclude that OXO-PE shows increased signs of degradation with time owing to the combination of abiotic and biotic processes, and its degradation is higher in the benthic than in the planktonic zone.

Pretreatment of Lignocellulosic Biomass with Cattle Rumen Fluid for Methane Production: Fate of Added Rumen Microbes and Indigenous Microbes of Methane Seed Sludge

The pretreatment of lignocellulosic substrates with cattle rumen fluid was successfully developed to increase methane production. In the present study, a 16S rRNA gene-targeted amplicon sequencing approach using the MiSeq platform was applied to elucidate the effects of the rumen fluid treatment on the microbial community structure in laboratory-scale batch methane fermenters. Methane production in fermenters fed rumen fluid-treated rapeseed (2,077.3 mL CH4 reactor-1 for a 6-h treatment) was markedly higher than that in fermenters fed untreated rapeseed (1,325.8 mL CH4 reactor-1). Microbial community profiling showed that the relative abundance of known lignocellulose-degrading bacteria corresponded to lignocellulose-degrading enzymatic activities. Some dominant indigenous cellulolytic and hemicellulolytic bacteria in seed sludge (e.g., Cellulosilyticum lentocellum and Ruminococcus flavefaciens) and rumen fluid (e.g., Butyrivibrio fibrisolvens and Prevotella ruminicola) became undetectable or markedly decreased in abundance in the fermenters fed rumen fluid-treated rapeseed, whereas some bacteria derived from seed sludge (e.g., Ruminofilibacter xylanolyticum) and rumen fluid (e.g., R. albus) remained detectable until the completion of methane production. Thus, several lignocellulose-degrading bacteria associated with rumen fluid proliferated in the fermenters, and may play an important role in the degradation of lignocellulosic compounds in the fermenter.

Bioinformatic strategies to address limitations of 16rRNA short-read amplicons from different sequencing platforms

Sequencing the 16S gene rRNA has become a popular method when identifying bacterial communities. However, recent studies address differences in the characterization based on sample preparation, sequencing platforms, and data analysis. In this work, we tested some of the available user-friendly protocols for data analysis with the reads obtained from the sequencing machines Illumina MiSeq and Ion Torrent Personal Genome Machine (PGM). We sought for the advantages and disadvantages of both platforms in terms of accuracy, detected species, and abundance, analyzing a staggered mock community. Four different pipelines were applied: QIIME 1.9.1 with default parameters, QIIME 1.9.1 with modified parameters and chimera removal, VSEARCH 2.3.4, and QIIME 2 v.2018.2. To address the limitations of species level detection we used species-classifier SPINGO. The optimal pipeline for PGM platform, was the use of QIIME 1.9.1 with default parameters (QIIME1), except when a study requires the detection of Bacteroides or other Bacteroidaceae members, in which QIIME1MOD (with chimera removal) seems to be a good alternative. For Illumina Miseq, VSEARCH strategy can be a good option. Our results also confirm that all the tested pipelines can be used for metagenomic analysis at family and genus level.

Characterization of complete mitochondrial genome of Pogonophryne albipinna (Perciformes: Artedidraconidae)

The complete mitochondrial genome of Pogonophryne albipinna was determined by the MiSeq platform, which was the first report in the family Artedidraconidae. The circular form of its mitochondrial genome was 17,086 bp, which contained the canonical eukaryotic 37 genes. The gene orders of P. albipinna was identical to the other icefish species, in which there was additional non-coding region and translocation of ND6 gene. Except for ATP6 gene and COI (GTG), 11 genes begin with the typical start codon, while incomplete stop codons (T– –) were identified in COII, ND4, and CytB. Phylogenetic tree with the currently known mitogenomes in suborder Notothenioidei showed that P. albipinna was located distinctly from those in Bathydraconidae and Nototheniidae forming a unique cluster as Artedidraconidae. The first complete mitochondrial genome of P. albipinna would be the fundamental data to understand the evolutional relationship of icefish species in the Antarctic Oceans.

Building a custom large-scale panel of novel microhaplotypes for forensic identification using MiSeq and Ion S5 massively parallel sequencing systems

A large number of new microhaplotype loci were identified in the human genome by applying a directed search with selection criteria emphasizing short haplotype length (<120 nucleotides) and maximum levels of polymorphism in the composite SNPs. From these searches, 107 autosomal microhaplotypes and 11 X chromosome microhaplotypes were selected, with well-spaced autosomal positions to ensure their independence in relationship tests. The 118 microhaplotypes were assembled into a single multiplex assay for the analysis of forensic DNA with massively parallel sequencing (MPS). A single AmpliSeq-adapted primer set was made for Illumina MiSeq and Thermo Fisher Ion S5 MPS platforms and the performance of the assay was comprehensively evaluated in both systems. Five microhaplotypes showed critical sequencing failures in both MPS platforms and were removed, while a further 13 required manual checks and the application of sequence quality thresholds in one or both systems to ensure the successful analysis of low-level DNA in these loci. The targeting of short microhaplotype spans during marker selection, with an average length of 51 nucleotides in the 118 loci, led to a high level of sensitivity for the panel when sequencing the very degraded DNA typically encountered in forensic casework and the identification of missing persons.

Comparative 16S rDNA metagenomics study of two samples of cassava peel heap from Nigeria and India

The microbiology of many cassava products and the wastes generated during the processing have been reported; however, majority of these reports used culture-dependent methods. This has resulted in a dearth of information on the bacterial diversity of cassava peels and peel heaps. Large amounts of cassava peels generated during the processing of cassava root are usually discharged on land or water as wastes and are allowed to rot in the open, especially in some developing countries. Culture-independent methods such as PCR-based amplification and sequencing of 16S rRNA genes, among others have been used in recent times to study the diversity of microbes in different environmental samples. In this study, bacterial isolates were screened for cellulase and xylanase enzyme activities on minimal agar and genomic DNA was isolated from cassava peel samples; metagenomics was carried out using MiSeq 2 × 300 with primers specific for V3-V4 bacterial region. Samples collected from Nigeria (AAG) had more species compared with samples from India (JHA) with 793 and 525 observed OTUs (operational taxonomic units), respectively. Five bacterial isolates from cassava peel-heap samples obtained from Ogbomoso, Nigeria showed no ability to produce cellulase enzyme, seven isolates from the Nigeria samples and three from Jorhat samples were positive for xylanase production; the highest amylase activity was shown by isolate AG18 (10,055 U/mL), while the lowest was recorded for isolate JA2 (2333 U/mL) with a significant difference observed in the amylase activities of isolates (p ≤ 0.05). Comparing the most abundant taxonomy for each of the samples at different taxonomic levels, the most abundant for sample AAG were phylum Firmicutes (42.11%), class Bacilli (41.27%), order Lactobacillales (33.11%), family Acetobacteraceae (31.30%), genus Acetobacter (30.02%) and unclassified species of Acetobacter (29.88%), while sample JHA had Actinobacteria (47.47%) as the highest phylum and class, order Actinomycetales (47.47%), family Brevibacteriaceae (46.97%), genus Brevibacterium (46.97%) and unclassified species of Brevibacterium (46.89%). This study provides an insight into the vast diversity of the bacteria associated with cassava peel heaps and the ability of some of the bacteria to produce selected extracellular enzymes.

Adapterama II: universal amplicon sequencing on Illumina platforms (TaggiMatrix)

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I, we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr_Demuxy.

Parameters Governing the Community Structure and Element Turnover in Kermadec Volcanic Ash and Hydrothermal Fluids as Monitored by Inorganic Electron Donor Consumption, Autotrophic CO2 Fixation and 16S Tags of the Transcriptome in Incubation Experiments

The microbial community composition and its functionality was assessed for hydrothermal fluids and volcanic ash sediments from Haungaroa and hydrothermal fluids from the Brothers volcano in the Kermadec island arc (New Zealand). The Haungaroa volcanic ash sediments were dominated by epsilonproteobacterial Sulfurovum sp. Ratios of electron donor consumption to CO₂ fixation from respective sediment incubations indicated that sulfide oxidation appeared to fuel autotrophic CO₂ fixation, coinciding with thermodynamic estimates predicting sulfide oxidation as the major energy source in the environment. Transcript analyses with the sulfide-supplemented sediment slurries demonstrated that Sulfurovum prevailed in the experiments as well. Hence, our sediment incubations appeared to simulate environmental conditions well suggesting that sulfide oxidation catalyzed by Sulfurovum members drive biomass synthesis in the volcanic ash sediments. For the Haungaroa fluids no inorganic electron donor and responsible microorganisms could be identified that clearly stimulated autotrophic CO₂ fixation. In the Brothers hydrothermal fluids Sulfurimonas (49%) and Hydrogenovibrio/Thiomicrospira (15%) species prevailed. Respective fluid incubations exhibited highest autotrophic CO₂ fixation if supplemented with iron(II) or hydrogen. Likewise catabolic energy calculations predicted primarily iron(II) but also hydrogen oxidation as major energy sources in the natural fluids. According to transcript analyses with material from the incubation experiments Thiomicrospira/Hydrogenovibrio species dominated, outcompeting Sulfurimonas. Given that experimental conditions likely only simulated environmental conditions that cause Thiomicrospira/Hydrogenovibrio but not Sulfurimonas to thrive, it remains unclear which environmental parameters determine Sulfurimonas’ dominance in the Brothers natural hydrothermal fluids.

Validation of NGS for mitochondrial DNA casework at the FBI Laboratory

As a first step towards integrating next generation sequencing (NGS) technology into the FBI Laboratory's operational casework, the PowerSeq™ CRM Nested System, an NGS-based mitochondrial DNA (mtDNA) control region assay, was developmentally and internally validated. The validation studies were conducted in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) Validation Guidelines for Forensic DNA Analysis Methods, and the FBI's Quality Assurance Standards (QAS) for Forensic DNA Testing Laboratories. The assay was shown to be highly reproducible, with variant frequencies across intra and inter-run replicates of the same sample differing, on average, by just 0.3% for substitutions and point heteroplasmies and 1.5% for insertions and deletions. The assay was also shown to be extremely sensitive, yielding complete control region sequence data from as few as 2000 copies of mtDNA. This is a more than 20-fold increase in sensitivity when compared to the FBI Laboratory's current Sanger sequencing-based protocols and, based on mtDNA quantitation values of samples routinely encountered in mtDNA casework, suggests that the percentage of questioned samples from which full control region data can be recovered will increase from our current 20% to approximately 90% success with NGS technology. In addition, the assay requires on average only 30% of the extract volume typically required to develop control region profiles from degraded samples via Sanger sequencing. Overall, these studies establish the reliability of the PowerSeq™ CRM Nested System for accurate mtDNA control region typing and can serve as a model for laboratories seeking to validate NGS protocols for forensic mtDNA analysis.

High-throughput amplicon sequencing datasets of the metacommunity DNA of the gut microbiota of naturally occurring and laboratory aquaculture green sea urchins Lytechinus variegatus

We present high-throughput amplicon sequence (HTS) datasets of the microbial metacommunity DNA of the gut tissue and the gut digesta of naturally occurring (n = 3) and laboratory aquaculture (n = 2) green sea urchins, Lytechinus variegatus. The HTS datasets were generated on an Illumina MiSeq by targeting the amplicons of the V4 region of the 16S rRNA gene. After the raw sequences were quality checked and filtered, 88% of the sequence reads were subjected to bioinformatics analyses to generate operation taxonomic units (OTUs), which were then verified for saturation by using rarefaction analysis at a 3% sequence variation. Further, the OTUs were randomly subsampled to the minimum sequence count values. Then, the FASTA-formatted representative sequences of the microbiota were assigned taxonomic identities through multiple databases using the SILVA ACT: Alignment, Classification and Tree Service (www.arb-silva.de/aligner). The HTS datasets of this metagenome can be accessed from the BioSample Submission Portal (https://www.ncbi.nlm.nih.gov/bioproject/) under the BioProject IDs PRJNA291441 and PRJNA326427.

Data on taxonomic annotation and diversity of 18S rRNA gene amplicon libraries derived from high throughput sequencing

This Data in Brief article is a supporting information for the research article entitled “Protistan community composition in anoxic sediments from three salinity-disparate Japanese lakes” by Kataoka and Kondo (2019) [1]. Summary of 18S rRNA gene sequences originated from anoxic sediment of three lakes in two seasons using high throughput sequencing techniques (MiSeq, Illumina) was shown in this data article. Supergroup-level taxonomy was compared between the SILVA search for SILVA database and BLASTn search for the PR2 database. Alpha diversity was calculated in each sample, and beta-diversity was calculated among the six amplicon libraries. Partial sequence length between the primer set of 574*f and 1132R Hugerth et al., 2015 was compared between the forward read and the combined read.

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.