# Independent recruitment of Igh alleles in V(D)J recombination.

Clara F Alves-Pereira, Raquel de Freitas, Telma Lopes, Rui Gardner, Filipa

Marta, Paulo Vieira, Vasco M Barreto

► To cite this version:

Clara F Alves-Pereira, Raquel de Freitas, Telma Lopes, Rui Gardner, Filipa Marta, et al.. Independent recruitment of Igh alleles in V(D)J recombination.. Nature Communications, 2014, 5, pp.5623. 10.1038/ncomms6623. pasteur-01164712

HAL Id: pasteur-01164712

https://pasteur.hal.science/pasteur-01164712

Submitted on 17 Jun 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

Received 14 Jul 2014 | Accepted 20 Oct 2014 | Published 17 Dec 2014

DOI: 10.1038/ncomms6623

OPEN

# Independent recruitment of *lgh* alleles in V(D)J recombination

Clara F. Alves-Pereira<fnref n="1" />, Raquel de Freitas<fnref n="1" />, Telma Lopes², Rui Gardner², Filipa Marta<fnref n="1" />, Paulo Vieira³,⁴ & Vasco M. Barreto<fnref n="1" />

How the vast majority of B cells express only one of the two alleles at their immunoglobulin loci remains a biological puzzle. Here, in mice reconstituted with a single haematopoietic stem cell, we demonstrate that each of the two immunoglobulin heavy chain (*Igh*) alleles has a similar probability to be the first to undergo $V_H$ to $DJ_H$ rearrangement. We also observe this similar probability in clones from multipotent and common lymphoid precursors. The extreme biases in the expression of the alleles that we find in more differentiated subsets are mostly due to constraints imposed by early rearrangements. Our data demonstrate that each of the two *Igh* alleles in a B cell behaves independently of the other, up to the moment when a successful rearrangement in one allele triggers a feedback mechanism that prevents further recombination.

<footnote n="1">¹Epigenetics and Soma Laboratory, Instituto Gulbenkian de Ciência, Rua da Quinta Grande, n° 6, 2780-156 Oeiras, Portugal. ²Flow Cytometry Laboratory UIC, Instituto Gulbenkian de Ciência, Rua da Quinta Grande, n° 6, 2780-156 Oeiras, Portugal. ³Unité Lymphopoièse, Institut Pasteur, 25, Rue du Dr Roux, 75724 Paris, France. ⁴INSERM U668, F-75015 Paris, France. Correspondence and requests for materials should be addressed to V.M.B. (email: vbarreto@igc.gulbenkian.pt).</footnote>
Each B cell expresses immunoglobulin heavy chains at its surface, typically encoded by only one of the two alleles¹. The question of why this phenomenon of allelic exclusion occurs has generated a number of non-exclusive and untested answers. For instance, it could be a solution to produce one or a limited set of antigen receptors per cell, as postulated by the clonal selection theory², or a safety mechanism to decrease the chances of detrimental translocations³, since antigen receptor genes undergo somatic recombination⁴. How allelic exclusion arises would seem to be more tractable; however, this is an ongoing question since the 1960s. Mostly because of the imprecise junction of the $V_H$, $D_H$ and $J_H$ gene segments that are found in the immunoglobulin heavy chain (Igh) locus, only a minor proportion of the rearrangements will be 'productive' or $V_H D_J^+$ (encoding a full-length chain), which intrinsically ensures some level of allelic exclusion. However, for the murine Igh locus, a minimal stochastic model of allelic exclusion based exclusively on the outcome of the rearrangement process⁵ is insufficient to account for the degree of exclusion observed at the surface of the cells⁶.

The $D_H$ to $J_H$ rearrangement step occurs in both Igh alleles⁷. This event has recently been proposed to mark the $DJ_H$ junction for secondary rearrangement, that is, the $V_H$ to $DJ_H$ stage⁸, and once a productive secondary rearrangement takes place in one of the alleles of the pro-B cell, expression of the pre-B-cell receptor (pre-BCR) at the surface signals a series of events leading to the shutdown of recombination at the second allele. This feedback mechanism, first suggested for the immunoglobulin kappa locus⁹, and subsequently elaborated for the Igh locus¹⁰, has been independently confirmed in a number of ways. Notably, a µH chain encoded by a transgene or a knocked-in allele inhibits $V_H$ to $DJ_H$ rearrangement at the endogenous Igh locus only if expressed at the surface, as part of the pre-BCR¹¹⁻¹³. In addition, in cells with genotypic inclusion (with both alleles rearranged productively) only one of the chains pairs efficiently with the surrogate light chain¹⁴ that also forms the pre-BCR and the frequency of cells with genotypic inclusion ($V_H D_J^+ / V_H D_J^+$) increases in the absence of the surrogate light chain¹⁵ or when the signalling from an intact pre-BCR is perturbed by mutations in one of its downstream targets¹⁶. More recently, two additional feedback mechanisms have been proposed, one triggered by an Ataxia Telangiectasia Mutated (ATM)-dependent sensing of the rearrangement at the first Igh allele³ and the other by a stable µH mRNA with a premature stop codon that cannot be translated into a µH chain¹⁷. Overall, feedback inhibition and the rearrangement of the second allele in its absence largely explain the proportion of $V_H D_J^+ / D_J^+$ and $V_H D_J^+ / V_H D_J^+ / $ IgM⁺ B cells¹⁸.

The events initially leading to $V_H$ to $DJ_H$ rearrangement that render the $V_H$ regions 'accessible' for rearrangement include the relocalization of the Igh locus away from the nuclear periphery in pro-B cells¹⁹ and long-range interactions mediated by looping of the Igh domains²⁰,²¹ that depend on regulatory regions, namely the one located in the $V_H$ to $D_H$ intergenic region²². However, it is unclear whether a specific mechanism ensures that each cell only tests one $V_H$ to $DJ_H$ rearrangement at a time. The acquisition of competence to undergo rearrangement by the Igh locus in the course of B-cell differentiation has been described at the molecular level (for example, ref. 23); however, typically these studies are carried out in bulk cultures, averaging out any existing differences between the alleles at the single-cell level. In turn, imaging-assisted studies (for example, refs 19–21,24) can discriminate between the alleles at the single-cell level, but measure features that are only correlates of $V_H(D_H)J_H$ accessibility, which are not always confirmed by different techniques²⁵, and provide mere snapshots of a process that is essentially dynamic.

Here we performed a clonal analysis of Igh allelic exclusion, and provide conclusive evidence that the allelic exclusion of immunoglobulin receptor loci differs from X-chromosome inactivation, as no stable epigenetic mark is propagated until pro-B cells start rearranging. More importantly, we also show that the Igh alleles rearrange independently from each other in the majority of clones derived from common lymphoid precursors.

## Results

**Clonal analysis tracks V(D)J recombination.** Clonal analysis is typically used to investigate the commitment of precursor cells to a given lineage or set of lineages. We have adapted this technique to test whether an Igh allele from a precursor cell is committed to recombine a $V_H$ segment before the other allele as the cell divides and enters the lymphoid lineage. As shown in Fig. 1, although in a bulk culture it is not possible to conclude whether the alleles from cells heterozygous for the Igh locus (Igh$^{a/b}$) are differently committed for rearrangement, in isolated clones the commitment and non-commitment scenarios produce sharply different outcomes. In the case of commitment, all the cells expressing the µ chain from the first allele to rearrange (in Fig. 1a, IgM$^a$) would have the second allele in the non-$V_H$-rearranged $DJ_H$ configuration and the remaining IgM⁺ cells (in Fig. 1a, IgM$^b$) would have the second allele in the $V_H DJ_H$ configuration. In contrast, if commitment does not occur (Fig. 1b), the IgM$^a$ and IgM$^b$ pools of cells would not differ in the status of rearrangement of the silenced allele, since each pool would have the same proportion of IgM⁺ cells with a unique $V_H$ to $DJ_H$ rearrangement versus cells that were rescued by the second $V_H$ to $DJ_H$ rearrangement. Thus, $D_H$ to $J_H$ and $V_H$ to $DJ_H$ junctions function as a natural reporter system, since the signature junction sequences somatically introduced in the Igh locus help to ascertain when rearrangement occurred in each clone as the cells divide.

**Both Igh alleles can rearrange first in clones derived from HSC.** The main feature of the Igh locus and other antigen receptor loci is V(D)J recombination; however, these genes are frequently included in a growing group of autosomal genes displaying what has been called random monoallelic expression²⁶. Unlike loci under parental origin-specific imprinting, in these genes roughly 50% of the cells express the maternal allele and the remaining 50% express the paternal allele. Since the parallel with X-chromosome inactivation has been made in several reports²⁷⁻²⁹ and haematopoietic stem cells (HSCs) show X-chromosome inactivation, we first asked whether in mice with single-HSC-derived haematopoietic systems, the emerging IgM$^a$ and IgM$^b$ populations would also show signs of a commitment that is already present in the HSC and that is then stably propagated.

Cells heterozygous for a pan-leukocitary marker polymorphism (CD45.1/45.2), the Igh allotypes (Igh$^{a/b}$) and the X-linked gene *Foxp3* and its *gfp*-fused knock-in (*Foxp3gfp/wt*) were single-cell-sorted from the HSC-long-term repopulation potential-enriched Lin⁻Rho⁻CD45intHoe⁻ side population (SP) pool of C57BL/6J female bone marrow (BMs; see Methods; Supplementary Fig. 1). Sublethally irradiated B6-Rag2⁻/⁻ was then single-cell-reconstituted and monitored over time to detect animals showing chimerism in the blood. The percentage of reconstituted animals was 4.8–22.9%. This low efficiency is in the range of what has been described³⁰, and further increases the likelihood of obtaining a collection mainly composed of animals with monoclonal haematopoietic systems. Animals with chimerism (N=15) were euthanized within 8–18 weeks from the single-cell transplant, analysed by flow cytometry for the presence of lymphoid and myeloid cells and in some cases used to produce secondary
<figure><img src="image_3.png" /></figure>

**Figure 1 | Rationale.** The commitment and no commitment models of *Igh* allelic choice predict a distinctive $V_H$ to $D_J_H$ rearrangement status in the two IgM<sup>a</sup> and IgM<sup>b</sup> fractions of B-cell clonal populations. (a) In a B-cell clone that is pre-determined to rearrange first a given allele (represented by a black locus with a white star), all the cells rearrange the marked allele first and the second allele (grey locus) remains non-rearranged, unless the previous rearrangement is non-productive ($V_HDJ_H^-$). The resulting clone contains two populations of B cells: one that expresses the product of the first rearrangement (a productive rearrangement, $V_HDJ_H^+$) and retains 50% of the alleles non-rearranged (IgM<sup>a</sup> population in the scheme) and another that expresses the product of the second rearrangement (IgM<sup>b</sup> population in the scheme) and has both alleles rearranged. (b) The alternative model (no commitment) assumes that both alleles are equally accessible, each with equal chance of rearranging first. This condition renders the IgM<sup>a</sup> and IgM<sup>b</sup> populations of a cell clone undistinguishable at the structure of the silent *Igh* allele, since both populations will have a similar frequency of silent alleles in the $V_HDJ_H$ and $DJ_H$ configurations. The present scheme assumes the ordered model of V(D)J rearrangement<sup>7</sup>, and represents the $V_H$ to $DJ_H$ rearrangement only (the D-to-J rearrangement is biallelic and has occurred before the developmental stage illustrated). In a sorted population of IgM<sup>+</sup> cells, the frequency of non-rearranged alleles can be measured by quantifying the $V_H-D_H$ intergenic fragment retention within the *Igh* locus.

recipients, attesting to their long-term repopulation potential (Table 1). The collection analysed for allelic exclusion comprised animals that were long-term and short-term reconstituted.

Female donor cells heterozygous both for the X-linked *Foxp3* and *Igh* loci were used to compare in the same experiment X-chromosome inactivation with *Igh* allelic exclusion. FOXP3 is a transcription factor expressed by regulatory T cells (Tregs). Owing to X-chromosome inactivation, in heterozygous animals for a knock-in green fluorescent protein (GFP fusion protein; *Foxp3*<sup>gfp/wt</sup>)<sup>31</sup> only 50% of Tregs are green (Fig. 2a). Control animals reconstituted with 1,000 HSC-enriched SP cells reproduced this balanced bimodal distribution of the GFP signal in Tregs (Fig. 2a, Table 1), which became biased when less cells were used in the reconstitution and strictly unimodal (0% GFP or 100% GFP) in monoclonal haematopoietic systems. We confirmed that each individual HSC has one inactivated X-chromosome and that inactivation is a stable epigenetic property, functioning here as an internal control for the monoclonality of the reconstitutions.

