HLA-A Locus is Associated With Sepsis and Septic Shock After Traumatic Injury.

OBJECTIVE: Determine whether variation in the HLA region is associated with the development of post-traumatic sepsis and septic shock. BACKGROUND: Sepsis-related deaths remain a major source of mortality after traumatic injury. Genetic characteristics may contribute to susceptibility to adverse outcomes including sepsis and septic shock. Recent advances in next-generation sequencing technology now allow comprehensive genotyping of the HLA region. METHODS: White adult trauma patients requiring more than 2 days of mechanical ventilation underwent HLA genotyping, and were followed for the development of sepsis and septic shock. Odds ratios (OR) for the associations between our outcomes and HLA variants were estimated, a correction for multiple comparisons was applied, and significant variants were included in regression models adjusting for potential confounders. RESULTS: A total of 1184 patients were included. Patients were severely injured (median injury severity score 33); 33% developed sepsis, 6% septic shock, and in-hospital mortality was 14%. An amino acid variant (156Q) within the HLA-A peptide-binding groove was associated with greater odds of sepsis [OR 1.50, (1.18-1.89)]. HLA-A∗02:01 was associated with lower odds of septic shock [OR 0.52, (0.32-0.82)]. These associations remained significant after adjusting for potential confounders. CONCLUSIONS: This is the first study to apply next-generation sequencing techniques to evaluate associations between immunogenetic factors and post-traumatic sepsis and septic shock. Associations with class I HLA variants are novel as they implicate adaptive immunity in post-traumatic sepsis. These findings are a step towards developing a panel of genetic markers assessing risk of infection-related complications as we move towards more personalized medicine.

HLA-DPB1 and HLA class I confer risk of and protection from narcolepsy.

Type 1 narcolepsy, a disorder caused by a lack of hypocretin (orexin), is so strongly associated with human leukocyte antigen (HLA) class II HLA-DQA1(∗)01:02-DQB1(∗)06:02 (DQ0602) that very few non-DQ0602 cases have been reported. A known triggering factor for narcolepsy is pandemic 2009 influenza H1N1, suggesting autoimmunity triggered by upper-airway infections. Additional effects of other HLA-DQ alleles have been reported consistently across multiple ethnic groups. Using over 3,000 case and 10,000 control individuals of European and Chinese background, we examined the effects of other HLA loci. After careful matching of HLA-DR and HLA-DQ in case and control individuals, we found strong protective effects of HLA-DPA1(∗)01:03-DPB1(∗)04:02 (DP0402; odds ratio [OR] = 0.51 [0.38-0.67], p = 1.01 × 10(-6)) and HLA-DPA1(∗)01:03-DPB1(∗)04:01 (DP0401; OR = 0.61 [0.47-0.80], p = 2.07 × 10(-4)) and predisposing effects of HLA-DPB1(∗)05:01 in Asians (OR = 1.76 [1.34-2.31], p = 4.71 × 10(-05)). Similar effects were found by conditional analysis controlling for HLA-DR and HLA-DQ with DP0402 (OR = 0.45 [0.38-0.55] p = 8.99 × 10(-17)) and DP0501 (OR = 1.38 [1.18-1.61], p = 7.11 × 10(-5)). HLA-class-II-independent associations with HLA-A(∗)11:01 (OR = 1.32 [1.13-1.54], p = 4.92 × 10(-4)), HLA-B(∗)35:03 (OR = 1.96 [1.41-2.70], p = 5.14 × 10(-5)), and HLA-B(∗)51:01 (OR = 1.49 [1.25-1.78], p = 1.09 × 10(-5)) were also seen across ethnic groups in the HLA class I region. These effects might reflect modulation of autoimmunity or indirect effects of HLA class I and HLA-DP alleles on response to viral infections such as that of influenza.

Digital microfluidic assay for protein detection.

Global studies of the human proteome have revealed a plethora of putative protein biomarkers. However, their application for early disease detection remains at a standstill without suitable methods to realize their utility in the clinical setting. There thus continues to be tremendous interest in developing new technology for sensitive protein detection that is both low in cost and carries a small footprint to be able to be used at the point of care. The current gold standard method for protein biomarker detection is the ELISA, which measures protein abundance using bulky fluorescent scanners that lack portability. Here, we present a digital microfluidic platform for protein biomarker detection that is low in cost compared with standard optical detection methods, without any compromise in sensitivity. This platform furthermore makes use of simple electronics, enabling its translation into a portable handheld device, and has been developed in a manner that can easily be adapted to assay different types of proteomic biomarkers. We demonstrate its utility in quantifying not only protein abundance, but also activity. Interleukin-6 abundance could be assayed from concentrations as low as 50 pM (an order of magnitude lower than that detectable by a comparable laboratory designed ELISA) using less than 5 µL of sample, and Abelson tyrosine kinase activity was detectable in samples containing 100 pM of kinase.

