Capture-based targeted sequencing using a T-cell control in myeloid malignancies and idiopathic cytopenias.

We report on a study of next-generation sequencing in 257 patients undergoing investigations for cytopenias. We sequenced bone marrow aspirates using a target enrichment panel comprising 82 genes and used T cells from paired blood as a control. One hundred and sixty patients had idiopathic cytopenias, 81 had myeloid malignancies and 16 had lymphoid malignancies or other diagnoses. Forty-seven of the 160 patients with idiopathic cytopenias had evidence of somatic pathogenic variants consistent with clonal cytopenias. Only 39 genes of the 82 tested were mutated in the 241 patients with either idiopathic cytopenias or myeloid neoplasms. We confirm that T cells can be used as a control to distinguish between germline and somatic variants. The use of paired analysis with a T-cell control significantly reduced the time molecular scientists spent reporting compared to unpaired analysis. We identified somatic variants of uncertain significance (VUS) in a higher proportion (24%) of patients with myeloid malignancies or clonal cytopenias compared to less than 2% of patients with non-clonal cytopenias. This suggests that somatic VUS are indicators of a clonal process. Lastly, we show that blood depleted of lymphocytes can be used in place of bone marrow as a source of material for sequencing.

DNA methylation-based classification of glioneuronal tumours synergises with histology and radiology to refine accurate molecular stratification.

AIMS: Glioneuronal tumours (GNTs) are poorly distinguished by their histology and lack robust diagnostic indicators. Previously, we showed that common GNTs comprise two molecularly distinct groups, correlating poorly with histology. To refine diagnosis, we constructed a methylation-based model for GNT classification, subsequently evaluating standards for molecular stratification by methylation, histology and radiology. METHODS: We comprehensively analysed methylation, radiology and histology for 83 GNT samples: a training cohort of 49, previously classified into molecularly defined groups by genomic profiles, plus a validation cohort of 34. We identified histological and radiological correlates to molecular classification and constructed a methylation-based support vector machine (SVM) model for prediction. Subsequently, we contrasted methylation, radiological and histological classifications in validation GNTs. RESULTS: By methylation clustering, all training and 23/34 validation GNTs segregated into two groups, the remaining 11 clustering alongside control cortex. Histological review identified prominent astrocytic/oligodendrocyte-like components, dysplastic neurons and a specific glioneuronal element as discriminators between groups. However, these were present in only a subset of tumours. Radiological review identified location, margin definition, enhancement and T2 FLAIR-rim sign as discriminators. When validation GNTs were classified by SVM, 22/23 classified correctly, comparing favourably against histology and radiology that resolved 17/22 and 15/21, respectively, where data were available for comparison. CONCLUSIONS: Diagnostic criteria inadequately reflect glioneuronal tumour biology, leaving a proportion unresolvable. In the largest cohort of molecularly defined glioneuronal tumours, we develop molecular, histological and radiological approaches for biologically meaningful classification and demonstrate almost all cases are resolvable, emphasising the importance of an integrated diagnostic approach.

Mutations in ALK signaling pathways conferring resistance to ALK inhibitor treatment lead to collateral vulnerabilities in neuroblastoma cells.

BACKGROUND: Development of resistance to targeted therapies has tempered initial optimism that precision oncology would improve poor outcomes for cancer patients. Resistance mechanisms, however, can also confer new resistance-specific vulnerabilities, termed collateral sensitivities. Here we investigated anaplastic lymphoma kinase (ALK) inhibitor resistance in neuroblastoma, a childhood cancer frequently affected by activating ALK alterations. METHODS: Genome-wide forward genetic CRISPR-Cas9 based screens were performed to identify genes associated with ALK inhibitor resistance in neuroblastoma cell lines. Furthermore, the neuroblastoma cell line NBLW-R was rendered resistant by continuous exposure to ALK inhibitors. Genes identified to be associated with ALK inhibitor resistance were further investigated by generating suitable cell line models. In addition, tumor and liquid biopsy samples of four patients with ALK-mutated neuroblastomas before ALK inhibitor treatment and during tumor progression under treatment were genomically profiled. RESULTS: Both genome-wide CRISPR-Cas9-based screens and preclinical spontaneous ALKi resistance models identified NF1 loss and activating NRASQ61K mutations to confer resistance to chemically diverse ALKi. Moreover, human neuroblastomas recurrently developed de novo loss of NF1 and activating RAS mutations after ALKi treatment, leading to therapy resistance. Pathway-specific perturbations confirmed that NF1 loss and activating RAS mutations lead to RAS-MAPK signaling even in the presence of ALKi. Intriguingly, NF1 loss rendered neuroblastoma cells hypersensitive to MEK inhibition. CONCLUSIONS: Our results provide a clinically relevant mechanistic model of ALKi resistance in neuroblastoma and highlight new clinically actionable collateral sensitivities in resistant cells.

