Comparison of the performance of MiSeq and NovaSeq in oral microbiome study

Objective

Next generation sequencing is commonly used to characterize the microbiome structure. MiSeq is most commonly used to analyze the microbiome due to its relatively long read length. Illumina also introduced the 250 × 2 chip for NovaSeq. The purpose of this study was to compare the performance of MiSeq and NovaSeq in the context of oral microbiome study.

Methods

Total read count, read quality score, relative bacterial abundance, community diversity, and correlation between two platforms were analyzed. Phylogenetic trees were analyzed for Streptococcus and periodontopathogens.

Results

NovaSeq produced significantly more read counts and assigned more operational taxonomic units (OTUs) compared to MiSeq. Community diversity was similar between MiSeq and NovaSeq. NovaSeq were able to detect more unique OTUs compared to MiSeq. When phylogenetic trees were constructed for Streptococcus and periodontopathogens, both platforms detected OTUs for most of the clades.

Conclusion

Taken together, while both MiSeq and NovaSeq platforms effectively characterize the oral microbiome, NovaSeq outperformed MiSeq in terms of read counts and detection of unique OTUs, highlighting its potential as a valuable tool for large scale oral microbiome studies.

Unraveling the shift in bacterial communities profile grown in sediments co-contaminated with chlorolignin waste of pulp-paper mill by metagenomics approach

Pulp-paper mills (PPMs) are known for consistently generating a wide variety of pollutants, that are often unidentified and highly resistant to environmental degradation. The current study aims to investigate the changes in the indigenous bacterial communities profile grown in the sediment co-contaminated with organic and inorganic pollutants discharged from the PPMs. The two sediment samples, designated PPS-1 and PPS-2, were collected from two different sites. Physico-chemical characterization of PPS-1 and PPS-2 revealed the presence of heavy metals (mg kg-1) like Cu (0.009-0.01), Ni (0.005-0.002), Mn (0.078-0.056), Cr (0.015-0.009), Pb (0.008-0.006), Zn (0.225-0.086), Fe (2.124-0.764), Al (3.477-22.277), and Ti (99.792-45.012) along with high content of chlorophenol, and lignin. The comparative analysis of organic pollutants in sediment samples using gas chromatography-mass spectrometry (GC-MS) revealed the presence of major highly refractory compounds, such as stigmasterol, β-sitosterol, hexadecanoic acid, octadecanoic acid; 2,4-di-tert-butylphenol; heptacosane; dimethyl phthalate; hexachlorobenzene; 1-decanol,2-hexyl; furane 2,5-dimethyl, etc in sediment samples which are reported as a potential toxic compounds. Simultaneously, high-throughput sequencing targeting the V3-V4 hypervariable region of the 16S rRNA genes, resulted in the identification of 1,249 and 1,345 operational taxonomic units (OTUs) derived from a total of 115,665 and 119,386 sequences read, in PPS-1 and PPS-2, respectively. Analysis of rarefaction curves indicated a diversity in OTU abundance between PPS-1 (1,249 OTUs) and PPS-2 (1,345 OTUs). Furthermore, taxonomic assignment of metagenomics sequence data showed that Proteobacteria (55.40%; 56.30%), Bacteoidetes (11.30%; 12.20%), and Planctomycetes (5.40%; 4.70%) were the most abundant phyla; Alphproteobacteria (20.50%; 23.50%), Betaproteobacteria (16.00%; 12.30%), and Gammaproteobacteria were the most recorded classes in PPS-1 and PPS-2, respectively. At the genus level, Thiobacillus (7.60%; 4.50%) was the most abundant genera grown in sediment samples. The results indicate significant differences in both the diversity and relative abundance of taxa in the bacterial communities associated with PPS-2 when compared to PPS-1. This study unveils key insights into contaminant characteristics and shifts in bacterial communities within contaminated environments. It highlights the potential for developing efficient bioremediation techniques to restore ecological balance in pulp-paper mill waste-polluted areas, stressing the importance of identifying a significant percentage of unclassified genera and species to explore novel genes.

Closed genome sequence of Clostridium botulinum type B1 strain isolated from an infant botulism case in the United States

We present the closed genome sequence of the Clostridium botulinum BT-22100019 strain isolated from the stool specimen of an infant diagnosed with botulism. With 4.33-Mb genome size and 28.0% G + C content, the bont/B1 gene encoded for botulinum neurotoxin serotype B was found on a 262 kb plasmid arranged in a ha+ orfx - cluster.

Bacteriomic Profiles of Rock-Dwelling Lichens from the Venezuelan Guiana Shield and the South African Highveld Plateau

Lichens are not only fungal-algal symbiotic associations but also matrices for association with bacteria, and the bacterial diversity linked to lichens has been receiving more attention in studies. This study compares the diversity and possible metabolism of lichen-associated bacteria from saxicolous foliose and fruticose taxa Alectoria, Canoparmelia, Crocodia, Menegazzia, Usnea, and Xanthoparmelia from the Venezuelan Guiana Shield and the South African Highveld Plateau. We used DNA extractions from the lichen thalli to amplify the eukaryotic 18S rRNA gene (rDNA) and the V3-V4 region of the bacterial 16S rDNA, of which amplicons were then Sanger- and MiSeq-sequenced, respectively. The V3-V4 sequences of the associated bacteria were grouped into operational taxonomic units (OTUs) ascribed to twelve bacterial phyla previously found in the rock tripe Umbilicaria lichens. The bacterial OTUs emphasized the uniqueness of each region, while, at the species and higher ranks, the regional microbiomes were shown to be somewhat similar. Nevertheless, regional biomarker OTUs were screened to predict relevant metabolic pathways, which implicated different regional metabolic features.

Genetic Diversity of Grapevine Virus A in Three Australian Vineyards Using Amplicon High Throughput Sequencing (Amplicon-HTS)

Shiraz disease (SD) is one of the most destructive viral diseases of grapevines in Australia and is known to cause significant economic loss to local growers. Grapevine virus A (GVA) was reported to be the key pathogen associated with this disease. This study aimed to better understand the diversity of GVA variants both within and between individual SD and grapevine leafroll disease (LRD) affected grapevines located at vineyards in South Australia. Amplicon high throughput sequencing (Amplicon-HTS) combined with median-joining networks (MJNs) was used to analyze the variability in specific gene regions of GVA variants. Several GVAII variant groups contain samples from both vineyards studied, suggesting that these GVAII variants were from a common origin. Variant groups analyzed by MJNs using the overall data set denote that there may be a possible relationship between variant groups of GVA and the geographical location of the grapevines.

Effects of error, chimera, bias, and GC content on the accuracy of amplicon sequencing

IMPORTANCE

Amplicon sequencing of targeted genes is the predominant approach to estimate the membership and structure of microbial communities. However, accurate reconstruction of community composition is difficult due to sequencing errors, and other methodological biases and effective approaches to overcome these challenges are essential. Using a mock community of 33 phylogenetically diverse strains, this study evaluated the effect of GC content on sequencing results and tested different approaches to improve overall sequencing accuracy while characterizing the pros and cons of popular amplicon sequence data processing approaches. The sequencing results from this study can serve as a benchmarking data set for future algorithmic improvements. Furthermore, the new insights on sequencing error, chimera formation, and GC bias from this study will help enhance the quality of amplicon sequencing studies and support the development of new data analysis approaches.

Conservation and diversity of the pollen microbiome of Pan-American maize using PacBio and MiSeq

Pollen is a vector for diversification, fitness-selection, and transmission of plant genetic material. The extent to which the pollen microbiome may contribute to host diversification is largely unknown, because pollen microbiome diversity within a plant species has not been reported, and studies have been limited to conventional short-read 16S rRNA gene sequencing (e.g., V4-MiSeq) which suffers from poor taxonomic resolution. Here we report the pollen microbiomes of 16 primitive and traditional accessions of maize (corn) selected by indigenous peoples across the Americas, along with the modern U.S. inbred B73. The maize pollen microbiome has not previously been reported. The pollen microbiomes were identified using full-length (FL) 16S rRNA gene PacBio SMRT sequencing compared to V4-MiSeq. The Pan-American maize pollen microbiome encompasses 765 taxa spanning 39 genera and 46 species, including known plant growth promoters, insect-obligates, plant pathogens, nitrogen-fixers and biocontrol agents. Eleven genera and 13 species composed the core microbiome. Of 765 taxa, 63% belonged to only four genera: 28% were Pantoea, 15% were Lactococcus, 11% were Pseudomonas, and 10% were Erwinia. Interestingly, of the 215 Pantoea taxa, 180 belonged to a single species, P. ananatis. Surprisingly, the diversity within P. ananatis ranged nearly 10-fold amongst the maize accessions analyzed (those with ≥3 replicates), despite being grown in a common field. The highest diversity within P. ananatis occurred in accessions that originated near the center of diversity of domesticated maize, with reduced diversity associated with the north-south migration of maize. This sub-species diversity was revealed by FL-PacBio but missed by V4-MiSeq. V4-MiSeq also mis-identified some dominant genera captured by FL-PacBio. The study, though limited to a single season and common field, provides initial evidence that pollen microbiomes reflect evolutionary and migratory relationships of their host plants.