Molecular indexing enables quantitative targeted RNA sequencing and reveals poor efficiencies in standard library preparations.

We present a simple molecular indexing method for quantitative targeted RNA sequencing, in which mRNAs of interest are selectively captured from complex cDNA libraries and sequenced to determine their absolute concentrations. cDNA fragments are individually labeled so that each molecule can be tracked from the original sample through the library preparation and sequencing process. Multiple copies of cDNA fragments of identical sequence become distinct through labeling, and replicate clones created during PCR amplification steps can be identified and assigned to their distinct parent molecules. Selective capture enables efficient use of sequencing for deep sampling and for the absolute quantitation of rare or transient transcripts that would otherwise escape detection by standard sequencing methods. We have also constructed a set of synthetic barcoded RNA molecules, which can be introduced as controls into the sample preparation mix and used to monitor the efficiency of library construction. The quantitative targeted sequencing revealed extremely low efficiency in standard library preparations, which were further confirmed by using synthetic barcoded RNA molecules. This finding shows that standard library preparation methods result in the loss of rare transcripts and highlights the need for monitoring library efficiency and for developing more efficient sample preparation methods.

Genomic responses in mouse models poorly mimic human inflammatory diseases.

A cornerstone of modern biomedical research is the use of mouse models to explore basic pathophysiological mechanisms, evaluate new therapeutic approaches, and make go or no-go decisions to carry new drug candidates forward into clinical trials. Systematic studies evaluating how well murine models mimic human inflammatory diseases are nonexistent. Here, we show that, although acute inflammatory stresses from different etiologies result in highly similar genomic responses in humans, the responses in corresponding mouse models correlate poorly with the human conditions and also, one another. Among genes changed significantly in humans, the murine orthologs are close to random in matching their human counterparts (e.g., R(2) between 0.0 and 0.1). In addition to improvements in the current animal model systems, our study supports higher priority for translational medical research to focus on the more complex human conditions rather than relying on mouse models to study human inflammatory diseases.

High-throughput, high-fidelity HLA genotyping with deep sequencing.

Human leukocyte antigen (HLA) genes are the most polymorphic in the human genome. They play a pivotal role in the immune response and have been implicated in numerous human pathologies, especially autoimmunity and infectious diseases. Despite their importance, however, they are rarely characterized comprehensively because of the prohibitive cost of standard technologies and the technical challenges of accurately discriminating between these highly related genes and their many allelles. Here we demonstrate a high-resolution, and cost-effective methodology to type HLA genes by sequencing, which combines the advantage of long-range amplification, the power of high-throughput sequencing platforms, and a unique genotyping algorithm. We calibrated our method for HLA-A, -B, -C, and -DRB1 genes with both reference cell lines and clinical samples and identified several previously undescribed alleles with mismatches, insertions, and deletions. We have further demonstrated the utility of this method in a clinical setting by typing five clinical samples in an Illumina MiSeq instrument with a 5-d turnaround. Overall, this technology has the capacity to deliver low-cost, high-throughput, and accurate HLA typing by multiplexing thousands of samples in a single sequencing run, which will enable comprehensive disease-association studies with large cohorts. Furthermore, this approach can also be extended to include other polymorphic genes.

Benchmarking outcomes in the critically injured burn patient.

OBJECTIVE: To determine and compare outcomes with accepted benchmarks in burn care at 6 academic burn centers. BACKGROUND: Since the 1960s, US morbidity and mortality rates have declined tremendously for burn patients, likely related to improvements in surgical and critical care treatment. We describe the baseline patient characteristics and well-defined outcomes for major burn injuries. METHODS: We followed 300 adults and 241 children from 2003 to 2009 through hospitalization, using standard operating procedures developed at study onset. We created an extensive database on patient and injury characteristics, anatomic and physiological derangement, clinical treatment, and outcomes. These data were compared with existing benchmarks in burn care. RESULTS: Study patients were critically injured, as demonstrated by mean % total body surface area (TBSA) (41.2 ± 18.3 for adults and 57.8 ± 18.2 for children) and presence of inhalation injury in 38% of the adults and 54.8% of the children. Mortality in adults was 14.1% for those younger than 55 years and 38.5% for those aged 55 years and older. Mortality in patients younger than 17 years was 7.9%. Overall, the multiple organ failure rate was 27%. When controlling for age and % TBSA, presence of inhalation injury continues to be significant. CONCLUSIONS: This study provides the current benchmark for major burn patients. Mortality rates, notwithstanding significant % TBSA and presence of inhalation injury, have significantly declined compared with previous benchmarks. Modern day surgical and medically intensive management has markedly improved to the point where we can expect patients younger than 55 years with severe burn injuries and inhalation injury to survive these devastating conditions.