High-throughput sequencing of peripheral blood for minimal residual disease monitoring in childhood precursor B-cell acute lymphoblastic leukemia: A prospective feasibility study.

BACKGROUND: Minimal residual disease (MRD) measured on end-of-induction bone marrow (BM) is the most important biomarker for guiding therapy in pediatric acute lymphoblastic leukemia (ALL). Due to limited sensitivity of current approaches, peripheral blood (PB) is not a reliable source for identifying patients needing treatment changes. We sought to determine if high-throughput sequencing (HTS) (next-generation sequencing) of rearranged immunoglobulin and T-cell receptor genes can overcome this and be used to measure MRD in PB. PROCEDURE: We employed a quantitative HTS approach to accurately measure MRD from one million cell equivalents of DNA from 17 PB samples collected at day 29 after induction therapy in patients with precursor B-cell ALL. We compared these results to the gold-standard real-time PCR result obtained from their paired BM samples, median follow-up 49 months. RESULTS: With the increased sensitivity, detecting up to one abnormal cell in a million normal cells, we were able to detect MRD in the PB by HTS in all those patients requiring treatment intensification (MRD ≥ 0.005% in BM). CONCLUSION: This is proof of principle that using the increased sensitivity of HTS, PB can be used to measure MRD and stratify children with ALL. The method is cost effective, rapid, accurate, and reproducible, with inherent advantages in children. Importantly, increasing the frequency testing by PB as opposed to intermittent BM sampling may allow extension of the dynamic range of MRD, giving a more complete picture of the kinetics of disease remission while improving relapse prediction and speed of detection.

Comprehensive molecular characterisation of epilepsy-associated glioneuronal tumours.

Glioneuronal tumours are an important cause of treatment-resistant epilepsy. Subtypes of tumour are often poorly discriminated by histological features and may be difficult to diagnose due to a lack of robust diagnostic tools. This is illustrated by marked variability in the reported frequencies across different epilepsy surgical series. To address this, we used DNA methylation arrays and RNA sequencing to assay the methylation and expression profiles within a large cohort of glioneuronal tumours. By adopting a class discovery approach, we were able to identify two distinct groups of glioneuronal tumour, which only partially corresponded to the existing histological classification. Furthermore, by additional molecular analyses, we were able to identify pathogenic mutations in BRAF and FGFR1, specific to each group, in a high proportion of cases. Finally, by interrogating our expression data, we were able to show that each molecular group possessed expression phenotypes suggesting different cellular differentiation: astrocytic in one group and oligodendroglial in the second. Informed by this, we were able to identify CCND1, CSPG4, and PDGFRA as immunohistochemical targets which could distinguish between molecular groups. Our data suggest that the current histological classification of glioneuronal tumours does not adequately represent their underlying biology. Instead, we show that there are two molecular groups within glioneuronal tumours. The first of these displays astrocytic differentiation and is driven by BRAF mutations, while the second displays oligodendroglial differentiation and is driven by FGFR1 mutations.

Deep sequencing reveals persistence of cell-associated mumps vaccine virus in chronic encephalitis.