Comparison of Nanopore and Synthesis-Based Next-Generation Sequencing Platforms for SARS-CoV-2 Variant Monitoring in Wastewater

Shortly after the beginning of the SARS-CoV-2 pandemic, many countries implemented sewage sentinel systems to monitor the circulation of the virus in the population. A fundamental part of these surveillance programs is the variant tracking through sequencing approaches to monitor and identify new variants or mutations that may be of importance. Two of the main sequencing platforms are Illumina and Oxford Nanopore Technologies. Here, we compare the performance of MiSeq (Illumina) and MinION (Oxford Nanopore Technologies), as well as two different data processing pipelines, to determine the effect they may have on the results. MiSeq showed higher sequencing coverage, lower error rate, and better capacity to detect and accurately estimate variant abundances than MinION R9.4.1 flow cell data. The use of different variant callers (LoFreq and iVar) and approaches to calculate the variant proportions had a remarkable impact on the results generated from wastewater samples. Freyja, coupled with iVar, may be more sensitive and accurate than LoFreq, especially with MinION data, but it comes at the cost of having a higher error rate. The analysis of MinION R10.4.1 flow cell data using Freyja combined with iVar narrows the gap with MiSeq performance in terms of read quality, accuracy, sensitivity, and number of detected mutations. Although MiSeq should still be considered as the standard method for SARS-CoV-2 variant tracking, MinION's versatility and rapid turnaround time may represent a clear advantage during the ongoing pandemic.

Evidence of the Involvement of Nematostoma parasiticum in Herpotrichia Needle Browning of Abies balsamea

Herpotrichia needle browning (HNB) is a disease that affects several species of fir trees in Europe and North America. HNB was first described by Hartig in 1884, who isolated a fungal pathogenic agent identified as responsible for the disease. This fungus was later named Herpotrichia parasitica but is currently named Nematostoma parasiticum. However, the identity of the pathogens causing HNB is regularly questioned and, to date, the true causal agent of this disease has not been definitely established. The present study aimed to identify the fungal populations present in needles of Christmas fir trees (Abies balsamea) and to correlate them with needle health status using robust molecular methods. PCR primers specific to N. parasiticum allowed detection of the presence of this fungus in DNA samples from symptomatic needles. Furthermore, high-throughput sequencing (Illumina MiSeq) clearly showed that N. parasiticum was associated with symptomatic needles. However, high-throughput sequencing results revealed that the presence of other species such as Sydowia polyspora and Rhizoctonia sp. may also correlate with the development of HNB. A diagnostic tool, based on quantitative PCR using a probe, was then developed to detect and quantify N. parasiticum in DNA samples. The efficacy of this molecular approach was validated through the detection of the pathogenic agent in symptomatic needle samples as well as in nonsymptomatic needles collected in trees affected by HNB. In contrast, N. parasiticum could not be found in needles from healthy trees. The present study argues for the importance of N. parasiticum in causing HNB symptoms.

Degradation of starch-based bioplastic bags in the pelagic and benthic zones of the Gulf of Oman

The Gulf of Oman is becoming increasingly polluted with plastics, hence bioplastics have been considered 'a substitute', although their biodegradability in marine environments has not been well investigated. Most research has been performed on cellulose-based bioplastics, whereas starch-based bioplastics have proven to be a suitable, but less researched, alternative. This study is the first of its kind designed to investigate the degradability of two different types of starch-based bioplastic bags, available in the market and labeled as "biodegradable", in the pelagic and benthic zones of one of the warmest marine environment in the world. Fourier-Transform Infrared Spectroscopy (FTIR) showed a clear reduction in the presence of OH, CH, and CO in the bioplastic bags after 5 weeks of immersion. Thermo-Gravimetric Analysis (TGA) indicated degradation of glycerol, starch, and polyethylene. The biofouling bacterial communities on bioplastic surfaces showed distinct grouping based on the immersion zone. Candidaatus saccharibacteria, Verrucomicrobiae, Acidimicrobiia and Planctomycetia sequences were only detectable on bioplastics in the pelagic zone, whereas Actinomyces, Pseudomonas, Sphingobium and Acinetobacter related sequences were only found on bioplastics in the benthic layer. We conclude that starch-based bioplastics are more readily degradable in the Gulf of Oman than conventional plastics, hence could serve as a better environmentally friendly alternative.

Occurrence, Distribution and Toxins of Benthic Cyanobacteria in German Lakes

Cyanobacteria are favored by climate change and global warming; however, to date, most research and monitoring programs have focused on planktic cyanobacteria. Benthic cyanobacteria blooms also increase and pose a risk to animal and human health; however, there is limited knowledge of their occurrence, distribution and the toxins involved, especially in relation to their planktic conspecifics. Therefore, we analyzed the benthic and planktic life forms of cyanobacterial communities in 34 lakes in Germany, including a monitoring of cyanotoxins. Community analyses were based on microscopic examination and Illumina sequencing of the 16S rRNA gene. The analyses of cyanotoxins were carried out using LC-MS/MS and ELISA. Observed benthic mats containing cyanobacteria consisted mainly of Nostocales and Oscillatoriales, being present in 35% of the lakes. Anatoxin was the most abundant cyanotoxin in the benthic samples, reaching maximum concentrations of 45,000 µg/L, whereas microcystin was the predominate cyanotoxin in the open-water samples, reaching concentrations of up to 18,000 µg/L. Based on the results, specific lakes at risk of toxic cyanobacteria could be identified. Our findings suggest that monitoring of benthic cyanobacteria and their toxins should receive greater attention, ideally complementing existing open-water sampling programs with little additional effort.

Evaluation of the Illumina iSeq whole genome sequencing system for enteric disease surveillance and outbreak detection

The Illumina iSeq low-capacity sequencing platform was evaluated for use in foodborne disease surveillance and outbreak detection. The platform produced high quality sequence data comparable to that of the Illumina MiSeq and was cost-effective with fast turn-around time in low sample volume environments.

Metagenomics and metabarcoding experimental choices and their impact on microbial community characterization in freshwater recirculating aquaculture systems

BACKGROUND

Microbial communities in recirculating aquaculture systems (RAS) play a role in system success, nutrient cycling, and water quality. Considering the increasing socio-economic role of fish farming, e.g., regarding food security, an in-depth understanding of aquaculture microbial communities is also relevant from a management perspective, especially regarding the growth, development, and welfare of the farmed animal. However, the current data on the composition of microbial communities within RAS is patchy, which is partly attributable to diverging method choices that render comparative analyses challenging. Therefore, there is a need for accurate, standardized, and user-friendly methods to study microbial communities in aquaculture systems.

RESULTS

We compared sequencing approach performances (3 types of 16S short amplicon sequencing, PacBio long-read amplicon sequencing, and amplification-free shotgun metagenomics) in the characterization of microbial communities in two commercial RAS fish farms. Results showed that 16S primer choice and amplicon length affect some values (e.g., diversity measures, number of assigned taxa or distinguishing ASVs) but have no impact on spatio-temporal patterns between sample types, farms and time points. This implies that 16S rRNA approaches are adequate for community studies. The long-read amplicons underperformed regarding the quantitative resolution of spatio-temporal patterns but were suited to identify functional services, e.g., nitrification cycling and the detection of pathogens. Finally, shotgun metagenomics extended the picture to fungi, viruses, and bacteriophages, opening avenues for exploring inter-domain interactions. All sequencing datasets agreed on major prokaryotic players, such as Actinobacteriota, Bacteroidota, Nitrospirota, and Proteobacteria.

CONCLUSION

The different sequencing approaches yielded overlapping and highly complementary results, with each contributing unique data not obtainable with the other approaches. We conclude that a tiered approach constitutes a strategy for obtaining the maximum amount of information on aquaculture microbial communities and can inform basic research on community evolution dynamics. For specific and/or applied questions, single-method approaches are more practical and cost-effective and could lead to better farm management practices.

Overview of microbial profiles in human hepatocellular carcinoma and adjacent nontumor tissues

BACKGROUND

