Lentiviral and targeted cellular barcoding reveals ongoing clonal dynamics of cell lines in vitro and in vivo
© Porter et al.; licensee BioMed Central Ltd. 2014
Received: 22 July 2013
Accepted: 30 May 2014
Published: 30 May 2014
Cell lines are often regarded as clonal, even though this simplifies what is known about mutagenesis, transformation and other processes that destabilize them over time. Monitoring these clonal dynamics is important for multiple areas of biomedical research, including stem cell and cancer biology. Tracking the contributions of individual cells to large populations, however, has been constrained by limitations in sensitivity and complexity.
We utilize cellular barcoding methods to simultaneously track the clonal contributions of tens of thousands of cells. We demonstrate that even with optimal culturing conditions, common cell lines including HeLa, K562 and HEK-293 T exhibit ongoing clonal dynamics. Starting a population with a single clone diminishes but does not eradicate this phenomenon. Next, we compare lentiviral and zinc-finger nuclease barcode insertion approaches, finding that the zinc-finger nuclease protocol surprisingly results in reduced clonal diversity. We also document the expected reduction in clonal complexity when cells are challenged with genotoxic stress. Finally, we demonstrate that xenografts maintain clonal diversity to a greater extent than in vitro culturing of the human non-small-cell lung cancer cell line HCC827.
We demonstrate the feasibility of tracking and quantifying the clonal dynamics of entire cell populations within multiple cultured cell lines. Our results suggest that cell heterogeneity should be considered in the design and interpretation of in vitro culture experiments. Aside from clonal cell lines, we propose that cellular barcoding could prove valuable in modeling the clonal behavior of heterogeneous cell populations over time, including tumor populations treated with chemotherapeutic agents.
Even under ideal growth conditions, cultured cells exhibit genetic heterogeneity. It is therefore valuable, although technically challenging, to track the behavior and interplay of clones within a cellular population. Furthermore, clonal dynamics play important roles in cancer and stem cell biology. We therefore aimed to develop a sensitive and quantitative method for tracking the clonal dynamics within populations of cells with minimal disruption to both individual cells and the population as a whole.
Early techniques, able to track one or a few clones, relied upon gross chromosomal markers [1, 2], heterozygous alleles [3, 4], or a rainbow of fluorescent markers . More recent methods have utilized viral integration to confer specific and theoretically unique heritable marks on a cell [6–9]. While these techniques greatly increase the number of clones that can be detected, the method is plagued by limitations in sensitivity and an inability to accurately measure the size of each clone, despite advances in detection [10–12]. To overcome these limitations, we decided to label cells with unique DNA barcodes, which can be recovered and sequenced to reveal the temporal and quantitative behavior of entire cell populations and also individual member clones.
The ability to track a limited subset of a cellular population with DNA barcodes has previously been demonstrated by several groups [13–17]. Here, we demonstrate the feasibility of monitoring entire cell populations using a barcode system that scales to many thousands or even a million individual clones. We also outline a novel non-viral barcoding method that targets barcodes to a single genomic locus through zinc-finger nuclease (ZFN)-induced homologous recombination and therefore avoids unpredictable viral insertional mutagenesis. With this more precise and scalable approach we are able to define the dynamics of an entire cell population rather than tracing the fates of only a few representative clones.
First, we validate the performance of our barcode method by tracking the in vitro dynamics of several common cell lines. We find that despite years in culture, common cell lines exhibit ongoing clonal instability. Next, we compare the clonal dynamics of cell populations barcoded by random insertion of a lentiviral vector versus targeted integration at a single genomic locus through homologous recombination and find that the nuclease-mediated insertion of the barcode sequence process itself results in dramatic changes in clonal representation. Finally, we measure the contributions of clones in primary xenograft tumors. By comparing the dynamics of the same population of clones in vitro and in vivo, we were able to show that the selective pressure that restricts clonal diversity is greater in culture than in a mouse xenograft. These findings add to our knowledge of in vitro and in vivo cellular behavior, and have important implications for the design and interpretation of experiments utilizing cultured cells.
Library validation and data analysis pipeline
To determine the complexity and distribution of the barcode library, as well as to determine the extent of error and bias introduced by sample preparation and sequencing, we independently PCR-amplified the plasmid barcode library for sequencing four separate times, and sequenced each amplified sample on an independent flow cell lane at a coverage of 400-fold.
The distribution of barcodes between the four replicates was nearly identical as well. Sequences were counted, and then the frequency of each barcode was calculated as a percentage of the whole. The mean percent frequency of each replicate was very similar to the expected for a library this size (Figure 2d). The median barcode frequencies of the four replicates were also very similar to one another, spanning 0.0066% to 0.0068% with a low standard deviation (expected median frequency in an unskewed population = 0.0068%) (Figure 2d). By comparing the frequencies of each barcode in each of the sample replicates, we were able to determine R2 values, which ranged from 0.989 to 0.996 (Additional file 1). From this, we were able to conclude that our method of PCR amplification, sequencing and analysis is highly reproducible and does not introduce significant error or bias.Our measure and quantification of bias within the replicate barcode library sequences are shown in Figure 2c-f. Figure 2c shows a histogram of barcode frequency distribution in this library across all four sequencing replicates. A completely normal distribution would result in a bell shaped curve. Figure 2e plots the percentage of barcodes against the percentage of sequences and an unbiased distribution would result in a 45-degree line (dotted line). In both of these figures the slight skewing of the original plasmid library is demonstrated by the deviation from a bell shaped curve in 2c and the deviation from the 45-degree line in 2e. We quantify the bias in Figure 2f by plotting the percentage of sequences that were accounted for by 10, 25, 50, and 75% of the most abundant barcodes. In the original plasmid library the top 10, 25, 50, and 75% most abundant barcodes account for approximately 27, 50, 77, and 93% of the sequences, respectively, thus providing a quantitative metric of bias in barcode representation. This slight skewing in the plasmid barcode library is most likely the result of its amplification through overnight growth in bacteria as part of its preparation.
Cellular barcode libraries and passaging experiments
K562 cellular barcode library passaging and results
For our first cellular barcode library and passaging experiments, we chose K562 cells, a common human leukemia cell line established in 1970 from a patient with chronic myelogenous leukemia . We found that in all three biological replicates, the number of barcodes detected in each population decreased over time and the clonal distribution within each population became more biased over time as the tails become larger than would be expected from a normal distribution (Figure 3c-g, Additional file 2). Each of the replicates also contained clones that came to constitute greater than 1% of the total population ('major clones'), but all clones constituted less than 10% of their respective population at PD 90. At each time point, the clones identified were categorized as rare (less than 0.0007% of the population), abundant (greater than 0.5% of the population), or average (all others) based on their individual contribution to the total number of cells in culture (Additional file 3).
In order to determine whether the clonal dynamics within the three populations were due to pre-existing cell-intrinsic factors, or if the populations underwent clonal selection after the split, we compared the identities of the major clones in each replicate. One clone (Figure 3g, yellow) was found in all three populations as a major clone, suggesting that factors intrinsic to this cell at the time it was marked caused its progeny to have a growth advantage over its neighbors. However, most of the other major clones within each replicate were unique to that population, suggesting that each clone’s growth advantage was gained after the clone was marked and the biological replicates had been separated, indicating ongoing clonal variation followed by selection during the course of the experiment. As the population doubling increased, the most abundant clones contributed to a larger and larger portion of the total population (Figure 3d,e). For example, at PD 0 the 10% most abundant clones accounted for 29% of the total cells in the culture, but by PD 90 the top 10% now accounted for almost 75% of the total cells in the population (Figure 3e). Importantly, the 10% most abundant clones at PD 90 were not the same as the top 10% at PD 0. Furthermore, the dominant clones identified at PD 90 were derived from clones in all three percentage contribution categories (rare, average and/or abundant) at PD 0 in all three biological replicate populations (Additional file 4). The distribution of clones widened, with greater percentages of clones showing up in the highest and lowest bins, indicating an increasing trend in high and low frequency clones (Figure 3c). Thus, these experiments demonstrate that K562 cells continue to display rapid clonal dynamics even under optimal culturing and passaging conditions.
K562 clonal cellular barcode library passaging and results
Since we observed ongoing clonal dynamics in our polyclonal K562 population, we hypothesized that this marked population of cells had developed significant heterogeneity over time from ongoing genetic and epigenetic changes that affected clonal fitness and dynamics. To test this hypothesis, we created a K562 line derived from a single cell, and repeated the barcoding experiment (as with the original K562 population). We found that although the rate of clone loss and diversification was slower, it still occurred (Additional files 5 and 6). There appears to be more overlap among the largest clones of the three biological replicates than seen with the polyclonal K562 cellular library, as well as a number of clones unique to each biological replicate, indicating ongoing clonal evolution). The slower but persistent changes observed in the population derived from a single cell are highlighted by the difference in percentage contribution of the top 10% most abundant clones identified. In the clonal K562 experiment, the top 10% of clones identified accounted for 32% of the population at PD 0 and 38% of the population at PD 90. This increase is dramatically less than that observed in the polyclonal K562 experiment wherein the top 10% most abundant barcodes accounted for 29% of the total sequences at PD 0 and 75% of the total sequences at PD 90 (compare Figure 3e with Additional file 5c).
Targeted barcode library in K562 cells
Clonal dynamics of HeLa and HEK-293 T-cell lines
In order to determine whether our findings of persistent and ongoing clonal dynamics in K562 cells were representative of other cell types, we marked and tracked the clonal dynamics of both the HeLa and HEK-293 T-cell lines. We created both cellular barcode libraries from the same lentiviral prep used in the K562 cell experiments, and passaged them in an identical manner. The results show that while relatively few clones were lost over 90 population doublings, we did see some skewing of the distribution of clones over time as well as development of major clones (Additional files 10, 11, 12 and 13). As with the original K562 experiments, we saw only a small number of major clones that recurred in different biological replicates, and a number of major clones that were unique to a single population. These results indicate that the HeLa and HEK-293 T-cell lines, as with K562 cells, show significant clonal dynamics even under ideal culture conditions.
The number of clones that contribute to 3T3 cell lines derived from mouse embryonic fibroblasts
The barcode system we describe here is applicable to a large number of biological questions, including quantifying the number and distribution of cells that contribute to downstream populations. To demonstrate this, we passaged barcode marked mouse embryonic fibroblasts in a 3T3 experiment  and found that a minimum of 0.7% of the fibroblasts transformed and contributed to the 3T3 population (data not shown).
Using the barcode marking system to compare clonal dynamics in vitroversus xenografts
We hypothesized that insertional mutagenesis caused by the barcode integration may have played a role in the growth advantage seen in this clone. We mapped the barcode insertion site of this clone to the second intron of MSRB3, a methionine sulfide reductase on human chromosome 12. Karyotype analysis of the HCC827 cells used in these experiments show the presence of three copies of chromosome 12. We therefore believe it unlikely that the integration of the barcode in this clone is the causal factor in its distinct advantage over the other clones in the population because this gene has no reported role in tumor cell proliferation and would not disrupt the coding region of the gene.
Quantifying changes in clonal representation using the Shannon-Weaver diversity index
The Shannon-Weaver diversity index is a powerful quantitative measure that accounts for both the number of different elements (in our case, cellular clones) and the relative representation of each element within the population (in our case, the relative abundance of each clone). It is broadly used in the ecology literature but applies very well to studies of clonal dynamics [29, 30]. In the Shannon-Weaver diversity index, a higher number shows that the population is more diverse and evenly represented while a lower number demonstrates a more restricted and more unequal population. In all of our experiments, the Shannon-Weaver diversity index decreased, usually quite dramatically over time (Additional file 16).
We have developed a system that genetically marks individual cells, allowing for the simple, simultaneous, and quantitative tracking of thousands of cells using a combination of barcode marking and high-throughput sequencing. In establishing and validating this method we have focused on a system in which we can track >12,000 different clones simultaneously, but have also extended this to develop barcode libraries of varying complexities, including libraries that consist of over one million different barcodes (data not shown). Just as with the >12,000 complexity barcode library, we confirmed the complexity of these larger libraries by sequencing and they are now being used in other work to study the dynamics of hematopoietic stem cell reconstitution in non-human primates. With these larger libraries even greater care must be taken at each step (the creation of lentivirus, the marking of cells, and so on) to maintain the complexity and not create bottlenecks that would artificially skew the clonal barcode representation.
We have validated the power of this system by using it to track the spontaneous clonal dynamics of cell lines grown in culture and after xenotransplantation. We found that cell lines that have been grown in culture for years continue to have ongoing, dramatic clonal fluctuations, with many clones disappearing over time and others outgrowing their neighbors. When we started the analysis from a single cell, these effects were reduced but not abrogated, demonstrating that the process is ongoing and not static even within a group of closely related cells. These results quantify and confirm what is understood, but often underappreciated, namely that the 'same' cell line in different laboratories or at different times in the same laboratory may contain markedly different clonal repertoires and thus may display strikingly different behavior. We designed our studies to put the populations under as little stress as possible and yet we still observed dynamic ongoing population changes. It is likely that under situations in which populations are stressed, even more distinct clonal dynamics will be observed. Moreover, this system should allow the study of diverse clonal repertoires of populations resulting from ongoing clonal dynamics when cells are placed under different environmental pressures. One great use of this system would be to study the clonal dynamics of chemotherapeutic resistance when cancer cells are exposed to different agents (both classical and 'targeted' therapies). Our studies were not designed to investigate the mechanism of the ongoing dynamics but the rapid pace of change suggests that epigenetic alterations may be an important aspect of the process, in addition to the background genetic changes that occur in all cells over time.
In demonstrating the utility of this powerful clone tracking strategy in multiple systems, two of our findings contradict general consensus in their respective fields. The first major surprise was that the degree of clonal skewing (more clonal dropout, higher proportion of clones becoming dominant within the population) was greater when we used ZFNs to target the barcode to a single 'safe harbor' site in the genome, compared with using lentiviral vectors to semi-randomly integrate the barcode into uncontrolled sites in the genome. This finding suggests that other cellular characteristics, such as pre-determined genetic and epigenetic factors, must be considered in any gene editing strategy. This finding also highlights this strategy as a sensitive method to potentially interrogate factors that might affect both genetic and epigenetic stability. The second surprise was finding that there is less clonal selection pressure when passaging human cancer cells through a xenograft mouse model than when culturing in a plastic dish. This finding highlights how this method can be applied to monitor the effects of both cell intrinsic and cell extrinsic factors in clonal diversity in complex biological systems. Thus, barcode marking is a useful system to interrogate different methods for their ability to maintain the complexity of a population of cells over time. Moreover, the results support the continued use of mouse xenografting as a method to maintain the clonal complexity of human tumors and thus may represent a better model system for studying cancer than in vitro cell culture. For example, this system could be used to study the clonal dynamics of primary tumor development, local tumor recurrence following resection, metastasis and the clonal dynamics of chemotherapy resistance.
Finally, there are several other potential important future uses of our cellular barcoding system. One use is as a control for any small interfering RNA or small hairpin RNA screens to establish the background rate of false positives created by clonal heterogeneity and dynamics of the cell population being challenged. A second use would be to use the system to better understand the dynamics of cancer, including the frequency of tumor initiating/propagating cells in a population and the clonality of recurrence and metastasis. While the power of barcoding has already been shown by studying reconstitution of the hematopoietic system after transplantation, our system is designed for scale and enables the tracking of a larger number of cells. We can track the efficiency and clonality of hematopoietic engraftment in larger, more clinically relevant model systems over time. Moreover, the system described here can be used to evaluate different interventions that might affect either of these factors in larger model systems. Lastly, this barcoding system could also be used to study the clonal dynamics of somatic cell re-programming in induced pluripotent stem cells.
We have developed a simple, cost-effective method to simultaneously track the fate of tens of thousands of cells using a DNA molecular barcode marking system. We believe this new method will be of broad use in a range of scientific fields, some of which we have outlined here, and is likely to lend insights into clonal population dynamics that are not possible using other available methods.
Materials and methods
The Illumina P5 adapter sequence was cloned upstream of UBC-GFP in pLGR7, a third generation lentiviral vector with self-inactivating long terminal repeats. Barcode oligos were synthesized by Integrated DNA Technologies (Coralville, IA, USA) and contained 20 randomized bases flanked by known anchor sequences and enzyme restriction sites. Approximately 1012 barcode oligos were annealed to a short (20 bp) oligo complementary to the 3’ end, followed by extension by Phusion polymerase to create double-stranded barcode fragments. These fragments were then cloned immediately downstream of the P5 adapter sequence within the vector. Pooled ligations were electroporated into Escherichia coli and expanded in SOC medium for 16 hours at 37°C with shaking. We plated 0.05% of this transformation on an antibiotic resistance plate and colonies were counted after overnight growth. The 1,063 colonies suggest the complexity of the remaining library to be approximately 20,000. Plasmid libraries were harvested and purified with an EndoFree Plasmid Maxi Kit (QIAGEN, Valencia, CA, USA).
Cell lines and cell culture
K562 and HCC827 cells were obtained from ATCC (Manassas, VA, USA) and maintained in RPMI 1640 (Hyclone, Logan, UT, USA) supplemented with 10% Bovine Growth Serum (Hyclone, Logan, UT, USA), 2 mM L-gluatmine, 100 units/ml penicillin, and 100 μg/ml streptomycin (Mediatech, Manassas, VA, USA). HEK-293T cells were obtained from ATCC and HeLa cells were a generous gift from Dr Alejandro Sweet-Cordero (Stanford University). Both were maintained in Dulbeco’s modified Eagle’s medium (Hyclone) supplemented with 10% Bovine Growth Serum (Hyclone), 2 mM L-gluatmine, 100 units/ml penicillin, and 100 μg/ml streptomycin (Mediatech).
Barcode lentivirus was produced in HEK-293 T cells by calcium phosphate transfection of 10 μg barcode plasmid vector along with 3 μg VSVG, 5 μg RRE, and 2.5 μg RSV/REV helper plasmids per 10 cm plate. At 18 hours post-transfection, fresh media with 4 mM caffeine was added to cells . Virus-containing supernatant was collected 24 hours later and concentrated by ultracentrifugation, aliquoted, and stored at -80°C.
Cellular barcode libraries and passaging experiments
For lentiviral barcoding, 2 × 106 cells were infected at low MOIs such that 5 to 10% of cells expressed GFP, in order to reduce the chances of any single cell being marked by more than one barcode. GFP-expressing cells were sorted from each population 4 days after transduction. These cells were expanded in culture for an additional 4 to 5 days, resulting in the PD 0 population. Aliquots of 3 × 105 PD 0 cells were used to start the parallel biological replicate populations A, B, and C for each cell line or condition. Additional PD 0 cells were harvested for genomic DNA for barcode analysis, and remaining cells were frozen and stored for future use.
Cells were passaged every 3 days by trypsinization if necessary for attached cells, mixing well and taking 3 × 105 live cells to a new well with fresh medium. Cells were monitored by flow cytometry for GFP expression as well as cell counts. At PD 30, 60, and 90, cells were harvested for genomic DNA as outlined below.
Targeted barcode insertion
ZFNs designed to target the human CCR5 locus have been described previously . The barcode targeting vector was created by cloning barcode oligos and the P5 Illumina adapter, similar to the lentiviral barcode library, downstream of the UBC-GFP within the homology arms of a CCR5 targeting vector described previously . Included outside the homology arms was an HSV-TK domain for negative selection.
K562 cells were nucleofected (Lonza, Basel, Switzerland) with 10 μg targeting vector and 1 μg of each ZFN plasmid using program T-016 and nucleofection buffer containing 100 mM KH2PO4, 15 mM NaHCO3, 12 mM MgCl2•6 H20, 8 mM adenosine 5′-triphosphate, 2 mM glucose, pH 7.4. Negative selection was performed with two pulses of 5 μM ganciclovir at days 5 and 11 post-transfection.
Eight-week-old female NU/NU athymic mice were obtained from Charles River Laboratories (Wilmington, MA, USA). 2 × 105 barcode-marked HCC827 cells in 200 μl phosphate-buffered saline were injected sub-cutaneously into the right flank of each mouse. Mice were maintained at Stanford University’s Research Animal Facility and monitored weekly until tumors were palpable and twice weekly thereafter. Animal procedures were approved by the International Animal Care and Use Committee.
DNA barcode extraction and sequencing
Genomic DNA from both cultured cells and tumors was harvested with the DNeasy Blood and Tissue Kit (QIAGEN). Barcodes were amplified by PCR using Phusion high-fidelity DNA polymerase (New England Biolabs, Beverly, MA, USA) with primers containing the adapters necessary for Illumina sequencing. Forward primers also contained 4 bp indexing tags to allow sample multiplexing. PCR bands were size-selected on 1.1% agarose gel and purified with a QIAquick Gel Extraction Kit (QIAGEN). DNA concentrations were checked by BioAnalyzer and high-throughput sequencing was performed using the Illumina GA II sequencer by the Stanford Functional Genomics Facility.
Barcode sequencing data analysis
We developed a computational tool to enable rapid readout of barcodes from raw Illumina FASTQ data. All methods and documentation are available via Github at . The software is open source and can be modified as needed to suit additional applications. The Illumina FASTQ data from these runs are available in the NIH Sequence Read Archive as project [SRA:SRP029299].
Monte Carlo simulation of the cell culture passaging experiments was performed and showed that, for experiments in which 3 × 105 cells were maintained at each passage, approximately 80 to 84% of clones should be present at PD 90 if the growth rate of all clones is equal. For experiments in which 2 × 106 cells were passaged, 100% of clones would be expected to persist at PD 90. Additional details of the simulations are available via Github at , including the histogram plots of various simulations. Briefly, we started the simulation with 300,000 cells barcoded with one of 14,000 unique identifiers. The cells were assumed to be evenly divided between the barcode identifiers.
After this, we perform the following steps repeatedly. First, the cell population is grown. The number of cells of each population is increased to reflect exponential growth of the population. We use a rate consistent with a doubling every 19 hours (this rate was determined experimentally by measuring the population doubling time of our K562 cells); this yields approximately 13-fold growth across 3 days (1 passage). Second, after 3 days, we randomly select 300,000 cells from the population (should be around 4,000,000 after 3 days of growth). Third, we report the number of unique barcodes present in the 300,000 selected cells.
The above is simulated for 25 passages (or 75 days or approximately 94 population doublings).
enhanced green fluorescent protein
multiplicity of infection
polymerase chain reaction
We thank Dr Jennifer Adair for her helpful comments on the work. We thank Dr Woodring Wright for helpful discussions in the developmental phase of the project. We thank Dr Stacey Wirt for her assistance with the mouse xenograft studies. We thank Gareth Highnam for his assistance with data analysis. We thank the Stanford Functional Genomics Facility headed by Dr John Coller for their help with the sequencing aspects of the project. MP was supported by the Laurie Krauss Lacob Faculty Scholar Award and we acknowledge her generous support. We thank the Alex’s Lemonade Stand Foundation for their support of SNP and this work through an Innovation Award to MHP.
- Ford CE, Hamerton JL, Barnes DW, Loutit JF: Cytological identification of radiation-chimaeras. Nature. 1956, 177: 452-454. 10.1038/177452a0.PubMedView ArticleGoogle Scholar
- Wu AM, Till JE, Siminovitch L, McCulloch EA: Cytological evidence for a relationship between normal hemotopoietic colony-forming cells and cells of the lymphoid system. J Exp Med. 1968, 127: 455-464. 10.1084/jem.127.3.455.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith LG, Weissman IL, Heimfeld S: Clonal analysis of hematopoietic stem-cell differentiation in vivo. Proc Natl Acad Sci U S A. 1991, 88: 2788-2792. 10.1073/pnas.88.7.2788.PubMedPubMed CentralView ArticleGoogle Scholar
- Harrison DE, Astle CM, Lerner C: Number and continuous proliferative pattern of transplanted primitive immunohematopoietic stem cells. Proc Natl Acad Sci U S A. 1988, 85: 822-826. 10.1073/pnas.85.3.822.PubMedPubMed CentralView ArticleGoogle Scholar
- Weber K, Thomaschewski M, Warlich M, Volz T, Cornils K, Niebuhr B, Tager M, Lutgehetmann M, Pollok JM, Stocking C, Dandri M, Benten D, Fehse B: RGB marking facilitates multicolor clonal cell tracking. Nat Med. 2011, 17: 504-509. 10.1038/nm.2338.PubMedView ArticleGoogle Scholar
- Stewart MH, Bendall SC, Levadoux-Martin M, Bhatia M: Clonal tracking of hESCs reveals differential contribution to functional assays. Nat Methods. 2010, 7: 917-922. 10.1038/nmeth.1519.PubMedView ArticleGoogle Scholar
- Ciuffi A, Ronen K, Brady T, Malani N, Wang G, Berry CC, Bushman FD: Methods for integration site distribution analyses in animal cell genomes. Methods. 2009, 47: 261-268. 10.1016/j.ymeth.2008.10.028.PubMedPubMed CentralView ArticleGoogle Scholar
- Harkey MA, Kaul R, Jacobs MA, Kurre P, Bovee D, Levy R, Blau CA: Multiarm high-throughput integration site detection: limitations of LAM-PCR technology and optimization for clonal analysis. Stem Cells Dev. 2007, 16: 381-392. 10.1089/scd.2007.0015.PubMedView ArticleGoogle Scholar
- Kim S, Kim N, Presson AP, An DS, Mao SH, Bonifacino AC, Donahue RE, Chow SA, Chen IS: High-throughput, sensitive quantification of repopulating hematopoietic stem cell clones. J Virol. 2010, 84: 11771-11780. 10.1128/JVI.01355-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu C, Jares A, Winkler T, Xie J, Metais JY, Dunbar CE: High efficiency restriction enzyme-free linear amplification-mediated polymerase chain reaction approach for tracking lentiviral integration sites does not abrogate retrieval bias. Hum Gene Ther. 2013, 24: 38-47.PubMedPubMed CentralView ArticleGoogle Scholar
- Cornils K, Lange C, Schambach A, Brugman MH, Nowak R, Lioznov M, Baum C, Fehse B: Stem cell marking with promotor-deprived self-inactivating retroviral vectors does not lead to induced clonal imbalance. Mol Ther. 2009, 17: 131-143. 10.1038/mt.2008.238.PubMedPubMed CentralView ArticleGoogle Scholar
- Brugman MH, Suerth JD, Rothe M, Suerbaum S, Schambach A, Modlich U, Kustikova O, Baum C: Evaluating a ligation-mediated PCR and pyrosequencing method for the detection of clonal contribution in polyclonal retrovirally transduced samples. Hum Gene Ther Methods. 2013, 24: 68-79. 10.1089/hgtb.2012.175.PubMedPubMed CentralView ArticleGoogle Scholar
- Lu R, Neff NF, Quake SR, Weissman IL: Tracking single hematopoietic stem cells in vivo using high-throughput sequencing in conjunction with viral genetic barcoding. Nat Biotechnol. 2011, 29: 928-933. 10.1038/nbt.1977.PubMedPubMed CentralView ArticleGoogle Scholar
- Gerrits A, Dykstra B, Kalmykowa OJ, Klauke K, Verovskaya E, Broekhuis MJ, de Haan G, Bystrykh LV: Cellular barcoding tool for clonal analysis in the hematopoietic system. Blood. 2010, 115: 2610-2618. 10.1182/blood-2009-06-229757.PubMedView ArticleGoogle Scholar
- Schepers K, Swart E, van Heijst JW, Gerlach C, Castrucci M, Sie D, Heimerikx M, Velds A, Kerkhoven RM, Arens R, Schumacher TN: Dissecting T cell lineage relationships by cellular barcoding. J Exp Med. 2008, 205: 2309-2318. 10.1084/jem.20072462.PubMedPubMed CentralView ArticleGoogle Scholar
- Verovskaya E, Broekhuis MJ, Zwart E, Ritsema M, van Os R, de Haan G, Bystrykh LV: Heterogeneity of young and aged murine hematopoietic stem cells revealed by quantitative clonal analysis using cellular barcoding. Blood. 2013, 122: 523-532. 10.1182/blood-2013-01-481135.PubMedView ArticleGoogle Scholar
- Gerlach C, Rohr JC, Perie L, van Rooij N, van Heijst JW, Velds A, Urbanus J, Naik SH, Jacobs H, Beltman JB, de Boer RJ, Schumacher TN: Heterogeneous differentiation patterns of individual CD8+ T cells. Science. 2013, 340: 635-639. 10.1126/science.1235487.PubMedView ArticleGoogle Scholar
- ClusterSeq. [https://www.github.com/adaptivegenome/clusterseq/]
- Zhang B, Metharom P, Jullie H, Ellem KA, Cleghorn G, West MJ, Wei MQ: The significance of controlled conditions in lentiviral vector titration and in the use of multiplicity of infection (MOI) for predicting gene transfer events. Genet Vaccines Ther. 2004, 2: 6-10.1186/1479-0556-2-6.PubMedPubMed CentralView ArticleGoogle Scholar
- Lozzio CB, Lozzio BB: Human chronic myelogenous leukemia cell-line with positive Philadelphia chromosome. Blood. 1975, 45: 321-334.PubMedGoogle Scholar
- Montini E, Cesana D, Schmidt M, Sanvito F, Ponzoni M, Bartholomae C, Sergi Sergi L, Benedicenti F, Ambrosi A, Di Serio C, Doglioni C, von Kalle C, Naldini L: Hematopoietic stem cell gene transfer in a tumor-prone mouse model uncovers low genotoxicity of lentiviral vector integration. Nat Biotechnol. 2006, 24: 687-696. 10.1038/nbt1216.PubMedView ArticleGoogle Scholar
- Modlich U, Navarro S, Zychlinski D, Maetzig T, Knoess S, Brugman MH, Schambach A, Charrier S, Galy A, Thrasher AJ, Bueren J, Baum C: Insertional transformation of hematopoietic cells by self-inactivating lentiviral and gammaretroviral vectors. Mol Ther. 2009, 17: 1919-1928. 10.1038/mt.2009.179.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu R, Paxton WA, Choe S, Ceradini D, Martin SR, Horuk R, MacDonald ME, Stuhlmann H, Koup RA, Landau NR: Homozygous defect in HIV-1 coreceptor accounts for resistance of some multiply-exposed individuals to HIV-1 infection. Cell. 1996, 86: 367-377. 10.1016/S0092-8674(00)80110-5.PubMedView ArticleGoogle Scholar
- Voit RA, McMahon MA, Sawyer SL, Porteus MH: Generation of an HIV resistant T-cell line by targeted 'stacking' of restriction factors. Mol Ther. 2013, 21: 786-795. 10.1038/mt.2012.284.PubMedPubMed CentralView ArticleGoogle Scholar
- Mussolino C, Morbitzer R, Lutge F, Dannemann N, Lahaye T, Cathomen T: A novel TALE nuclease scaffold enables high genome editing activity in combination with low toxicity. Nucleic Acids Res. 2011, 39: 9283-9293. 10.1093/nar/gkr597.PubMedPubMed CentralView ArticleGoogle Scholar
- Perez EE, Wang J, Miller JC, Jouvenot Y, Kim KA, Liu O, Wang N, Lee G, Bartsevich VV, Lee YL, Guschin DY, Rupniewski I, Waite AJ, Carpenito C, Carroll RG, Orange JS, Urnov FD, Rebar EJ, Ando D, Gregory PD, Riley JL, Holmes MC, June CH: Establishment of HIV-1 resistance in CD4+ T cells by genome editing using zinc-finger nucleases. Nat Biotechnol. 2008, 26: 808-816. 10.1038/nbt1410.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu J: Preparation, culture, and immortalization of mouse embryonic fibroblasts. Curr Protoc Mol Biol. 2005, Chapter 28:Unit 28 21Google Scholar
- Amann J, Kalyankrishna S, Massion PP, Ohm JE, Girard L, Shigematsu H, Peyton M, Juroske D, Huang Y, Stuart Salmon J, Kim YH, Pollack JR, Yanagisawa K, Gazdar A, Minna JD, Kurie JM, Carbone DP: Aberrant epidermal growth factor receptor signaling and enhanced sensitivity to EGFR inhibitors in lung cancer. Cancer Res. 2005, 65: 226-235.PubMedGoogle Scholar
- Martins FC, De S, Almendro V, Gonen M, Park SY, Blum JL, Herlihy W, Ethington G, Schnitt SJ, Tung N, Garber JE, Fetten K, Michor F, Polyak K: Evolutionary pathways in BRCA1-associated breast tumors. Cancer Discov. 2012, 2: 503-511. 10.1158/2159-8290.CD-11-0325.PubMedPubMed CentralView ArticleGoogle Scholar
- Maley CC, Galipeau PC, Finley JC, Wongsurawat VJ, Li X, Sanchez CA, Paulson TG, Blount PL, Risques RA, Rabinovitch PS, Reid BJ: Genetic clonal diversity predicts progression to esophageal adenocarcinoma. Nat Genet. 2006, 38: 468-473. 10.1038/ng1768.PubMedView ArticleGoogle Scholar
- Ellis BL, Potts PR, Porteus MH: Creating higher titer lentivirus with caffeine. Hum Gene Ther. 2011, 22: 93-100. 10.1089/hum.2010.068.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.