Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “amplicon sequencing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

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

ABSTRACT Targeted amplicon sequencing is widely used in microbial ecology studies. However, sequencing artifacts and amplification biases are of great concern. To identify sources of these artifacts, a systematic analysis was performed using mock communities comprised of 16S rRNA genes from 33 bacterial strains. Our results indicated that while sequencing errors were generally isolated to low-abundance operational taxonomic units, chimeric sequences were a major source of artifacts. Singleton and doubleton sequences were primarily chimeras. Formation of chimeric sequences was significantly correlated with the GC content of the targeted sequences. Low-GC-content mock community members exhibited lower rates of chimeric sequence formation. GC content also had a large impact on sequence recovery. The quantitative capacity was notably limited, with substantial recovery variations and weak correlation between anticipated and observed strain abundances. The mock community strains with higher GC content had higher recovery rates than strains with lower GC content. Amplification bias was also observed due to the differences in primer affinity. A two-step PCR strategy reduced the number of chimeric sequences by half. In addition, comparative analyses based on the mock communities showed that several widely used sequence processing pipelines/methods, including DADA2, Deblur, UCLUST, UNOISE, and UPARSE, had different advantages and disadvantages in artifact removal and rare species detection. These results are important for improving sequencing quality and reliability and developing new algorithms to process targeted amplicon sequences. 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.

59 BASIC BIOLOGICAL SCIENCES↗

Multiplex PCR Targeted Amplicon Sequencing (MTA-Seq): Simple, Flexible, and Versatile SNP Genotyping by Highly Multiplexed PCR Amplicon Sequencing

Next-generation sequencing (NGS) technologies have enabled genome re-sequencing for exploring genome-wide polymorphisms among individuals, as well as targeted re-sequencing for the rapid and simultaneous detection of polymorphisms in genes associated with various biological functions. Therefore, a simple and robust method for targeted re-sequencing should facilitate genotyping in a wide range of biological fields. In this study, we developed a simple, custom, targeted re-sequencing method, designated “multiplex PCR targeted amplicon sequencing (MTA-seq),” and applied it to the genotyping of the model grass Brachypodium distachyon. To assess the practical usability of MTA-seq, we applied it to the genotyping of genome-wide single-nucleotide polymorphisms (SNPs) identified in natural accessions (Bd1-1, Bd3-1, Bd21-3, Bd30-1, Koz-1, Koz-3, and Koz-4) by comparing the re-sequencing data with that of reference accession Bd21. Examination of SNP-genotyping accuracy in 443 amplicons from eight parental accessions and an F1 progeny derived by crossing of Bd21 and Bd3-1 revealed that ~95% of the SNPs were correctly called. The assessment suggested that the method provided an efficient framework for accurate and robust SNP genotyping. The method described here enables easy design of custom target SNP-marker panels in various organisms, facilitating a wide range of high-throughput genetic applications, such as genetic mapping, population analysis, molecular breeding, and genomic diagnostics.

59 BASIC BIOLOGICAL SCIENCES↗

Amplicon Sequencing Reveals Microbiological Signatures in Spent Nuclear Fuel Storage Basins

Water quality is an important determinant for the structural integrity of alloy cladded fuels and assemblies during long-term wet storage. Detailed characterization of a water filled storage basin for spent nuclear reactor fuel was performed following the formation and proliferation of an amorphous white flocculent. White precipitant was sampled throughout the storage basin for chemical and spectroscopic characterization, and eDNA was extracted for pyrosequencing of bacterial rRNA gene diversity. Accordingly, spectroscopic analyses indicated the precipitant to be primarily amorphous to crystalline aluminum (oxy) hydroxides with minor associated elemental components including Fe, Si, Ti, and U. High levels of dissolved carbon were co-localized with the precipitant relative to bulk water. Bacterial densities were highly variable between sampling locations and with depth; cell numbers (log scale) ranged from 5.6 to 4.89 cells / mL. Bacterial diversity that was physically associated with the aluminum (oxy) hydroxide complexes exceeded an estimated 4,000 OTUs / amplicon library (3% cutoff) and the greatest percent majority of sequences were aligned to the families Burkholderiales (23%), Nitrospiraceae (23%), Hyphomicrobiaceae (17%), and Comamonadaceae (6%). We surmise that episodic changes in the physical and chemical properties of the basin contribute to the polymerization of aluminum (oxy) hydroxides, which in turn can chemisorb nutrients, carbon ligands and bacterial cells from the surrounding bulk aqueous phase. As such, these precipitants should establish favorable microhabitats for bacterial colonization and growth. Comparative analyses of 16S rRNA gene amplicon libraries across diverse environmental landscapes were performed and microbiological signatures unique to the spent nuclear fuel storage basin environment were revealed. These insights could spur the development of tractable bioindicators that are specific of and diagnostic for water quality at discrete locations and finer scales of resolution, marking an important contribution for improved water quality and management of spent nuclear fuel storage facilities.