Routine childhood vaccination against measles, mumps and rubella has virtually abolished virus-related morbidity and mortality. Notwithstanding this, we describe here devastating neurological complications associated with the detection of live-attenuated mumps virus Jeryl Lynn (MuVJL5) in the brain of a child who had undergone successful allogeneic transplantation for severe combined immunodeficiency (SCID). This is the first confirmed report of MuVJL5 associated with chronic encephalitis and highlights the need to exclude immunodeficient individuals from immunisation with live-attenuated vaccines. The diagnosis was only possible by deep sequencing of the brain biopsy. Sequence comparison of the vaccine batch to the MuVJL5 isolated from brain identified biased hypermutation, particularly in the matrix gene, similar to those found in measles from cases of SSPE. The findings provide unique insights into the pathogenesis of paramyxovirus brain infections.

Cell cycle status in AML blast cells from peripheral blood, bone marrow aspirates and trephines and implications for biological studies and treatment.

Using immunohistochemistry and flow cytometry to define phases of the cell cycle, this study shows that a high proportion of acute myeloid leukaemia (AML) blasts obtained from trephine biopsies are cycling, whereas >95% of peripheral blood-derived blasts are arrested in G1 . Results obtained from bone marrow aspirates are more similar to those from blood rather than from trephine biopsies. These differences were confirmed by gene expression profiling in a patient with high count AML. This has implications for cell cycle and other biological studies using aspirates rather than trephine biopsies and for the use of cell mobilising agents before chemotherapy.

The kinetochore protein, CENPF, is mutated in human ciliopathy and microcephaly phenotypes.

BACKGROUND: Mutations in microtubule-regulating genes are associated with disorders of neuronal migration and microcephaly. Regulation of centriole length has been shown to underlie the pathogenesis of certain ciliopathy phenotypes. Using a next-generation sequencing approach, we identified mutations in a novel centriolar disease gene in a kindred with an embryonic lethal ciliopathy phenotype and in a patient with primary microcephaly. METHODS AND RESULTS: Whole exome sequencing data from a non-consanguineous Caucasian kindred exhibiting mid-gestation lethality and ciliopathic malformations revealed two novel non-synonymous variants in CENPF, a microtubule-regulating gene. All four affected fetuses showed segregation for two mutated alleles [IVS5-2A>C, predicted to abolish the consensus splice-acceptor site from exon 6; c.1744G>T, p.E582X]. In a second unrelated patient exhibiting microcephaly, we identified two CENPF mutations [c.1744G>T, p.E582X; c.8692 C>T, p.R2898X] by whole exome sequencing. We found that CENP-F colocalised with Ninein at the subdistal appendages of the mother centriole in mouse inner medullary collecting duct cells. Intraflagellar transport protein-88 (IFT-88) colocalised with CENP-F along the ciliary axonemes of renal epithelial cells in age-matched control human fetuses but did not in truncated cilia of mutant CENPF kidneys. Pairwise co-immunoprecipitation assays of mitotic and serum-starved HEKT293 cells confirmed that IFT88 precipitates with endogenous CENP-F. CONCLUSIONS: Our data identify CENPF as a new centriolar disease gene implicated in severe human ciliopathy and microcephaly related phenotypes. CENP-F has a novel putative function in ciliogenesis and cortical neurogenesis.

Astrovirus VA1/HMO-C: an increasingly recognized neurotropic pathogen in immunocompromised patients.

BACKGROUND: An 18-month-old boy developed encephalopathy, for which extensive investigation failed to identify an etiology, 6 weeks after stem cell transplant. To exclude a potential infectious cause, we performed high-throughput RNA sequencing on brain biopsy. METHODS: RNA-Seq was performed on an Illumina Miseq, generating 20 million paired-end reads. Nonhost data were checked for similarity to known organisms using BLASTx. The full viral genome was sequenced by primer walking. RESULTS: We identified an astrovirus, HAstV-VA1/HMO-C-UK1(a), which was highly divergent from human astrovirus (HAstV 1-8) genotypes, but closely related to VA1/HMO-C astroviruses, including one recovered from a case of fatal encephalitis in an immunosuppressed child. The virus was detected in stool and serum, with highest levels in brain and cerebrospinal fluid (CSF). Immunohistochemistry of the brain biopsy showed positive neuronal staining. A survey of 680 stool and 349 CSF samples identified a related virus in the stool of another immunosuppressed child. CONCLUSIONS: The discovery of HAstV-VA1/HMO-C-UK1(a) as the cause of encephalitis in this case provides further evidence that VA1/HMO-C viruses, unlike HAstV 1-8, are neuropathic, particularly in immunocompromised patients, and should be considered in the differential diagnosis of encephalopathy. With a turnaround from sample receipt to result of <1 week, we confirm that RNA-Seq presents a valuable diagnostic tool in unexplained encephalitis.