Single cell mutational analysis of PIK3CA in circulating tumor cells and metastases in breast cancer reveals heterogeneity, discordance, and mutation persistence in cultured disseminated tumor cells from bone marrow.

BACKGROUND: Therapeutic decisions in cancer are generally guided by molecular biomarkers or, for some newer therapeutics, primary tumor genotype. However, because biomarkers or genotypes may change as new metastases emerge, circulating tumor cells (CTCs) from blood are being investigated for a role in guiding real-time drug selection during disease progression, expecting that CTCs will comprehensively represent the full spectrum of genomic changes in metastases. However, information is limited regarding mutational heterogeneity among CTCs and metastases in breast cancer as discerned by single cell analysis. The presence of disseminated tumor cells (DTCs) in bone marrow also carry prognostic significance in breast cancer, but with variability between CTC and DTC detection. Here we analyze a series of single tumor cells, CTCs, and DTCs for PIK3CA mutations and report CTC and corresponding metastatic genotypes. METHODS: We used the MagSweeper, an immunomagnetic separation device, to capture live single tumor cells from breast cancer patients' primary and metastatic tissues, blood, and bone marrow. Single cells were screened for mutations in exons 9 and 20 of the PIK3CA gene. Captured DTCs grown in cell culture were also sequenced for PIK3CA mutations. RESULTS: Among 242 individual tumor cells isolated from 17 patients and tested for mutations, 48 mutated tumor cells were identified in three patients. Single cell analyses revealed mutational heterogeneity among CTCs and tumor cells in tissues. In a patient followed serially, there was mutational discordance between CTCs, DTCs, and metastases, and among CTCs isolated at different time points. DTCs from this patient propagated in vitro contained a PIK3CA mutation, which was maintained despite morphological changes during 21 days of cell culture. CONCLUSIONS: Single cell analysis of CTCs can demonstrate genotypic heterogeneity, changes over time, and discordance from DTCs and distant metastases. We present a cautionary case showing that CTCs from any single blood draw do not always reflect metastatic genotype, and that CTC and DTC analyses may provide independent clinical information. Isolated DTCs remain viable and can be propagated in culture while maintaining their original mutational status, potentially serving as a future resource for investigating new drug therapies.

High-quality DNA sequence capture of 524 disease candidate genes.

The accurate and complete selection of candidate genomic regions from a DNA sample before sequencing is critical in molecular diagnostics. Several recently developed technologies await substantial improvements in performance, cost, and multiplex sample processing. Here we present the utility of long padlock probes (LPPs) for targeted exon capture followed by array-based sequencing. We found that on average 92% of 5,471 exons from 524 nuclear-encoded mitochondrial genes were successfully amplified from genomic DNA from 63 individuals. Only 144 exons did not amplify in any sample due to high GC content. One LPP was sufficient to capture sequences from <100-500 bp in length and only a single-tube capture reaction and one microarray was required per sample. Our approach was highly reproducible and quick (<8 h) and detected DNA variants at high accuracy (false discovery rate 1%, false negative rate 3%) on the basis of known sample SNPs and Sanger sequence verification. In a patient with clinical and biochemical presentation of ornithine transcarbamylase (OTC) deficiency, we identified copy-number differences in the OTC gene at exon-level resolution. This shows the ability of LPPs to accurately preserve a sample's genome information and provides a cost-effective strategy to identify both single nucleotide changes and structural variants in targeted resequencing.

Human transcriptome array for high-throughput clinical studies.

