Home - News ( United States | United Kingdom | Italy | Germany ) - Football scores

Showing content from https://grunwaldlab.github.io/demulticoder/ below:

Simultaneous Analysis of Multiplexed Metabarcodes • demulticoder

The demulticoder package is a Cutadapt and DADA2 wrapper package for metabarcoding analyses. The main commands and outputs are intuitive and comprehensive, which helps to account for the complex and iterative nature of metabarcoding analyses.

Here is a brief schematic of the general workflow:

• It automates the use of DADA2 to analyze data derived from multiple metabarcodes.
  
• It reduces the number of manual input steps
  
• Handles analysis of two metabarcodes multiplexed into the same sequencing batch
  
• Analyze different types of metabarcodes simultaneously
  
• Reproducible workflows for oomycetes
• Supported metabarcodes: 16S rDNA, ITS1, rps10, and up to two additional metabarcodes

Dependencies:
First install Cutadapt program following the instructions here: https://cutadapt.readthedocs.io/en/stable/installation.html
Let’s locate where the cutadapt executable is. You must do this from a Terminal window:

#If you installed with pip or pipx, or homebrew, run this command from a Terminal window
which cutadapt
cutadapt --version

If you followed the cutadapt installation instructions to create a conda environment called cutadapt (change to whatever you named your environment), to install it in, open up a Terminal window and type these commands:

#Run commands from a Terminal window
conda activate cutadapt
which cutadapt
cutadapt --version

Second, make sure the following R packages are installed:

• demulticoder (Available through CRAN)
  
• DADA2 (Latest version is 3.20)
  • To install, follow these instructions: https://www.bioconductor.org/packages/release/bioc/html/dada2.html
    
• phyloseq
  • To install: https://www.bioconductor.org/packages/release/bioc/html/phyloseq.html
    
• metacoder (Available through CRAN)

To install demulticoder from CRAN:

To instead install the development version for the newest features or bug fixes:

Let’s now make sure all packages are loaded (following installation):

1. Set-up input directory and files

To demonstrate how to use the package, we have a small test data set that comes loaded with the package. This data set will be used in the workflow example below.

Already loaded in the test data set directory are the following files:

• PE short read amplification data
  • Files: S1_R1.fastq.gz, S1_R2.fastq.gz, S2_R1.fastq.gz, S2_R1.fastq.gz
    
  • The files must end in either R1.fastq.gz , or R2.fastq.gz and each sample must have both R1 and R2 files.
• metadata.csv
  • New row for each unique sample
  • Samples entered twice if samples contain two pooled metabolites, as in the test data template
• primerinfo_params.csv
  • New row for each unique metabarcode (applicable options: ‘rps10’, ‘its’, ‘r16S’, ‘other1’, ‘other2’) and associated primer sequence
  • Optional Cutadapt and DADA2 parameters
• Taxonomy databases
  • UNITE fungal database (abridged version)
  • completed

See Documentation for how to format databases and input files.

For more details on each step, check out the Getting Started tab on the package website

2. Prepare reads

output<-prepare_reads(
  data_directory = system.file("extdata", package = "demulticoder"), # This allows us to use the test directory located within the package
  output_directory = tempdir(), # OR change to you preferred location on your local computer (Example: "~/demulticoder_test")
  overwrite_existing = TRUE)

3. Cut and trim reads User must install cutadapt on their local machine and append the path to the executable.

cut_trim(
  output,
  cutadapt_path="/usr/bin/cutadapt", # Change to the location on your computer. (Example: "/usr/bin/cutadapt")
  overwrite_existing = TRUE) 

4. Make ASV abundance matrix

5. Assign taxonomy

assign_tax(
  output,
  asv_abund_matrix,
  overwrite_existing = TRUE)

6. Convert ASV matrix to taxmap and phyloseq objects

The package was developed by Martha Sudermann, Zachary Foster, Samantha Dawson, Hung Phan, Jeff Chang, and Niklaus Grünwald

Sudermann MA, Foster ZSL, Dawson SCL, Phan H, Fieland VJ, Martin FN, Chang JH, Grünwald NJ. Demulticoder: An R Package for the Simultaneous Analysis of Multiplexed Metabarcodes. Phytopathology. 2025. 10.1094/PHYTO-02-25-0043-FI. Epub ahead of print.

This project was funded by the 2020 CDFA Specialty Block Grant Program Project 45.

RetroSearch is an open source project built by @garambo | Open a GitHub Issue

Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo

HTML: 3.2 | Encoding: UTF-8 | Version: 0.7.4