Capturing needles in haystacks: a comparison of B-cell receptor sequencing methods.

BACKGROUND: Deep-sequencing methods are rapidly developing in the field of B-cell receptor (BCR) and T-cell receptor (TCR) diversity. These promise to revolutionise our understanding of adaptive immune dynamics, identify novel antibodies, and allow monitoring of minimal residual disease. However, different methods for BCR and TCR enrichment and amplification have been proposed. Here we perform the first systematic comparison between different methods of enrichment, amplification and sequencing for generating BCR and TCR repertoires using large sample numbers. RESULTS: Resampling from the same RNA or cDNA pool results in highly correlated and reproducible repertoires, but resampling low frequency clones leads to stochastic variance. Repertoires generated by different sequencing methods (454 Roche and Illumina MiSeq) and amplification methods (multiplex PCR, 5' Rapid amplification of cDNA ends (5'RACE), and RNA-capture) are highly correlated, and resulting IgHV gene frequencies between the different methods were not significantly different. Read length has an impact on captured repertoire structure, and ultimately full-length BCR sequences are most informative for repertoire analysis as diversity outside of the CDR is very useful for phylogenetic analysis. Additionally, we show RNA-based BCR repertoires are more informative than using DNA. CONCLUSIONS: Repertoires generated by different sequencing and amplification methods are consistent, but we show that read lengths, depths and error profiles should be considered in experimental design, and multiple sampling approaches could be employed to minimise stochastic sampling variation. This detailed investigation of immune repertoire sequencing methods is essential for informing basic and clinical research.

Risk HLA-DQA1 and PLA(2)R1 alleles in idiopathic membranous nephropathy.

BACKGROUND: Idiopathic membranous nephropathy is a major cause of the nephrotic syndrome in adults, but its etiologic basis is not fully understood. We investigated the genetic basis of biopsy-proven cases of idiopathic membranous nephropathy in a white population. METHODS: We performed independent genomewide association studies of single-nucleotide polymorphisms (SNPs) in patients with idiopathic membranous nephropathy from three populations of white ancestry (75 French, 146 Dutch, and 335 British patients). The patients were compared with racially matched control subjects; population stratification and quality controls were carried out according to standard criteria. Associations were calculated by means of a chi-square basic allele test; the threshold for significance was adjusted for multiple comparisons (with the Bonferroni method). RESULTS: In a joint analysis of data from the 556 patients studied (398 men), we identified significant alleles at two genomic loci associated with idiopathic membranous nephropathy. Chromosome 2q24 contains the gene encoding M-type phospholipase A(2) receptor (PLA(2)R1) (SNP rs4664308, P=8.6×10(-29)), previously shown to be the target of an autoimmune response. Chromosome 6p21 contains the gene encoding HLA complex class II HLA-DQ alpha chain 1 (HLA-DQA1) (SNP rs2187668, P=8.0×10(-93)). The association with HLA-DQA1 was significant in all three populations (P=1.8×10(-9), P=5.6×10(-27), and P=5.2×10(-36) in the French, Dutch, and British groups, respectively). The odds ratio for idiopathic membranous nephropathy with homozygosity for both risk alleles was 78.5 (95% confidence interval, 34.6 to 178.2). CONCLUSIONS: An HLA-DQA1 allele on chromosome 6p21 is most closely associated with idiopathic membranous nephropathy in persons of white ancestry. This allele may facilitate an autoimmune response against targets such as variants of PLA2R1. Our findings suggest a basis for understanding this disease and illuminate how adaptive immunity is regulated by HLA.