We then analysed *Igh* allelic exclusion in the single HSC-reconstituted animals. We consistently found roughly equal numbers of IgM<sup>a</sup> and IgM<sup>b</sup> B cells per animal (Fig. 2a and Supplementary Table 2), similar to what is found in polyclonal haematopoietic systems<sup>6</sup>. This observation is in sharp contrast to the X-chromosome inactivation data and suggests the lack of stable allelic commitment in the *Igh*; however, as explained above, the definitive read-out relies on the examination of the rearrangement status of the non-expressed allele. Thus, a quantitative PCR assay was designed to measure the fraction of alleles that retain the $V_H-D_H$ intergenic region, which is excised as an episome during $V_H$ to $DJ_H$ recombination and is diluted out as the pre-BCR-expressing cell divides and progresses in development (Fig. 2b). The assay reports the expected differences in the retention of the $V_H-D_H$ intergenic fragment

in $Rag2^{-/-}$ pro-B cells (100%) and splenic B cells from wild-type (WT) animals (roughly 60% of cells in $V_HDJ_H/DJ_H$ configuration and 40% in the $V_HDJ_H/V_HDJ_H$ configuration<sup>32</sup>, which equals 30% of retention), $Jht^{+/-}$ IgM<sup>+</sup> cells (the knockout allele cannot rearrange and therefore all IgM<sup>+</sup> cells are in the $V_HDJ_H$/germline configuration<sup>33</sup>, which equals 50% intergenic fragment retention) and $Atm^{-/-}$ (which have an intergenic fragment retention below that of WT animals<sup>3</sup>). If allelic commitment exists before rearrangement, then the isolated IgM<sup>a</sup> and IgM<sup>b</sup> subpopulations of a given clone must differ drastically in the amount of the $V_H-D_H$ intergenic fragment retained (50 versus 0%); however, we did not detect any difference in the four animals analysed (Fig. 2c). Similarly, if allelic pre-commitment exists, the cells would carry a productive rearrangement in the pre-committed allele (either IgM<sup>a</sup> or IgM<sup>b</sup>) together with no $V_H$ rearrangement in the other, while the cells expressing the other allele would have roughly the same number of productive and non-productive rearrangements. Since the $V_HDJ_H^+/V_HDJ_H^-$ ratio may change depending on the $V_H$ family, we matched the type of families selected for analysis per isolated fraction, focusing on the most prominent families: J558 and 7183. We did not find any striking difference in the frequency of $V_HDJ_H^+$ and $V_HDJ_H^-$ rearrangements of IgM<sup>a</sup> and IgM<sup>b</sup>-isolated splenic B cells from the monoclonally reconstituted animals (Fig. 2d). Taken together, these data conclusively show that, contrary to X-chromosome inactivation, there is no stable propagation of an allelic commitment to rearrange first the $V_H$ segment in the clones that emerge from a single HSC and differentiate *in vivo* into pro-B cells.

**IgM<sup>a/b</sup> expression becomes skewed in clones derived from CLP.** We next asked whether *Igh* allelic commitment for $V_H$ to $DJ_H$ rearrangement would emerge at some stage during the lineage
differentiation to pro-B cells, by applying the previous clonal approach to progressively more differentiated precursor cells. Given the low expansion potential of more differentiated precursors compared with the HSC, we did not reconstitute mice with single B-cell precursors, and used instead a cell culture system. Single cells from pools enriched for HSCs, multipotent progenitors (MPPs), common lymphoid precursors (CLPs) or Pro-B cells were plated under stimulation conditions that are routinely used to mimic the development that occurs in the BM³⁴. As expected, the B-cell-plating efficiency increased from HSCs to CLPs (see Methods; Supplementary Fig. 4), while the

time required to detect surface IgM expression decreased. In Fig. 3a, we plot the percentage of IgMª-expressing cells (relative to total IgM⁺) per clone growing from seeded HSCs, MPPs and CLPs. As expected, within each group the distribution of IgMª- and IgMᵇ-expressing cells is similar. However, the dispersion of each of these percentages is much more pronounced in the clones emerging from CLPs than in the clones resulting from the HSCs and the MMP progenitors (CLP versus MPP, P = 4.5 × 10⁻¹² by F-test; CLP versus HSC, P = 1.1 × 10⁻⁵ by F-test). Notably, the collection of the CLP-derived clones includes examples of clones completely skewed to IgMª or IgMᵇ (Fig. 3b). Given the

<table><caption>**Table 1 | Single HSC-reconstituted animals and representative controls.**</caption><thead><tr><th>Sample (no.)</th><th>No. injected cells</th><th>Graft age (weeks)</th><th>(%) CD45.1<sup>+</sup> spleen</th><th>(%) CD45.1<sup>+</sup> LN</th><th>(%) CD45.1<sup>+</sup> BM</th><th>(%) CD45.1<sup>+</sup> thymus</th><th>FOXP3<sup>-</sup> GFP<sup>-</sup></th><th>FOXP3<sup>-</sup> GFP<sup>+</sup></th></tr></thead><tbody><tr><td>1.</td><td>4.1_2.7M</td><td>1</td><td>10.8</td><td>26.6</td><td>14.6</td><td>1.6</td><td>2.9</td><td>0.0</td><td>100.0</td></tr><tr><td>2.</td><td>4.2_2.7M</td><td>1</td><td>10.8</td><td>18.6</td><td>32.5</td><td>2.4</td><td>1.8</td><td>0.0</td><td>100.0</td></tr><tr><td>3.</td><td>5.1_3.5M</td><td>1</td><td>14.0</td><td>71.0</td><td>27.8</td><td>73.6</td><td>—</td><td>0.6</td><td>99.4</td></tr><tr><td>4.</td><td>5.2_3.5M</td><td>1</td><td>14.0</td><td>24.9</td><td>89.9</td><td>10.9</td><td>46.4</td><td>99.6</td><td>0.4</td></tr><tr><td>5.</td><td>8.1_3.7M</td><td>1</td><td>14.8</td><td>59.9</td><td>69.5</td><td>2.3</td><td>0.1</td><td>98.3</td><td>1.7</td></tr><tr><td>6.</td><td>4.6_4.0M</td><td>1</td><td>16.0</td><td>41.6</td><td>12.4</td><td>2.3</td><td>3.0</td><td>0.0</td><td>100.0</td></tr><tr><td>7.</td><td>4.3_4.5M</td><td>1</td><td>18.0</td><td>15.4</td><td>1.6</td><td>1.3</td><td>0.1</td><td>—</td><td>—</td></tr><tr><td>8.</td><td>4.5_4.5M</td><td>1</td><td>18.0</td><td>23.9</td><td>72.3</td><td>1.6</td><td>3.7</td><td>100.0</td><td>0.0</td></tr><tr><td>9.</td><td>8.2_5.4M</td><td>1</td><td>21.6</td><td>39.3</td><td>3.1</td><td>17.8</td><td>2.6</td><td>0.0</td><td>100.0</td></tr><tr><td>10.</td><td>8.3_5.4M</td><td>1</td><td>21.6</td><td>35.7</td><td>—</td><td>3.2</td><td>6.3</td><td>—</td><td>—</td></tr><tr><td>11.</td><td>8.4_5.4M</td><td>1</td><td>21.6</td><td>32.7</td><td>55.1</td><td>4.1</td><td>32.6</td><td>0.6</td><td>99.4</td></tr><tr><td>12.</td><td>8.5_5.4M</td><td>1</td><td>21.6</td><td>—</td><td>59.8</td><td>6.2</td><td>18.8</td><td>99.6</td><td>0.4</td></tr><tr><td>13.</td><td>8.6_5.4M</td><td>1</td><td>21.6</td><td>19.4</td><td>—</td><td>7.2</td><td>7.6</td><td>—</td><td>—</td></tr><tr><td>14.</td><td>8.7_5.4M</td><td>1</td><td>21.6</td><td>37.5</td><td>—</td><td>6.4</td><td>19.2</td><td>—</td><td>—</td></tr><tr><td>15.</td><td>5.4_6.9M</td><td>1</td><td>27.4</td><td>76.7</td><td>82.1</td><td>13.6</td><td>8.6</td><td>99.8</td><td>0.2</td></tr><tr><td>16.</td><td>4.1_7.9M_2nd_a</td><td>1</td><td>31.6</td><td>1.6</td><td>2.1</td><td>0.6</td><td>30.8</td><td>0.0</td><td>100.0</td></tr><tr><td>17.</td><td>4.1_7.9M_2nd_b</td><td>1</td><td>31.6</td><td>8.1</td><td>50.7</td><td>1.4</td><td>1.4</td><td>0.6</td><td>99.4</td></tr><tr><td>18.</td><td>5.1_7.4M_2nd_a</td><td>1</td><td>29.6</td><td>37.1</td><td>68.1</td><td>2.6</td><td>6.1</td><td>0.8</td><td>99.2</td></tr><tr><td>19.</td><td>c4.7_2.2M</td><td>6</td><td>8.8</td><td>—</td><td>73.1</td><td>—</td><td>12.5</td><td>99.5</td><td>0.5</td></tr><tr><td>20.</td><td>c4.4_4.5M</td><td>6</td><td>18.0</td><td>51.5</td><td>25.0</td><td>13.0</td><td>0.2</td><td>100.0</td><td>0.0</td></tr><tr><td>21.</td><td>c4.12_7.9M</td><td>12</td><td>31.6</td><td>18.2</td><td>22.3</td><td>0.3</td><td>0.3</td><td>—</td><td>—</td></tr><tr><td>22.</td><td>c8.49_3.7M</td><td>&lt;50</td><td>14.8</td><td>78.9</td><td>77.6</td><td>14.3</td><td>72.8</td><td>97.4</td><td>2.6</td></tr><tr><td>22.</td><td>c4.6_2.2M</td><td>50</td><td>8.8</td><td>27.0</td><td>78.9</td><td>—</td><td>57.3</td><td>6.5</td><td>93.5</td></tr><tr><td>23.</td><td>c8.50_3.7M</td><td>50</td><td>14.8</td><td>25.7</td><td>82.9</td><td>12.1</td><td>0.2</td><td>13.4</td><td>86.6</td></tr><tr><td>24.</td><td>c5.1000_2.7M</td><td>1,000</td><td>10.8</td><td>69.7</td><td>19.4</td><td>13.9</td><td>69.6</td><td>42.4</td><td>57.6</td></tr><tr><td>25.</td><td>B6.Foxp3<sup>wt/wt</sup></td><td>GFP<sup>-/-</sup> Control</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>99.9</td><td>0.1</td></tr><tr><td>26.</td><td>B6.Foxp3<sup>gfp/gfp</sup></td><td>GFP<sup>+/-</sup> Control</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.2</td><td>99.8</td></tr><tr><td>27.</td><td>B6.Foxp3<sup>gfp/wt</sup></td><td>Donor Control</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>53.5</td><td>46.5</td></tr></tbody></table>

BM, bone marrow; CWBM, clonal whole bone marrow; GFP, green fluorescent protein; HSC, haematopoietic stem cell; LN, lymph node.

Entries 1-15 correspond to single-cell-injected mice, 16-18 are secondary recipients of 5,000 WBM cells of clone no. 4.1_2.7 M (line 1) and no. 5.1_3.5 M (line 3), lines 19-24 correspond to animals injected with more than one cell, and animals 25-27 are non-reconstituted controls. Graft age is the time in weeks at which the sample was collected after the injection of the single cell. Frequency of CD45.1⁺ cells in the spleen, LN, BM and thymus is presented if available. The last two columns (Foxp3 - GFP⁺⁻) show the frequency of LN CD45.1⁺ Foxp3⁺ cells that were staining positive or negative for Foxp3⁻ -GFP⁺ by flow cytometry (representative histograms are shown in Fig. 2a). The animal no. 8.1_3.7 M was excluded from further analysis because a small population of GFP⁺ cells were observed in addition to GFP⁻ -Foxp3⁺ cells. On the basis of these data, we estimate that the percentage of false monoclonal reconstitutions within the selected pool of animals with unimodal distributions of GFP in the regulatory T-cell population is below 5%.

