Provides immunoglobulin (Ig) sequence lineage reconstruction, diversity profiling, and amino acid property analysis.
Alakazam is part of the Immcantation analysis framework for Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) and provides a set of tools to investigate lymphocyte receptor clonal lineages, diversity, gene usage, and other repertoire level properties, with a focus on high-throughput immunoglobulin (Ig) sequencing.
Alakazam serves five main purposes:
General:
seqDist() wherein distance was not properly calculated in
some sequences containing gap characters.getAAMatrix() return matrix.General:
readChangeoDb() to wrap data.table::fread() instead of
utils::read.table() if the input file is not compressed.testSeqEqual(), getSeqDistance() and getSeqMatrix() to C++ to
improve performance of collapseDuplicates() and other dependent functions.testSeqEqual(), getSeqDistance() and getSeqMatrix() to
seqEqual(), seqDist() and pairwiseDist(), respectively.pairwiseEqual() which creates a logical sequence distance matrix;
TRUE if sequences are identical, FALSE if not, excluding Ns and gaps.X in
translateDNA().collapseDuplicates() wherein the input data type sanity check
would cause the vignette to fail to build under R 3.3.ExampleDb.gz file with a larger, more clonal, ExampleDb
data object.ExampleTrees with a larger set of trees.multiggplot() to gridPlot().Amino Acid Analysis:
normalize=FALSE for charge calculations to be more consistent
with previously published repertoire sequencing results.Diversity Analysis:
progress argument to rarefyDiversity() and testDiversity() to
enable the (previously default) progress bar.estimateAbundance() were the function would fail if there
was only a single input sequence per group.data and summary slots of DiversityTest to
uppercase for consistency with other tools.plot to plotDiversityCurve for DiversityCurve
objects.Gene Usage:
sortGenes() function to sort V(D)J genes by name or locus position.clone argument to countGenes() to allow restriction of gene
abundance to one gene per clone.Topology Analysis:
General:
base::nchar().General:
Amino Acid Analysis:
aliphatic() function were not being
passed through the ellipsis argument of aminoAcidProperties().aminoAcidProperties().AA_TRANS to ABBREV_AA.Diversity:
rarefyDiversity()
output.Lineage:
ExampleTrees data with example output from buildPhylipLineage().General:
getDNADistMatrix() and getAADistMatrix() to getDNAMatrix and
getAAMatrix(), respectively.getSeqMatrix() which calculates a pairwise distance matrix for a set
of sequences.multiggplot() function for performing multiple panel plots.Amino Acid Analysis:
gravy(), bulk(), aliphatic(), polar(),
charge(), countPatterns() and aminoAcidProperties().Annotation:
getSegment(), getAllele(), getGene() and getFamily(). May be
disabled by providing the argument strip_d=FALSE.countGenes() to tabulate V(D)J allele, gene and family usage.Diversity:
countClones(), estimateAbundance() and plotAbundance().resampleDiversity() to rarefyDiversity() and changed many of
the internals. Bootstrapping is now performed on an inferred complete
relative abundance distribution.rarefyDiversity() and testDiversity().rarefyDiversity()
and testDiversity() are now calculated using the mean and standard
deviation of the bootstrap realizations, rather than the median and
upper/lower quantiles.plotDiversityCurve().Initial public release.
General:
citation("alakazam") command.Lineage:
buildPhylipLineage().Lineage:
buildPhylipLineage() would hang on R 3.2 due to R change
request PR#15508.Prerelease for review.