Genetic complexity in hypertrophic cardiomyopathy revealed by high-throughput sequencing.

BACKGROUND: Clinical interpretation of the large number of rare variants identified by high throughput sequencing (HTS) technologies is challenging. The aim of this study was to explore the clinical implications of a HTS strategy for patients with hypertrophic cardiomyopathy (HCM) using a targeted HTS methodology and workflow developed for patients with a range of inherited cardiovascular diseases. By comparing the sequencing results with published findings and with sequence data from a large-scale exome sequencing screen of UK individuals, we sought to quantify the strength of the evidence supporting causality for detected candidate variants. METHODS AND RESULTS: 223 unrelated patients with HCM (46±15 years at diagnosis, 74% males) were studied. In order to analyse coding, intronic and regulatory regions of 41 cardiovascular genes, we used solution-based sequence capture followed by massive parallel resequencing on Illumina GAIIx. Average read-depth in the 2.1 Mb target region was 120. Rare (frequency<0.5%) non-synonymous, loss-of-function and splice-site variants were defined as candidates. Excluding titin, we identified 152 distinct candidate variants in sarcomeric or associated genes (89 novel) in 143 patients (64%). Four sarcomeric genes (MYH7, MYBPC3, TNNI3, TNNT2) showed an excess of rare single non-synonymous single-nucleotide polymorphisms (nsSNPs) in cases compared to controls. The estimated probability that a nsSNP in these genes is pathogenic varied between 57% and near certainty depending on the location. We detected an additional 94 candidate variants (73 novel) in desmosomal, and ion-channel genes in 96 patients (43%). CONCLUSIONS: This study provides the first large-scale quantitative analysis of the prevalence of sarcomere protein gene variants in patients with HCM using HTS technology. Inclusion of other genes implicated in inherited cardiac disease identifies a large number of non-synonymous rare variants of unknown clinical significance.

Albuminuria is associated with too few glomeruli and too much testosterone.

Normally, the glomerular filtration barrier almost completely excludes circulating albumin from entering the urine. Genetic variation and both pre- and postnatal environmental factors may affect albuminuria in humans. Here we determine whether glomerular gene expression in mouse strains with naturally occurring variations in albuminuria would allow identification of proteins deregulated in relatively 'leaky' glomeruli. Albuminuria increased in female B6 to male B6 to female FVB/N to male FVB/N mice, whereas the number of glomeruli/kidney was the exact opposite. Testosterone administration led to increased albuminuria in female B6 but not female FVB/N mice. A common set of 39 genes, many expressed in podocytes, were significantly differentially expressed in each of the four comparisons: male versus female B6 mice, male versus female FVB/N mice, male FVB/N versus male B6 mice, and female FVB/N versus female B6 mice. The transcripts encoded proteins involved in oxidation/reduction reactions, ion transport, and enzymes involved in detoxification. These proteins may represent novel biomarkers and even therapeutic targets for early kidney and cardiovascular disease.

Epilepsy, ataxia, sensorineural deafness, tubulopathy, and KCNJ10 mutations.

BACKGROUND: Five children from two consanguineous families presented with epilepsy beginning in infancy and severe ataxia, moderate sensorineural deafness, and a renal salt-losing tubulopathy with normotensive hypokalemic metabolic alkalosis. We investigated the genetic basis of this autosomal recessive disease, which we call the EAST syndrome (the presence of epilepsy, ataxia, sensorineural deafness, and tubulopathy). METHODS: Whole-genome linkage analysis was performed in the four affected children in one of the families. Newly identified mutations in a potassium-channel gene were evaluated with the use of a heterologous expression system. Protein expression and function were further investigated in genetically modified mice. RESULTS: Linkage analysis identified a single significant locus on chromosome 1q23.2 with a lod score of 4.98. This region contained the KCNJ10 gene, which encodes a potassium channel expressed in the brain, inner ear, and kidney. Sequencing of this candidate gene revealed homozygous missense mutations in affected persons in both families. These mutations, when expressed heterologously in xenopus oocytes, caused significant and specific decreases in potassium currents. Mice with Kcnj10 deletions became dehydrated, with definitive evidence of renal salt wasting. CONCLUSIONS: Mutations in KCNJ10 cause a specific disorder, consisting of epilepsy, ataxia, sensorineural deafness, and tubulopathy. Our findings indicate that KCNJ10 plays a major role in renal salt handling and, hence, possibly also in blood-pressure maintenance and its regulation.

