chimera_removal_vs: Search for a list of sequence in an object to remove chimera...

In MiscMetabar: Miscellaneous Functions for Metabarcoding Analysis

Search for a list of sequence in an object to remove chimera taxa using vsearch

Description

Use the VSEARCH software.

Usage

chimera_removal_vs(object, type = "Discard_only_chim", clean_pq = FALSE, ...)

Arguments

object	(required) A phyloseq-class object or one of dada, derep, data.frame or list coercible to sequences table using the function dada2::makeSequenceTable()
type	(default "Discard_only_chim"). The type define the type of filtering. "Discard_only_chim" will only discard taxa classify as chimera by vsearch "Select_only_non_chim_seqlen_filtered" will only select taxa classify as non-chimera by vsearch(after filtering taxa based on their sequence length by the parameter min_seq_length from the chimera_detection_vs() function) "Select_only_chim" will only select taxa classify as chimera by vsearch (after filtering taxa based on their sequence length by the parameter min_seq_length from the chimera_detection_vs() function)
clean_pq	(logical; default FALSE) If TRUE, return the phyloseq object after cleaning using the default parameter of clean_pq() function.
...	Additional arguments passed on to chimera_detection_vs() function

Details

This function is mainly a wrapper of the work of others. Please cite vsearch.

Value

• I/ a sequences tables if object is of class dada, derep, data.frame or list.

• II/ a phyloseq object without (or with if type = 'Select_only_chim') chimeric taxa

Author(s)

Adrien Taudière

See Also

chimera_detection_vs(), dada2::removeBimeraDenovo()

Examples

data_fungi_nochim <- chimera_removal_vs(data_fungi)

## Not run: 
# Adding a chimeric sequence for the example
data_fungi_with_chim <- data_fungi
data_fungi_with_chim@refseq["ASV1710"] <- Biostrings::xscat(
  Biostrings::subseq(data_fungi_with_chim@refseq[1], start = 1, end = 150),
  Biostrings::subseq(data_fungi_with_chim@refseq[4], start = 151, end = 300)
)
data_fungi_nochim <- chimera_removal_vs(data_fungi)

# Higher value of abskew parameter is less stringent
data_fungi_nochim_16 <- chimera_removal_vs(data_fungi,
  abskew = 16, min_seq_length = 10
)

# Potential Chimeric ASVs detected by vsearch
chim_asv <- taxa_names(data_fungi_with_chim)[!taxa_names(data_fungi_with_chim)
%in% taxa_names(data_fungi_nochim)]
"ASV1710" %in% chim_asv
track_wkflow(list(data_fungi_with_chim, data_fungi_nochim))

data_fungi_nochim2 <-
  chimera_removal_vs(data_fungi, type = "Select_only_non_chim_seqlen_filtered")
data_fungi_chimera <-
  chimera_removal_vs(data_fungi, type = "Select_only_chim")

## End(Not run)

MiscMetabar documentation built on June 8, 2026, 5:07 p.m.