R Tools for Inferring New Immunoglobulin Alleles from Rep-Seq Data

Infers the V genotype of an individual from immunoglobulin (Ig) repertoire-sequencing (Rep-Seq) data, including detection of any novel alleles. This information is then used to correct existing V allele calls from among the sample sequences.


High-throughput sequencing of B cell immunoglobulin receptors is providing unprecedented insight into adaptive immunity. A key step in analyzing these data involves assignment of the germline V, D and J gene segment alleles that comprise each immunoglobulin sequence by matching them against a database of known V(D)J alleles. However, this process will fail for sequences that utilize previously undetected alleles, whose frequency in the population is unclear.

TIgGER is a computational method that significantly improves V(D)J allele assignments by first determining the complete set of gene segments carried by an individual (including novel alleles) from V(D)J-rearrange sequences. TIgGER can then infer a subject's genotype from these sequences, and use this genotype to correct the initial V(D)J allele assignments.

The application of TIgGER continues to identify a surprisingly high frequency of novel alleles in humans, highlighting the critical need for this approach. (TIgGER, however, can and has been used with data from other species.)

  • Detecting novel alleles
  • Inferring a subject's genotype
  • Correcting preliminary allele calls
  • A table of sequences from a single individual, with columns containing the following:
    • V(D)J-rearranged nucleotide sequence (in IMGT-gapped format)
    • Preliminary V allele calls
    • Preliminary J allele calls
    • Length of the junction region
  • Germline Ig sequences in IMGT-gapped fasta format (e.g., as those downloaded from IMGT/GENE-DB)

The former can be created through the use of IMGT/HighV-QUEST and Change-O.

For help, questions, or suggestions, please contact Daniel Gadala-Maria or use the issue tracker.

News

  • Bugfix following recent update of alakazam (0.2.5) to import selectively.
  • Removed unneeded dependency on shazam package (not needed as of 0.2.5.999).
  • More updates to work with the latest version of dplyr (0.5.0).
  • Bugfix in findNovelAlleles when allele passed germline_min but not min_seqs.
  • Fixed vignette typo and updated findUnmutatedCalls man page.
  • Updated code to work with the latest version of dplyr (0.5.0).
  • Fixed a bug werein findNovelAlleles() was not running in parallel, even when nproc > 1.
  • Changed default to nproc=1 in findNovelAlleles().
  • Initial CRAN release.

Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.

install.packages("tigger")

0.2.10 by Jason Vander Heiden, 3 months ago


http://tigger.readthedocs.io


Report a bug at https://bitbucket.org/kleinstein/tigger/issues


Browse source code at https://github.com/cran/tigger


Authors: Daniel Gadala-Maria [aut], Jason Vander Heiden [ctb, cre], Steven Kleinstein [aut, cph]


Documentation:   PDF Manual  


CC BY-SA 4.0 license


Imports alakazam, tidyr, dplyr, doParallel, foreach, graphics, grid, iterators, lazyeval, parallel, stats

Depends on ggplot2

Suggests knitr, testthat


See at CRAN