ChIP-on-chip analysis reveals angiopoietin 2 (Ang2, ANGPT2) as a novel target of steroidogenic factor-1 (SF-1, NR5A1) in the human adrenal gland.

The nuclear receptor steroidogenic factor-1 (SF-1, NR5A1) is a key regulator of adrenal and gonadal biology. Disruption of SF-1 can lead to disorders of adrenal development, while increased SF-1 dosage has been associated with adrenocortical tumorigenesis. We aimed to identify a novel subset of SF-1 target genes in the adrenal by using chromatin immunoprecipitation (ChIP) microarrays (ChIP-on-chip) combined with systems analysis. SF-1 ChIP-on-chip was performed in NCI-H295R human adrenocortical cells using promoter tiling arrays, leading to the identification of 445 gene loci where SF-1-binding regions were located from 10 kb upstream to 3 kb downstream of a transcriptional start. Network analysis of genes identified as putative SF-1 targets revealed enrichment for angiogenic process networks. A 1.1-kb SF-1-binding region was identified in the angiopoietin 2 (Ang2, ANGPT2) promoter in a highly repetitive region, and SF-1-dependent activation was confirmed in luciferase assays. Angiogenesis is paramount in adrenal development and tumorigenesis, but until now a direct link between SF-1 and vascular remodeling has not been established. We have identified Ang2 as a potentially important novel target of SF-1 in the adrenal gland, indicating that regulation of angiogenesis might be an important additional mechanism by which SF-1 exerts its actions in the adrenal gland.

Global analysis of gene expression in NGF-deprived sympathetic neurons identifies molecular pathways associated with cell death.

BACKGROUND: Developing sympathetic neurons depend on nerve growth factor (NGF) for survival and die by apoptosis after NGF withdrawal. This process requires de novo gene expression but only a small number of genes induced by NGF deprivation have been identified so far, either by a candidate gene approach or in mRNA differential display experiments. This is partly because it is difficult to obtain large numbers of sympathetic neurons for in vitro studies. Here, we describe for the first time, how advances in gene microarray technology have allowed us to investigate the expression of all known genes in sympathetic neurons cultured in the presence and absence of NGF. RESULTS: We have used Affymetrix Exon arrays to study the pattern of expression of all known genes in NGF-deprived sympathetic neurons. We identified 415 up- and 813 down-regulated genes, including most of the genes previously known to be regulated in this system. NGF withdrawal activates the mixed lineage kinase (MLK)-c-Jun N-terminal kinase (JNK)-c-Jun pathway which is required for NGF deprivation-induced death. By including a mixed lineage kinase (MLK) inhibitor, CEP-11004, in our experimental design we identified which of the genes induced after NGF withdrawal are potential targets of the MLK-JNK-c-Jun pathway. A detailed Gene Ontology and functional enrichment analysis also identified genetic pathways that are highly enriched and overrepresented amongst the genes expressed after NGF withdrawal. Five genes not previously studied in sympathetic neurons - trib3, ddit3, txnip, ndrg1 and mxi1 - were validated by real time-PCR. The proteins encoded by these genes also increased in level after NGF withdrawal and this increase was prevented by CEP-11004, suggesting that these genes are potential targets of the MLK-JNK-c-Jun pathway. CONCLUSIONS: The sympathetic neuron model is one of the best studied models of neuronal apoptosis. Overall, our microarray data gives a comprehensive overview of, and provides new information about, signalling pathways and transcription factors that are regulated by NGF withdrawal.

Coordinated oncogenic transformation and inhibition of host immune responses by the PAX3-FKHR fusion oncoprotein.