aluminum↗

Diazotroph Community Characterization via a High-Throughput nifH Amplicon Sequencing and Analysis Pipeline

The dinitrogenase reductase gene ( nifH ) is the most widely established molecular marker for the study of nitrogen-fixing prokaryotes in nature. A large number of PCR primer sets have been developed for nifH amplification, and the effective deployment of these approaches should be guided by a rapid, easy-to-use analysis protocol. Bioinformatic analysis of marker gene sequences also requires considerable expertise. In this study, we advance the state of the art for nifH analysis by evaluating nifH primer set performance, developing an improved amplicon sequencing workflow, and implementing a user-friendly bioinformatics pipeline. Furthermore, the developed amplicon sequencing workflow is a three-stage PCR-based approach that uses established technologies for incorporating sample-specific barcode sequences and sequencing adapters. Based on our primer evaluation, we recommend the Ando primer set be used with a modified annealing temperature of 58°C, as this approach captured the largest diversity of nifH templates, including paralog cluster IV/V sequences. To improve nifH sequence analysis, we developed a computational pipeline which infers taxonomy and optionally filters out paralog sequences. In addition, we employed an empirical model to derive optimal operational taxonomic unit (OTU) cutoffs for the nifH gene at the species, genus, and family levels. A comprehensive workflow script named TaxADivA (TAXonomy Assignment and DIVersity Assessment) is provided to ease processing and analysis of nifH amplicons. Our approach is then validated through characterization of diazotroph communities across environmental gradients in beach sands impacted by the Deepwater Horizon oil spill in the Gulf of Mexico, in a peat moss-dominated wetland, and in various plant compartments of a sugarcane field. IMPORTANCE: Nitrogen availability often limits ecosystem productivity, and nitrogen fixation, exclusive to prokaryotes, comprises a major source of nitrogen input that sustains food webs. The nifH gene, which codes for the iron protein of the nitrogenase enzyme, is the most widely established molecular marker for the study of nitrogen-fixing microorganisms (diazotrophs) in nature. In this study, a flexible sequencing/analysis pipeline, named TaxADivA, was developed for nifH amplicons produced by Illumina paired-end sequencing, and it enables an inference of taxonomy, performs clustering, and produces output in formats that may be used by programs that facilitate data exploration and analysis. Diazotroph diversity and community composition are linked to ecosystem functioning, and our results advance the phylogenetic characterization of diazotroph communities by providing empirically derived nifH similarity cutoffs for species, genus, and family levels. Thus, the utility of our pipeline is validated for diazotroph communities in a variety of ecosystems, including contaminated beach sands, peatland ecosystems, living plant tissues, and rhizosphere soil.

59 BASIC BIOLOGICAL SCIENCES↗

Intensity of sample processing methods impacts wastewater SARS-CoV-2 whole genome amplicon sequencing outcomes

