1. Introduction
Xenotransplantation of human cells into immunodeficient mice is an essential component of cell therapy development [
1,
2,
3]. By determining the distribution and retention of cell products after administration to mice, important questions can be addressed, such as whether the delivered cells actually reach the target site(s), whether they engraft in sufficient numbers to produce the desired effect, whether they are distributed to unwanted non-target tissues, and how long they persist in the host organism. Answering these questions is critical to uncovering the underlying mechanism(s) of action, exploring interventions to improve cell engraftment ratios and evaluating the biosafety profile of cell-based therapy strategies as required by the regulatory authorities [
3,
4,
5].
Bioanalysis studies evaluating the distribution, persistence and clearance of cell therapy products place special demands on the analytical method. First, the method must be highly specific and sensitive to detect and quantify very small populations of human cells that are vastly outnumbered by the mouse cells surrounding them. Secondly, because transplanted cells can be distributed in an inhomogeneous manner within organs [
6,
7], the method must ensure that its accuracy is not biased by sampling errors due to non-random cell distribution. Thirdly, the method should ideally not require the cells to be pre-labeled, since this may affect their viability and/or functionality [
8,
9,
10] and would therefore require the comparability of the labeled cells to the (unlabeled) cells intended for human use to be established [
10,
11,
12]. Given these requirements, quantitative polymerase chain reaction (qPCR) has evolved as a method of choice in bioanalysis of cell therapy products. qPCR assays enable specific and extremely sensitive tracking and absolute quantification of the donor cells via detection of species-specific DNA sequences within the whole host organ. This makes qPCR suitable for rapid and convenient systematic quantification of unlabeled donor cells even at very low cell numbers in a broad range of host organs and tissues [
4,
5,
13,
14].
However, the performance of a qPCR assay depends on a number of parameters, including but not limited to tissue sampling and nucleic acid extraction methods, choice of primers and probes, selection of reagents and reaction conditions, determination of the quantification cycle (Cq) values, and matrix effects [
15,
16]. Therefore, any qPCR method must be validated to ensure the reliability of the data that will be obtained. Remarkably, however, despite the growing number of qPCR applications in regulated bioanalysis, and even though qPCR is one of the officially recognized methods for determining the pharmacokinetics of cell therapy products [
11], regulatory guidance on the use of qPCR is limited [
17,
18,
19,
20,
21]. The International Council for Harmonisation of Technical Requirements for Pharmaceuticals for Human Use (ICH) Guideline M10 on Bioanalytical Method Validation and Study Sample Analysis [
22,
23], which recently superseded the guidelines on bioanalytical method validation released by the European Medicines Agency (EMA) [
24] and the US Food and Drug Administration (FDA) [
25], focuses on pharmacokinetic methods that are suitable for classical small-molecule drugs and large-molecule biologics, such as chromatographic methods and ligand-binding assays, but do not address the specific requirements for proper validation of PCR assays [
17,
18,
19].
Another challenge in bioanalytical qPCR studies is that the different organs and tissues of the host organism can affect amplification efficiency and amplicon detection to varying degrees, and ignoring such matrix effects can easily lead to an underestimation of the true number of target cells in a sample. Matrix effects have been studied predominantly in environmental microbiology, microbial food safety and forensic analyses, where the amounts of target nucleic acids are often extremely small and the matrices particularly diverse and challenging [
26,
27,
28,
29]. In contrast, with the exception of forensically and microbiologically relevant body fluids and secretions as well as food safety-relevant matrices such as muscle tissue and milk, the effects of mammalian matrices on qPCR results have been assessed and discussed only in a very limited manner and only for a limited number of tissue types [
13,
30,
31,
32,
33,
34].
With particular attention to these challenges, we describe here the development and validation of a qPCR assay for reliable detection and quantification of human cells in mouse tissues and blood. The validation included a systematic assessment of the matrix effects of 14 different mouse tissues including blood. This assay enabled the generation of preclinical biodistribution data acceptable to regulatory authorities [
35], which were required for the approval of a medicinal product based on skin-derived ABCB5
+ mesenchymal stromal cells (MSCs) [
36,
37] to be tested in clinical trials. The insights presented here may also be relevant for a range of other scientific contexts and purposes.
2. Materials and Methods
2.1. Assay design
A TaqMan® qPCR assay was developed for detection and quantification of human ABCB5
+ MSCs in mouse tissues by detection of a DNA sequence of human α-satellite DNA [
38]. As an internal control of efficient DNA extraction from mouse tissues, in mouse tissue homogenates also a mouse-specific DNA sequence of the prostaglandin E receptor 2 (
PTGER2) gene [
39] was detected. The assays were run in a certified GLP-compliant test facility (Accelero Bioanalytics, Berlin, Germany). Reporting follows the MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments) Guidelines where applicable [
40].
2.2. Primers and Probes
Human-specific DNA was detected by amplification of a sequence of the α-satellite DNA on chromosome 17 (HSSATA17, GenBank Acc. No. M13882), using forward primer GGGATAATTTCAGCTGACTAAACAG, reverse primer AAACGTCCACTTGCAGATTCTA and 6-carboxyfluorescein-labeled probe CACGTTTGAAACACTCTTTTTGCA carrying the Black Hole Quencher® BHQ®-1. Mouse-specific DNA was detected by amplification of a mouse-specific DNA fragment of the
PTGER2 gene using forward primer TACCTGCAGCTGTACGCCAC, reverse primer GCCAGGAGAATGAGGTGGTC and carboxytetramethylrhodamine-labeled probe CCTGCTGCTTATCGTGGCTG carrying BHQ®-2. The specificity of these sequences has been confirmed previously [
38,
39]. All primers and probes were supplied by Microsynth (Balgach, Switzerland). Primer and probe lyophilizates were dissolved in DNase-free water to prepare 100-µM stock solutions, from which 18-µM (primers) and 5-µM (probes) working solutions were prepared.
2.3. Human ABCB5+ MSCs
Human ABCB5
+ MSCs were derived from skin samples taken from human subjects undergoing abdominoplasties or other medical interventions providing leftover skin tissue. Skin sampling was performed in accordance with the German Act on Organ and Tissue Donation, Removal and Transplantation after written informed consent was obtained from each donor. Cell production was carried out in an EU-GMP grade A cabinet in a grade B clean room facility under laminar air flow according to a validated GMP-compliant manufacturing protocol as described previously [
36]. In brief, skin tissue was enzymatically digested, and cells were expanded as an unsegregated culture by serial passaging upon adherence selection in an in-house MSC-favoring culture medium. From these cultures, ABCB5
+ MSCs were isolated by antibody-coupled magnetic bead sorting using a mouse anti-human monoclonal antibody directed against the extracellular loop 3 of the ABCB5 molecule [
41] (Maine Biotechnology Services, Portland, Maine; GMP-compliant purification: Bibitec, Bielefeld, Germany). After enzymatic detachment of the beads from the cell surface, the isolated ABCB5
+ MSCs were cryo-preserved in CryoStor® CS10 freeze medium (BioLife Solution, Bothell, WA) containing 10% dimethyl sulfoxide and stored in the vapor phase of liquid nitrogen.
For spiking, ABCB5+ MSCs suspensions were thawed at 37 °C in a thermal mixer. To remove cell debris and free, degraded DNA, the cell suspension was washed with 1× PBS. After centrifugation for 5 min at 500×g, the cell pellet was taken up in PBS to adjust the cell concentration to approximately 2000 cells/µl. Effective cell concentration as assessed by cell counting under a light microscope using a hemocytometer was 2008 cells/µl. The cell suspension was aliquoted and stored at -80 to -60 °C.
2.4. Mouse tissue sampling
Animal breeding, care, necropsy and tissue collection were conducted by a specialized contract research organization (Preclinics, Potsdam, Germany). SCID/beige mice (21-23 weeks old, 5 male, 5 female) were anesthetized with isoflurane and whole blood was collected by cardiac puncture into EDTA collection tubes. The animals were then sacrificed by an overdose of xylazine and the following tissues were collected: mouse skin (neck region), thigh muscle (M. quadriceps femoris), lymph nodes (cervical, axial and inguinal lymph nodes pooled), liver, spleen, lung, brain, bone (femur) including marrow, kidneys, thymus, thyroid and ovaries/testes.
To avoid contamination with human DNA and cross-contamination between animals, necropsy and tissue collection were performed under a laminar airflow workbench. All work areas were disinfected before work commenced. Personnel wore disposable lab coats, hoods and face masks, and two pairs of gloves. The gloves were disinfected before work commenced and the outer pair of gloves was changed after each animal. A separate autoclaved dissection set was used for every four animals (one cage) and was disinfected after each animal. Tissue pads for necropsy were changed after each animal. Tissues were collected in DNAse-free collection tubes.
2.5. Preparation of reference standards
Reference standards were prepared and samples were lysed and DNA extracted using the NucleoSpin® 96 Tissue kit (Macherey-Nagel, Düren, Germany) according to the manufacturer’s instructions, following the protocols described below.
2.5.1. Calibration standards and quality control standards
Seven calibration standards (range: 125 to 20,000 ABCB5+ MSCs) and five quality control standards (range: 125 to 15,000 ABCB5+ MSCs) were prepared by serial dilution of ABCB5+ MSC suspension in Tris-EDTA buffer. Tris-EDTA buffer without cells was used as blank sample. Twenty microliters of each cell suspension or blank sample was added to 180 µl lysis buffer T1, followed by the addition of 25 µl proteinase K. Samples were lysed at 56 °C for 15-30 min. The lysates were stored at -80 to -60 °C until DNA was eluted using 60 µl of pre-heated (70 °C) elution buffer BE.
2.5.2. Tissue quality control standards
Mouse tissues were taken up in lysis buffer T1 (amount as required to adjust the intended tissue concentration, see
Table S1) into homogenization tubes filled with ceramic beads (for soft tissues) or steel beads (for bone tissue; addition of lysis buffer T1 only after two “dry” homogenization cycles without lysis buffer), and homogenized in the Precellys Evolution homogenizer (Bertin Technologies, Frankfurt, Germany) at 6000 rpm and room temperature for 20 s per cycle with at least 30 s pause between cycles. If more than two cycles were required for complete homogenization (see
Table S1 for the total number of cycles for the different tissues), the samples were cooled between each cycle. Homogenates (225 µl each) were spiked with 25 µl cell suspension containing 0, 125, 625 or 5000 ABCB5
+ MSCs. The spiked homogenates were homogenized for a further cycle, and then 25 µl proteinase K was added. Samples were lysed at 56 °C for 15-30 min.
For preparation of blood quality control standards, 100 µl EDTA blood was spiked with 20 µl cell suspension (containing 0, 125, 625 or 5000 ABCB5+ MSCs), filled up to 400 µl with PBS, and then 25 µl (assays 1 and 2) or 50 µl (assay 3) proteinase K and 400 µl binding buffer BQ1 was added. The samples were lysed at room temperature for 5 min (assays 1 and 2) or 30 min (assay 3) followed by 70 °C for 15 min.
All lysates were stored at -80 to -60 °C until DNA was eluted using 60 µl of pre-heated (70 °C) elution buffer BE.
All steps related to tissue transfer, cutting, splitting, lysis, DNA extraction and transfer of the eluate to the qPCR plate were performed under a laminar airflow workbench using sterile, disposable equipment. Personnel wore two pairs of sterile gloves; the outer pair and the tools used to split and transfer the samples were changed between each sample. Where tissue samples had to be sectioned, sterile DNA-free 24-well culture plates were used as a sectioning surface.
2.5.3. Freeze-thaw stability
To test for freeze-thaw stability of the extracted DNA, DNA eluted from quality control standards was divided into two aliquots. One set of aliquots was stored at 2 to 8 °C and the other set at -25 to -15 °C for at least 12 hours prior to analysis.
2.6. Amplification
The master mix for a singleplex PCR reaction consisted of 5 µl GoTaq Probe qPCR Master Mix (Promega, Madison, WI, USA) added with 0.5 µl each of DNase-free water, 18-µM forward primer, 18-µM reverse primer and 5-µM probe. The master mix (7 µl) was mixed with 3 µl of a 10-fold dilution (in DNase-free water) of template DNA, resulting in a reaction volume of 10 µl. For no-template control, DNAse-free water was used. Amplifications were run on an Applied BiosystemsTM ViiATM 7 Dx Real-Time PCR instrument (Thermo Fisher Scientific, Langenselbold, Germany). The cycling program consisted of an initial denaturation step of 10 min at 95 °C followed by 45 cycles of 15 seconds at 95 °C and 1 minute at 60 °C. For detection of human DNA, samples were assayed in triplicate; except for samples spiked with lower cell numbers (250 and below), which were assayed in sextuplicate. For detection of mouse DNA, samples were assayed in monoplicate.
2.7. Assay validation and acceptance criteria
The assay was validated in accordance with the general requirements for bioanalytical method validation set out in the EMA Guideline on Bioanalytical Method Validation [
24], recently superseded by the ICH guideline M10 on Bioanalytical Method Validation and Study Sample Analysis [
22], evaluating the parameters linearity, intra-assay and inter-assay accuracy and precision, specificity, freeze-thaw stability and tissue matrix effects against predefined acceptance criteria (
Table 1).
2.7.1. Linearity
Calibration curves were generated by plotting the Cq number against the logarithm of the cell numbers in the calibration standards. A linear calibration function:
was fitted by least-square regression. Linearity was assumed if the correlation coefficient r
2 was ≥ 0.95. The equation of the best-fit line was used to back-calculate the numbers of cell equivalents. Assay efficiency (E) was calculated as:
2.7.2. Accuracy and precision
Intra-assay and inter-assay accuracy (expressed as percent bias of the calculated cell number from the nominal cell number; acceptance criterium: within ± 40%) and precision [expressed as percent coefficient of variation (CV) between a series of measurements of the same sample; acceptance criterium: ≤ 40%] were determined on the calibration standards and quality control standards assayed in triplicate or sextuplicate (as specified above). Three independent runs were performed on three different days. For assay validation, at least 75% of the calibration standard samples and at least 67% of the quality control standards had to meet the acceptance criteria for accuracy and precision.
2.7.3. Specificity
To demonstrate specificity of the assay, the no-template controls were required to give either no amplification signal or a Cq value unequivocally distinguishable from the lower limit of quantification (LLOQ). To confirm that the assay specifically quantifies human cells even in the presence of mouse cells, tissue blanks, i.e. mouse tissue samples not spiked with human cells, were also analyzed.
2.7.4. Tissue matrix effects
The effects of mouse tissues on the DNA extraction and assay performance were assessed by determining the spike recovery in the tissue quality control standards (mouse tissues spiked with ABCB5+ MSCs). Three independent runs were performed on three different days. Spike recovery rates (ratio of measured to the nominal cell counts, expressed as a percentage of the nominal cell count) were used to calculate matrix factors (defined as the reciprocal of the percent recovery rate) for each tissue analyzed.
4. Discussion
Although qPCR-based assays have emerged as an important bioanalytical method for assessing the pharmacokinetics of human cell-based medicinal products [
4,
5,
13,
19], the regulatory guidelines on the validation of bioanalytical methods released by the EMA, FDA and ICH [
22,
23,
24,
25] focus on methods suitable for small- and large-molecule drugs such as chromatographic and ligand-binding assays. While the basic concepts and parameters of method validation described in these guidelines can be adapted to cell quantification by qPCR, in the absence of specific regulatory recommendations including definitive acceptance criteria for a validated qPCR assay, researchers must rely on published evidence from bioanalytical scientists as well as recently issued best practice recommendations [
18,
42,
43] and white papers from scientific networks [
19,
20,
21,
44].
The validation presented here followed the validation parameters set out in the EMA Guideline on Bioanalytical Method Validation [
24], recently superseded by the ICH Guideline M10 on Bioanalytical Method Validation and Study Sample Analysis [
22], using predefined acceptance criteria (
Table 1). The acceptance criteria for accuracy and precision, which were set to within ± 40% and ≤ 40%, respectively, are within the range of those recently recommended by the European Bioanalysis Forum [
20]. The data obtained show that human MSCs can be detected and quantified with acceptable linearity, accuracy and precision within the range of 125 (LLOQ) and 20,000 (ULOQ) cells/200 µl lysate (
Table 1).
When developing and validating a bioanalytical assay, it is essential to match the setup to the actual study in which the assay will be used [
22,
23]. Besides factors such as the intended mouse strain(s) or the amount of available sample material, an important aspect of bioanalytical cell detection by qPCR is the fact that many components in biological matrices can bias the cell quantification results. Such components include molecules that can impair DNA extraction or, after being co-extracted with the DNA, can affect amplification by disturbing annealing of the primers or inhibiting the DNA polymerase, or interfere with amplicon detection by quenching fluorescence or interacting with the fluorophore [
13,
16,
27,
29]. In contrast to environmental, food or forensic qPCR, there is only limited information on substances that may interfere with qPCR-based cell detection in the field of bioanalysis. Known inhibitory molecules present in bioanalytically relevant tissues and body fluids include molecules present in skin, muscle and bone such as melanin, myoglobin, collagen and calcium ions, as well as various blood constituents including added anticoagulants (
Table 6). However, the issue of matrix effects on cell quantification by qPCR must be considered in all tissue types of the body. The yield and quality of extracted genomic DNA can vary widely depending on the physical and biochemical nature of each tissue [
13,
30], and the tissue from which the DNA was extracted can have a significant effect on the efficiency, accuracy and precision of the qPCR assay [
33]. Therefore, current regulatory guidelines [
22,
23] and best-practice recommendations [
18,
19,
20,
43,
44] advise to assess potential matrix effects by determining the recovery of target DNA spiked into each tissue of interest during assay validation, with recovery rates in a wide range between 30% and 100% considered to be expected [
18,
43].
Our experiments on a wide variety of tissues revealed an even wider range of recovery rates between the various tissue types, with two out of the 14 tissues showing recovery rates below 30% (one of which, testes, with a recovery rate of 29% just missing the range) and three out of the 14 tissues showing recovery rates above 100% (
Figure 1). The lowest recovery was achieved in blood samples, where only 11% of spiked cells were detected. In this respect, it is important to note that blood is generally considered to be a particularly challenging matrix [
51,
52,
55]. PCR mixtures based on
Taq DNA polymerases have been reported to be inhibited in the presence of 1% [
56] or even 0.004% (v/v) EDTA whole blood [
57], and whole blood components co-extracted with the DNA can cause several negative effects such as loss of amplifiable target DNA, reduction of amplification efficiency and quenching of fluorescence [
50,
51,
52,
53,
54]. On the other hand, thorough purification of DNA extracts can result in significant loss of DNA, which can also reduce recovery rates to as low as 10% [
29]. Modifications to the extraction protocol, such as increasing the volume of proteinase K and the incubation time, have been reported to increase the yield of amplifiable DNA from blood samples [
58], but did not improve the recovery rates of human ABCB5
+ MSCs spiked into mouse blood samples in the present study (
Table S4).
Overall, the data from the spike-and-recovery experiments demonstrate that a qPCR assay for bioanalytical studies runs the risk of substantially underestimating or even overestimating cell numbers, if the potential matrix effects due to physical and biochemical differences between the various tissues of the host organisms are not taken into account. Interestingly, a trend towards higher recovery rates at lower spiked cell concentrations was observed in almost all tissues in our assay (
Figure 1). Although we have no causal explanation for this observation, this trend suggests that the risk of underestimation decreases towards the lower limit of the validated quantification range.
In any case, while regulatory authorities require the determination of tissue matrix effects as part of bioanalytical method validation [
22,
25], they do not provide guidance on how to handle the results. An ideal, thorough assay optimization to minimize the impact of the different tissue matrices on cell quantification results would require an elaborate, costly and animal-intensive program to evaluate the amount, integrity and purity of the DNA extracted in different ways and/or to assess the amplification efficiency for, e.g., different buffer compositions and/or added facilitators [
16,
27,
29]. Such complex programs, which would have to be performed and validated separately for each tissue of interest, would be beyond the resources of a research group or cell therapy-developing company, especially as bioanalytical studies require the analysis of a wide range of different tissues. Instead, for reasons of feasibility in the sense of a fit-for-purpose approach [
22,
23], we suggest correcting for the effects of the different tissue matrices on cell quantification in the subsequent actual sample measurements by multiplying the back-calculated cell count by a matrix-specific factor, representing the reciprocal of the percent recovery rate (
Table 5).
An important question in bioanalytical studies is whether the cells that are detected are actually alive and, as such, potentially active, or whether they are not [
59]. Researchers need to note that, unlike other cell detection methods such as flow cytometry and optical imaging [
4,
60], qPCR by itself cannot distinguish between nucleic acids extracted from live cells, cell fragments, or cell corpses engulfed by local macrophages [
5]. However,
in vivo, genomic DNA is rapidly degraded upon cell death as an intrinsic part of the apoptotic program and/or by lysosomal DNAses of cells that have phagocytosed the apoptotic or necrotic cell corpses [
61,
62]. Therefore, it is widely accepted that inadvertent quantification of DNA isolated from dead cells is rather unlikely and, in the context of human cell xenotransplantation, the presence of dead human cells or residual human genomic DNA would not be expected to significantly bias the quantification of live human cells [
34,
59,
63,
64]. As with the distinction between living and dead cells, qPCR is also not able to distinguish between proliferating and non-proliferating cells. For cells that are intended for use as cell therapy products, however, any potential for unwanted proliferation must be excluded [
5,
65]. Therefore, in bioanalytical studies assessing the biosafety of a cell-based therapy, cell quantification by qPCR needs to be complemented by an appropriate method to assess the proliferative activity of the detected cells, e.g., immunohistochemical double-staining of tissue slides with a human-specific antibody and an antibody against a proliferation marker such as Ki67 [
35].
Author Contributions
Conceptualization: H.M.S., J.E. and M.A.K., Validation, H.M.S. and M.A.K., Formal Analysis, H.M.S. and E.N.R., Data Curation, H.M.S. and A.N., Writing – Original Draft Preparation, E.N.-R., Writing, Review & Editing, H.M.S., A.N., J.E., C.G., M.H.F. and M.A.K., Visualization, E.N.-R., Supervision, C.G., M.H.F. and M.A.K., Project Administration, H.M.S.