A 6.9 million-feature oligonucleotide array of the human transcriptome [Glue Grant human transcriptome (GG-H array)] has been developed for high-throughput and cost-effective analyses in clinical studies. This array allows comprehensive examination of gene expression and genome-wide identification of alternative splicing as well as detection of coding SNPs and noncoding transcripts. The performance of the array was examined and compared with mRNA sequencing (RNA-Seq) results over multiple independent replicates of liver and muscle samples. Compared with RNA-Seq of 46 million uniquely mappable reads per replicate, the GG-H array is highly reproducible in estimating gene and exon abundance. Although both platforms detect similar expression changes at the gene level, the GG-H array is more sensitive at the exon level. Deeper sequencing is required to adequately cover low-abundance transcripts. The array has been implemented in a multicenter clinical program and has generated high-quality, reproducible data. Considering the clinical trial requirements of cost, sample availability, and throughput, the GG-H array has a wide range of applications. An emerging approach for large-scale clinical genomic studies is to first use RNA-Seq to the sufficient depth for the discovery of transcriptome elements relevant to the disease process followed by high-throughput and reliable screening of these elements on thousands of patient samples using custom-designed arrays.

High-throughput VDJ sequencing for quantification of minimal residual disease in chronic lymphocytic leukemia and immune reconstitution assessment.

The primary cause of poor outcome following allogeneic hematopoietic cell transplantation (HCT) for chronic lymphocytic leukemia (CLL) is disease recurrence. Detection of increasing minimal residual disease (MRD) following HCT may permit early intervention to prevent clinical relapse; however, MRD quantification remains an uncommon diagnostic test because of logistical and financial barriers to widespread use. Here we describe a method for quantifying CLL MRD using widely available consensus primers for amplification of all Ig heavy chain (IGH) genes in a mixture of peripheral blood mononuclear cells, followed by high-throughput sequencing (HTS) for disease-specific IGH sequence quantification. To achieve accurate MRD quantification, we developed a systematic bioinformatic methodology to aggregate cancer clone sequence variants arising from systematic and random artifacts occurring during IGH-HTS. We then compared the sensitivity of IGH-HTS, flow cytometry, and allele-specific oligonucleotide PCR for MRD quantification in 28 samples collected from 6 CLL patients following allogeneic HCT. Using amplimer libraries generated with consensus primers from patient blood samples, we demonstrate the sensitivity of IGH-HTS with 454 pyrosequencing to be 10(-5), with a high correlation between quantification by allele-specific oligonucleotide PCR and IGH-HTS (r = 0.85). From the same dataset used to quantify MRD, IGH-HTS also allowed us to profile IGH repertoire reconstitution after HCT-information not provided by the other MRD methods. IGH-HTS using consensus primers will broaden the availability of MRD quantification in CLL and other B cell malignancies, and this approach has potential for quantitative evaluation of immune diversification following transplant and nontransplant therapies.

Development of a genomic metric that can be rapidly used to predict clinical outcome in severely injured trauma patients.

OBJECTIVE: Many patients have complicated recoveries following severe trauma due to the development of organ injury. Physiological and anatomical prognosticators have had limited success in predicting clinical trajectories. We report on the development and retrospective validation of a simple genomic composite score that can be rapidly used to predict clinical outcomes. DESIGN: Retrospective cohort study. SETTING: Multi-institutional level 1 trauma centers. PATIENTS: Data were collected from 167 severely traumatized (injury severity score >15) adult (18-55 yr) patients. METHODS: Microarray-derived genomic data obtained from 167 severely traumatized patients over 28 days were assessed for differences in messenger RNA abundance among individuals with different clinical trajectories. Once a set of genes was identified based on differences in expression over the entire study period, messenger RNA abundance from these subjects obtained in the first 24 hours was analyzed in a blinded fashion using a rapid multiplex platform, and genomic data reduced to a single metric. RESULTS: From the existing genomic dataset, we identified 63 genes whose leukocyte expression differed between an uncomplicated and complicated clinical outcome over 28 days. Using a multiplex approach that can quantitate messenger RNA abundance in less than 12 hours, we reassessed total messenger RNA abundance from the first 24 hours after trauma and reduced the genomic data to a single composite score using the difference from reference. This composite score showed good discriminatory capacity to distinguish patients with a complicated outcome (area under a receiver-operator curve, 0.811; p <0.001). This was significantly better than the predictive power of either Acute Physiology and Chronic Health Evaluation II or new injury severity score scoring systems. CONCLUSIONS: A rapid genomic composite score obtained in the first 24 hours after trauma can retrospectively identify trauma patients who are likely to develop complicated clinical trajectories. A novel platform is described in which this genomic score can be obtained within 12 hours of blood collection, making it available for clinical decision making.