Intratumoral microbial communities have been recently discovered to exist in a variety of cancers and have been found to be intricately involved in tumour progression. Therefore, investigating the profiles and functions of intratumoral microbial distribution in hepatocellular carcinoma (HCC) is imperative.

METHODS

To verify the presence of microorganisms in HCC, we performed fluorescence in situ hybridization (FISH) using HCC tissues and conducted MiSeq using 99 HCC and paracancerous tissues to identify the key microorganisms and changes in metabolic pathways affecting HCC progression through a variety of bioinformatics methods.

RESULTS

Microbial diversity was significantly higher in HCC tissues than in adjacent tissues. The abundances of microorganisms such as Enterobacteriaceae, Fusobacterium and Neisseria were significantly increased in HCC tissues, while the abundances of certain antitumour bacteria such as Pseudomonas were decreased. Processes such as fatty acid and lipid synthesis were significantly enhanced in the microbiota in HCC tissues, which may be a key factor through which intratumoral microbes influence tumour progression. There were considerable differences in the microbes and their functions within tumour tissue collected from patients with different clinical features.

CONCLUSION

We comprehensively evaluated the intratumoral microbial atlas of HCC tissue and preliminarily explored the mechanism of the effects of the microbial community involving changes in lipid metabolism and effects on HCC progression, which lays the foundation for further research in this field.

Fungal and Bacterial Diversity in the Tuber magnatum Ecosystem and Microbiome

Fungi belonging to the genus Tuber produce edible ascocarps known as truffles. Tuber magnatum Picco may be the most appreciated truffle species given its peculiar aroma. While its life cycle is not yet fully elucidated, some studies demonstrated an active role of microorganisms. The main goal of this study was to determine how the T. magnatum microbiome varies across space and time. To address this, we characterized microbial communities associated with T. magnatum through high-throughput amplicon sequencing of internal transcribed spacer (ITS) and 16S rDNAs in three productive natural sites in Italy across 2 years. At each site, four truffles were sampled as well as the soil underneath and at 40, 100, and 200 cm from the harvesting points, to assess for microbial variation between substrates, years, and sites. A statistically significant site-related effect on microbial communities was identified, whereas only the prokaryotic community was significantly affected by the distance of soil from the truffle. Significant differences between sampling years were also found, demonstrating a possible relation among rainfall precipitation and Firmicutes and Actinobacteria. Thirty-six bacterial OTUs in truffles and 11 bacterial OTUs in soils beneath truffles were identified as indicator taxa. As shown for other truffle species, the dominance of Bradyrhizobium, Rhizobium, and Ensifer spp. within the truffle fruiting body suggests an evolutionary adaptation of this microorganism to the genus Tuber. The present work offers novel and relevant insights into the microbial ecology of T. magnatum ecosystems and fruiting bodies. The function and role of these bacteria in the truffle microbiome and life cycle are in need of further investigation.

Emodin modulates gut microbial community and triggers intestinal immunity

BACKGROUND

The gut microbiota (GM) plays an important role in human health and is being investigated as a possible target for new therapies. Although there are many studies showing that emodin can improve host health, emodin-GM studies are scarce. Here, the effects of emodin on the GM were investigated in vitro and in vivo.

RESULTS

In vitro single bacteria cultivation showed that emodin stimulated the growth of beneficial bacteria Akkermansia, Clostridium, Roseburia, and Ruminococcus but inhibited major gut enterotypes (Bacteroides and Prevotella). Microbial community analysis from a synthetic gut microbiome model through co-culture indicated the consistent GM change by emodin. Interestingly, emodin stimulated Clostridium and Ruminococcus (which are related to Roseburia and Faecalibacterium) in a mice experiment and induced anti-inflammatory immune cells, which may correlate with its impact on specific gut bacteria.

CONCLUSION

Emodin (i) showed similar GM changes in monoculture, co-culture, and in an in vivo mice experiment and (ii) simulated regulatory T-cell immune responses in vivo. This suggest that emodin may be used to modulate the GM and improve health. © 2022 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Microbial community and predictive functionalities associated with the marine sediment of Coastal Gujarat

Marine sediments are complex ecosystems where structures and functions constantly change due to natural and anthropogenic influences. In this investigation, a comprehensive and comparative analysis of the bacterial communities and their functional potential of the pristine and polluted marine sediments were carried out using MiSeq. The phylum Proteobacteria was dominant in all study sites. Other phyla were Actinobacteria, Bacteroidetes, Planctomycetes, Acidobacteria, Chloroflexi, Nitrospirae, Cyanobacteria, Verrucomicrobia, Tenericutes, and Chlorobi. Interestingly, about 50% of genera belong to the unclassified categories. The key genera were identified as Acinetobacter, Bacillus, Pseudomona, Idiomarina, Thalassospira, and Marinobacter, Halomonas, Planctomyces, Psychrobacter, and Vogesella. PICRUSt analysis revealed that major functions are associated with the metabolism category. Additionally, metabolism related to amino acids, carbohydrates, energy generation, xenobiotics degradation, nitrogen, sulfate, and methane were prominent. Similarly, the predicted metabolisms by COG and KEGG were observed in the microbial communities of the marine sediments. To date, a comprehensive description of the microbial life with metabolic potential in these study sites has not been investigated. This study therefore significantly adds to our understanding of the microbiome and its functional attributes of marine sediments.

Comparable outcomes from long and short read random sequencing of total RNA for detection of pathogens in chicken respiratory samples

Co-infections of avian species with different RNA viruses and pathogenic bacteria are often misdiagnosed or incompletely characterized using targeted diagnostic methods, which could affect the accurate management of clinical disease. A non-targeted sequencing approach with rapid and precise characterization of pathogens should help respiratory disease management by providing a comprehensive view of the causes of disease. Long-read portable sequencers have significant potential advantages over established short-read sequencers due to portability, speed, and lower cost. The applicability of short reads random sequencing for direct detection of pathogens in clinical poultry samples has been previously demonstrated. Here we demonstrate the feasibility of long read random sequencing approaches to identify disease agents in clinical samples. Experimental oropharyngeal swab samples (n = 12) from chickens infected with infectious bronchitis virus (IBV), avian influenza virus (AIV) and Mycoplasma synoviae (MS) and field-collected clinical oropharyngeal swab samples (n = 11) from Kenyan live bird markets previously testing positive for Newcastle disease virus (NDV) were randomly sequenced on the MinION platform and results validated by comparing to real time PCR and short read random sequencing in the Illumina MiSeq platform. In the swabs from experimental infections, each of three agents in every RT-qPCR-positive sample (Ct range 19-34) was detectable within 1 h on the MinION platform, except for AIV one agent in one sample (Ct = 36.21). Nine of 12 IBV-positive samples were assigned genotypes within 1 h, as were five of 11 AIV-positive samples. MinION relative abundances of the test agent (AIV, IBV and MS) were highly correlated with RT-qPCR Ct values (R range-0.82 to-0.98). In field-collected clinical swab samples, NDV (Ct range 12-37) was detected in all eleven samples within 1 h of MinION sequencing, with 10 of 11 samples accurately genotyped within 1 h. All NDV-positive field samples were found to be co-infected with one or more additional respiratory agents. These results demonstrate that MinION sequencing can provide rapid, and sensitive non-targeted detection and genetic characterization of co-existing respiratory pathogens in clinical samples with similar performance to the Illumina MiSeq.

Method comparison for Japanese encephalitis virus detection in samples collected from the Indo-Pacific region

Introduction

Japanese encephalitis virus (JEV) is a mosquito-borne viral pathogen, which is becoming a growing public health concern throughout the Indo-Pacific. Five genotypes of JEV have been identified. Current vaccines are based on genotype III and provide a high degree of protection for four of the five known genotypes.

Methods

RT-PCR, Magpix, Twist Biosciences Comprehensive Viral Research Panel (CVRP), and SISPA methods were used to detect JEV from mosquito samples collected in South Korea during 2021. These methods were compared to determine which method would be most effective for biosurveillance in the Indo-Pacific region.

Results

Our data showed that RT-PCR, Twist CVRP, and SISPA methods were all able to detect JEV genotype I, however, the proprietary Magpix panel was only able to detect JEV genotype III. Use of minION sequencing for pathogen detection in arthropod samples will require further method development.

Conclusion

Biosurveillance of vectorborne pathogens remains an area of concern throughout the Indo-Pacific. RT-PCR was the most cost effective method used in the study, but TWIST CVRP allows for the identification of over 3,100 viral genomes. Further research and comparisons will be conducted to ensure optimal methods are used for large scale biosurveillance.

Comparison of iSeq and MiSeq as the two platforms for 16S rRNA sequencing in the study of the gut of rat microbiome

Abstract