**Figure 2 | Direct test for a stable *Igh* epigenetic commitment in HSC that dictates which allele undergoes *V<sub>H</sub>* to *DJ<sub>H</sub>* rearrangement first in the B lineage.** (a) The expression of FOXP3GFP in lymph node (LN)-derived regulatory T cells is bimodal in the non-clonal populations, and unimodal in mice reconstituted with a single HSC (top histograms), but the expression of the *Igh* alleles is balanced in the splenic B cells of both the non-clonal and single-cell-reconstituted mice (bottom contour plots). Top histograms represent GFP expression profiles of CD45.1⁺CD4⁺FOXP3⁺-gated LN cells, and bottom contour plots show the CD45.1⁺CD19⁺-gated splenic B-cell surface IgMª/IgMᵇ profiles of the same clonal and non-clonal samples. (b) Genomic quantification with qPCR of *in vivo*-differentiated B-cell clones. The retention of the *V<sub>H</sub>*-*D<sub>H</sub>* intergenic fragment in IgMª and IgMᵇ-sorted splenic populations from each of the single-HSC clones—the same samples shown in **a**—is similar to each other, and in the same range of IgM⁺-sorted cells from non-clonal populations. The qPCR on DNA from IgM⁺ cells purified from *Jht*⁺⁻ mice (that can only rearrange one of the *Igh* alleles) and *Atm*⁻⁻ mice (which have a higher frequency of *V<sub>H</sub>*-*DJ<sub>H</sub>*/*V<sub>H</sub>*-*DJ<sub>H</sub>* cells than the WT³,²²) produced the expected results. Error bars are the s.e.m. of three independent runs of the same samples in different input amounts. (c) The reduction in the amount of the *V<sub>H</sub>*-*D<sub>H</sub>* intergenic fragment from the fractions A and B of pro-B cells to immature B cells shows that the assay is reporting the changes in the structure of the *Igh* locus throughout development. (d) Comparison of the frequency of productive and non-productive *V<sub>H</sub>*-*DJ<sub>H</sub>* rearrangements detected in the splenic IgMª and IgMᵇ-sorted pools from clone no. 5.1_3.5 M (P = 1.00, Fisher's exact test) and no.5.2_3.5 M (P = 0.55, Fisher's exact test).
heterogeneous nature of the CLP populations, we decided to further dissect them to identify the earliest population containing progenitors with potential for the highly biased expression of

IgMª or IgMᵇ. The surface marker Ly6d (Fig. 3c) identifies one population (Ly6d⁻) that conserves the full lymphoid potential and one population (Ly6d⁺) that is to date the earliest B-cell

<figure><img src="image_4.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 3** | IgM<sup>a/b</sup> expression in CLP-derived clones. As the developmental stage of the isolated precursor cell approaches the stage when recombination typically occurs, the distribution of IgM<sup>a</sup>/IgM<sup>b</sup>-expressing cells in the resulting B-cell clones differentiated *in vitro* becomes biased. The graphs show *P* values from an F-test. (a) IgM<sup>a</sup> percentage (IgM<sup>a</sup>/[IgM<sup>a</sup> + IgM<sup>b</sup>]) obtained from flow cytometric analysis of B-cell clones differentiated in culture from HSCs, MPPs and CLPs; each dot represents a clone, open and closed dots distinguish two independent assays. (b) Representative examples of flow cytometric profiles of CLP-derived clones presenting extreme bias in the expression of IgM<sup>a</sup> or IgM<sup>b</sup> (clone no. 054 and no. 003) and a mild biased/balanced (clone no. 039). (c) Left: density plot shows the Lin<sup>-</sup>IL7R<sup>+</sup>-gated BM population, staining cKit<sup>int</sup> and Sca-1<sup>int</sup>, thus identified as the CLP population. Right: histogram shows the Ly6d-staining profile in the CLP-gated population. (d) Percentages of IgM<sup>a</sup>-expressing cells in clones derived from CLPs partitioned according to the presence of the Ly6d marker<sup>35</sup>. The percentage of IgM-expressing cell clones emerging per CLP.Ly6d<sup>-</sup> plated cell was 2%, while Ly6d<sup>+</sup> cells render IgM<sup>+</sup> B cell clones at least 10 × more frequently. (e) Distribution of the percentages of heavy chain alleles (IgM<sup>a</sup>/[IgM<sup>a</sup> + IgM<sup>b</sup>]) and light chain alleles (hCκ/[hCκ + mCκ]) in CLP Ly6d<sup>+</sup>-derived clones from the same experiment. (f) CLP Ly6d<sup>+</sup> clone no. 010 flow cytometric profile for the immunoglobulin heavy chain alleles (IgM<sup>a</sup>/IgM<sup>b</sup>) and light chain alleles (hCκ/mCκ).</figcaption></figure>
progenitor³⁵. We could confirm that, under conditions promoting B-cell differentiation, the B-cell-plating efficiency of the Lyd6⁻ CLP-derived clones was at least 10 × lower than in Ly6d⁺ CLP-derived clones (in one experiment, 15 wells presented colonies out of 960 single-cell Lyd6⁻ CLP-plated cells, while 68 wells presented colonies out of 384 single-cell-plated Lyd6⁺ CLP-plated cells). We observed that the highly skewed clones are more frequent in the Ly6d⁺ pool (Fig. 3d), although the difference is not statistically significant. We have further differentiated B-cell clones from whole B220⁺CD43⁺ pro-B cells and from the fraction B and found that, although the lower expansion potential of pro-B cells originated smaller clones, CLP- and B220⁺CD43⁺-derived clones have a similar frequency of skewed clones (Supplementary Fig. 4c). Fraction B presented an even lower expansion potential, which precluded further molecular analysis. We have also compared *Igh* and *Igk* loci with respect to the allelic bias observed in *in vitro* CLP-derived clones from mice heterozygous for both *Igh* and light chain immunoglobulin loci (B6.Cg-*Igh*<sup>a/b</sup>, *hCk/mCk*). The *Igh* locus typically rearranges before the kappa light chain and usually four to six cell divisions occur between the rearrangement of these immunoglobulin genes³⁶,³⁷. As observed in Fig. 3e,f, the bias observed for the *Igh* is not reproduced for the *Igk* (*P* = 0.00046 by F-test) presumably because the number of cells rearranging the *Igk* alleles independently is much higher than the number of cells rearranging the *Igh* alleles within the clone. The balanced ratio of the Igk-expressing cells also suggests that the *Igk* locus is not pre-committed in the CLP stage of B-cell development. Taken together, these data show that both the probability of finding clones with allelic bias and the degree of bias increase as the stage at which the cells are seeded gets closer to the time of V(D)J rearrangement.

**Independent recruitment of *Igh* alleles in CLP-derived clones.**

Different scenarios could explain the distribution of the IgMª and IgMᵇ percentages in CLP-derived clones (Fig. 3a). In the first scenario, the clones would be synchronized for V(D)J rearrangement and in each clone there would be a wave of IgMª (or IgMᵇ)-expressing cells followed by a partially overlapping second wave of IgMᵇ (or IgMª) cells. In this scenario, consistent with an asymmetric allele choice for rearrangement within each clone, the time point at which the clones were analysed would only be a snapshot, failing to capture the ongoing dynamics in the cultures. We have sampled the cultures at different time points and found no clear evidence for the existence of such waves in the majority of clones (Supplementary Fig. 5). The remaining scenarios can be distinguished by the status of rearrangement of the silent allele in the IgMª and IgMᵇ populations of each clone. This status was estimated with the real-time quantitative PCR (qPCR) of the *V_H-D_H* intergenic fragment and also by sequencing *V_HDJ_H* and *DJ_H* rearrangements from the sorted populations, taking into account polymorphisms downstream of *J_H1*, *J_H2* and *J_H3* segments that enable us to assign rearrangements to the BALB/c (*Igh*ª) or C57BL/6 (*Igh*ᵇ) allotypes. We first examined the CLP-derived clones with both the IgMª and IgMᵇ pools clearly discernible and balanced with the help of flow cytometry (Fig. 4a,b). In the genome of these IgM⁺ populations, the percentage of *V_H-D_H* intergenic fragment retention when compared with the Rag2⁻⁻ genome was typically between 40% (similar to the mean retention observed in non-clonal *in vitro*-derived IgM⁺ cells, Fig. 4a) and 50% (the percentage in the genome of non-clonal *in vitro*-derived IgM⁺ cells of *Jht*⁺⁻ mice). More importantly, within the same clone, the IgMª and IgMᵇ genomes retain approximately the same percentage of the *V_H-D_H* intergenic fragment, which means that each allele can undergo *V_H* to *DJ_H* rearrangement first (Figs 1 and 4a).

Furthermore, we sequenced *DJ_H* and *V_HDJ_H* rearrangements and found diverse *V_H-D_H* junctions as well as diverse *D_H-J_H* junctions, and a limited number of these *D_H-J_H* junctions contributing to the diverse *V_HDJ_H* rearrangements per clone. This shows that the Ly6d⁺ single-cell-sorted cells are devoid of *D_H-J_H* rearrangements before culture, and suggests that in some clones we are testing subclonal expansions and tracking the choice to rearrange shortly before the decision to undergo *V_H* to *DJ_H* rearrangement. Representatively, in clone Xp40, the same two *DJ_H* junctions (one originating from the *Igh*ᵇ allele and the other from the *Igħ*ª allele) prevail in the genome of the IgMª, IgMᵇ and IgM⁻ pools (Table 2). The *Igħ*ᵇ *DJ_H* junction is shown to be predominantly contributing to the diverse *V_HDJ_H⁺* rearrangements of the IgMᵇ⁺ pool and in IgMª pool the same *DJ_H* junction is part of non-productive (*V_HDJ_H⁻*) and the same junction was detected in the germline-*DJ_H* configuration, while the dominant *Igħ*ª *DJ_H* junction was detected as part of productive (*V_HDJ_H⁺*) rearrangements in the IgMª cell pool and in non-productive and *DJ_H* rearrangements in the IgMᵇ pool. Different *DJ_H* junctions are detected in all three fractions of the clone but were not detected as part of *V_HDJ_H* rearrangements. From the analysis of rearrangements we conclude that the culture system mimics the features of V(D)J recombination *in vivo* and that considerable intraclonal diversity is generated in our clones, extending previous observations made in human B-cell lines³⁸. More importantly, we observed that any of the alleles can undergo *V_H* to *DJ_H* rearrangement first within each clone.

**Founder effects in highly skewed clones.** The most striking behaviour of the CLP-derived clones is the extreme bias displayed by some clones (Figs 3a,b and 5). This pattern could emerge from a pre-determined choice to rearrange a given allele first, if the culture conditions do not favour the second attempt to rearrange *V_H* to *DJ_H*. The comparison between the amount of *V_H-D_H* intergenic fragment in IgM⁺ non-clonal cell populations differentiated *in vivo* and *in vitro* shows that in culture conditions there is indeed a minor increase in the retention of the *V_H-D_H* intergenic fragment (41.6 ± 5.8% versus 34.1 ± 4.9%), but a sizable proportion of cells have two rearrangements (Supplementary Fig. 5). Thus, it is difficult to explain the extreme skewing observed in these clones solely by attributing a fixed choice in the allele to undergo *V_H* to *DJ_H* rearrangement. Next, the status of rearrangement of the expressed and silent alleles in the clones was investigated. In the 18 IgM⁺ populations of the CLP-derived clones with both IgMª and IgMᵇ clearly distinguishable (Fig. 4a), no example was found of a population with a *V_H-D_H* intergenic fragment below 30%. However, in the 10 IgM⁺ populations from CLP-derived clones with extreme biases, three populations have a percentage of *V_H-D_H* intergenic fragment below 30% (*P* = 0.037, Fisher's exact test). Interestingly, two of those clones were essentially depleted of the *V_H-D_H* intergenic fragment (Fig. 5a), showing that these populations are mostly composed of cells that also have the silenced allele in the *V_HDJ_H* configuration. This observation suggests that the extreme bias in this subset of clones results from a rearrangement occurring early on in the culture, that supposedly is non-productive and permanently silences the allele. In the case of the populations without a *V_H-D_H* intergenic fragment, such rearrangement could be a *V_HDJ_H⁻* or a *V_HDJ_H⁺* encoding a μ chain that cannot form the pre-BCR¹⁴, in the case of the populations with a sizable *V_H-D_H* intergenic fragment, such rearrangements could be *DJ_H* junctions using the reading frame 3 that frequently contains premature stop codons³⁹. For the CLP-derived B-cell clone Xp29, which does not retain the *V_H-D_H* intergenic fragment (Fig. 5a), sequencing data support the hypothesis of a non-productive rearrangement-imposed founder
<figure><img src="image_9.png" /><figcaption>**Figure 4** | Rearrangement status of the silent *Igh* allele in CLP-derived clones. The status of rearrangement of the silent allele in IgMª and IgMᵇ-sorted cells from *in vitro* CLP-derived B-cell clones shows independent allele recruitment for rearrangement. (a) *V*<sub>H</sub>-*D*<sub>H</sub> intergenic fragment retention quantification data are shown in two blocks: in black, a battery of control profiles: *Rag2*<sup>-/-</sup> and *Jht*<sup>-/-</sup> CD19<sup>+</sup> cells, that do not rearrange the *Igh* locus, *Jht*<sup>+/−</sup> as in Fig. 2c, the PD31 cell line, which lacks the *V*<sub>H</sub>-*D*<sub>H</sub> intergenic fragment<sup>55</sup>, and additional WT controls, all differentiated *in vitro*. In grey, data from pairs of IgMª- and IgMᵇ-sorted samples from mildly biased/balanced clones. Sample Yn4 is a CLP-Ly6d<sup>-</sup>-derived clone, and the remaining samples are CLP-Ly6d<sup>+</sup>-derived. Error bars denote s.e.m. (b) Flow cytometric profiles of the WT bulk control and clones represented in the graph. (c) Top panel: schematic representation of the diversity of rearrangements obtained by sequencing *V*<sub>H</sub>*DJH* and *DJH* rearrangements of IgMᵇ, IgMª and IgM<sup>-</sup> populations of the clone Xp40 (see Table 2 for a detailed description). Legend: long brush strokes represent D-J joints and short ones the *V*<sub>H</sub>-*D*<sub>H</sub> joints. Sequences of *Igh*ª, *Igh*ᵇ and unknown allotype are represented by cold, warm and greyish colours, respectively. A long stroke not followed by a short stroke represents a *DJH* rearrangement. A dot represents a non-productive rearrangement. Different colours in longer strokes represent diverse D-to-J joints as different colours in shorter strokes represent diverse *V*<sub>H</sub> to *DJH* joints. The same colours represent the same *DJH* joint and identical *VDJH* joints are not represented. Left lower panel: clone Xp40 flow cytometric profile as in b. Right lower panel: map showing polymorphisms in the *JH* cluster of *Igh*ª (Balb/c) and *Igh*ᵇ (C57BL/6) allotypes, the *JH* segment position and relative positions of 3'-*JH*1', -*JH*2 and -*JH*3 primers (arrowheads, *JH*4 primer position not shown in map).</figcaption></figure>
<table><caption>Table 2 | Sequencing V<sub>H</sub>DJ<sub>H</sub> and DJ<sub>H</sub> junctions of clone Xp40.</caption><thead><tr><th>ID*</th><th>V<sub>H</sub>DJ<sub>H</sub> (CDR3)/DJ<sub>H</sub> junction (5'-3')</th><th>V<sub>H</sub></th><th>D<sub>H</sub></th><th>RF‡</th><th>J<sub>H</sub></th><th>A§</th><th>F</th><th>R||</th></tr></thead><tbody><tr><td colspan="9"><strong>Xp40, IgM<sup>b+</sup></strong></td></tr><tr><td colspan="9"><strong>V<sub>H</sub>DJ<sub>H</sub></strong></td></tr><tr><td>2.1.</td><td>tgt gca aga <u>aga gat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>J558.26.116</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.2.</td><td>tgt gca aga <u>ggg agg gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>V23</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.3.</td><td>tgt gcc aga <u>agc ctc tat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>Q52.2.4</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.4.</td><td>tgt gcc aga <u>aat tgg gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>Q52.9.29</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.5.</td><td>tgt gca aga <u>agg ggc tat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>36-60.8.74</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.6.</td><td>tgt gca aga <u>ttt gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>36-60.8.74</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.7.</td><td>tgt aca aga <u>ggc tat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>7183.14.25</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.8.</td><td>tgt gca aga <u>ggt cgg gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>7183.21b</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.9.</td><td>tgc aca <u>ggg cta gat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>J606.1.79</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.10.</td><td>tgc aca <u>ggc ctg gat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>J606.1.79</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>11.1.</td><td>tgt gca aga <u>gaa gac ccc ctc ccg</u> ttt gct tac tgg</td><td>J558.84.190</td><td>?</td><td>—</td><td>3</td><td>b</td><td>p</td><td>2</td></tr><tr><td>1.1.</td><td>tgt gcc aga <u>aat act acg g.t</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>Q52.a2.4</td><td>FL16.2</td><td>2</td><td>3</td><td>a</td><td>np</td><td>2</td></tr><tr><td>1.2.</td><td><u>tgac gca cat tac tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>7183.a2.3</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>np</td><td>2</td></tr><tr><td>1.3.</td><td>tgt acc agg <u>gca cta cgg.</u> <u>tt cc</u> gcc tgg ttt gct tac tgg</td><td>J606.a6.127</td><td>FL16.2</td><td>3</td><td>3</td><td>a</td><td>np</td><td>2</td></tr><tr><td colspan="9"><strong>DJ<sub>H</sub></strong></td></tr><tr><td>1.</td><td>tt cat tac tac <u>ggt ccc gcc</u> tgg ttt gct tac tgg</td><td>—</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>—</td><td>8</td></tr><tr><td>3.</td><td><u>tt tat tac tac ggt agt agc</u> tat gct atg gac tac tgg</td><td>—</td><td>FL16.1</td><td>1</td><td>4</td><td>?</td><td>—</td><td></td></tr><tr><td>4.</td><td>tt tat tac tac <u>ggt agt agc</u> tat tac tat gct atg gac tac tgg</td><td>—</td><td>FL16.1</td><td>1</td><td>4</td><td>?</td><td>—</td><td></td></tr><tr><td>5.</td><td><u>tt tta cta cgg tag tag</u> cat tac tat gct atg gac tac tgg</td><td>—</td><td>FL16.1</td><td>3</td><td>4</td><td>?</td><td>stp</td><td></td></tr><tr><td colspan="9"><strong>Xp40, IgM<sup>b+</sup></strong></td></tr><tr><td colspan="9"><strong>V<sub>H</sub>DJ<sub>H</sub></strong></td></tr><tr><td>2.11.</td><td>tgt gat cta <u>tgga aa..</u> <u>ac</u> tac ttt gac tac tgg</td><td>J558.21pg.111</td><td>SP2.9</td><td>3</td><td>2</td><td>b</td><td>np</td><td></td></tr><tr><td>2.12.</td><td>tgt gca <u>tgga ggc ggg gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>J558.1.85</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>np</td><td></td></tr><tr><td>2.13.</td><td>tgt gca aga <u>agc tga aa..</u> <u>ac</u> tac ttt gac tac tgg</td><td>J558.26.116</td><td>SP2.9</td><td>3</td><td>2</td><td>b</td><td>np</td><td></td></tr><tr><td>2.14.</td><td>tgt gca aga <u>tgga cta tga aa..</u> <u>ac</u> tac ttt gac tac tgg</td><td>J558.72.173</td><td>SP2.9</td><td>3</td><td>2</td><td>b</td><td>np</td><td></td></tr><tr><td>2.15.</td><td>tgt gca aga <u>tat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>X24.1pg.45</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td>2</td></tr><tr><td>2.16.</td><td>tgt gca aga <u>aag gga</u> <u>aac</u> tac ttt gac tac tgg</td><td>X24.2.50</td><td>?</td><td>—</td><td>2</td><td>b</td><td>p</td><td>2</td></tr><tr><td>2.17.</td><td>tgt gca aga <u>ctatga aa..</u> <u>ac</u> tac ttt gac tac tgg</td><td>7183.a2.3</td><td>SP2.9</td><td>3</td><td>2</td><td>b</td><td>np</td><td></td></tr><tr><td>2.18.</td><td>tgt gca aga <u>gat gtg aaa..</u> <u>ac</u> tac ttt gac tac tgg</td><td>7183.21b</td><td>SP2.9</td><td>2</td><td>2</td><td>b</td><td>np</td><td></td></tr><tr><td>1.4.</td><td>tgt gca agg <u>aac tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>J558.47</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td></td></tr><tr><td>1.5.</td><td>tgt gca aga <u>tcg gga ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>124</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td></td></tr><tr><td>1.6.</td><td>tgt gca aga <u>ttcc ctc tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>J558.7</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td></td></tr><tr><td>1.7.</td><td>tgt gca aga <u>gaa ggg gac</u> <u>ttcc gcc</u> tgg ttt gct tac tgg</td><td>S107.a3.106</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td>2</td></tr><tr><td>12.1</td><td>tgt gca aga <u>ttcc tac tat gat ggt</u> <u>ttcc gct</u> tac tat gct atg gac tac tgg</td><td>J558.84.190</td><td>SP2.9</td><td>1</td><td>4</td><td>?</td><td>p</td><td>4</td></tr><tr><td>13.1</td><td>tgt gca aga <u>cat agt aac.</u> t gct atg gac tac tgg</td><td>7183.2.3</td><td>SP2.x</td><td>1</td><td>4</td><td>?</td><td>np</td><td>2</td></tr><tr><td colspan="9"><strong>DJ<sub>H</sub></strong></td></tr><tr><td>2.</td><td><u>tc tat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>—</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td></td><td></td></tr><tr><td>14.</td><td><u>tt tat tac tac ggt agt agc</u> <u>ccg</u> ttt gct tac tgg</td><td>—</td><td>FL16.1</td><td>1</td><td>3</td><td>b</td><td></td><td></td></tr><tr><td>15.</td><td><u>tt tat tac tac ggt agt agc</u> ttt gct tac tgg</td><td>—</td><td>FL16.1</td><td>1</td><td>3</td><td>b</td><td></td><td></td></tr><tr><td>16.</td><td><u>ttt att act acg gta gta</u> <u>gct gcc</u> tgg ttt gct tac tgg</td><td>—</td><td>FL16.1</td><td>2</td><td>3</td><td>b</td><td></td><td></td></tr><tr><td>17.</td><td><u>c cta ctatag tag</u> <u>ggc gcc</u> tgg ttt gct tac tgg</td><td>—</td><td>SP2.x</td><td>3</td><td>3</td><td>b</td><td>stp</td><td></td></tr><tr><td>18.</td><td><u>t tta cta cgg</u> tac tat gct atg gac tac tgg</td><td>—</td><td>FL16.1</td><td>3</td><td>4</td><td>?</td><td></td><td></td></tr><tr><td>19.</td><td><u>ttt att act acg gta gta</u> <u>gca tat</u> tac tat gct atg gac tac tgg</td><td>—</td><td>FL16.1</td><td>2</td><td>4</td><td>?</td><td></td><td></td></tr><tr><td>20.</td><td><u>tt tat tac tac ggt agt agc</u> tac tat gct atg gac tac tgg</td><td>—</td><td>FL16.1</td><td>1</td><td>4</td><td>?</td><td></td><td></td></tr><tr><td colspan="9"><strong>Xp40, IgM<sup>-</sup></strong></td></tr><tr><td colspan="9"><strong>V<sub>H</sub>DJ<sub>H</sub></strong></td></tr><tr><td>1.4.</td><td>tgt gca agg <u>aac tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>J558.17</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td>2</td></tr><tr><td>1.8.</td><td>tgt aca aga <u>gcg att cat tac tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>J558.40</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td></td></tr><tr><td>1.9.</td><td>tgt gcc aga <u>gag cat tac tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>Q52.a27.79</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td>3</td></tr><tr><td>1.10.</td><td>tgt gcc agt <u>cac tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>Q52.a2.4</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td></td></tr><tr><td>1.11.</td><td>tgt gca agg <u>ttc att act acg gtc</u> <u>cc</u> gcc tgg ttt gct tac tgg</td><td>36-60.a6.114</td><td>FL16.2</td><td>2</td><td>3</td><td>a</td><td>np</td><td></td></tr><tr><td>1.12.</td><td>tgt gca aga <u>ttcc cac tac ggl[c]</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>36-60.a9.121</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td></td></tr><tr><td>1.13.</td><td>tgt gca <u>atc cat tac tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>36-60.a9.121</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td>2</td></tr><tr><td>1.14.</td><td>tgt gtc aga <u>cgaggg</u> <u>ttc tac ggt</u> <u>ccc gcc</u> ggg ttt gct tac tgg</td><td>7183.a19.31</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td>2</td></tr><tr><td>1.15.</td><td>tgt gca aga <u>ta cta cgg tcc</u> <u>c.. gcc</u> tgg ttt gct tac tgg</td><td>81X</td><td>FL16.2</td><td>3</td><td>3</td><td>a</td><td>np</td><td></td></tr><tr><td>1.16.</td><td>tgt acg agg <u>ttc ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>7183.a33.55</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td></td></tr><tr><td>1.17.</td><td>tgt gca <u>aaa tca att</u> <u>ttc tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>X24.a1.84</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td>3</td></tr><tr><td>1.18.</td><td>tgt gca aga <u>ggg aat tac ggt</u> <u>ccc gcc</u> tgg ttt gct tac tgg</td><td>X24.a1.84</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td>p</td><td></td></tr><tr><td>1.19.</td><td>tgt atg aga <u>ta cta cgg tcc</u> <u>c.. gcc</u> tgg ttt gct tac tgg</td><td>11.a2.92</td><td>FL16.2</td><td>3</td><td>3</td><td>a</td><td>np</td><td>3</td></tr><tr><td>1.20.</td><td>tgt atg aga <u>tgg gga cta cgg</u> <u>ttc c</u> <u>c.. gcc</u> tgg ttt gct tac tgg</td><td>11.a2.92</td><td>FL16.2</td><td>3</td><td>3</td><td>a</td><td>np</td><td></td></tr><tr><td>2.19.</td><td>tgt gca <u>att tat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>J558.26.116</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.20.</td><td>tgt aca aga <u>ttc tat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>J558.3.90</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.21.</td><td>tgt gca <u>ta tat gaa</u> <u>aac</u> tac ttt gac tac tgg</td><td>J558.78.182</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.22.</td><td>tgt gcc aga <u>atg aaa .ac</u> <u>ac</u> tac ttt gac tac tgg</td><td>Q52.2.4</td><td>SP2.9</td><td>2</td><td>2</td><td>b</td><td>np</td><td>4</td></tr></tbody></table>
<table><caption>Table 2 (Continued)</caption><thead><tr><th>ID<sup>*</sup></th><th>V<sub>H</sub>DJ<sub>H</sub> (CDR3)/DJ<sub>H</sub> junction (5'-3')<sup>†</sup></th><th>V<sub>H</sub></th><th>D<sub>H</sub></th><th>RF<sup>‡</sup></th><th>J<sub>H</sub></th><th>A<sup>§</sup></th><th>F</th><th>R<sup>||</sup></th></tr></thead><tbody><tr><td>2.23.</td><td>tgt gca aga tat <b>ggg</b> atc tat gaa <b>aac</b> tac ttt gac tac tgg</td><td>36-60.8.74</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.24.</td><td>tgt gca aga <b>tcc</b> <b>tcc</b> ttc tat gaa <b>aac</b> tac ttt gac tac tgg</td><td>36-60.8.74</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.25.</td><td>tgt gca agg <b>cac</b> tat gaa <b>aac</b> tac ttt gac tac tgg</td><td>7183.7.10</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.26.</td><td>tgt acc agg <b>gaa</b> <b>aac</b> tac ttt gac tac tgg</td><td>J606.4.82</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.27.</td><td>tgc aca <b>ggc</b> <b>ccc</b> ctc tat gaa <b>aac</b> tac ttt gac tac tgg</td><td>J606.1.79</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td></td></tr><tr><td>2.28.</td><td>tgt acc agc tat gaa <b>aac</b> tac ttt gac tac tgg</td><td>J606.4.82</td><td>SP2.9</td><td>1</td><td>2</td><td>b</td><td>p</td><td>2</td></tr><tr><td>DJ<sub>H</sub></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1.</td><td>tt cat tac tac <b>ggt</b> <b>ccc</b> gcc tgg ttt gct tac tgg</td><td>—</td><td>FL16.2</td><td>1</td><td>3</td><td>a</td><td></td><td>6</td></tr><tr><td>6.</td><td>tt tat tac tac <b>ggt</b> agt agt <b>ccc</b> tgg ttt gct tac tgg</td><td>—</td><td>FL16.1</td><td>1</td><td>3</td><td>b</td><td></td><td>4</td></tr><tr><td>7.</td><td>t tta tta cta <b>cga</b> <b>tag</b> agg ttt gct tac tgg</td><td>—</td><td>FL16.1</td><td>3</td><td>3</td><td>b</td><td>stp</td><td></td></tr><tr><td>8.</td><td>ttt att act acg <b>gta</b> <b>gta</b> <b>ta</b> c tat gcc atg gac tac tgg</td><td>—</td><td>FL16.1</td><td>2</td><td>4</td><td>?</td><td></td><td>2</td></tr><tr><td>9.</td><td>t tta tta cta <b>cggtag</b> <b>tag</b> cta <b>ta</b> tat tac tat gct atg gac tac tgg</td><td>—</td><td>FL16.1</td><td>3</td><td>4</td><td>?</td><td>stp</td><td></td></tr><tr><td>10.</td><td>t tta tta cta <b>cgg</b> <b>tag</b> <b>tag</b> ctt gct atg gac tac tgg</td><td>—</td><td>FL16.1</td><td>3</td><td>4</td><td>?</td><td>stp</td><td></td></tr></tbody></table>

A, allotype; F, functionality prediction; p, productive; np, non-productive; stp, stop codon in DJ<sub>H</sub> conformation.

*Each line preceded by a number shows the CDR3 sequences or the D-J<sub>H</sub> junctions of B-cell clone-sorted populations (IgM<sup>a</sup> + or IgM<sup>b</sup> + genomic fractions) sequenced from amplicons produced using a battery of 5' primers specific for the mouse V<sub>H</sub> families or DJ<sub>H</sub> gene segments and 3' J<sub>H</sub>-specific primers (see Methods). Numbers are used as sequence identification tags. The first number before the first dot identifies a unique DJ<sub>H</sub> junction, and the second number (between dots) identifies a unique V<sub>H</sub>DJ<sub>H</sub>. In the table, a total of 20 different DJ<sub>H</sub> junctions are described, but only five were also detected within V<sub>H</sub>DJ<sub>H</sub> rearrangements ('1', '2', '11', '12' and '13'). Of these DJ<sub>H</sub> junctions, two are predominant among all fractions of the clone: an Ig<sup>h</sup>-derived DJ<sub>H</sub>2 joint and an Ig<sup>h</sup>-derived DJ<sub>H</sub>3 junction. These recombined independently with multiple V<sub>H</sub>-gene segments distributed along the locus: DJ<sub>H</sub> sequence number '1', was embedded in N = 20 V<sub>H</sub>DJ<sub>H</sub> sequences (1.1.-1.2.0.), and DJ<sub>H</sub> sequence number '2', that was embedded in N = 28 V<sub>H</sub>DJ<sub>H</sub> sequences (2.1.-2.28.). These two particular DJ<sub>H</sub> sequences were detected in all three populations of the clone, with the V<sub>H</sub>DJ<sub>H</sub> configuration presumably corresponding to the expressed allele, and V<sub>H</sub>DJ<sub>H</sub> and DJ<sub>H</sub> corresponding to the non-expressed one. The existence of many independent V<sub>H</sub>DJ<sub>H</sub> sequences of the non-expressed allele in both IgM<sup>+</sup> fractions is indicative of a recruitment of V<sub>H</sub>-gene segments for V<sub>H</sub> to DJ<sub>H</sub> rearrangement on each allele that is independent of the other allele. The diversity of V<sub>H</sub>DJ<sub>H</sub> rearrangements sharing the same pair of DJ<sub>H</sub> joints, and a number of independent DJ<sub>H</sub> joints not detected in V<sub>H</sub>DJ<sub>H</sub> rearrangements in the IgM<sup>+</sup>, IgM<sup>b</sup> and IgM<sup>-</sup> populations, further suggests a subclonal expansion within the clone.

†The analysis was initially performed by running the IMGT/HIGHV-QUEST<sup>56</sup>, and subsequently the sequences were manually curated. The CDR3 was defined as the starting cysteine at position 104 and ending at the JH 5' invariant tryptophan. Underlined sequences are the nucleotides attributed to DJ<sub>H</sub> sequences; N and P additions are in bold; frame shifts are represented by dots; stop codons are in italics; parenthesis are used to highlight probable PCR errors. This notation applies to all sequencing data tables.

‡F, DH-reading frame.

§a, Ig<sup>h</sup>-coding for IgM<sup>a</sup> and 'b' designates Ig<sup>h</sup>-coding for IgM<sup>b</sup>.

||R, the number of times the same sequence was detected (by Sanger sequencing) each time a band was cloned: typically, four colonies were screened per cloned band, and in some relevant cases where no diversity was found, more colonies were sequenced from the same cloned bands (usually in batches of four).

effect: the V<sub>H</sub>DJ<sub>H</sub> rearrangements detected are diverse for the productive V<sub>H</sub> to DJ<sub>H</sub> junctions, but limited concerning both the DJ<sub>H</sub> to J<sub>H</sub> junctions (N=2 detected), and the V<sub>H</sub> to DJ<sub>H</sub> non-productive junction (unique sequence detected, several times; Fig. 5c and Supplementary Table 3). Thus, this clone appears to have made a first or very few non-productive V<sub>H</sub> to DJ<sub>H</sub> rearrangements early on in the culture, by chance in the same allele, which was followed by a number of cell divisions and then an extensive rearrangement of the second allele. Clones 13-33, Xp39, no. 5.3 and no. 4.2, which also show absolute biases (Fig. 5b) retain an intermediate amount of V<sub>H</sub>-DJ<sub>H</sub> intergenic fragment (Fig. 5a). Sequencing data for clone no. 5.3 (showing a stop codon in the DJ<sub>H</sub> joint of one of the alleles) and no. 4.2 (with a DJ<sub>H</sub> joint that lacks a highly conserved codon) also confirm the founder effects imposed by rearrangement (Fig. 5c and Supplementary Tables 4 and 5). In clones 13-2, 13-18 Yp62 and Yp4, the silent allele is mostly in the DJ<sub>H</sub> configuration (Fig. 5a), which suggests that one of the alleles may have a lower probability to rearrange. Sequencing analysis of the IgM<sup>b</sup>-skewed clone 13-18 indeed shows that the IgM<sup>+</sup> population is composed of cells that have rearranged V<sub>H</sub> to DJ<sub>H</sub> only in the Ig<sup>h</sup><sup>b</sup> allele, leaving Ig<sup>h</sup><sup>a</sup> in the DJ<sub>H</sub> configuration (Fig. 5c and Supplementary Table 6). Interestingly, the IgM<sup>-</sup> population sequencing analysis shows V<sub>H</sub>DJ<sub>H</sub> rearrangements in both alleles, with Ig<sup>h</sup><sup>a</sup> rearrangements that are more frequently non-productive (seven non-productive with frequent stop codons out of eight sequences) than the Ig<sup>h</sup><sup>b</sup> ones (one non-productive out of 10 sequences). These results do not exclude that, for a few subclonal expansions, one of the alleles acquires an increased probability of rearrangement. Overall, we conclude that a minority of the CLP-derived clones with extreme skewing may have a relative bias to rearrange predominantly one of the alleles and that the majority of CLP-derived clones with extreme biases is explained by founder effects imposed by rearrangement that introduce a genetic constraint in the ability to rearrange a given allele in the precursor cells from which the IgM<sup>+</sup> cells emerge.

## Discussion

Here we provide the first conclusive evidence that the allelic exclusion of immunoglobulin receptor loci differs from X-chromosome inactivation, as no stable epigenetic mark is propagated until pro-B cells start rearranging. More importantly, we also show that the Ig<sup>h</sup> alleles rearrange independently from each other in the majority of clones derived from common lymphoid precursors.

Roughly, the number of cells expressing the maternal or the paternal Ig<sup>h</sup> allele is the same. Similarly, in female mammals each cell inactivates the paternal or the maternal X-chromosome with identical probability. This common feature is frequently emphasized, and it has been hypothesized that autosomal genes under monoallelic expression of one or the other allele are flanked by regions with a high density of long interspersed nuclear elements<sup>40</sup>, the same elements that play a key role in the generation and propagation of the silencing epigenetic marks along the X-chromosome<sup>41</sup>. Furthermore, it is known that the alleles replicate asynchronously in parent-of-origin imprinting<sup>42</sup>, other cases of monoallelic expression<sup>28,43</sup> and X-linked inactivation<sup>44</sup>. In particular, antigen receptor loci were found to establish this pattern early in development, around the time X-inactivation is also established, and then proceed to maintain it<sup>28</sup>. This led to the proposal that the allele to replicate first retains its epigenetic mark as the cell divides and undergoes differentiation, so that in clonally related developing B cells, it is always the same allele that undergoes V<sub>H</sub> to DJ<sub>H</sub> rearrangement first<sup>28</sup>. The analysis of haematopoietic systems resulting from a single female HSC is a direct test of this hypothesis because the HSC population is heterogeneous for the inactivated X-chromosome, and if the Ig<sup>h</sup> locus carries an epigenetic mark that is transmitted to the progeny as stably as the X-chromosome marks, V<sub>H</sub> to DJ<sub>H</sub> rearrangement would occur first in the same allele in all, or the vast majority, of the pro-B cells. As expected<sup>45</sup>, X-chromosome inactivation is stably transmitted during all cell divisions that occur as a HSC differentiates into a regulatory T
<figure><img src="image_10.png" /><figcaption>**Figure 5** | Rearrangement status of the silent *Igh* allele in highly skewed CLP-derived clones. (a) Quantitative PCR of the V<sub>H</sub>-D intergenic DNA fragment in IgM<sup>+</sup> fractions of the *in vitro*-derived B-cell clones expressing only one allele (grey bars) and *in vitro*-derived controls. Error bars denote s.e.m. (b) CD19<sup>+</sup>-gated clones' dot plots showing IgM<sup>a</sup>/IgM<sup>b</sup> flow cytometric profiles for some of the clones in **a**. (c) Top panels' dot plots show the remaining clones in **a**, not shown in **b**; shaded quadrants in the plots show the fractions that are represented in the tree-like scheme of the lower panels. Each scheme represents the diversity of V(D)J rearrangements obtained by sequencing *VDJ<sub>H</sub>* and *DJ<sub>H</sub>* rearrangements of the clone fractions indicated, as described in Fig. 4. From the represented clones, only clone 13–18 profile (for the IgM<sup>+</sup> fraction) is compatible with a model of allelic preference. Clone no. 5.3 has a stop codon in the *Igh<sup>a</sup>*-derived *DJ<sub>H</sub>* joint that prevents IgM<sup>a</sup>-productive rearrangements; clone no. 4.2 contains a prohibitive rearrangement in the *Igh<sup>b</sup>*-derived *DJ<sub>H</sub>* joint (see Supplementary Table 6); clone Xp29 contains a single non-productive *V<sub>H</sub>DJ<sub>H</sub>* rearrangement in the *Igh<sup>a</sup>*-derived allele as inferred by phenotype.</figcaption></figure>

cell; in contrast, from the status of rearrangement of the silent allele in IgM<sup>a</sup>- and IgM<sup>b</sup>-sorted populations we have found that in cells clonally related to a single HSC, either the maternal or the paternal alleles can rearrange the V<sub>H</sub> first with similar frequency, thus showing that there is no stable epigenetic mark established

early in development that defines a temporal allelic hierarchy for V<sub>H</sub> to DJ<sub>H</sub> rearrangement. Farago et al.<sup>46</sup> have reported the same conclusion for the *Igκ* locus. In their study, multiple HSCs carrying distinguishable *Igκ* alleles (hCκ/mCκ) were transduced with lentiviruses (that integrate the genome of each cell in a
unique manner, thus identifying its clonal progeny), and then used in bulk to reconstitute mice. This technique was instrumental in the past to establish whether cell lineages share a precursor cell⁴⁷; however, we argue that it provides inconclusive results for allelic commitment. Since it is well established that 40–50% of B cells rearrange both *Igκ* alleles⁴⁸, the progeny of each precursor cell is inevitably present in both the hCk⁺ and mCk⁺ populations and it follows that the pattern of lentiviral integrations in the isolated hCk⁺ and mCk⁺ populations is expected to be roughly the same, with or without commitment. Instead, from the *in vivo* clonal analysis here described based on the analysis of the status of rearrangement of the silent allele, we conclude that the parallel between the allelic exclusion of *Igh* and X-chromosome inactivation is misleading. It thus appears that either allele from the daughter pro-B cells of a given HSC can undergo *VH* to *DIH* first, and that asynchronous replication is not a good predictor of the fate of the alleles, in V(D)J recombination and other contexts. Indeed, it has been shown that replication timing is stable after X-inactivation but switches between alleles before the onset of the inactivation and is independent from the subsequent choice to inactivate the X-chromosome⁴⁹.

The data from CLP-derived clones revealed that both *Igh* alleles have the potential to undergo *VH* to *DIH* first, until shortly before the onset of rearrangement. It cannot be completely excluded that the accessibility of the *Igh* locus in culture may differ from the *in vivo* conditions, particularly at the level of the initiation of the *Igh* locus accessibility, nuclear localization of the alleles and overall rearrangement kinetics. However, our culture conditions are consensually accepted as mimicking the B-cell development that occurs *in vivo* in the bone marrow³⁴, allelic exclusion was observed in the flow cytometric profiles of the clones, *in vitro* experiments performed months apart in different laboratories with cytokine supernatants from different batches produced similar results and a direct comparison of rearrangement efficiencies *in vitro* and *in vivo* did not reveal statistically significant differences (Supplementary Fig. 6). Strikingly, most of the clones with an extreme bias in the expression of IgMª or IgMᵇ could be explained as a consequence of rearrangement dead ends, that is, those that permanently prevent one allele to be productively expressed at the membrane, giving the other the chance to undergo multiple productive rearrangements within the clone. This explanation of a bias determined by a limiting number of rearranging cells is compatible with the more balanced IgMª/IgMᵇ ratios in clones derived from HSCs and MPPs, in which rearrangement starts relatively late in the cultures, after the initial cell has expanded substantially. It also explains the slightly more pronounced IgMª/IgMᵇ bias in Ly6d⁺-derived clones compared with the clones derived from the precursor Ly6d⁻, since the former start rearranging in culture earlier than the latter. Finally, it makes sense that, in CLP-derived clones, the allele rearrangement bias observed for the *Igh* alleles is not reproduced by the *Igκ* alleles, which typically rearrange after the substantial expansion of cells expressing the pre-BCR³⁷,⁵⁰. This last result is the only one addressing the *Igκ* locus in this work and does not reproduce the results presented in ref. 46. Further analysis is needed to address this issue thoroughly.

From these results, it can be concluded that the mechanism of allelic exclusion of the *Igh* locus differs from the one claimed for the *Igκ* locus⁴⁶. In that study, it was proposed that stable epigenetic marks acquired around the time of lymphoid differentiation pre-determine a temporal hierarchy for rearrangement; however, the status of rearrangement of the silenced allele in the clones analysed was not reported, which limits the interpretation of the data. It is nevertheless possible that the *Igh* and the *Igκ* loci behave differently, particularly given the organization of the *VH* and *JH* elements that allows successive

rearrangements in the *Igκ* alleles and the existence of *D*<sub>H</sub> elements in the *Igh* locus. Indeed, the *Igh* locus has more structural similarities with the *Tcrβ* locus than with the *Igκ* locus, namely the ordered rearrangement. Of note, the *Tcrβ* alleles were found to independently associate with pericentromeric heterochromatin⁵¹, an interaction that inhibits rearrangement, and this observation is compatible with our conclusions for the *Igh* locus. Furthermore, this model is also the simplest explanation. Unlike the ‘pre-determination’ model, in which the strict temporal hierarchy for recombination implies an *ad hoc* mechanism to render the second allele active for recombination, in the model of independent recruitment the second rearrangement will eventually occur unless the first rearrangement triggers the feedback mechanism that shuts down recombination or the cell dies. The idea of independent recruitment is also compatible with the existence of a proportion of cells with ‘synchronous’ rearrangements, that is, cells in which the second *VH* to *DIH* rearrangement occurs in the time window before the feedback mechanism triggered by the formation of the pre-BCR blocks further rearrangement of the *Igh* locus (Supplementary Fig. 7). It is not known what proportion of cells rearrange the alleles ‘synchronously’. However, in the absence of complementary feedback mechanisms mediated by *ATM*³ and the μ chain messenger RNA¹⁷, which respond to the first rearrangement faster than the feedback triggered by the pre-BCR, there is a measurable, although small, increase in the frequency of allelically included cells. We have confirmed this increase in genomic allelic inclusion for *Atm*⁻/⁻ cells in our assay (Fig. 2c). Finally, this model of independent rearrangement is not incompatible with the existence of moderate and transient allelic biases in the probability to undergo *VH* to *DIH* rearrangement, if each allele behaves independently from the other while becoming accessible to rearrangement, but not from its previous copy in the mother cell. Thus, the simplest interpretation of our data is that in the vast majority of clones the two alleles gain recombination competence in a sufficiently synchronized way for the V(D)J mechanism to reveal its stochastic component. Any case of an initial asymmetry that may persist for a number of divisions and lead to biases in the recombination does not contradict the model. Indeed, we may have identified such possible bias in only one of the clones studied, and these rare cases are expected from a truly independent opening of the alleles (Supplementary Fig. 7).

For over three decades, the explanations of allelic exclusion have essentially differed on the weight given to the stochastic and the regulated events. Using *in vivo* and *in vitro* clonal approaches that are based on the outcome of the rearrangement process in genetically unmanipulated murine cells, we found that an independent recruitment of the *Igh* alleles underlies the mechanism of allelic exclusion of the *Igh* locus. This intrinsic stochastic feature establishes the default mechanism that ensures two rearrangement attempts per cell, unless the first rearrangement triggers the definitive feedback that shuts down further rearrangement.

## Methods

**Mice.** All mice were bred and maintained at the specific pathogen-free animal facilities of the Instituto Gulbenkian de Ciência (Oeiras, Portugal). BALB/c, C57BL/6j, B6.Rag2⁻/⁻, B6.SJL-Ptprca (the pan-leukocytic marker gene Ptprca is also known as CD45.1 or Ly5.1; the C57BL/6j strain carries the Ptprcb allele, also known as CD45.2 or Ly5.2) and B6.Cg-Ighª were originally received from The Jackson Laboratory (Bar Harbor, ME, USA). Mice with Foxp3gfp knock-in³¹ (B6.Foxp3gfp, courtesy of Dr Alexander Rudensky—Memorial Sloan-Kettering Cancer Center, New York, USA) were backcrossed with the C57BL/6j strain at our animal facility. Animals were bred to generate female heterozygous B6.Cg-CD45.1/45.2Igha/b Foxp3gfp/wt and used in cell transfer experiments. Female 8–12 week-old-mice were
used as donors. For the *in vitro* cell culture system (BALB/c x C57BL/6)F1, B6.Cg-*Igha/b* or B6;129.*Igha/b*Igk*<sup>h/m</sup> mice heterozygous for the *Igh* and/or the *Igk* loci were used (B6;129.*Igk*<sup>h</sup> mice<sup>52</sup> were courtesy of Dr Michel Nussenzweig (Rockefeller University, NY, USA). Mice with a targeted deletion for the immunoglobulin heavy chain joining region (B6.Igh-Jht<sup>-/-</sup> mice<sup>33</sup>) were courtesy of Dr Jocelyne Demengeot-IGC, Oeiras, Portugal) or ATM knockout mice (B6;129.*Atm*<sup>-/-</sup> mice<sup>53</sup>, courtesy of Dr Almudena Ramiro-CNIO, Madrid, Spain). This research project was reviewed and approved by the Ethics Committee of the Instituto Gulbenkian de Ciência, and by the Portuguese National Entity that regulates the use of laboratory animals (Direcção Geral de Alimentação e Veterinária. All experiments conducted on animals followed the Portuguese (Decreto-Lei n° 113/2013) and European (Directive 2010/63/EU) legislations concerning housing, husbandry and animal welfare.

**Cell suspensions.** All cells were suspended in Hanks' balanced-salt solution (HBSS, Gibco) supplemented with 5% fetal bovine serum (FBS, Life Technologies). BMs were flushed and single-cell-suspended from the femora and tibiae with a 26-gauge needle of 1-ml syringe, filtered and rinsed with HBSS 5% FBS. Thymic spleens and lymph nodes were dissected and single-cell suspensions obtained using a 70-µm nylon mesh. Blood was collected from the facial vein in EDTA. Erythrocytes from blood and spleen samples were lysed using an NH₄Cl hypotonic solution at a 9:1 NH₄Cl:cell suspension ratio (no. 07850, STEMCELL Technologies), for 10 min and immediately rinsed at least twice and resuspended in HBSS with 5% FBS.

**SP staining.** BM cell suspensions were incubated with the FcBlock reagent (anti-CD16/32, in Rat serum—no. 13551, STEMCELL Technologies), labelled with the biotinylated antibodies anti-CD45R/B220, anti-CD19, anti-CD11b/Bac-1, anti-Ly6G/GR-1, anti-pan-NK, anti-Ly-76/TER and anti-CD3 (full list of antibodies and dyes in Supplementary Table 1), and incubated with Streptavidin microbeads (Miltenyi Biotec), for negative selection of lineage-positive cells by immunomagnetic separation using a MACS column (Miltenyi Biotec). Cells were further incubated in serum-free media (StemSpam SFEM, no. 09650, StemCell Technologies) with 20 µg ml⁻¹ Rhodamine 123 (no. R302, Invitrogen) for 30 min in a water bath at 37 °C (in the dark), cooled, rinsed with HBSS, 5% FBS and then stained with Hoechst 33342 (no. H3570, Invitrogen) at a concentration of 2.5 µg ml⁻¹ diluted in the serum-free media as before, for 90 min at 37 °C, in the dark. Cells were then labelled in HBSS with 5% FBS for 30 min on ice, with Streptavidin-APC and anti-CD45-PE and 1 µg ml⁻¹ Propidium Iodide (no. P3566, Invitrogen). For the long-term reconstituting SP cells gating adjustment, an aliquot of the sample was incubated with 25 mg ml⁻¹ of Verapamil (no. 676777, Calbiochem) before staining with the dyes. Further stainings with different HSC-specific marker combinations were tested, confirming the enrichment of HSCs within the SP gates.

**Cell sorting and *in vivo* single-cell repopulation assay.** Single Lin⁻ Rho⁻ CD45intHoe⁻ long-term reconstituting SP cells of adult female mice in the C57BL/6J background, heterozygous for CD45.1/45.2, *Foxp3*<sup>gfp</sup>/x and *Igha/b*, stained as above described, were sorted using the single-cell deposition unit of a MoFlo (Beckman Coulter) or FACS Aria (Becton Dickinson) fluorescence-activated cell sorter machine into the individual wells of Terasaki plates (no. 452256, MicroWell 60-well MiniTray, Nunc Brand, Thermo Fisher Scientific Inc.) preloaded with 15 µl of serum-free media (StemSpam SFEM, StemCell Technologies). Sorted IgMª and IgMᵇ populations had purities above 95% (for most samples the purity was assessed with flow cytometry immediately after the sorting and for a few samples the purity was assessed with flow cytometry and/or semiquantitative PCR analysis of an allele-specific polymorphism in the VH-D intergenic region). To check for the purity of the single-cell sorting used in the *in vivo* experiments, each well was examined in a 4 °C room using an inverted microscope and the wells containing single cells were used in reconstitutions (Supplementary Fig. 1).

Eight- to twelve-week-old recipient females of the B6.Rag2⁻/⁻ background received sublethal whole-body γ-irradiation with 600 cGy (Gammacell 2000 Mølsgaard Medical), 4–8 h before an intravenous retro-orbital injection with one or more test cells, as indicated. Recipient animals were analysed routinely every 2 weeks for up to 10 weeks post reconstitution for the presence of chimeric cells in the peripheral blood by staining with anti-CD45.1 and anti-CD45.2 antibodies (Supplementary Table 1). Clonally repopulated animals were killed and cells from BM, spleen, thymus and lymph nodes (LNs) were analysed for the presence of terminally differentiated B cells, T cells and myeloid cells of donor origin, as well as progenitor donor cells (antibodies used are described in Supplementary Table 1). Animals selected for subsequent analysis showed in Tregs *Foxp3*<sup>gfp</sup> knock-in gene expression from a single X-chromosome. Splenic CD19<sup>+</sup> IgMª<sup>+</sup> and IgMᵇ<sup>+</sup> cells were then sorted on a FACS Aria instrument, their DNA extracted and the recombination status of the *Igh* gene was analysed as described below.

All flow cytometric data collection was performed on a FACS Calibur (Becton Dickinson), a FACS Aria (Becton Dickinson) or a CyAN ADP (Dako Cytomation, Beckman Coulter). The analysis was performed using the FlowJo software (Tree Star Inc.).

**In vitro B-cell clonal differentiation from precursors.** BM cell suspensions enriched for haematopoietic progenitors (Lin⁻) were treated as described in ref. 34. Briefly, HSC, MPP and CLP cells were sorted using a FACS Aria instrument after staining with anti-CD117(c-Kit), anti-Ly6A/E (Sca-1), anti-CD127(IL7Alpha) and anti-CD135(FLK2; Supplementary Table 1). Single cells were directly sorted into flat-bottom 96-well plates (Techno Plastic Products) previously seeded with non-irradiated 2.5 × 10³ OP-9 cells per well, in complete culture media (OPTI-MEM, Invitrogen), supplemented with 10% FBS (Life Technologies) penicillin (50 U ml⁻¹), streptomycin (50 µg ml⁻¹), β-mercaptoethanol (50 µM) and Gentamycin (50 µg ml⁻¹) and cultured in saturating amounts of the following cytokines: interleukin-7 (IL-7), c-Kit ligand and Flt3 ligand. Supernatants of myeloma cell lines transfected with the DNA encoding the above-mentioned cytokines were kindly provided by Dr Ana Cumano, Institut Pasteur, Paris, France.

After 12–14, 17–20 or 34 days of culture of CLPs, MPPs or HSCs, respectively, individual plates were analysed for single colonies. Wells with single colonies were stained with anti-CD19, anti-IgMª and anti-IgMᵇ (Supplementary Table 1) analysed with flow cytometry and each IgM⁺ allelic fraction (and in some cases also the IgM⁻ fraction) was sorted in bulk using a FACS Aria cell sorter instrument as described above.

**DNA samples.** Genomic DNA (gDNA) from *in vivo*-differentiated clonal CD19⁺IgMª and CD19⁺IgMᵇ pools of at least 10⁴ sorted cells was isolated by proteinase K (Sigma) digestion, phenol-chloroform-extracted, precipitated by ethanol and resuspended in TE buffer. The set of *in vitro*-differentiated sorted samples included also CD19⁺IgMᵇ populations and the numbers of sorted cells were invariably smaller than those of the *ex vivo* samples (at least 2,000 cells); for these *in vitro*-derived samples, gDNA was obtained by incubation with cell digestion medium containing 400 ng µl⁻¹ proteinase K and 17 µM sodium dodecyl sulfate (SDS, Sigma) in PBS (Invitrogen) for 1 h at 50 °C, followed by enzyme heat inactivation for 30 min at 95 °C. gDNA thus obtained was diluted to homogeneous cell-equivalents per µl, as determined by the number of cells sorted, and used directly in PCR reactions, as described below. Treatment of control samples was matched to that of each set of samples. As controls for *in vivo*-differentiated clones, *ex vivo* splenocytes of (C57BL/6*BALB/c)F1 WT animals and mice of the *Rag2*<sup>-/-</sup>, *Jht*<sup>-/-</sup>, *Jht*<sup>+</sup>/*- and *Atm*<sup>-/-</sup> genotype were sorted and treated as the test samples. qDNA from the cells lines 18–81 and PD31 (cell lines kindly provided by Dr Michele Goodhardt, Institut Universitaire d'Hématologie, Paris, France and Dr Patricia Cortes, Icahn Medical Institute, New York, USA) were used as negative controls for the presence of *V_H-D_H* intergenic fragment.

**Real-time quantitative PCR.** The frequency of *DJH* only rearranged alleles in the IgMª and IgMᵇ populations of the same B-cell clones was determined with real-time qPCR, with the test and reference sequences localized in the recombining and non-recombining regions of the *Igh* locus, as shown in Supplementary Fig. 2a. The test primer pairs, VDI4FW_5'-GCAGGAAGACCACAGGAGC-3' and VDI4RE_5'-AGCCACTGTAGATGAGGGGT-3', amplify a ~100-bp fragment herein termed 'V_H-D intergenic fragment' downstream of the last *V_H-gene* segment of the *Igh* locus and upstream of the first *D_H-gene* segment. The reference primer pairs, CmuFW_5'-GTTCTGTGCCTCCGTCTAGC-3' and CmuRE_5'-CTTGGGGTGACTGAGCATTTGC-3', also amplify a ~100-bp fragment located in the µ-constant region of the *Igh* locus (*Cµ*). Each primer pair was included in a different reaction mixture and runs were performed in triplicates or quadruplates using the 7900HT Fast Real-Time PCR System (Applied Biosystems), in 384-well plates with optical tape sealers (Applied Biosystems). Reactions were performed in 10 µl final volume containing 1 × the Power SYBR Green PCR Master Mix (Applied Biosystems) and 0.2 µM of each primer. Thermocycling conditions were 95 °C for 10 min, followed by 45 cycles of 95 °C for 15 s and 60 °C for 60 s. The melting curve of the finished PCR reaction was systematically used to confirm that a single amplicon was amplified.

Calibration curves were performed using six standard dilutions, in triplicates, of B6.Jht⁺⁻ gDNA from *ex vivo*-sorted IgM⁺ cells, ranging from 17 cell equivalents up to 54.240 cell equivalents per reaction (Supplementary Fig. 2b), and test sample concentrations used per run ranged from 200 to 700 cell equivalents per run. No template control reactions were performed using water or a diluted digestion medium solution without template. Beyond replicas, independent repetition of each sample quantification was performed whenever the sample amount was sufficient, in diverse dilutions, but given the low amount of DNA in the majority of samples it was not possible to repeat the PCR on a different day. The SDS software v.2.2 (Applied Biosystems) was used to calculate the quantification cycle (Cq) value. Real-time quantification data were averaged between replicates and relative amount of DNA was calculated with the *Cµ* constant fragment quantification as the reference, using the delta Cq method adjusted for PCR efficiency⁵⁴.

**Igh rearrangement PCR and sequencing.** For the amplification of *DJH* and *V_HDJH* rearrangements, we used the primer collection and the two-step nested PCR conditions previously described¹⁸. One thousand cell-equivalents per run were used for the first round multiplex PCR of *V_HDJH*-rearrangement amplifications in a 100-µl final master mix volume. For the nested round of PCR,
alternative reverse primers were *JH*1-1, *JH*2- and *JH*3-specific (JH1A-5′-TGTGGCA GATGCCTGACATGGGG-3′, JH2A-5′-TCAGCCAGGGCTCCCAATGACC C-3′, and JH3A-5′-CAGCAGGCAGAGAGTCCCTGACCC-3′), allowing allotype discrimination upon sequencing (Supplementary Fig. 3). For sequencing, PCR bands were cloned in the pGEM-T Easy VECTOR SYSTEM (Promega) and a number of clones were analysed. Typically, four clones were sequenced per band. Given both the high number of PCR cycles used and the cloning step before sequencing the *V*<sub>H</sub>*D*<sub>H</sub> amplicons, it was necessary to evaluate the impact in the analysis of PCR/sequencing-introduced errors. This was estimated using the formula (f<sub>mut</sub> × *P*<sub>stop</sub> + *f*<sub>indel</sub> − *f*<sub>mut</sub> × *f*<sub>stop</sub> × *f*<sub>indel</sub>) × L/1,000, where *f*<sub>mut</sub> and *f*<sub>indel</sub> stand for the number of PCR/sequencing-introduced point mutations (2.6) and the number of PCR/sequencing-introduced gaps/insertions (0.18) per 1,000 nucleotides, respectively, estimated from alignments of the germline regions (a total of 16,020 bp) downstream of the *JH* elements that we have sequenced, *p*<sub>stop</sub> is the probability of producing a stop codon from a PCR-introduced point mutation (0.071) taking in consideration the codon usage of the *V*<sub>H</sub>*D*<sub>H</sub> sequences obtained in this study and *L* is the average size of the amplicons. If the entire *V*<sub>H</sub> and CDR3 regions (*L* = 269 bp) of the amplicon are considered, the percentage of productive *V*<sub>H</sub>*D*<sub>H</sub> rearrangements that become non-productive because of PCR errors is 8.9%. This percentage decreases to 2.0% if only a sequence region including the entire CDR3 is used (*L* = 60 bp), which is feasible since the PCR errors in the *V*<sub>H</sub> region can be aligned with the known *V*<sub>H</sub> germline sequences from antigen receptor databases. In any case, the errors introduced by PCR do not affect our conclusions.

## References

1. Cebra, J. J., Colberg, J. E. & Dray, S. Rabbit lymphoid cells differentiated with respect to alpha-, gamma-, and mu- heavy polypeptide chains and to allotypic markers Aal1 and Aa2. *J. Exp. Med.* **123**, 547–558 (1966).

2. Burnet, F. M. A modification of Jerne's theory of antibody production using the concept of clonal selection. *Aust. J. Sci.* **20**, 67–68 (1957).

3. Hewitt, S. L. RAG-1 and ATM coordinate monoallelic recombination and nuclear positioning of immunoglobulin loci. *Nat. Immunol.* **10**, 655–664 (2009).

4. Tonegawa, S. Somatic generation of antibody diversity. *Nature* **302**, 575–581 (1983).

5. Coleclough, C., Perry, R. P., Karjalainen, K. & Weigert, M. Abrant rearrangements contribute significantly to the allelic exclusion of immunoglobulin gene expression. *Nature* **290**, 372–378 (1981).

6. Barreto, V. & Cumano, A. Frequency and characterization of phenotypic Ig heavy chain allelically included IgM-expressing B cells in mice. *J. Immunol.* **164**, 893–899 (2000).

7. Alt, F. W. *et al.* Ordered rearrangement of immunoglobulin heavy chain variable region segments. *EMBO J.* **3**, 1209–1219 (1984).

8. Selimyan, R. *et al.* Localized DNA demethylation at recombination intermediates during immunoglobulin heavy chain gene assembly. *PLoS Biol.* **11**, e1001475–e1001475 (2013).

9. Rose, S. M. S., Kuehl, W. M. W. & Smith, G. P. Cloned MPC 11 myeloma cells express two kappa genes: a gene for a complete light chain and a gene for a constant region polypeptide. *Cell* **12**, 453–462 (1977).

10. Alt, F. F., Rosenberg, N. N., Lewis, S. S., Thomas, E. E. & Baltimore, D. D. Organization and reorganization of immunoglobulin genes in A-MULV-transformed cells: rearrangement of heavy but not light chain genes. *Cell* **27**, 381–390 (1981).

11. Weaver, D., Costantini, F., Imanishi-Kari, T. & Baltimore, D. A transgenic immunoglobulin mu gene prevents rearrangement of endogenous genes. *Cell* **42**, 117–127 (1985).

12. Nussenzweig, M. C. *et al.* Allelic exclusion in transgenic mice that express the membrane form of immunoglobulin mu. *Science* **236**, 816–819 (1987).

13. Kitamura, D. & Rajewsky, K. Targeted disruption of mu chain membrane exon causes loss of heavy-chain allelic exclusion. *Nature* **356**, 154–156 (1992).

14. Boekel, E. T., Melchers, F. & Rolink, A. G. Precursor B cells showing H chain allelic inclusion display allelic exclusion at the level of pre-B cell receptor surface expression. *Immunity* **8**, 199–207 (1998).

15. Löffert, D., Ehlich, A., Muller, W. & Rajewsky, K. Surrogate light chain expression is required to establish immunoglobulin heavy chain allelic exclusion during early B cell development. *Immunity* **4**, 133–144 (1996).

16. Hauser, J., Grundström, C. & Grundström, T. Allelic exclusion of IgH through inhibition of E2A in a VDJ recombination complex. *J. Immunol.* **192**, 2460–2470 (2014).

17. Lutz, J. *et al.* Pro-B cells sense productive immunoglobulin heavy chain rearrangement irrespective of polypeptide production. *Proc. Natl Acad. Sci. USA* **108**, 10644–10649 (2011).

18. Ehlich, A., Martin, V., Muller, W. & Rajewsky, K. Analysis of the B-cell progenitor compartment at the level of single cells. *Curr. Biol.* **4**, 573–583 (1994).

19. Kosak, S. T. *et al.* Subnuclear compartmentalization of immunoglobulin loci during lymphocyte development. *Science* **296**, 158–162 (2002).

20. Sayegh, C. E., Sayegh, C., Jhunjhunwala, S., Riblet, R. & Murre, C. Visualization of looping involving the immunoglobulin heavy-chain locus in developing B cells. *Genes Dev.* **19**, 322–327 (2005).

21. Roldán, E. *et al.* Locus 'decontraction' and centromeric recruitment contribute to allelic exclusion of the immunoglobulin heavy-chain gene. *Nat. Immunol.* **6**, 31–41 (2005).

22. Guo, C. C. *et al.* CTCF-binding elements mediate control of V(D)J recombination. *Nature* **477**, 424–430 (2011).

23. Fuxa, M. *et al.* Pax5 induces V-to-DJ rearrangements and locus contraction of the immunoglobulin heavy-chain gene. *Genes Dev.* **18**, 411–422 (2004).

24. Skok, J. A. *et al.* Nonequivalent nuclear location of immunoglobulin alleles in B lymphocytes. *Nat. Immunol.* **2**, 848–854 (2001).

25. Holwerda, S. J. B. *et al.* Allelic exclusion of the immunoglobulin heavy chain locus is independent of its nuclear localization in mature B cells. *Nucleic Acids Res.* **41**, 6905–6916 (2013).

26. Chess, A. Random and non-random monoallelic expression. *Neuropsychopharmacology* **38**, 55–61 (2012).

27. Singh, N. *et al.* Coordination of the random asynchronous replication of autosomal loci. *Nat. Genet.* **33**, 339–341 (2003).

28. Mostoslavsky, R. *et al.* Asynchronous replication and allelic exclusion in the immune system. *Nature* **414**, 221–225 (2001).

29. Pereira, J. P., Girard, R., Chaby, R., Cumano, A. & Vieira, P. Monoallelic expression of the murine gene encoding Toll-like receptor 4. *Nat. Immunol.* **4**, 464–470 (2003).

30. Osawa, M., Hanada, K., Hamada, H. & Nakauchi, H. Long-term lymphohematopoietic reconstitution by a single CD34-low/negative hematopoietic stem cell. *Science* **273**, 242–245 (1996).

31. Fontenot, J. D. *et al.* Regulatory T cell lineage specification by the forkhead transcription factor foxp3. *Immunity* **22**, 329–341 (2005).

32. Vettermann, C. & Schlissel, M. S. Allelic exclusion of immunoglobulin genes: models and mechanisms. *Immunol. Rev.* **237**, 22–42 (2010).

33. Gu, H. H., Zou, Y. R. Y. & Rajewsky, K. K. Independent control of immunoglobulin switch recombination at individual switch regions evidenced through Cre-loxP-mediated gene targeting. *Cell* **73**, 1155–1164 (1993).

34. Vieira, P. & Cumano, A. Differentiation of B lymphocytes from hematopoietic stem cells. *Methods Mol. Biol.* **271**, 67–76 (2004).

35. Inlay, M. A. *et al.* Ly6d marks the earliest stage of B-cell specification and identifies the forkpoint between B-cell and T-cell development. *Genes Dev.* **23**, 2376–2381 (2009).

36. Decker, D. J., Boyle, N. E., Koziol, J. A. & Klinman, N. R. The expression of the Ig H chain repertoire in developing bone marrow B lineage cells. *J. Immunol.* **146**, 350–361 (1991).

37. Hess, J. *et al.* Induction of pre-B cell proliferation after de novo synthesis of the pre-B cell receptor. *Proc. Natl Acad. Sci. USA* **98**, 1745–1750 (2001).

38. Wang, Y.-H. *et al.* V(D)J recombinatorial repertoire diversification during intraclonal pro-B to B-cell differentiation. *Blood* **101**, 1030–1037 (2003).

39. Gu, H., Kitamura, D. & Rajewsky, K. B cell development regulated by gene rearrangement: arrest of maturation by membrane-bound Dμ protein and selection of DHElement reading frames. *Cell* **65**, 47–54 (1991).

40. Allen, E. *et al.* High concentrations of long interspersed nuclear element sequence distinguish monoallelically expressed genes. *Proc. Natl Acad. Sci. USA* **100**, 9940–9945 (2003).

41. Chow, J. C. *et al.* LINE-1 activity in facultative heterochromatin formation during X chromosome inactivation. *Cell* **141**, 956–969 (2010).

42. Kitsberg, D. *et al.* Allele-specific replication timing of imprinted gene regions. *Nature* **364**, 459–463 (1993).

43. Chess, A., Simon, I., Cedar, H. & Axel, R. Allelic inactivation regulates olfactory receptor gene expression. *Cell* **78**, 823–834 (1994).

44. Takagi, N. Differentiation of X chromosomes in early female mouse embryos. *Exp. Cell Res.* **86**, 127–135 (1974).

45. Lebon, J. M., Tam, P. P., Singer-Sam, J., Riggs, A. D. & Tan, S. S. Mouse endogenous X-linked genes do not show lineage-specific delayed inactivation during development. *Genet. Res.* **65**, 223–227 (1995).

46. Farago, M. *et al.* Clonal allelic predetermination of immunoglobulin-κ rearrangement. *Nature* **490**, 561–565 (2012).

47. Keller, G., Paige, C., Gilboa, E. & Wagner, E. F. Expression of a foreign gene in myeloid and lymphoid cells derived from multipotent haematopoietic precursors. *Nature* **318**, 149–154 (1985).

48. Mostoslavsky, R. & Alt, F. W. The lingering enigma of the allelic exclusion mechanism. *Cell* **118**, 539–544 (2004).

49. Gribnau, J. X chromosome choice occurs independently of asynchronous replication timing. *J. Cell Biol.* **168**, 365–373 (2005).

50. Yamagami, T., Boekel, ten, E., Andersson, J., Rolink, A. & Melchers, F. Frequencies of multiple IgL chain gene rearrangements in single normal or kappaL chain-deficient B lineage cells. *Immunity* **11**, 317–327 (1999).
51. Schlimgen, R. J., Reddy, K. L., Singh, H. & Krangel, M. S. Initiation of allelic exclusion by stochastic interaction of Tcrb alleles with repressive nuclear compartments. *Nat. Immunol.* **9**, 802–809 (2008).

52. Casellas, R. R. et al. Contribution of receptor editing to the antibody repertoire. *Science* **291**, 1541–1544 (2001).

53. Elson, A. et al. Pleiotropic defects in ataxia-telangiectasia protein-deficient mice. *Proc. Natl Acad. Sci. USA* **93**, 13084–13089 (1996).

54. Pfaffl, M. W. A new mathematical model for relative quantification in real-time RT-PCR. *Nucleic Acids Res.* **29**, e45 (2001).

55. Hesse, J. E., Lieber, M. R., Gellert, M. & Mizuuchi, K. Extrachromosomal DNA substrates in pre-B cells undergo inversion or deletion at immunoglobulin V-(D)-J joining signals. *Cell* **49**, 775–783 (1987).

56. Alamyar, E., Duroux, P., Lefranc, M. P. & Giudicelli, V. IMGT(®) tools for the nucleotide analysis of immunoglobulin (IG) and T cell receptor (TR) V-(D)-J repertoires, polymorphisms, and IG mutations: IMGT/V-QUEST and IMGT/HighV-QUEST for NGS. *Methods Mol. Biol.* **882**, 569–604 (2012).

## Acknowledgements

We thank Marie Bonnet, Thiago Carvalho, Ana Cumano, Revati Masilamani, Gillian Wu and all the members of our laboratory for the critical reading of the manuscript. C.F.A.P. is supported by an FCT fellowship (SFRH/BD/33217/2007). This work was supported by an FCT grant (PTDC/BIA-BCM/108020/2008).

## Author contributions

C.F.A.P. and V.M.B. designed the experiments, analysed the data, discussed the results and wrote the manuscript. C.F.A.P. performed all the experiments. R.de.F and F.M. provided technical assistance. T.L. and R.G. gave support concerning flow cytometry. P.V. advised and gave support to all experiments. V.M.B. supervised all the work.

## Additional information

**Supplementary Information** accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Pereira, C. F. A. et al. Independent recruitment of Igh alleles in V(D)J recombination. *Nat. Commun.* **5**:5623 doi: 10.1038/ncomms6623 (2014).

<figure><img src="image_11.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/