Identification of rare DNA variants in mitochondrial disorders with improved array-based sequencing.

A common goal in the discovery of rare functional DNA variants via medical resequencing is to incur a relatively lower proportion of false positive base-calls. We developed a novel statistical method for resequencing arrays (SRMA, sequence robust multi-array analysis) to increase the accuracy of detecting rare variants and reduce the costs in subsequent sequence verifications required in medical applications. SRMA includes single and multi-array analysis and accounts for technical variables as well as the possibility of both low- and high-frequency genomic variation. The confidence of each base-call was ranked using two quality measures. In comparison to Sanger capillary sequencing, we achieved a false discovery rate of 2% (false positive rate 1.2 × 10⁻5, false negative rate 5%), which is similar to automated second-generation sequencing technologies. Applied to the analysis of 39 nuclear candidate genes in disorders of mitochondrial DNA (mtDNA) maintenance, we confirmed mutations in the DNA polymerase gamma POLG in positive control cases, and identified novel rare variants in previously undiagnosed cases in the mitochondrial topoisomerase TOP1MT, the mismatch repair enzyme MUTYH, and the apurinic-apyrimidinic endonuclease APEX2. Some patients carried rare heterozygous variants in several functionally interacting genes, which could indicate synergistic genetic effects in these clinically similar disorders.

High-throughput method for analyzing methylation of CpGs in targeted genomic regions.

A unique microarray-based method for determining the extent of DNA methylation has been developed. It relies on a selective enrichment of the regions to be assayed by target amplification by capture and ligation (mTACL). The assay is quantitatively accurate, relatively precise, and lends itself to high-throughput determination using nanogram amounts of DNA. The measurements using mTACLs are highly reproducible and in excellent agreement with those obtained by sequencing (r = 0.94). In the present work, the methylation status of >145,000 CpGs from 5,472 promoters in 221 samples was measured. The methylation levels of nearby CpGs are correlated, but the correlation falls off dramatically over several hundred base pairs. In some instances, nearby CpGs have very different levels of methylation. Comparison of normal and tumor samples indicates that in tumors, the promoter regions of genes involved in differentiation and signaling are preferentially hypermethylated, whereas those of housekeeping genes remain hypomethylated. mTACL is a platform for profiling the state of methylation of a large number of CpG in many samples in a cost-effective fashion, and is capable of scaling to much larger numbers of CpGs than those collected here.

Analysis of factorial time-course microarrays with application to a clinical study of burn injury.

Time-course microarray experiments are capable of capturing dynamic gene expression profiles. It is important to study how these dynamic profiles depend on the multiple factors that characterize the experimental condition under which the time course is observed. Analytic methods are needed to simultaneously handle the time course and factorial structure in the data. We developed a method to evaluate factor effects by pooling information across the time course while accounting for multiple testing and nonnormality of the microarray data. The method effectively extracts gene-specific response features and models their dependency on the experimental factors. Both longitudinal and cross-sectional time-course data can be handled by our approach. The method was used to analyze the impact of age on the temporal gene response to burn injury in a large-scale clinical study. Our analysis reveals that 21% of the genes responsive to burn are age-specific, among which expressions of mitochondria and immunoglobulin genes are differentially perturbed in pediatric and adult patients by burn injury. These new findings in the body's response to burn injury between children and adults support further investigations of therapeutic options targeting specific age groups. The methodology proposed here has been implemented in R package "TANOVA" and submitted to the Comprehensive R Archive Network at http://www.r-project.org/. It is also available for download at http://gluegrant1.stanford.edu/TANOVA/.

Benchmarking outcomes in the critically injured trauma patient and the effect of implementing standard operating procedures.

OBJECTIVE: To determine and compare outcomes with accepted benchmarks in trauma care at 7 academic level I trauma centers in which patients were treated on the basis of a series of standard operating procedures (SOPs). BACKGROUND: Injury remains the leading cause of death for those younger than 45 years. This study describes the baseline patient characteristics and well-defined outcomes of persons hospitalized in the United States for severe blunt trauma. METHODS: We followed 1637 trauma patients from 2003 to 2009 up to 28 hospital days using SOPs developed at the onset of the study. An extensive database on patient and injury characteristics, clinical treatment, and outcomes was created. These data were compared with existing trauma benchmarks. RESULTS: The study patients were critically injured and were in shock. SOP compliance improved 10% to 40% during the study period. Multiple organ failure and mortality rates were 34.8% and 16.7%, respectively. Time to recovery, defined as the time until the patient was free of organ failure for at least 2 consecutive days, was developed as a new outcome measure. There was a reduction in mortality rate in the cohort during the study that cannot be explained by changes in the patient population. CONCLUSIONS: This study provides the current benchmark and the overall positive effect of implementing SOPs for severely injured patients. Over the course of the study, there were improvements in morbidity and mortality rates and increasing compliance with SOPs. Mortality was surprisingly low, given the degree of injury, and improved over the duration of the study, which correlated with improved SOP compliance.