Wastewater SARS-CoV-2 surveillance has been deployed since the beginning of the COVID-19 pandemic to monitor the dynamics in virus burden in local communities. Genomic surveillance of SARS-CoV-2 in wastewater, particularly efforts aimed at whole genome sequencing for variant tracking and identification, are still challenging due to low target concentration, complex microbial and chemical background, and lack of robust nucleic acid recovery experimental procedures. The intrinsic sample limitations are inherent to wastewater and are thus unavoidable. Here, we use a statistical approach that couples correlation analyses to a random forest-based machine learning algorithm to evaluate potentially important factors associated with wastewater SARS-CoV-2 whole genome amplicon sequencing outcomes, with a specific focus on the breadth of genome coverage. We collected 182 composite and grab wastewater samples from the Chicago area between November 2020 to October 2021. Samples were processed using a mixture of processing methods reflecting different homogenization intensities (HA + Zymo beads, HA + glass beads, and Nanotrap), and were sequenced using one of the two library preparation kits (the Illumina COVIDseq kit and the QIAseq DIRECT kit). Technical factors evaluated using statistical and machine learning approaches include sample types, certain sample intrinsic features, and processing and sequencing methods. The results suggested that sample processing methods could be a predominant factor affecting sequencing outcomes, and library preparation kits was considered a minor factor. Finally, a synthetic SARS-CoV-2 RNA spike-in experiment was performed to validate the impact from processing methods and suggested that the intensity of the processing methods could lead to different RNA fragmentation

60 APPLIED LIFE SCIENCES↗

16S rRNA Amplicon sequence data and community assembly process file

Soils along a salinity gradient from a coastal forested floodplain (Beaver Creek in Western Washington) were analyzed for their microbial community composition and community assembly processes. These files contain processed 16S rRNA amplicon sequences into OTU table, the corresponding tree file, metadata, community assembly process analysis and results, microbial phyla abundances, and codes evaluating BNTI indices and community assembly process estimates.

54 ENVIRONMENTAL SCIENCES↗

Metatranscriptomics and Amplicon Sequencing Reveal Mutualisms in Seagrass Microbiomes

Terrestrial plants benefit from many well-understood mutualistic relationships with root and leaf-associated microbiomes, but relatively little is known about these relationships for seagrass and other aquatic plants. We used 16S rRNA gene amplicon sequencing and metatranscriptomics to assess potential mutualisms between microorganisms and the seagrasses Zostera marina and Zostera japonica collected from mixed beds in Netarts Bay, OR, United States. The phylogenetic composition of leaf-, root-, and water column-associated bacterial communities were strikingly different, but these communities were not significantly different between plant species. Many taxa present on leaves were related to organisms capable of consuming the common plant metabolic waste product methanol, and of producing agarases, which can limit the growth of epiphytic algae. Taxa present on roots were related to organisms capable of oxidizing toxic sulfur compounds and of fixing nitrogen. Metatranscriptomic sequencing identified expression of genes involved in all of these microbial metabolic processes at levels greater than typical water column bacterioplankton, and also identified expression of genes involved in denitrification and in bacterial synthesis of the plant growth hormone indole-3-acetate. These results provide the first evidence using metatranscriptomics that seagrass microbiomes carry out a broad range of functions that may benefit their hosts, and imply that microbe–plant mutualisms support the health and growth of aquatic plants.

59 BASIC BIOLOGICAL SCIENCES↗

Using high-throughput amplicon sequencing to determine diet of generalist lady beetles in agricultural landscapes

Determining feeding relationships is central to understanding biological control potential in the field. However, methods to differentiate actual (or realized) feeding relationships from potential feeding relationships is lacking especially for small, generalist predators such as lady beetles. In this study, we used DNA metabarcoding approaches to characterize actual feeding relationships of lady beetles (Coccinellidae) in the field and validated our methods with a lab study. We first asked whether high-throughput amplicon sequencing (HTS) can characterize diets of lady beetles ranging from monotypic diets to diverse diet mixtures in the lab. We then examined whether diet composition and breadth of lady beetles collected from different habitat types in southern WI varied between monocultures of soybean and corn, diverse tallgrass prairie, and urban habitats. Lastly, we asked whether different body or tissue types (partial-body versus whole-body specimens) would change the likelihood of prey detection for both studies. In our controlled lab study, we found that HTS can accurately assess diet composition and diet breadth for lady beetle populations, but at the individual level, HTS has limitation for individuals feeding on more than three species of prey at any given time. In our field study, we documented lower prey richness and diet breadth in corn compared to soybean and grassland, and greater prey DNA in regurgitants than in whole-body specimens. Finally, we found that reduced diet diversity was associated with an increased prevalence of intraguild predation, but that habitat and prey diversity are not necessarily correlated. The prey detection rates (proportion of specimens with prey DNA) in our study were comparable to other studies (55% in the lab study, 33% in field study) and varied with diet composition and habitat type. Furthermore, the relatively low detection rates suggest that many consumer individuals would need to be assayed to fully assess diet diversity, especially in diverse systems.

