The somatic recombination of V, D, and J gene segments in B-cells introduces a great deal of diversity, and divergence from reference segments. Many recent studies of antibodies focus on the population of antibody transcripts that show which V, D, and J gene segments have been favored for a particular antigen, a repertoire. To properly describe the antibody repertoire, each antibody must be labeled by its constituting V, D, and J gene segment, a task made difficult by somatic recombination and hypermutation events. While previous approaches to repertoire analysis were based on sequential alignments, we describe a new de Bruijn graph–based algorithm to perform VDJ labeling and benchmark its performance.
Get full access to this article
View all access options for this article.
References
1.
AnglyF.E., WillnerD., RohwerF., et al.2012. Grinder: A versatile amplicon and shotgun sequence simulator. Nucleic Acids Res. 40, e94–e94.
2.
ArnaoutR., LeeW., CahillP., et al.2011. High-resolution description of antibody heavy-chain repertoires in humans. PLoS One, 6, e22365.
3.
BasuM., HegdeM.V., and ModakM.J.1983. Synthesis of compositionally unique DNA by terminal deoxynucleotidyl transferase. Biochem Biophys Res Commun, 111, 1105–1112.
4.
BrochetX., LefrancM., and GiudicelliV.2008. IMGT/V-QUEST: The highly customized and integrated system for IG and TR standardized VJ and VDJ sequence analysis. Nucleic Acids Res. 36, W503–W508.
5.
ChenW., PrabakaranP., ZhuZ., et al.2012. Identification of cross-reactive IgG antibodies from an acute HIV-1-infected patient using phage display and high-throughput sequencing technologies. Exp. Mol. Pathol., 93.
6.
ClarkL.A., GanesanS., PappS., and van VlijmenH.W.2006. Trends in antibody sequence changes during the somatic hypermutation process. J. Immunol., 177, 333–340.
7.
CompeauP.E., PevznerP.A., and TeslerG.2011. How to apply de Bruijn graphs to genome assembly. Nat. Biotechnol., 29, 987–991.
8.
DesiderioS.V., YancopoulosG.D., PaskindM., et al.1984. Insertion of N regions into heavy-chain genes is correlated with expression of terminal deoxytransferase in B cells. Nature, 311, 752–755.
9.
DörnerT., FosterS.J., FarnerN.L., and LipskyP.E.1998. Somatic hypermutation of human immunoglobulin heavy chain genes: Targeting of RGYW motifs on both DNA strands. Eur. J. Immunol., 28, 3384–3396.
10.
FowlkesE.B., and MallowsC.L.1983. A method for comparing two hierarchical clusterings. J. Am. Stat. Assoc., 78, 553–569.
11.
GaétaB.A., MalmingH.R., JacksonK.J., et al.2007. iHMMune-align: Hidden Markov model-based alignment and identification of germline genes in rearranged immunoglobulin gene sequences. Bioinformatics, 23, 1580–1587.
12.
HalemanoK., GuoK., HeilmanK.J., et al.2014. Immunoglobulin somatic hypermutation by apobec3/rfv3 during retroviral infection. Proc. Natl. Acad. Sci. USA, 111, 7759–7764.
13.
IqbalZ., CaccamoM., TurnerI., et al.2012. De novo assembly and genotyping of variants using colored de Bruijn graphs. Nat. Genet., 44, 226–232.
14.
JaccardP.1908. Nouvelles recherches sur la distribution orale. Bull. Soc. Vaudoise Sci. Nat., 44, 223–370.
15.
JacksonK.J., BoydS., GaétaB.A., and CollinsA.M.2010. Benchmarking the performance of human antibody gene alignment utilities using a 454 sequence dataset. Bioinformatics, 26, 3129–3130.
16.
JacksonK.J., GaetaB., SewellW., and CollinsA.M.2004. Exonuclease activity and P nucleotide addition in the generation of the expressed immunoglobulin repertoire. BMC Immunol. 5, 19.
17.
JiangN., HeJ., WeinsteinJ.A., et al.2013. Lineage structure of the human antibody repertoire in response to inuenza vaccination. Sci. Transl. Med. 5, 171ra19.
18.
JiangN., WeinsteinJ., PenlandL., et al.2011. Determinism and stochasticity during maturation of the zebrafish antibody repertoire. Proc. Natl. Acad. Sci. USA, 108, 5348–5353.
19.
Ohm-LaursenL., NielsenM., LarsenS.R., and BaringtonT.2006. No evidence for the use of DIR, D–D fusions, chromosome 15 open reading frames or VH replacement in the peripheral repertoire was found on application of an improved algorithm, JointML, to 6329 human immunoglobulin H rearrangements. Immunology, 119, 265–277.
20.
PevznerP.A., TangH., and TeslerG.2004. De novo repeat classification and fragment assembly. Genome Res. 14, 1786–1796.
21.
RobinsonJ., HalliwellJ.A., McWilliamH., et al.2013. The IMGT/HLA database. Nucleic Acids Res. 41, D1222–D1227.
22.
RogozinI.B., and KolchanovN.A.1992. Somatic hypermutagenesis in immunoglobulin genes: II. Infuence of neighbouring base sequences on mutagenesis. Biochim. Biophys. Acta, 1171, 11–18.
23.
SafonovaY., BonissoneS., KurpilyanskyE., et al.2015. IgRepertoireConstructor: A novel algorithm for antibody repertoire construction and immunoproteogenomics analysis. Bioinformatics, 31, i53–i61.
24.
Souto-CarneiroM.M., LongoN.S., RussD.E., et al.2004. Characterization of the human Ig heavy chain antigen binding complementarity determining region 3 using a newly developed software algorithm, JOINSOLVER. J. Immunol., 172, 6790–6802.
25.
VolpeJ.M., CowellL.G., and KeplerT.B.2006. SoDA: Implementation of a 3D alignment algorithm for inference of antigen receptor recombinations. Bioinformatics, 22, 438–444.
26.
WangX., WuD., ZhengS., et al.2008. Ab-origin: An enhanced tool to identify the sourcing gene segments in germline for rearranged antibodies. BMC Bioinform. 9, S20.
27.
WeinsteinJ., JiangN., WhiteR., et al.2009. High-throughput sequencing of the zebrafish antibody repertoire. Science, 324, 807–810.
28.
WineY., BoutzD.R., LavinderJ.J., et al.2013. Molecular deconvolution of the monoclonal antibodies that comprise the polyclonal serum response. Proc. Natl. Acad. Sci. USA, 110, 2993–2998.
29.
YeJ., MaN., MaddenT.L., and OstellJ.M.2013. IgBLAST: An immunoglobulin variable domain sequence analysis tool. Nucleic Acids Res, 41, W34–W40.
30.
ZerbinoD.R., and BirneyE.2008. Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 18, 821–829.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