Amplicon-based next-generation sequencing (NGS) of the 16S ribosomal RNA (16S) regions is a culture-free method used to identify and analyze Procaryota occurring within a given sample. The prokaryotic 16S rRNA gene contains conserved regions and nine variable regions (V1-V9) frequently used for phylogenetic classification of genus or species in diverse microbial populations. This work compares the accuracy and efficacy of two platforms, iSeq and MiSeq from Illumina, used in sequencing 16S rRNA. The most important similarities and differences of 16S microbiome sequencing in 20 fecal rat samples were described. Genetic libraries were prepared according to 16S Metagenomic Sequencing Library Preparation (Illumina) for the V3 and V4 regions of the 16S. The species richness obtained using iSeq technology was lower compared to MiSeq. At the second taxonomy level (L2), the abundance of taxa was comparable for both platforms. At the L7, the taxa abundance was significantly different, and the number of taxa was higher for the MiSeq. The alpha diversity was lower for iSeq than for MiSeq, starting from the order to the species level. The beta diversity estimation revealed statistically significant differences in microbiota diversity starting from the class level to the species level in samples sequenced on two investigated platforms. This work disclosed that the iSeq platform could be used to evaluate the bacterial profile of the samples to characterize the overall profile. The MiSeq System seems to be better for a detailed analysis of the differences in the microbiota composition.

Key points

• iSeq platform allows to shorten the sequencing time three times compared to the MiSeq.

• iSeq can only be used for an initial and quick microbiome assessment.

• MiSeq is better for a detailed analysis of the differences in the microbiota composition.

Supplementary information

The online version contains supplementary material available at 10.1007/s00253-022-12251-z.

Exploring bacterial communities through metagenomics during bioremediation of polycyclic aromatic hydrocarbons from contaminated sediments

The goal of this study was to evaluate the degradation effectiveness of PAHs degrading bacteria at the mesocosm level, including Stenotrophomonas maltophilia (SC), mixed culture (MC), and enriched native microflora (EC) at the mesocosm level. Maximum degradation was found in the mesocosm MC (26.67 %), followed by SC (25.08 %) and EC (18.25 %) after 60 days. Thus, mixed culture and Stenotrophomonas maltophilia could be a game changer in the PAHs bioremediation at the chronically contaminated sites. MiSeq sequencing has revealed dominancy of γ-Proteobacteria, α-Proteobacteria, β-Proteobacteria at class level and Sphingomonadales, oceanospirillales, Rhodothermales at Order level. Families Alcanivoracaceae, Alteromonadaceae, Nocardiaceae, Rhodospirillaceae and genus Stenotrophomonas, Alcanivorax, Methylophaga, Fluviicola and Rhodoplanes were considerably increased which play key role in the PAHs degradation. Dominant bacterial communities have revealed resilience community to enable potential PAHs degradation process in all the mesocosms. To the best our knowledge this is the first ever attempt in PAHs biodegradation study conducted at the mesocosm level mimicking natural environmental conditions. Consequently, this study could be a benchmark against which future progress studies for the policy makers and stakeholders to design appropriate bioremediation study for the historically PAHs polluted contaminate sites.

Development and characterization of microsatellite primers for Triops granarius (Branchiopoda: Notostraca) using MiSeq technology

Background

Next-generation sequencing technology has allowed for the rapid development of microsatellites, neutral polymorphic markers that can be used for the analysis of population structure.

Methods and Results

In this study, we performed whole-genome sequencing using the Illumina MiSeq system and de novo assembly to design microsatellite primers for Triops granarius populations in Qatar. The developed microsatellites are suitable for future studies of genetic structuring among geographically isolated freshwater pools. A total of 23 different primer pairs produced typical microsatellite results, with each pair successfully amplified in up to 40 individuals. Only five of the loci produced a significant departure from Hardy-Weinberg equilibrium.

Conclusions

Some of the underlying mechanisms regarding the few loci that deviated from HWE may be further investigated to determine the source of deviation. As T. granarius is the most widely distributed species of the family, the development of these molecular markers would be useful for conducting population genetics and biogeographical studies broadly.

Shifts in the intestinal microflora of meat rabbits in response to glucocorticoids

BACKGROUND

As major stress hormones, glucocorticoids can directly or indirectly affect the intestinal microflora, although few studies have focused on changes in the composition of the intestinal microflora. In this study, rabbits were randomly divided into two groups: gavage administration with saline, and the same doses of dexamethasone (1 mg kg^-1 ). After 7 days, the microbial diversity of the jejunum contents was analysed.

RESULTS

The gut microflora richness and diversity had no significant difference between the two groups. The proportions of Firmicutes and Bacteroidetes were the most abundant in the jejunum of meat rabbits. Dexamethasone injection led to a change in the structure of the gut microflora composition, and we found that there were six biomarkers with linear discriminant analysis score >4 (Firmicutes, Caproiciproducens, Clostridiales, Clostridia, Psychrobacter, and Psychrobacter faecalis), moreover, the results of this study provide new insight into alleviating the effects of stress on meat rabbits.

CONCLUSION

It was concluded that glucocorticoids caused changes in the composition of intestinal microflora. © 2022 Society of Chemical Industry.

Insights into the bacterial community structure of marine mucilage by metabarcoding

This prospective study was aimed to explore the bacterial diversity of marine mucilage developed in the Marmara Sea and the North Aegean Sea by metabarcoding. For this purpose, mucilage samples were collected from five different sampling locations, and the bacterial community structure was analyzed by 16S rRNA gene amplicon sequencing. The results highlighted a diverse bacterial community dominated by Proteobacteria and Bacteroidetes species. A negative and significant correlation between pH level and Campylobacterales, Clostridiales, and Vibronales abundances was detected, while a strong positive correlation was determined between total phosphorus (TP) and Campylobacterales. Results revealed that the bacterial community in the mucilage samples was predominated by particle-attached species preferring high-nutrient concentrations. This is the first study evaluating the bacterial diversity in a mucilage outbreak using a metabarcoding approach. Its results may contribute to this growing area of research and provide a database for further studies.

16S rRNA Based Profiling of Bacterial Communities Colonizing Bakery-Production Environments

Conventional culture-based techniques are largely inadequate in elucidating the microbiota contained in an environment, due to low recovery within a complex bacterial community. This limitation has been mitigated by the use of next-generation sequencing (NGS)-based approaches thereby facilitating the identification and classification of both culturable and uncultivable microorganisms. Amplicon targeted NGS methods, such as 16S ribosomal RNA (16S rRNA) and shotgun metagenomics, are increasingly being applied in various settings such as in food production environments to decipher the microbial consortium therein. Even though multiple food matrices/food production environments have been studied, the low-moisture environment associated with bakery food production remains to be investigated. To address this knowledge gap, in this study, we investigated the microbiome associated with two bakery production sites (designated as A and B) located in Ireland using 16S rRNA-amplicon-based sequencing. Amplicons corresponding to a hypervariable region contained within the 16S rRNA gene were amplified from DNA samples purified from environmental swabs and ingredients collected at both sites at various stages (preparation, production, postproduction, and storage) across the bakery production chain, over three seasons (winter, spring, and summer). These amplicons were sequenced, and data were analyzed using the mothur pipeline and visualized using MicrobiomeAnalyst and a series of R packages. The top seven bacterial phyla identified at both sites were composed of Proteobacteria, Firmicutes, Actinobacteria, Bacteroidetes, Deinococcus-Thermus, Patescibacteria, and Verrucomicrobia. In addition, the phyla Tenericutes (Mycoplasmatota) and Acidobacteria were observed only in samples taken at site B. Different bacterial compositions were identified at each stage of production. These same bacteria were also found to be present in the final processed food suggesting the influence of the environment on the food matrix. This study is the first demonstration of the utility of 16S rRNA amplicon-based sequencing to describe the microbiota associated with bakery processing environments.

Speeding up the detection of invasive bivalve species using environmental DNA: a Nanopore and Illumina sequencing comparison

Traditional detection of aquatic invasive species via morphological identification is often time-consuming and can require a high level of taxonomic expertise, leading to delayed mitigation responses. Environmental DNA (eDNA) detection approaches of multiple species using Illumina-based sequencing technology have been used to overcome these hindrances, but sample processing is often lengthy. More recently, portable nanopore sequencing technology has become available, which has the potential to make molecular detection of invasive species more widely accessible and substantially decrease sample turnaround times. However, nanopore-sequenced reads have a much higher error rate than those produced by Illumina platforms, which has so far hindered the adoption of this technology. We provide a detailed laboratory protocol and bioinformatic tools (msi package) to increase the reliability of nanopore sequencing to detect invasive species, and we test its application using invasive bivalves while comparing it with Illumina-based sequencing. We sampled water from sites with pre-existing bivalve occurrence and abundance data, and contrasting bivalve communities, in Italy and Portugal. Samples were extracted, amplified, and sequenced by the two platforms. The mean agreement between sequencing methods was 69% and the difference between methods was non-significant. The lack of detections of some species at some sites could be explained by their known low abundances. This is the first reported use of MinION to detect aquatic invasive species from eDNA samples.