59 BASIC BIOLOGICAL SCIENCES↗

Amplicon Sequencing Assessment to Measure Microbial Community Response from Heavy Metal Contaminated Soils in Savannah River site, Tims Branch watershed

Long-term presence of heavy metal contamination in soils is correlated with changes in microbial community structures and can lead to the reduction of indigenous species. Tolerance to soil heavy metal toxicity varies among different microbial communities. It remains unclear whether long term effects caused by metal perturbations are associated with soil microbial community dysregulation in Tims Branch watershed. Research Questions: Does heavy metal contamination in Tims Branch soils decrease microbial species diversity significantly compared to non-contaminated soils? Does the species richness decrease in contaminated samples? between samples? Is diversity in low, mid and high contamination soils significantly different compared to control samples? Hypothesis: The relative abundance and diversity of species are significantly altered in soils which are contaminated with heavy metals in the Tims Branch watershed. Objectives: Use amplicon sequencing technology along with bioinformatics to compare bacterial communities between four sites in Tims Branch watershed. Evaluate the percent relative abundance and diversity across the site. Conclusions: All bacterial communities exposed to different heavy metal concentrations were dominated by four major groups (Proteobacteria, Acidobacteria, Chloroflexi and Verrucomicrobia). The dominant phyla Proteobacteria, Acidobacteria and Chloroflexi accounted for 61 % of the relative abundance. Soil bacterial α-diversity, expressed as observed species richness metric and Shannon's diversity index, was the highest in sample location four (high contamination). The lowest observed was at the control location. Species richness increased when concentrations of heavy metals were localized, inferring that the metals can alter microbial community structure. PERMANOVA analysis between contaminated groups tested significant using the Adonis (F (4,29) =1.9712; p < 0.001). There is a 20 percent variation explained between groups (R2 = 0.214). This analysis proves to be beneficial in detecting microbial communities altered by contaminated soils. Moreover, certain concentrations of heavy metals can alter microbial community structures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Characterization of prokaryotic communities in Puerto Rican caves using 16S rDNA amplicon sequencing

The cave ecosystems host microbial communities adapted to extreme environments. This study utilized 16S rDNA to investigate the prokaryotic diversity across seven caves in Puerto Rico’s northern limestone karst belt. Microbial profiling revealed distinct subterranean communities, enhancing our understanding of cave microbiology and potential applications in environmental conservation and microbial research.

16S↗

evSeq: Cost-Effective Amplicon Sequencing of Every Variant in a Protein Library

Widespread availability of protein sequence-fitness data would revolutionize both our biochemical understanding of proteins and our ability to engineer them. Unfortunately, even though thousands of protein variants are generated and evaluated for fitness during a typical protein engineering campaign, most are never sequenced, leaving a wealth of potential sequence-fitness information untapped. Primarily, this is because sequencing is unnecessary for many protein engineering strategies; the added cost and effort of sequencing is thus unjustified. It also results from the fact that, even though many lower cost sequencing strategies have been developed, they often require at least some sequencing or computational resources, both of which can be barriers to access. In this work, we present every variant sequencing (evSeq), a method and collection of tools/standardized components for sequencing a variable region within every variant gene produced during a protein engineering campaign at a cost of cents per variant. evSeq was designed to democratize low-cost sequencing for protein engineers and, indeed, anyone interested in engineering biological systems. Execution of its wet-lab component is simple, requires no sequencing experience to perform, relies only on resources and services typically available to biology labs, and slots neatly into existing protein engineering workflows. Analysis of evSeq data is likewise made simple by its accompanying software (found at github.com/fhalab/evSeq, documentation at fhalab.github.io/evSeq), which can be run on a personal laptop and was designed to be accessible to users with no computational experience. Here, low-cost and easy to use, evSeq makes collection of extensive protein variant sequence-fitness data practical.