Antibodies specifically target AML antigen NuSAP1 after allogeneic bone marrow transplantation.

Identifying the targets of immune response after allogeneic hematopoietic cell transplantation (HCT) promises to provide relevant immune therapy candidate proteins. We used protein microarrays to serologically identify nucleolar and spindle-associated protein 1 (NuSAP1) and chromatin assembly factor 1, subunit B (p60; CHAF1b) as targets of new antibody responses that developed after allogeneic HCT. Western blots and enzyme-linked immunosorbent assays (ELISA) validated their post-HCT recognition and enabled ELISA testing of 120 other patients with various malignancies who underwent allo-HCT. CHAF1b-specific antibodies were predominantly detected in patients with acute myeloid leukemia (AML), whereas NuSAP1-specific antibodies were exclusively detected in patients with AML 1 year after transplantation (P < .001). Complete genomic exon sequencing failed to identify a nonsynonymous single nucleotide polymorphism (SNP) for NuSAP1 and CHAF1b between the donor and recipient cells. Expression profiles and reverse transcriptase-polymerase chain reaction (RT-PCR) showed NuSAP1 was predominately expressed in the bone marrow CD34(+)CD90(+) hematopoietic stem cells, leukemic cell lines, and B lymphoblasts compared with other tissues or cells. Thus, NuSAP1 is recognized as an immunogenic antigen in 65% of patients with AML following allogeneic HCT and suggests a tumor antigen role.

Loss of the plant DEAD-box protein ISE1 leads to defective mitochondria and increased cell-to-cell transport via plasmodesmata.

Plants have intercellular channels, plasmodesmata (PD), that span the cell wall to enable cell-to-cell transport of micro- and macromolecules. We identified an Arabidopsis thaliana embryo lethal mutant increased size exclusion limit 1 (ise1) that results in increased PD-mediated transport of fluorescent tracers. The ise1 mutants have a higher frequency of branched and twinned PD than wild-type embryos. Silencing of ISE1 in mature Nicotiana benthamiana leaves also leads to increased PD transport, as monitored by intercellular movement of a GFP fusion to the tobacco mosaic virus movement protein. ISE1 encodes a putative plant-specific DEAD-box RNA helicase that localizes specifically to mitochondria. The N-terminal 100 aa of ISE1 specify mitochondrial targeting. Mitochondrial metabolism is compromised severely in ise1 mutant embryos, because their mitochondrial proton gradient is disrupted and reactive oxygen species production is increased. Although mitochondria are essential for numerous cell-autonomous functions, the present studies demonstrate that mitochondrial function also regulates the critical cell non-cell-autonomous function of PD.

Isolating highly enriched populations of circulating epithelial cells and other rare cells from blood using a magnetic sweeper device.

The enumeration of rare circulating epithelial cells (CEpCs) in the peripheral blood of metastatic cancer patients has shown promise for improved cancer prognosis. Moving beyond enumeration, molecular analysis of CEpCs may provide candidate surrogate endpoints to diagnose, treat, and monitor malignancy directly from the blood samples. Thorough molecular analysis of CEpCs requires the development of new sample preparation methods that yield easily accessible and purified CEpCs for downstream biochemical assays. Here, we describe a new immunomagnetic cell separator, the MagSweeper, which gently enriches target cells and eliminates cells that are not bound to magnetic particles. The isolated cells are easily accessible and can be extracted individually based on their physical characteristics to deplete any cells nonspecifically bound to beads. We have shown that our device can process 9 mL of blood per hour and captures >50% of CEpCs as measured in spiking experiments. We have shown that the separation process does not perturb the gene expression of rare cells. To determine the efficiency of our platform in isolating CEpCs from patients, we have isolated CEpCs from all 47 tubes of 9-mL blood samples collected from 17 women with metastatic breast cancer. In contrast, we could not find any circulating epithelial cells in samples from 5 healthy donors. The isolated CEpCs are all stored individually for further molecular analysis.