Exploring statistical weight estimates for mitochondrial DNA matches involving heteroplasmy

Massively parallel sequencing (MPS) of mitochondrial (mt) DNA allows forensic laboratories to report heteroplasmy on a routine basis. Statistical approaches will be needed to determine the relative frequency of observing an mtDNA haplotype when including the presence of a heteroplasmic site. Here, we examined 1301 control region (CR) sequences, collected from individuals in four major population groups (European, African, Asian, and Latino), and covering 24 geographically distributed haplogroups, to assess the rates of point heteroplasmy (PHP) on an individual and nucleotide position (np) basis. With a minor allele frequency (MAF) threshold of 2%, the data was similar across population groups, with an overall PHP rate of 37.7%, and the majority of heteroplasmic individuals (77.3%) having only one site of heteroplasmy. The majority (75.2%) of identified PHPs had an MAF of 2-10%, and were observed at 12.6% of the nps across the CR. Both the broad and phylogenetic testing suggested that in many cases the low number of observations of heteroplasmy at any one np results in a lack of statistical association. The posterior frequency estimates, which skew conservative to a degree depending on the sample size in a given haplogroup, had a mean of 0.152 (SD 0.134) and ranged from 0.031 to 0.83. As expected, posterior frequency estimates decreased in accordance with 1/n as the sample size (n) increased. This provides a proposed conservative statistical framework for assessing haplotype/heteroplasmy matches when applying an MPS technique in forensic cases and will allow for continual refinement as more data is generated, both within the CR and across the mitochondrial genome.

Different Urban Forest Tree Species Affect the Assembly of the Soil Bacterial and Fungal Community

The selection of tree species used for the afforestation of urban forests is very important for maintaining the urban ecosystem, while soil microbe is one of the driving factors of material cycling in the urban forest ecosystem and for health of forests. In this study, the characteristics of surface soil bacterial and fungal community structure in four urban forests (primarily composed of Fraxinus mandshurica (Fm), Quercus mongolica (Qm), Pinus sylvestris var. mongolica (Ps), and Pinus tabulaeformis var. Mukdensis (Pt) as the main dominant tree species, respectively) were investigated by high-throughput sequencing. Our results showed that the alpha diversity of the soil microbial community in the Fm urban forest was the highest, while the lowest was in the Ps urban forest. In the bacterial community, Proteobacteria was the most predominant phylum in soils from Fm, Ps, and Pt urban forests. The most relatively abundant phylum of the Qm urban forest soil was Acidobacteria. The relative abundances of the bacterial communities at the genus level in the soil of four urban forests were significantly different. The soil bacterial communities in Ps and Pt urban forests were more similar, and Qm and Fm were also more similar. In the fungal community, Basidiomycota was the most predominant phylum in soils from Qm, Ps, and Pt urban forests. The phylum with the greatest relative abundance in the Fm urban forest soil was Ascomycota. There were differences in the fungal community between Qm, Fm, Ps, and Pt urban forests. Soil microbial community composition was affected by environmental factors: soil bacterial and fungal community compositions were significantly related to soil electrical conductivity (EC), alkali hydrolysable nitrogen (AHN), total nitrogen (TN), and total phosphorus (TP). In conclusion, the soil microbial community structure was related to both forest's tree species and soil properties.

Estimating individual mtDNA haplotypes in mixed DNA samples by combining MinION and MiSeq

We tried to estimate individual mtDNA haplotypes in mixed DNA samples by combining MinION and MiSeq. The BAM files produced by MiSeq were viewed using Integrative Genomics Viewer (IGV) to verify mixed bases. By sorting the reads according to base type for each mixed base, partial haplotypes were determined. Then, the BAM files produced by MinKNOW were viewed using IGV. To determine haplotypes with IGV, only mixed bases determined by MiSeq were used as target bases. By sorting the reads according to base type for each target base, each contributor's haplotype was estimated. In mixed samples from two contributors, even a haplotype with a minor contribution of 5% could be distinguished from the haplotype of the major contributor. In mixed samples of three contributors (mixture ratios of 1:1:1 and 4:2:1), each haplotype could also be distinguished. Sequences of C-stretches were determined very inaccurately in the MinION analysis. Although the analysis method was simple, each haplotype was correctly detected in all mixed samples with two or three contributors in various mixture ratios by combining MinION and MiSeq. This should be useful for identifying contributors to mixed samples.

Deep sequencing across multiple host species tests pine-endophyte specificity

PREMISE

Foliar fungal endophytes vary in their distributions across landscapes or plant host taxa, indicative of specialized ecologies and host specific adaptations. Accounts of specialization, however, depend on the taxonomic breadth and geographic range of the host plants included in each study. A broad region-scale study or deep sampling of diverse potential host species still remains relatively rare but is becoming increasingly possible with high-throughput sequencing.

METHODS

Amplicon sequencing was used to rapidly identify the fungal endophytic community among six pine (Pinus, Pinaceae) species co-occurring across northeastern United States and to test for site and host specialization. We focused on the endophytic genus Lophodermium (Rhytismataceae), whose species members are thought to specialize on different pine species, to test if amplicon sequencing could rapidly verify previously implied or discover new patterns of host specificity.

RESULTS

While amplicon sequencing could analyze more samples at greater depths and recover greater numbers of unique Lophodermium taxa than when endophyte communities were surveyed with traditional culturing methods, patterns of specialization were not better supported. This may be because amplicon sequencing can indiscriminately capture non-host specific organisms found incidentally from plant tissues or because we have overestimated host-specificity in the past with biased culturing techniques.

CONCLUSIONS

Amplicon sequencing can quickly identify patterns of host specificity by allowing large-scale surveys but has limitations in quantifying the level of intimacy of these relationships.

Variations in bacterial and archaeal community structure and diversity along the soil profiles of a peatland in Southwest China

As bacteria and archaea are key components in the ecosystem, information on their dynamics in soil profiles is important for understanding the biogeochemical cycles in peatlands. However, little is known about the vertical distribution patterns of bacteria and archaea in the Bitahai peatland, or about their relationships with soil chemical properties. Here, bacterial and archaeal abundance, diversity, and composition of the Bitahai peatlands at 0-100 cm soil depths were analyzed by sequencing of 16S rRNA genes (Illumina, MiSeq). Soil pH, total C, N, and P concentrations and stoichiometric ratios were also estimated. The results revealed that total C and total N contents, as well as C:P and N:P ratios, significantly increased with increasing peatland soil depths, while total P decreased. The top three dominant phyla were Proteobacteria (39.64%), Acidobacteria (12.93%), and Chloroflexi (12.81%) in bacterial communities, and were Crenarchaeota (58.67%), Thaumarchaeota (14.34%), and Euryarchaeota (10.82%) in archaeal communities in the Bitahai peatland, respectively. The total relative abundance of methanogenic groups and ammonia-oxidizing microorganisms all significantly decreased with soil depth. Both bacterial and archaeal diversities were significantly affected by the soil depth. Soil C, N, and P concentrations and stoichiometric ratios markedly impacted the community structure and diversity in archaea, but not in bacteria. Therefore, these results highlighted that the microbial community structure and diversity depended on soil depth for the Bitahai peatlands, and the factors affecting bacteria and archaea in the Bitahai peatlands were different.

An experimental strategy for preparing circular ssDNA virus genomes for next-generation sequencing

The ability of begomoviruses to evolve rapidly threatens many crops and underscores the importance of detecting these viruses quickly and to understand their genome diversity. This study presents an improved protocol for the enhanced amplification and enrichment of begomovirus DNA for use in next generation sequencing of the viral genomes. An enhanced rolling circle amplification (RCA) method using EquiPhi29 polymerase was combined with size selection to generate a cost-effective, short-read sequencing method. This improved short-read sequencing produced at least 50 % of the reads mapping to the target viral reference genomes, African cassava mosaic virus and East African cassava mosaic virus. This study provided other insights into common misconceptions about RCA and lessons that could be learned from the sequencing of single-stranded DNA virus genomes. This protocol can be used to examine the viral DNA as it moves from host to vector, thus producing valuable information for viral DNA population studies, and would likely work well with other circular Rep-encoding ssDNA viruses (CRESS) DNA viruses.