59 BASIC BIOLOGICAL SCIENCES↗

16S and ITS Amplicon Sequencing Fastq files and metadata from PARCHED Panama Tropical Forest soils, 2019-2020,

Model projections predict tropical forests will experience longer periods of drought and more intense precipitation cycles under a changing climate. Such transitions have implications for structure-function relationships within microbial communities. We examine how chronic drying might reshape prokaryotic and fungal communities across four lowland forests in Panama with a wide variation in mean annual precipitation and soil fertility. Four sites were established across a 1000 mm span in mean annual precipitation (2335 to 3300 mm). We expected microbial communities at sites with lower MAP to be less sensitive to chronic drying than sites with higher MAP; while fungal communities to be more resistant to disturbance than prokaryotes. At each location, partial throughfall exclusion structures were established over 10 x 10 m plots to reduce direct precipitation input. Raw demultiplexed sequences (bacteria, archaea, fungal) from soil samples taken from PARCHED Panama Tropical Forest throughfall exclusion experiments. Files that contain 16S are sequences from prokaryotes, ITS indicates sequences from fungi. Compressed fastq files are contained in Field_PARCHED_ITS_fastq_2020.zip, Field_PARCHED_ITS_fastq_2019.zip, Field_PARCHED_16S_fastq_2020.zip, Field_PARCHED_16S_fastq_2019.zip.There are four sites across the isthmus of Panama, each with 4 control and 4 exclusion plots. Throughfall exclusion shelters were built to intercept 50% of throughfall precipitation that hits the soil. Samples were taken on May 2019 and Jan 2020 from 0-10 cm and 10-20 cm depth approximately 9 and 18 months after shelter installation. Metadata and sample IDs for fastq files for 2019 sampling are within Field_PARCHED_Metadata_2019_ITS.csv and Field_PARCHED_Metadata_2019_16S.csv. The metadata for both 16S and ITS fastq files for January 2020 is included in Field_PARCHED_Metadata_2020.csv. No data processing or QA/QC was done on the raw data. Data processing example provided in R notebook file.

54 ENVIRONMENTAL SCIENCES↗

Carbon dioxide (CO2) flux, amplicon sequencing and liquid chromatography–mass spectrometry (LC-MS) data for pyrogenic organic matter (PYOM) amended soil incubations in lab, Madison, WI, 2021-22

This dataset comprises soil-pyrogenic organic matter (PyOM) incubation experiments conducted in a laboratory setting. The focus was on measuring carbon dioxide (CO2) flux emitted in the headspace of incubation jars over a one-month period. The aim was to estimate carbon mineralization from two different carbon fractions in PyOM, as well as carbon mineralization from soil organic carbon (SOC). Bacterial and fungal community profiles were tracked at various time points during the incubation to observe changes in the overall community and specific responders to PyOM. LC-MS analysis of soil samples was also performed to monitor changes in the soil's chemical composition.The provided dataset includes processed CO2 flux measurements obtained from the Picarro-multiplexer measurement setup for each sample incubation jar. This data is available in the "CO2Flux.csv" file. Additionally, the dataset includes CO2 flux data partitioned to estimate mean carbon mineralization from the two PyOM carbon fractions, bulk PyOM, and SOC. The cumulative carbon mineralization data is provided in "cml_respired.csv," while the rate of carbon mineralization data is in "rate_respired.csv." The file "cml_respired_lastcyc.csv" contains cumulative carbon mineralization data for the PyOM carbon fractions, bulk PyOM, and SOC throughout the entire incubation period.For the LC-MS analysis, the raw data can be found in the "20220721_NZ_PeakHeights_FinalAnalysis.csv" file, while sample metadata is provided in "20220721_NZpilot_Metadata.csv". Processed LC-MS data is available in the file "20220721_NZ_FINAL_msdat_sub3.csv". Lastly, the dataset includes the relative abundance data for bacterial genera that exhibited a significant positive response to the addition of PyOM produced at 350 degrees Celsius, which can be found in "resp_genera_350.csv".

54 ENVIRONMENTAL SCIENCES↗