Tumors have evolved elaborate mechanisms for evading immune detection, such as production of immunoinhibitory cytokines and down-regulation of major histocompatibility complex (MHC) expression. We have studied PAX3-FKHR as an example of an oncogenic fusion protein associated with an aggressive metastatic cancer. We show that PAX3-FKHR alters expression of genes that are normally regulated by Janus kinase/signal transducer and activator of transcription (STAT) signaling pathways. This occurs as a result of a specific interaction between PAX3-FKHR and the STAT3 transcription factor, which results in a dramatic reduction in tumor MHC expression, and an alteration in local cytokine concentrations to inhibit surrounding inflammatory cells and immune detection. Collectively, these data show that an oncogenic transcription factor can promote tumor growth and tissue invasion while inhibiting local inflammatory and immune responses. This is the first time that an immunomodulatory role has been described for an oncogenic fusion protein.

Biologic predictors of extension of oligoarticular juvenile idiopathic arthritis as determined from synovial fluid cellular composition and gene expression.

OBJECTIVE: To identify biomarkers in the first synovial fluid (SF) aspirate obtained from children with oligoarticular juvenile idiopathic arthritis (JIA), which could be used to identify children whose disease is likely to extend to a more severe phenotype. METHODS: Patients with recent-onset oligoarticular JIA were identified and grouped according to those whose mild disease persisted (persistent disease) or those whose disease would extend from a mild to more severe phenotype (extended-to-be disease) at 1 year after diagnosis. Flow cytometry was used to delineate differences in the mononuclear cell populations between the first blood sample and first SF aspirate from the same patient and between outcome (persistent versus extended-to-be) groups. Proportions of lymphocytes in the joint were modeled on chemotaxis of lymphocytes to CCL5, using Transwell migration assays. Levels of CCL5 in the SF were quantified by enzyme-linked immunosorbent assay. RNA profiles of SF mononuclear cells were compared between groups using the Affymetrix GeneChip hybridization protocol and hierarchical clustering analyses. RESULTS: Compared with peripheral blood mononuclear cells, SF mononuclear cells displayed an expansion of CD8+ T cells, reduced proportion of B cells, and expansion of CD16- natural killer cells. The lower CD4:CD8 ratio in the SF was recapitulated in vitro by the observed migration of blood T cells in response to CCL5. Synovial CCL5 levels were higher in children whose disease extended to a more severe phenotype. The CD4:CD8 ratio in the SF was significantly lower in patients with extended-to-be oligoarticular JIA (0.57 compared with 0.90 in the persistent disease group, difference 0.33, 95% confidence interval 0.04-0.62; P = 0.009). Gene expression profiling revealed that 344 genes were >1.5-fold differentially expressed between outcome groups (P < 0.05), and these included genes associated with inflammation and macrophage differentiation, which showed increased levels in patients with extended disease at 1 year, and genes associated with immune regulation, which showed increased levels in patients with persistent disease at 1 year. CONCLUSION: Analyses of the proportions of synovial lymphocytes, levels of CCL5, and differential gene expression yielded potential biomarkers with which to predict the likelihood of extension of oligoarticular JIA to a more severe disease phenotype.

Chemotherapy induces canalization of cell state in childhood B-cell precursor acute lymphoblastic leukemia.

Comparison of intratumor genetic heterogeneity in cancer at diagnosis and relapse suggests that chemotherapy induces bottleneck selection of subclonal genotypes. However, evolutionary events subsequent to chemotherapy could also explain changes in clonal dominance seen at relapse. We, therefore, investigated the mechanisms of selection in childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) during induction chemotherapy where maximal cytoreduction occurs. To distinguish stochastic versus deterministic events, individual leukemias were transplanted into multiple xenografts and chemotherapy administered. Analyses of the immediate post-treatment leukemic residuum at single-cell resolution revealed that chemotherapy has little impact on genetic heterogeneity. Rather, it acts on extensive, previously unappreciated, transcriptional and epigenetic heterogeneity in BCP-ALL, dramatically reducing the spectrum of cell states represented, leaving a genetically polyclonal but phenotypically uniform population with hallmark signatures relating to developmental stage, cell cycle and metabolism. Hence, canalization of cell state accounts for a significant component of bottleneck selection during induction chemotherapy.