# Recovering probabilities for nucleotide trimming processes for T cell receptor TRA and TRG V-J junctions analyzed with IMGT tools

- Kevin Bleakley
^{1, 2, 3}Email author, - Marie-Paule Lefranc
^{4}and - Gérard Biau
^{5}

**9**:408

https://doi.org/10.1186/1471-2105-9-408

© Bleakley et al; licensee BioMed Central Ltd. 2008

**Received: **06 February 2008

**Accepted: **02 October 2008

**Published: **02 October 2008

## Abstract

### Background

Nucleotides are trimmed from the ends of variable (V), diversity (D) and joining (J) genes during immunoglobulin (IG) and T cell receptor (TR) rearrangements in B cells and T cells of the immune system. This trimming is followed by addition of nucleotides at random, forming the N regions (N for nucleotides) of the V-J and V-D-J junctions. These processes are crucial for creating diversity in the immune response since the number of trimmed nucleotides and the number of added nucleotides vary in each B or T cell. IMGT^{®} sequence analysis tools, IMGT/V-QUEST and IMGT/JunctionAnalysis, are able to provide detailed and accurate analysis of the final observed junction nucleotide sequences (tool "output"). However, as trimmed nucleotides can potentially be replaced by identical N region nucleotides during the process, the observed "output" represents a *biased* estimate of the "true trimming process."

### Results

A probabilistic approach based on an analysis of the standardized tool "output" is proposed to infer the probability distribution of the "true trimmming process" and to provide plausible biological hypotheses explaining this process. We collated a benchmark dataset of TR alpha (TRA) and TR gamma (TRG) V-J rearranged sequences and junctions analysed with IMGT/V-QUEST and IMGT/JunctionAnalysis, the nucleotide sequence analysis tools from IMGT^{®}, the international ImMunoGeneTics information system^{®}, http://imgt.cines.fr. The standardized description of the tool output is based on the IMGT-ONTOLOGY axioms and concepts. We propose a simple first-order model that attempts to transform the observed "output" probability distribution into an estimate closer to the "true trimming process" probability distribution. We use this estimate to test the hypothesis that Poisson processes are involved in trimming. This hypothesis was not rejected at standard confidence levels for three of the four trimming processes: TRAV, TRAJ and TRGV.

### Conclusion

By using trimming of rearranged TR genes as a benchmark, we show that a probabilistic approach, applied to IMGT^{®} standardized tool "outputs" opens the way to plausible hypotheses on the events involved in the "true trimming process" and eventually to an exact quantification of trimming itself. With increasing high-throughput of standardized immunogenetics data, similar probabilistic approaches will improve understanding of processes so far only characterized by the "output" of standardized tools.

## Keywords

## Background

The diversity of the chains of immunoglobulins (IG) or antibodies and T cell receptors (TR) depends on several mechanisms [1–10]: first, combinatorial diversity, which is a consequence of the number of variable (V), diversity (D) and joining (J) genes in the IG and TR loci [9, 10], second, exonuclease trimming of V, D and J nucleotides and third, addition at random of nucleotides at the V-J and V-D-J junction (N region diversity).