A novel high-throughput sequencing approach reveals the presence of a new virus infecting Rosa: rosa ilarvirus-1 (RIV-1)

Roses are one of the most valuable ornamental flowering shrubs grown worldwide. Despite the widespread of rose viruses and their impact on cultivation, they have not been studied in detail in the United Kingdom (UK) since the 1980's. As part of a survey of rose viruses entering the UK, 35 samples were collected at Heathrow Airport (London, UK) and were tested by RT-qPCR for different common rose viruses. Of the 35 samples tested using RT-qPCR for prunus necrotic ringspot virus (PNRSV; genus Ilarvirus), 10 were positive. Confirmatory testing was performed using RT-PCR with both PNRSV-specific and ilarvirus-generic primers, and diverse results were obtained: One sample was exclusively positive when using the ilarvirus-generic primers, and subsequent sequencing of the RT-PCR product revealed homology to other ilarviruses but not PNRSV. Further work to characterise the virus was performed using high throughput sequencing, both the MinION Flongle and Illumina MiSeq. The sequencing confirmed the presence of a new virus within group 2 of the genus Ilarvirus and we propose the name "rosa ilarvirus-1″ (RIV-1). Here, we describe the identification of a novel virus using the low-cost Flongle flow cell and discuss its potential as a front-line diagnostic tool.

First Genomic Resource of the Columbia Lance Nematode Hoplolaimus columbus

The Columbia lance nematode Hoplolaimus columbus has been reported frequently from North America due to its negative impact on agricultural production. In this study, for the first time, we sequenced the whole genome of a female specimen by using whole-genome amplification and Illumina MiSeq. Data were de novo assembled to form scaffolds of 205.75 Mbp consisting of 118,374 contigs. The largest scaffold was 636,881 bp. Benchmarking Universal Single-Copy Orthologs completeness was 66.6% (eukaryotic dataset), and >8,000 unique genes were predicted by GeneMark-ES. In total, 61,855 protein sequences were predicted by AUGUSTUS, and 10,085 of them were annotated by PANNZER2 with at least one function. These data will provide valuable resources for studies focusing on pathogenicity and phylogenomics of plant-parasitic nematodes.

Whole-Genome Sequencing of SARS-CoV-2: Assessment of the Ion Torrent AmpliSeq Panel and Comparison with the Illumina MiSeq ARTIC Protocol

Fast and effective methods are needed for sequencing of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome to track genetic mutations and to identify new and emerging variants during the ongoing pandemic. The objectives were to assess the performance of the SARS-CoV-2 AmpliSeq research panel and S5 plug-in analysis tools for whole-genome sequencing analysis of SARS-CoV-2 and to compare the results with those obtained with the MiSeq-based ARTIC analysis pipeline, using metrics such as depth, coverage, and concordance of single-nucleotide variant (SNV) calls. A total of 191 clinical specimens and a single cultured isolate were extracted and sequenced with AmpliSeq technology and analysis tools. Of the 191 clinical specimens, 83 (with threshold cycle [C_T] values of 15.58 to 32.54) were also sequenced using an Illumina MiSeq-based method with the ARTIC analysis pipeline, for direct comparison. A total of 176 of the 191 clinical specimens sequenced on the S5XL system and prepared using the SARS-CoV-2 research panel had nearly complete coverage (>98%) of the viral genome, with an average depth of 5,031×. Similar coverage levels (>98%) were observed for 81/83 primary specimens that were sequenced with both methods tested. The sample with the lowest viral load (C_T value of 32.54) achieved 89% coverage using the MiSeq method and failed to sequence with the AmpliSeq method. Consensus sequences produced by each method were identical for 81/82 samples in areas of equal coverage, with a single difference present in one sample. The AmpliSeq approach is as effective as the Illumina-based method using ARTIC v3 amplification for sequencing SARS-CoV-2 directly from patient specimens across a range of viral loads (C_T values of 15.56 to 32.54 [median, 22.18]). The AmpliSeq workflow is very easily automated with the Ion Chef and S5 instruments and requires less training and experience with next-generation sequencing sample preparation than the Illumina workflow.

Impact of Soil Microbes and Oxygen Availability on Bacterial Community Structure of Decomposing Poultry Carcasses

The impact of soil with an intact microbial community and oxygen availability on moisture content, soil pH, and bacterial communities during decomposition of poultry carcasses was investigated. Poultry carcasses were decomposed in soil with or without a microbial community, under aerobic or anaerobic conditions. The samples collected in each microcosm burial set-up were analyzed by targeted 16S rRNA amplicon sequencing and Amplicon sequence variants (ASV) method. Our results showed that moisture was high in the burial set-ups under anaerobic conditions and pH was high in the burial set-ups under aerobic conditions. Meanwhile, the Chao1 and Shannon index significantly differed between the different burial set-ups and across different time points. In addition, bacterial taxa composition during the early period of decomposition differed from that of the late period. A total of 23 phyla, 901 genera, and 1992 species were identified. Firmicutes was the most dominant phyla in all burial set-ups throughout the decomposition. At day 60, Pseudogracilibacillus was dominant in the burial set-ups under aerobic conditions, while Lentibacillus dominated in the burial set-ups under anaerobic conditions. This study demonstrated that the soil microbial community and availability of oxygen significantly affected the changes in moisture content, pH, and bacterial composition during the decomposition process.

Feasting, not fasting: winter diets of cave hibernating bats in the United States

Temperate bat species use extended torpor to conserve energy when ambient temperatures are low and food resources are scarce. Previous research suggests that migratory bat species and species known to roost in thermally unstable locations, such as those that roost in trees, are more likely to remain active during winter. However, hibernating colonies of cave roosting bats in the southeastern United States may also be active and emerge from caves throughout the hibernation period. We report what bats are eating during these bouts of winter activity. We captured 2,044 bats of 10 species that emerged from six hibernacula over the course of 5 winters (October–April 2012/2013, 2013/2014, 2015/2016, 2016/2017, and 2017/2018). Using Next Generation sequencing of DNA from 284 fecal samples, we determined bats consumed at least 14 Orders of insect prey while active. Dietary composition did not vary among bat species; however, we did record variation in the dominant prey items represented in species’ diets. We recorded Lepidoptera in the diet of 72.2% of individual Corynorhinus rafinesquii and 67.4% of individual Lasiurus borealis. Diptera were recorded in 32.4% of Myotis leibii, 37.4% of M. lucifugus, 35.5% of M. sodalis and 68.8% of Perimyotis subflavus. Our study is the first to use molecular genetic techniques to identify the winter diet of North American hibernating bats. The information from this study is integral to managing the landscape around bat hibernacula for insect prey, particularly in areas where hibernating bat populations are threatened by white-nose syndrome.

Phytate and Microbial Suspension Amendments Increased Soybean Growth and Shifted Microbial Community Structure

Phytate represents an organic pool of phosphorus in soil that requires hydrolysis by phytase enzymes produced by microorganisms prior to its bioavailability by plants. We tested the ability of a microbial suspension made from an old growth maple forest's undisturbed soil to mineralize phytate in a greenhouse trial on soybean plants inoculated or non-inoculated with the suspension. MiSeq Amplicon sequencing targeting bacterial 16S rRNA gene and fungal ITS was performed to assess microbial community changes following treatments. Our results showed that soybean nodulation and shoot dry weight biomass increased when phytate was applied to the nutrient-poor substrate mixture. Bacterial and fungal diversities of the root and rhizosphere biotopes were relatively resilient following inoculation by microbial suspension; however, bacterial community structure was significantly influenced. Interestingly, four arbuscular mycorrhizal fungi (AMF) were identified as indicator species, including Glomus sp., Claroideoglomus etunicatum, Funneliformis mosseae and an unidentified AMF taxon. We also observed that an ericoid mycorrhizal taxon Sebacina sp. and three Trichoderma spp. were among indicator species. Non-pathogenic Planctobacteria members highly dominated the bacterial community as core and hub taxa for over 80% of all bacterial datasets in root and rhizosphere biotopes. Overall, our study documented that inoculation with a microbial suspension and phytate amendment improved soybean plant growth.

Impact of Dietary Protein Content on Soil Bacterial and Fungal Communities in a Rice-Crab Co-culture System

Although co-culture of paddy fields with aquatic animals is lucrative, the ecological impacts of high-protein content entering the agricultural soil via animal pellet feed and feces have not been well studied. Moreover, the effects of dietary protein on soils and soil microbes remain unclear. To elucidate this, we examined soil bacterial and fungal community composition and temporal changes in paddy fields subjected to different protein-content diets via 16S/18S rRNA gene amplicon sequencing analysis with a high-throughput next-generation sequencer. MiSeq sequencing revealed that protein content significantly impacted fungal community structure. High-protein diets reduced bacterial community diversity and relative abundance in both July and October. The phylum-level bacterial taxonomic composition was not affected by diet treatment, while in fungi, a major phylum-level shift was evident. Hierarchically clustered analysis showed that high-protein diets significantly reduced the relative abundance of Brevundimonas in both July and October. Saprotrophic macrofungal diversity was negatively related to dietary protein content. Considering microbial community structure and environmental factors, ca. 15% protein content is appropriate for the rice-crab co-culture system that we studied.

Long-term Impact of Gold and Platinum on Microbial Diversity in Australian Soils

The effects of platinum (Pt) and gold (Au) and on the soil bacterial community was evaluated in four different Australian soil types (acidic Burn Grounds (BGR), organic matter-rich Fox Lane, high silt/metal Pinpinio (PPN), and alkali Minnipa (MNP) spiked with either Pt or Au at 1, 25, and 100 mg kg^-1 using a next-generation sequencing approach (amplicon-based, MiSeq). Soil type and metal concentrations were observed to be key drivers of Pt and Au effects on soil microbial community structure. Different trends were therefore observed in the response of the bacterial community to Pt and Au amendments; however in each soil type, Pt and Au amendment caused a detectable shift in community structure that in most samples was positively correlated with increasing metal concentrations. New dominant groups were only observed in BGR and PPN soils at 100 mg kg^-1 (Kazan-3B-28 and Verrucomicrobia groups (BGR, Pt) and Firmicutes and Caldithrix groups (PPN, Pt) and WS2 (BGR, Au). The effects of Pt on soil microbial diversity were largely adverse at 100 mg kg^-1 and were pronounced in acidic, basic, and metal/silt-rich soils. However, this effect was concentration-related; Au appeared to be more toxic to soil bacterial communities than Pt at 25 mg kg^-1 but Pt was more toxic at 100 mg kg^-1. More bacterial groups such as those belonging to Burkholderiales/Burkholderiaceae, Alicyclobacillaceae, Rubrobacteraceae, Cytophagaceae, Oxalobacteraceae were selectively enriched by Pt compared to Au (Sphingomonadaceae and Rhodospirillaceae) amendments irrespective of soil type. The research outcomes have important implications in the management (remediation) of Pt- and Au-contaminated environments.

Human Mitochondrial Control Region and mtGenome: Design and Forensic Validation of NGS Multiplexes, Sequencing and Analytical Software

Forensic mitochondrial DNA (mtDNA) analysis conducted using next-generation sequencing (NGS), also known as massively parallel sequencing (MPS), as compared to Sanger-type sequencing brings modern advantages, such as deep coverage per base (herein referred to as read depth per base pair (bp)), simultaneous sequencing of multiple samples (libraries) and increased operational efficiencies. This report describes the design and developmental validation, according to forensic quality assurance standards, of end-to-end workflows for two multiplexes, comprised of ForenSeq mtDNA control region and mtDNA whole-genome kits the MiSeq FGx^TM instrument and ForenSeq universal analysis software (UAS) 2.0/2.1. Polymerase chain reaction (PCR) enrichment and a tiled amplicon approach target small, overlapping amplicons (60–150 bp and 60–209 bp for the control region and mtGenome, respectively). The system provides convenient access to data files that can be used outside of the UAS if desired. Studies assessed a range of environmental and situational variables, including but not limited to buccal samples, rootless hairs, dental and skeletal remains, concordance of control region typing between the two multiplexes and as compared to orthogonal data, assorted sensitivity studies, two-person DNA mixtures and PCR-based performance testing. Limitations of the system and implementation considerations are discussed. Data indicated that the two mtDNA multiplexes, MiSeq FGx and ForenSeq software, meet or exceed forensic DNA quality assurance (QA) guidelines with robust, reproducible performance on samples of various quantities and qualities.

Characterization of the complete mitochondrial genome of the Northern Mud Gudgeon, Ophiocara porocephala (Perciformes: Eleotridae) with phylogenetic implications

The first mitochondrial genome of Ophiocara porocephala was determined by the combination of next-generation sequencing (NGS) and Sanger sequencing methods. A complete circular mitogenome of O. porocephala (16,529 bp) consisted of 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and two non-coding regions, including a control region (D-loop) and a light strand origin of replication (O_L). Two start codons (ATG and GTG) and four stop codons (TAG, TAA, TA-, and T-) were used in all the PCGs. Except for ND6 and eight transfer RNAs (tRNAs), all the other genes were encoded in the heavy strand. Based on phylogenetic analysis, O. porocephala formed a clade with three other species in the subfamily Butinae, while the other 10 made a subfamily Eleotrinae clade.

Ramadan Fasting Leads to Shifts in Human Gut Microbiota Structured by Dietary Composition

The structure and diversity of human gut microbiota are directly related to diet, though less is known about the influences of ethnicity and diet-related behaviors, such as fasting (intermittent caloric restriction). In this study, we investigated whether fasting for Ramadan altered the microbiota in Chinese and Pakistani individuals. Using high-throughput 16S rRNA gene sequencing and self-reported dietary intake surveys, we determined that both the microbiota and dietary composition were significantly different with little overlap between ethnic groups. Principal Coordinate Analyses (PCoA) comparison of samples collected from both groups before and after fasting showed partial separation of microbiota related to fasting in the Pakistani group, but not in the Chinese group. Measurement of alpha diversity showed that Ramadan fasting significantly altered the coverage and ACE indices among Chinese subjects, but otherwise incurred no changes among either group. Specifically, Prevotella and Faecalibacterium drove predominance of Bacteroidetes and Firmicutes in the Pakistani group, while Bacteroides (phylum Bacteroidetes) were the most prevalent among Chinese participants both before and after fasting. We observed significant enrichment of some specific taxa and depletion of others in individuals of both populations, suggesting that fasting could affect beta diversity. Notably, Dorea, Klebsiella, and Faecalibacterium were more abundant in the Chinese group after fasting, while Sutterella, Parabacteroides, and Alistipes were significantly enriched after fasting in the Pakistani group. Evaluation of the combined groups showed that genera Coprococcus, Clostridium_XlV, and Lachnospiracea were all significantly decreased after fasting. Analysis of food intake and macronutrient energy sources showed that fat-derived energy was positively associated with Oscillibacter and Prevotella, but negatively associated with Bacteroides. In addition, the consumption of sweets was significantly positively correlated with the prevalence of Akkermansia. Our study indicated that diet was the most significant influence on microbiota, and correlated with ethnic groups, while fasting led to enrichment of specific bacterial taxa in some individuals. Given the dearth of understanding about the impacts of fasting on microbiota, our results provide valuable inroads for future study aimed at novel, personalized, behavior-based treatments targeting specific gut microbes for prevention or treatment of digestive disorders.

Bacterial diversity changes in agricultural soils influenced by poultry litter fertilization

Poultry litter is widely applied as agricultural fertilizer and can affect the soil microbiome through nutrient overload and antibiotic contamination. In this study, we assessed changes in soil bacterial diversity using high-throughput sequencing approaches. Four samples in triplicate were studied: soils with short- and long-term fertilization by poultry litter (S1 = 10 months and S2 = 30 years, respectively), a soil inside a poultry shed (S3), and a forest soil used as control (S0). Samples S0, S1, and S2 revealed a relatively high richness, with confirmed operational taxonomic units (OTUs) in the three replicates of each sample ranging from 1243 to 1279, while richness in S3 was about three times lower (466). The most abundant phyla were Proteobacteria, Bacteroidetes, and Actinobacteria. Acidobacteria, Planctomycetes, and Verrucomicrobia were also abundant but highly diminished in S3, while Firmicutes was less abundant in S0. Changes in bacterial communities were very evident at the genera level. The genera Gaiella, Rhodoplanes, Solirubacter, and Sphingomonas were predominant in S0 but strongly decreased in the other soils. Pedobacter and Devosia were the most abundant in S1 and were diminished in S2, while Herbiconiux, Brevundimonas, Proteiniphilum, and Petrimonas were abundant in S2. The most abundant genera in S3 were Deinococcus, Truepera, Rhodanobacter, and Castellaniella. A predictive analysis of the metabolic functions with Tax4Fun2 software suggested the potential presence of enzymes associated with antibiotic resistance as well as with denitrification pathways, indicating that the S3 soil is a potential source of nitrous oxide, a powerful greenhouse gas.

Assessment of Illumina® Human mtDNA Genome assay: workflow evaluation with development of analysis and interpretation guidelines

Mitochondrial DNA (mtDNA) is a small but significant part of the human genome, whose applicability potential has gradually increased with the advent of massively parallel sequencing (MPS) technology. Knowledge of the particular workflow, equipment, and reagents used, along with extensive usage of negative controls to monitor all preparation steps constitute the prerequisites for confident reporting of results. In this study, we performed an assessment of Illumina® Human mtDNA Genome assay on MiSeq FGx™ instrument. Through analysis of several types of negative controls, as well as mtDNA positive controls, we established thresholds for data analysis and interpretation, consisting of several components: minimum read depth (220 reads), minimum quality score (41), percentage of minor allele sufficient for analysis (3.0%), percentage of minor allele sufficient for interpretation (6.0%), and percentage of major allele sufficient for homoplasmic variant call (97.0%). Based on these criteria, we defined internal guidelines for analysis and interpretation of mtDNA results obtained by MPS. Our study shows that the whole mtDNA assay on MiSeq FGx™ produces repeatable and reproducible results, independent of the analyst, which are also concordant with Sanger-type sequencing results for mtDNA control region, as well as with MPS results produced by NextSeq®. Overall, established thresholds and interpretation guidelines were successfully applied for the sequencing of complete mitochondrial genomes from high-quality samples. The underlying principles and proposed methodology on the definition of internal laboratory guidelines for analysis and interpretation of MPS results may be applicable to similar MPS workflows, e.g. targeting good-quality samples in forensic genetics and molecular diagnostics.

Blood meal source and mixed blood-feeding influence gut bacterial community composition in Aedes aegypti

BACKGROUND

The guts of blood-sucking insects host a community of bacteria that can shift dramatically in response to biotic and abiotic factors. Identifying the key factors structuring these microbial communities has important ecological and epidemiological implications.

METHODS

We used the yellow fever mosquito, Aedes aegypti, to investigate the impact of mixed blood meals on gut microbiota of vector mosquitoes. Adult females were experimentally fed on sugar or blood from chicken, rabbit or a mixture of chicken and rabbit blood, and their gut microbiota were characterized using 16S rRNA gene amplification and MiSeq sequencing.

RESULTS

The gut bacterial communities of mosquitoes fed on the three blood meal treatments clustered separately, suggesting that host species identity and mixed blood-feeding are key determinants of gut bacterial community composition in mosquitoes. Mixed blood meal had a synergistic effect on both operational taxonomic unit (OTU) richness and the Shannon diversity index, suggesting that mixed blood-feeding can offset the nutritional deficit of blood meals from certain host species. The microbial communities observed in this study were distinct from those identified from similarly fed Ae. aegypti from our previous study.

CONCLUSIONS

These findings demonstrate that vector host-feeding preferences can influence gut microbial composition and diversity, which could potentially impact pathogen acquisition and transmission by the vector. The results also demonstrate that different microenvironmental conditions within the laboratory may play an important role in structuring the microbial communities of independently reared mosquito colonies.

Treatments of porcine fecal samples affect high-throughput virome sequencing results

The porcine enteric microbiota is currently extensively studied, taking advantage of developments in high-throughput sequencing technologies. However, the viral part of the microbiota, the virome, is being lightly explored, and the impact of the pretreatments used before sequencing the viruses is barely considered. In this study, the impacts of filtration, RNase and DNase treatments on virus reads recovery and diversity after sequencing on a MiSeq platform were assessed on fecal samples individually taken at <3, 5, 12 and 20 weeks from two piglets. None of the four pretreatment series affected the virus read averages or influenced diversity, but the samples with the higher proportion of reads corresponding to an entry in the "nt" database were those receiving the least number of pretreatments. The enzymatic pretreatments affected the detection of the single-stranded RNA viruses of Aichivirus C, porcine astrovirus, Sapovirus and posavirus, which is worrisome, as these viruses can be involved in swine diarrhea. If enzymatic pretreatments are used when sequencing using a high-throughput method, it may impact single-stranded RNA virus recovery, but not the overall virome diversity. Therefore, filtrated samples may be the better option, reducing the amount of bacterial genetic material while preserving the virus reads.

Evaluation of a custom GeneRead™ massively parallel sequencing assay with 210 ancestry informative SNPs using the Ion S5™ and MiSeq platforms

A custom GeneRead DNAseq SNP panel with 210 markers was evaluated using the Ion S5 and MiSeq sequencing platforms. Sensitivity, PCR cycle number, and the use of half volume of reagents for target enrichment and library preparation were tested. Furthermore, genotype concordance between results obtained with the different sequencing platforms and with known profiles generated using other sequencing assays was analysed. The GeneRead DNASeq SNP assay gave reproducible results with an input of 200 pg DNA on both platforms. A total of 204 loci were successfully sequenced. Three loci failed completely in the PCR amplification, and three additional loci displayed frequent locus drop-outs due to low read depth or high heterozygote imbalance. Overall, the read depth across the loci was more well-balanced with the MiSeq, while the heterozygote balance was less variable with the Ion S5. Noise levels were low on both platforms (median< 0.2 %). Two simple criteria for genotyping were applied: A minimum threshold of 45 reads and an acceptable heterozygote balance range of 0.3-3.0. Complete concordance between platforms was observed except for three genotypes in one of the poorly performing loci, rs1470637. This locus had relatively low read depths on both platforms, skewed heterozygote balance, and frequent locus drop-outs. There was also full genotype concordance between the results from the GeneRead assay and known profiles generated with the QIAseq and Ion AmpliSeq assays. The few discordant results were either due to locus drop-outs in the poorly performing loci or allele drop-outs in the QIAseq assay. Profiles with a minimum of 179 SNPs were obtained from four challenging case work samples (blood swabs, bone, or blood from a corpse). Overall, the GeneRead DNASeq assay showed considerable potential and could provide a reliable method for SNP genotyping in cases involving identification of individuals, prediction of phenotypic traits, and ancestry inference.

Characterization of Fecal Microbiota with Clinical Specimen Using Long-Read and Short-Read Sequencing Platform

Accurate and rapid identification of microbiotic communities using 16S ribosomal (r)RNA sequencing is a critical task for expanding medical and clinical applications. Next-generation sequencing (NGS) is widely considered a practical approach for direct application to communities without the need for in vitro culturing. In this report, a comparative evaluation of short-read (Illumina) and long-read (Oxford Nanopore Technologies (ONT)) platforms toward 16S rRNA sequencing with the same batch of total genomic DNA extracted from fecal samples is presented. Different 16S gene regions were amplified, bar-coded, and sequenced using the Illumina MiSeq and ONT MinION sequencers and corresponding kits. Mapping of the sequenced amplicon using MinION to the entire 16S rRNA gene was analyzed with the cloud-based EPI2ME algorithm. V3–V4 reads generated using MiSeq were aligned by applying the CLC genomics workbench. More than 90% of sequenced reads generated using distinct sequencers were accurately classified at the genus or species level. The misclassification of sequenced reads at the species level between the two approaches was less substantial as expected. Taken together, the comparative results demonstrate that MinION sequencing platform coupled with the corresponding algorithm could function as a practicable strategy in classifying bacterial community to the species level.

Porcine reproductive and respiratory syndrome virus whole-genome sequencing efficacy with field clinical samples using a poly(A)-tail viral genome purification method

The genomic surveillance of porcine reproductive and respiratory syndrome virus (PRRSV) is based on sequencing of the ORF5 gene of the virus, which covers only 4% of the entire viral genome. It is expected that PRRSV whole-genome sequencing (WGS) will improve PRRSV genomic data and allow better understanding of clinical discrepancies observed in the field when using ORF5 sequencing. Our main objective was to implement an efficient method for WGS of PRRSV from clinical samples. The viral genome was purified using a poly(A)-tail viral genome purification method and sequenced using Illumina technology. We tested 149 PRRSV-positive samples: 80 sera, 33 lungs, 33 pools of tissues, 2 oral fluids, and 1 processing fluid (i.e., castration liquid). Overall, WGS of 67.1% of PRRSV-positive cases was successful. The viral load, in particular for tissues, had a major impact on the PRRSV WGS success rate. Serum was the most efficient type of sample to conduct PRRSV WGS poly(A)-tail assays, with a success rate of 76.3%, and this result can be explained by improved sequencing reads dispersion matching throughout the entire viral genome. WGS was unsuccessful for all pools of tissue and lung samples with Cq values > 26.5, whereas it could still be successful with sera at Cq ≤ 34.1. Evaluation of results of highly qualified personnel confirmed that laboratory skills could affect PRRSV WGS efficiency. Oral fluid samples seem very promising and merit further investigation because, with only 2 samples of low viral load (Cq = 28.8, 32.8), PRRSV WGS was successful.