These processes together create a huge diversity in V-J and V-D-J junctions as exemplified by the rearranged IG and TR sequences from IMGT/LIGM-DB [11]. In addition, rearranged V-J and V-D-J genes from IG (but not those from TR) are specifically submitted to the mechanism of somatic hypermutations [9] (IMGT Education, Tutorials, http://imgt.cines.fr). The number of different antigen receptors (IG and TR) per individual is estimated to be 2 × 10^{12} in humans and the only limiting factor seems to be the number of B cells (for the IG) and T cells (for the TR) which is genetically programmed in a given species.

Trimming by exonuclease occurs at the ends of the 3'V-REGION and 5'J-REGION [12] (IMGT labels from the DESCRIPTION axiom of IMGT-ONTOLOGY are in capital letters [13, 14]) and at both ends of the D-REGION, present in the IG heavy (IGH), TR beta (TRB) and TR delta (TRD) loci [9, 10]. Little is known about the mechanisms that regulate trimming of V, D and J genes during V-J and V-D-J rearrangement. Given the importance of trimming in the creation of the vast diversity of V-J and V-D-J junctions, it is of great interest to better understand this process.

Based on the IMGT-ONTOLOGY axioms and concepts of classification (IMGT gene names) [9, 10, 15, 16], description (IMGT labels) [17, 18] and numerotation (IMGT concepts for numbering, in particular, IMGT unique numbering for V, C and G domains) [19–21], on-line tools have been developed by IMGT^{®}, the international ImMunoGeneTics information system^{®}, http://imgt.cines.fr[22], for the standardized analysis of immunogenetics data.

*number of trimmed nucleotides*and

*N-REGION length*(number of added nucleotides) can be obtained [12]. However, these numbers represent what is observed in the final "output" but do not necessarily represent the extent of the "true" trimming or nucleotide addition processes. Indeed, randomly trimmed nucleotides can be replaced by

*identical*randomly added N region nucleotides. As a consequence, the number of trimmed V or J nucleotides (represented by the dots in Figure 1) will sometimes be underestimated.

There is therefore a need to quantify this bias if we want to investigate the underlying processes. The goal of the present article is to explore this possibility using TRA and TRG trimming processes, where only V and J genes are involved [10].

Our strategy is the following: given an IMGT/V-QUEST+JCTA standardized output, we aim to calculate the probabilities of all possible trimming events that are consistent with this output. Then, using many such outputs, we aim to probabilistically transform the set of tool "output" data into a representation of the "true trimming process" (i.e., the amount of trimming that *actually* occurred). This probabilistic framework appears naturally by first taking the "output" dataset and simply calculating the empirical probability that the tool "output" shows that 0,1,2... nucleotides were trimmed. Then, understanding how the tool works, we aim to "correct" these empirical probabilities with respect to the tool's biases. A comprehensive introduction to probability distributions (empirical, true) can be found in [26, 27] and a simple introduction to Bernoulli and Poisson distributions is included in Supplementary Data [see Additional file 1].

A first-order model is presented in Results, along with statistical tests on the transformed probability distributions. A proof of the first-order model and a proposed second-order model (also with proof) can be found in Supplementary Data [see Additional file 1].

## Results and discussion

### A first-order model

As potentially more nucleotides are trimmed in the "true process" than appear to have been trimmed according to the tool "output," we would like to transform the "output" data into "true process" data.

A factor also to take into consideration are the quantities of data at zero (except for TRGJ), which do not match the relatively smooth form of the tool "output" data distributions (see Figures 2 and 3). This may be evidence of a two-step process: either the trimming process is activated, or not. If activated, it follows some as yet unknown law. If not, no trimming occurs. Obviously, if the unknown law also takes the value zero, the fraction of data that takes the value zero would then have two sources (either the first process is not activated, or is activated and the second process gives the value zero). Thankfully, as will be shown under the following first-order model, probabilistically transforming the "output" distribution towards the "true process" distribution (under the hypotheses of the model) does not cause further complications. Indeed, the transformed masses (i.e., fractions of the total number of data found at each possible data value) *above zero* do not depend on the original fraction *at zero*. This means that performing maximum likelihood estimation of the parameters of a two-step process is well-defined on the transformed data.

#### Recovering an estimation of the true process probability distribution

Here we introduce a mathematical result that allows us to recover an estimation of the true process probability distribution of the number of trimmed V nucleotides. This result is almost (but not entirely) valid for the true process probability distribution of the number of trimmed J nucleotides. The potential problem is that IMGT/V-QUEST+JCTA selects the J gene after the V gene (see Methods and [25] for more details), thus there is a non-zero chance that 5'J-REGION nucleotides will accidentally be included in the V gene prediction when there has been no N region nucleotide addition. After reanalyzing the data, we found that in the TRAV-TRAJ dataset, this happened at most 3 times and thus was rare enough to be ignored. However, for the TRGV-TRGJ data, this potentially happened quite often, so estimated probability distribution results for the TRGJ trimming process must be used with caution.

Let ℙ {*B* = *k*} mean 'the probability that *k* 3'V-REGION nucleotides are trimmed under the (unknown) true trimming process distribution *f*_{
B
}'. We want to estimate this for *k* ≥ 0. Let ℙ {*F* = *i*} mean 'the probability that *i* nucleotides appear *k* to have been trimmed.' That is, the random variable *F* represents the 3'V-REGION trimming distribution of the tool "output." We do not know the distribution *f*_{
F
}of *F* exactly, but through our datasets we have an empirical estimate of it.

*f*

_{ F }to estimate

*f*

_{ B }. To begin, Theorem 1 [see Additional file 1] shows that under some simple hypotheses (the 'first-order' model), there is an explicit link between the law of the observed 3'V-REGION tool "output" trimming distribution and the "true" (or more correctly, "bias-corrected": technically, it is "true" only if the hypotheses of the first-order model hold in general) process distribution. Indeed, for any

*k*≥ 1 we find:

*k*= 0 we find:

We call this the (4/3, 1/3) rule. Supposing the first-order hypotheses are correct, we would have for example that the bias-corrected probability that 5 V nucleotides were trimmed is equal to (4/3) the probability the tool "output" gives 5 trimmed nucleotides minus (1/3) the probability it gives 6 trimmed nucleotides. We see indeed that under these hypotheses, transformed fractions of data at each data value *above zero* do not depend on the original fraction of data *at zero*.

We remark that it is unlikely that the probabilities of appearance of A, C, G and T nucleotides in the N region are equal (= 1/4, as is assumed in the first-order model), nor in the 3'V-REGION or 5'J-REGION. A second-order model, giving much more freedom to possible A, C, G and T frequencies (each frequency taking some value between 1/6 and 1/3) can be found in Supplementary Data [see Additional file 1]. In brief, we find that the first-order model approximates well the more general second-order model. Thus for simplicity, the first-order result can be used in the place of the second-order result to form hypotheses on trimming processes.

### Testing the transformed V and J trimming distributions

Under the hypotheses of the first-order model, we transformed the TRA and TRG tool "output" data following the law *f*_{
F
}into probability distributions following the law *f*_{
B
}.

*p*unknown) followed by a Poisson process (parameter

*λ*unknown) if the Bernoulli process gave a success. This meant a density function of:

Maximum likelihood was then performed in order to simultaneously estimate the parameters *p* and *λ*, this being necessary to subsequently test the hypothesis that we are dealing with a two-step Bernoulli-Poisson process having parameters *p* and *λ*.

Given data *x*_{1}, *x*_{2},..., *x*_{
n
}, it is easy to show that maximum likelihood estimation gives the equations *g*(*λ*) = (1 - exp(-*λ*))*C* - *mλ* = 0 and *p* = *m*/*n*(1 - exp(-*λ*)) to be solved, where *m* is the number of *x*_{
i
}> 0 and *C* the sum of the values of the *x*_{
i
}> 0. As *m* and *C* are thus constants given any dataset, we see that resolving *g*(*λ*) = 0 for *λ* then allows us to solve for *p* in the second equation. Upon performing the first-order transformation, we found (*m*, *C*) = (517/3, 708), (580/3, 3286/3), (152, 1682/3), (670/3, 4238/3) for the TRAV, TRAJ, TRGV and TRGJ datasets, respectively.

To see that *g*(*λ*) = 0 has a unique solution (and thus *p* also) here, we first remark that for each of these *m*, *C* > 0, lim_{λ→0}*g'*(*λ*) > 0 and *g''*(*λ*) < 0 for *λ* > 0, lim_{λ→∞}*g'*(*λ*) = -*m* < 0, and *g'*(*λ*) is a continuous function for *λ* > 0. Thus, by the intermediate value theorem, there exists at least one *λ* > 0 such that *g'*(*λ*) = 0, and since *g''*(*λ*) < 0 for *λ* > 0, there is in fact a unique solution, which can be easily found numerically for each given *m*, *C* > 0. Indeed, we find (*p*, *λ*) = (0.83, 4.04), (0.92, 5.65), (0.71, 3.59), (1, 6.31) for the TRAV, TRAJ, TRGV and TRGJ datasets, respectively.

*f*in each of the four cases. We tested the four empirical distributions against the theoretical Bernoulli-Poisson distribution

*f*using Pearson's χ

^{2}test. The null hypothesis ${\mathscr{H}}_{0}$ is that the distribution follows

*f*with parameters (

*p*,

*λ*). In order to keep within the assumptions of the test, the data were re-binned into

*n*= 8, 10, 8 and 9 bins for the TRAV, TRAJ, TRGV and TRGJ trimming distributions, respectively. As shown in [28], since the parameters (

*p*,

*λ*) were initially estimated using maximum likelihood, the degree of freedom lies somewhere between

*n*- 1 -

*r*and

*n*- 1, where

*r*is the number of parameters estimated using maximum likelihood. We have thus that

*r*= 2.

We found χ^{2} = 7.97, 11.93, 7.27 and 31.62 for the TRAV, TRAJ, TRGV and TRGJ trimming distributions, respectively. For TRAV, TRAJ and TRGV, we find that at all standard values of statistical significance (*p* = 0.05, 0.01, 0.005), the null hypothesis is not rejected, and thus it is plausible that the empirical results follow a Bernoulli-Poisson-type law. However, for TRGJ, the null hypothesis is rejected at all of the same values of statistical significance. Thus, as it stands, the Bernoulli-Poisson law hypothesis would seem unlikely for the TRGJ trimming process.

## Conclusion

Exploiting standardized "output" datasets of IMGT/V-QUEST+JCTA, we have shown how to recover, under several hypotheses, a representation of the probability distributions of the "true" (or "bias-corrected") TRAV, TRAJ, TRGV and TRGJ trimming processes.

We proceeded by constructing a simple first-order model, known as the (4/3, 1/3) rule, followed by a second-order model [see Additional file 1] which had more general hypotheses. It it clear that the first-order model is a good approximation to the second-order model. We then showed that a kind of two-step Bernoulli-Poisson distribution could plausibly explain the transformed TRAV, TRAJ and TRGV trimming distributions.

We remark that for the TRA and TRG data available to us, the first-order model is "close" to the original IMGT/V-QUEST output data. This is partially due to the relatively smoothly varying data distributions being only slightly modified by performing the operation 4/3 ℙ {*F* = *k*} - 1/3 ℙ {*F* = *k* + 1} (this would not necessarily be true for more irregular probability distributions). An implication of this, for biologists, is that when hypothesis testing on TRA and TRG data sets, as long as the data is relatively smoothly varying from one value to the next, there should be no problem using the IMGT/V-QUEST+JCTA output data, without transformation. Indeed, for our 4 data sets, the same hypothesis tests gave the same statistical result both on the IMGT/V-QUEST+JCTA output data as well as the first-order transformed data.

The statistical analysis of TR and IG junction sequences is a very young field due to the need of having large, clean datasets, unthought-of until recently. Since processes such as the trimming process examined in this article are very little understood from a physical point of view (i.e., what is the exact series of events? By which enzyme is trimming performed? How is exonuclease activity controlled [29]?), we see this work as opening a window to making hypotheses about the very nature of these physical processes and eventually improve our understanding of the complex molecular mechanisms of V-(D)-J recombination [30–33]. IMGT^{®} standardized criteria will eventually enable dealing with datasets numbering in the thousands or millions, impossible to deal with by hand. Under this framework of much larger datasets, we hope the present work will inspire improved models that eventually allow a series of specific, testable hypotheses to be made.

## Methods

### Datasets

T cell receptor (TR) genes were chosen for their absence of somatic hypermutation (in contrast to the IG) [9, 10]. Among TR, the TRA and TRG rearrangements were selected because these loci have only two types of rearranging genes, V and J, in contrast to the TRB and TRD rearrangements which also have D genes [9]. The TRA dataset consisted of 212 human rearranged TRAV-TRAJ junction sequences, selected after alignment and analysis by the integrated IMGT/V-QUEST+JCTA software [23–25] and for which the output was agreed upon by experts (any sequence with potential but not yet confirmed allelic polymorphisms or with some unusual characteristics in the 3'V-REGION or 5'J-REGION was not included in the dataset). This same dataset was used in [12] to perform some preliminary statistical analyses.

An identical methodology was used to collate a dataset of 220 human rearranged TRGV-TRGJ junction sequences. Figures 2 and 3 show the IMGT/V-QUEST+JCTA output for the 'number of trimmed V (and J) nucleotides' for TRA and TRG, respectively.

### Junction analysis

The methodology for the detailed analysis of the junction is described in [25]. Briefly, IMGT/JunctionAnalysis [25] uses the 3'V-REGION of the 'best' aligned germline V gene and allele identified by IMGT/V-QUEST [23, 24] to analyse the junction and delimit the 3' end of the 3'V-REGION in the analysed sequence (checking as far as possible in the 3' direction until encountering a nucleotide that is different from the germline 3'V-REGION, as by default no mutation is allowed for TR). Then, IMGT/JunctionAnalysis uses the 5'J-REGION of the 'best' aligned germline J gene and allele identified by IMGT/V-QUEST to delimit the 5' end of the 5'J-REGION in the analysed sequence (checking as far as possible in the 5' direction until encountering a nucleotide that is different from the germline 5'J-REGION, as by default no mutation is allowed for TR). The remaining nucleotides between the post-trimming 3'V-REGION and post-trimming 5'J-REGION nucleotides are denoted the N region (or if no trimming has occurred, short nucleotide sequences known as the P3'V-REGION or P5'J-REGION may be present [34, 35]). The variables used for statistical analyses of TRA V-J junctions are described in [12]. The same variables were used for the TRG V-J junctions.

## Declarations

### Acknowledgements

We would like to thank the referees for numerous useful remarks that helped to improve the article. We are grateful to Yan Wu for data analysis, Véronique Giudicelli, Xavier Brochet, François Ehrenmann and Patrice Duroux for their contribution to upgrading the IMGT/V-QUEST software. We thank Gérard Lefranc for fruitful discussion and the IMGT^{®} team for its constant motivation and expertise. KB is the recipient of a doctoral grant from the Ministère de l'Enseignement Supérieur et de la Recherche (MESR) Université Montpellier 2. IMGT^{®} is a registered Centre National de la Recherche Scientifique (CNRS) mark. IMGT^{®} is a National Bioinformatics RIO platform since 2001 (CNRS, INSERM, CEA, INRA) and a National Bioinformatics IBiSA platform since 2007. IMGT^{®} was funded in part by the BIOMED1 (BIOCT930038), Biotechnology BIOTECH2 (BIO4CT960037), 5th PCRDT Quality of Life and Management of Living Resources (QLG2-2000-01287) programmes of the European Union and received subventions from the Réseau National des Génopoles (RNG), Génopole-Montpellier-Languedoc-Roussillon. IMGT^{®} is currently supported by the CNRS, the MESR (Université Montpellier 2 Plan Pluri-Formation), the Région Languedoc-Roussillon (Grand Plateau Technique pour la Recherche), Agence Nationale de la Recherche (ANR BIOSYS06-135457) and the ImmunoGrid project (IST-2004-028069) of the 6th framework programme of the European Union.

## Authors’ Affiliations

## References

- Brack C, Hirama M, Lenhard-Schuller R, Tonegawa S: A complete immunoglobulin gene is created by somatic recombination.
*Cell*1978, 15: 1–14.View ArticlePubMedGoogle Scholar - Sakano H, Hüppi K, Heinrich G, Tonegawa S: Sequences at the somatic recombination sites of immunoglobulin light-chain genes.
*Nature*1979, 280: 288–294.View ArticlePubMedGoogle Scholar - Weigert M, Perry R, Kelley D, Hunkapiller T, Schilling J, Hood L: The joining of V and J gene segments creates antibody diversity.
*Nature*1980, 283: 497–499.View ArticlePubMedGoogle Scholar - Early P, Huang H, Davis M, Calame K, Hood L: An immunoglobulin heavy chain variable region gene is generated from three segments of DNA: VH, D and JH.
*Cell*1980, 19: 981–992.View ArticlePubMedGoogle Scholar - Alt F, Baltimore D: Joining of immunoglobulin heavy chain gene segments: implications from a chromosome with evidence of three D-JH fusions.
*Proc Natl Acad Sci USA*1982, 79: 4118–4122.PubMed CentralView ArticlePubMedGoogle Scholar - Tonegawa S: somatic generation of antibody diversity.
*Nature*1983, 302: 575–581.View ArticlePubMedGoogle Scholar - Okazaki K, Davis D, Sakano H: T cell receptor beta gene sequences in the circular DNA of thymocyte nuclei: direct evidence for intramolecular DNA deletion in V-D-J joining.
*Cell*1987, 49: 477–485.View ArticlePubMedGoogle Scholar - Toda M, Fujimoto S, Iwasato T, Takeshita S, Tezuka K, Ohbayahshi T, Yamagishi H: Structure of extrachromosomal circular DNAs excised from T-cell antigen receptor alpha and delta-chain loci.
*J Mol Biol*1988, 202: 219–231.View ArticlePubMedGoogle Scholar - Lefranc MP, Lefranc G:
*The Immunoglobulin FactsBook*. London, UK: Academic Press; 2001:458.Google Scholar - Lefranc MP, Lefranc G:
*The T cell receptor FactsBook*. London, UK: Academic Press; 2001:398.Google Scholar - Giudicelli V, Ginestoux C, Folch G, Jabado-Michaloud J, Chaume D, Lefranc MP: IMGT/LIGM-DB, the IMGT
^{®}comprehensive database of immunoglobulin and T cell receptor nucleotide sequences.*Nucleic Acids Research*2006, 34: D781-D784.PubMed CentralView ArticlePubMedGoogle Scholar - Bleakley K, Giudicelli V, Wu Y, Lefranc MP, Biau G: IMGT standardization for statistical analyses of T cell receptor junctions: The TRAV-TRAJ example.
*In Silico Biol*2006, 6: 573–588.PubMedGoogle Scholar - Giudicelli V, Lefranc MP: Ontology for Immunogenetics: IMGT-ONTOLOGY.
*Bioinformatics*1999, 15: 1047–1054.View ArticlePubMedGoogle Scholar - Duroux P, Kaas Q, Brochet X, Lane J, Ginestoux C, Lefranc MP, Giudicelli V: IMGT-Kaleidescope, the Formal IMGT-ONTOLOGY paradigm.
*Biochimie*2008, 90: 570–583.View ArticlePubMedGoogle Scholar - Lefranc MP: WHO-IUIS Nomenclature Subcommittee for Immunoglobulins and T cell receptors report.
*Immunogenetics*2007, 59: 899–902.View ArticlePubMedGoogle Scholar - Lefranc MP: WHO-IUIS Nomenclature Subcommittee for Immunoglobulins and T cell receptors report August 13th International Congress of Immunobiology, Rio de Janeiro, Brazil.
*Dev Comp Immunol*2008, 32: 461–463.View ArticlePubMedGoogle Scholar - Lefranc MP, Giudicelli V, Ginestoux C, Bosc N, Folch G, Guiraudou D, Jabado-Michaloud J, Magris S, Scaviner D, Thouvenin V, Combres K, Girod D, Jeanjean S, Protat C, Yousfi Monod M, Duprat E, Kaas Q, Pommié C, Chaume D, Lefranc G: IMGT-ONTOLOGY for Immunogenetics and Immunoinformatics.
*In Silico Biol*2004, 4: 17–29.PubMedGoogle Scholar - Lefranc MP, Clément O, Kaas Q, Duprat E, Chastellan P, Coelho I, Combres K, Ginestoux C, Giudicelli V, Chaume D, Lefranc G: IMGT-Choreography for Immunogenetics and Immunoinformatics.
*In Silico Biol*2005, 5: 45–60.PubMedGoogle Scholar - Lefranc MP, Pommié C, Ruiz M, Giudicelli V, Foulquier E, Truong L, Thouvenin-Contet V, Lefranc G: IMGT unique numbering for immunoglobulin and T cell receptor variable domains and Ig superfamily V-like domains.
*Dev Comp Immunol*2003, 27: 55–77.View ArticlePubMedGoogle Scholar - Lefranc MP, Pommié C, Kaas Q, Duprat E, Bosc N, Guiraudou D, Jean C, Ruiz M, Da Piedade I, Rouard M,
*et al*.: IMGT unique numbering for immunoglobulin and T cell receptor constant domains and Ig superfamily C-like domains.*Dev Comp Immunol*2005, 29: 185–203.View ArticlePubMedGoogle Scholar - Lefranc MP, Duprat E, Kaas Q, Tranne M, Thiriot A, Lefranc G: IMGT unique numbering for MHC groove G-DOMAIN and MHC superfamily (MhcSF) G-LIKE-DOMAIN.
*Dev Comp Immunol*2005, 29: 917–938.View ArticlePubMedGoogle Scholar - Lefranc MP, Giudicelli V, Kaas Q, Duprat E, Jabado-Michaloud J, Scaviner D, Ginestoux C, Clément O, Chaume D, Lefranc G: IMGT
^{®}, the international ImMunoGeneTics information system.*Nucleic Acids Res*2005, 33: D593-D597.PubMed CentralView ArticlePubMedGoogle Scholar - Giudicelli V, Chaume D, Lefranc MP: IMGT/V-QUEST, an integrated software for immunoglobulin and T cell receptor V-J and V-D-J rearrangement analysis.
*Nucleic Acids Res*2004, 32: W435-W440.PubMed CentralView ArticlePubMedGoogle Scholar - Brochet X, Lefranc MP, Giudicelli V: IMGT/V-QUEST: the highly customized and integrated system for IG and TR standardized V-J and V-D-J sequence analysis.
*Nucleic Acids Research*2008, 36: W503–508.PubMed CentralView ArticlePubMedGoogle Scholar - Yousfi Monod M, Giudicelli V, Chaume D, Lefranc MP: IMGT/JunctionAnalysis: the first tool for the analysis of the immunoglobulin and T cell receptor complex V-J and V-D-J JUNCTIONs.
*Bioinformatics*2004, 20: I379-I385.View ArticlePubMedGoogle Scholar - Shao J:
*Mathematical Statistics*. NY, USA: Springer; 2003.View ArticleGoogle Scholar - Silverman B:
*Density Estimation for Statistics and Data Analysis*. FL, USA: Chapman and Hall/CRC; 1992.Google Scholar - Chernoff H, Lehmann EL: The use of maximum likelihood estimates in χ
^{2}tests for goodness of fit.*Ann Math Stat*1954, 25: 579–586.View ArticleGoogle Scholar - Souto-Carneiro MM, Fritsch R, Sepúlveda N, Lagareiro MJ, Morgado N, Longo NS, Lipsky PE: The NF-kappaB canonical pathway is involved in the control of the exonucleolytic processing of coding ends during V(D)J recombination.
*J Immunol*2008, 180: 1040–1049.View ArticlePubMedGoogle Scholar - Market E, Papavasiliou FN: V(D)J recombination and the evolution of the adaptive immune system.
*PLoS Biol*2003, 1: E16.PubMed CentralView ArticlePubMedGoogle Scholar - Schatz DG, Spanopoulou E: Biochemistry of V(D)J recombination.
*Curr Top Microbiol Immunol*2005, 290: 49–85.PubMedGoogle Scholar - Lu H, Schwarz K, Lieber MR: Extent to which hairpin opening by the Artemis:DNA-PKcs complex can contribute to junctional diversity in V(D)J recombination.
*Nucl Acids Res*2007, 35: 6917–6923.PubMed CentralView ArticlePubMedGoogle Scholar - Raval P, Kriatchko AN, Kumar S, Swanson PC: Evidence for Ku70/Ku80 association with full-length RAG1.
*Nucl Acids Res*2008, 36: 2060–2072.PubMed CentralView ArticlePubMedGoogle Scholar - Lafaille JJ, DeCloux A, Bonneville M, Takagaki Y, Tonegawa S: Junctional sequences of T cell receptor gamma delta genes: implications for gamma delta T cell lineages and for a novel intermediate of V-(D)-J joining.
*Cell*1989, 59(5):859–870.View ArticlePubMedGoogle Scholar - Lewis SM: P nucleotide insertions and the resolution of hairpin DNA structures in mammalian cells.
*Proc Natl Acad Sci USA*1994, 91: 1332–1336.PubMed CentralView ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.