Standardization of Workflow and Flow Cytometry Panels for Quantitative Expression Profiling of Surface Antigens on Blood Leukocyte Subsets: An HCDM CDMaps Initiative.

Background: The Human Cell Differentiation Molecules (HCDM) organizes Human Leukocyte Differentiation Antigen (HLDA) workshops to test and name clusters of antibodies that react with a specific antigen. These cluster of differentiation (CD) markers have provided the scientific community with validated antibody clones, consistent naming of targets and reproducible identification of leukocyte subsets. Still, quantitative CD marker expression profiles and benchmarking of reagents at the single-cell level are currently lacking. Objective: To develop a flow cytometric procedure for quantitative expression profiling of surface antigens on blood leukocyte subsets that is standardized across multiple research laboratories. Methods: A high content framework to evaluate the titration and reactivity of Phycoerythrin (PE)-conjugated monoclonal antibodies (mAbs) was created. Two flow cytometry panels were designed: an innate cell tube for granulocytes, dendritic cells, monocytes, NK cells and innate lymphoid cells (12-color) and an adaptive lymphocyte tube for naive and memory B and T cells, including TCRγδ+, regulatory-T and follicular helper T cells (11-color). The potential of these 2 panels was demonstrated via expression profiling of selected CD markers detected by PE-conjugated antibodies and evaluated using 561 nm excitation. Results: Using automated data annotation and dried backbone reagents, we reached a robust workflow amenable to processing hundreds of measurements in each experiment in a 96-well plate format. The immunophenotyping panels enabled discrimination of 27 leukocyte subsets and quantitative detection of the expression of PE-conjugated CD markers of interest that could quantify protein expression above 400 units of antibody binding capacity. Expression profiling of 4 selected CD markers (CD11b, CD31, CD38, CD40) showed high reproducibility across centers, as well as the capacity to benchmark unique clones directed toward the same CD3 antigen. Conclusion: We optimized a procedure for quantitative expression profiling of surface antigens on blood leukocyte subsets. The workflow, bioinformatics pipeline and optimized flow panels enable the following: 1) mapping the expression patterns of HLDA-approved mAb clones to CD markers; 2) benchmarking new antibody clones to established CD markers; 3) defining new clusters of differentiation in future HLDA workshops.

Disease progression of 213 patients hospitalized with Covid-19 in the Czech Republic in March-October 2020: An exploratory analysis.

We collected a multi-centric retrospective dataset of patients (N = 213) who were admitted to ten hospitals in Czech Republic and tested positive for SARS-CoV-2 during the early phases of the pandemic in March-October 2020. The dataset contains baseline patient characteristics, breathing support required, pharmacological treatment received and multiple markers on daily resolution. Patients in the dataset were treated with hydroxychloroquine (N = 108), azithromycin (N = 72), favipiravir (N = 9), convalescent plasma (N = 7), dexamethasone (N = 4) and remdesivir (N = 3), often in combination. To explore association between treatments and patient outcomes we performed multiverse analysis, observing how the conclusions change between defensible choices of statistical model, predictors included in the model and other analytical degrees of freedom. Weak evidence to constrain the potential efficacy of azithromycin and favipiravir can be extracted from the data. Additionally, we performed external validation of several proposed prognostic models for Covid-19 severity showing that they mostly perform unsatisfactorily on our dataset.

A novel class of ZNF384 aberrations in acute leukemia.

Fusion of the ZNF384 gene as the 3' partner to several different 5' partner genes occurs recurrently in B-cell precursor acute lymphoblastic and mixed phenotype B/myeloid leukemia. These canonical fusions (ZNF384r) contain the complete ZNF384 coding sequence and are associated with a specific gene expression signature. Cases with this signature, but without canonical ZNF384 fusions (ZNF384r-like cases), have been described previously. Although some have been shown to harbor ZNF362 fusions, the primary aberrations remain unknown in a major proportion. We studied 3 patients with the ZNF384r signature and unknown primary genetic background and identified a previously unknown class of genetic aberration affecting the last exon of ZNF384 and resulting in disruption of the C-terminal portion of the ZNF384 protein. Importantly, in 2 cases, the ZNF384 aberration, indel, was missed during the bioinformatic analysis but revealed by the manual, targeted reanalysis. Two cases with the novel aberrations had a mixed (B/myeloid) immunophenotype commonly associated with canonical ZNF384 fusions. In conclusion, we present leukemia cases with a novel class of ZNF384 aberrations that phenocopy leukemia with ZNF384r. Therefore, we show that part of the so-called ZNF384r-like cases represent the same genetic subtype as leukemia with canonical ZNF384 fusions.

Nucleoporin TPR Affects C2C12 Myogenic Differentiation via Regulation of Myh4 Expression.

The nuclear pore complex (NPC) has emerged as a hub for the transcriptional regulation of a subset of genes, and this type of regulation plays an important role during differentiation. Nucleoporin TPR forms the nuclear basket of the NPC and is crucial for the enrichment of open chromatin around NPCs. TPR has been implicated in the regulation of transcription; however, the role of TPR in gene expression and cell differentiation has not been described. Here we show that depletion of TPR results in an aberrant morphology of murine proliferating C2C12 myoblasts (MBs) and differentiated C2C12 myotubes (MTs). The ChIP-Seq data revealed that TPR binds to genes linked to muscle formation and function, such as myosin heavy chain (Myh4), myocyte enhancer factor 2C (Mef2C) and a majority of olfactory receptor (Olfr) genes. We further show that TPR, possibly via lysine-specific demethylase 1 (LSD1), promotes the expression of Myh4 and Olfr376, but not Mef2C. This provides a novel insight into the mechanism of myogenesis; however, more evidence is needed to fully elucidate the mechanism by which TPR affects specific myogenic genes.

DUX4r, ZNF384r and PAX5-P80R mutated B-cell precursor acute lymphoblastic leukemia frequently undergo monocytic switch.

Recently, we described B-cell precursor acute lymphoblastic leukemia (BCP-ALL) subtype with early switch to the monocytic lineage and loss of the B-cell immunophenotype, including CD19 expression. Thus far, the genetic background has remained unknown. Among 726 children consecutively diagnosed with BCP-ALL, 8% patients experienced switch detectable by flow cytometry (FC). Using exome and RNA sequencing, switch was found to positively correlate with three different genetic subtypes: PAX5-P80R mutation (5 cases with switch out of 5), rearranged DUX4 (DUX4r; 30 cases of 41) and rearranged ZNF384 (ZNF384r; 4 cases of 10). Expression profiles or phenotypic patterns correlated with genotypes, but within each genotype they could not identify cases who subsequently switched. If switching was not taken into account, the B-cell-oriented FC assessment underestimated the minimal residual disease level. For patients with PAX5-P80R, a discordance between FC-determined and PCR-determined MRD was found on day 15, resulting from a rapid loss of the B-cell phenotype. Discordance on day 33 was observed in all the DUX4r, PAX5-P80R and ZNF384r subtypes. Importantly, despite the substantial phenotypic changes, possibly even challenging the appropriateness of BCP-ALL therapy, the monocytic switch was not associated with a higher incidence of relapse and poorer prognosis in patients undergoing standard ALL treatment.

M2-like macrophages dictate clinically relevant immunosuppression in metastatic ovarian cancer.

BACKGROUND: The immunological microenvironment of primary high-grade serous carcinomas (HGSCs) has a major impact on disease outcome. Conversely, little is known on the microenvironment of metastatic HGSCs and its potential influence on patient survival. Here, we explore the clinical relevance of the immunological configuration of HGSC metastases. METHODS: RNA sequencing was employed on 24 paired primary tumor microenvironment (P-TME) and metastatic tumor microenvironment (M-TME) chemotherapy-naive HGSC samples. Immunohistochemistry was used to evaluate infiltration by CD8+ T cells, CD20+ B cells, DC-LAMP+ (lysosomal-associated membrane protein 3) dendritic cells (DCs), NKp46+ (natural killer) cells and CD68+CD163+ M2-like tumor-associated macrophages (TAMs), abundance of PD-1+ (programmed cell death 1), LAG-3+ (lymphocyte-activating gene 3) cells, and PD-L1 (programmed death ligand 1) expression in 80 samples. Flow cytometry was used for functional assessments on freshly resected HGSC samples. RESULTS: 1468 genes were differentially expressed in the P-TME versus M-TME of HGSCs, the latter displaying signatures of extracellular matrix remodeling and immune infiltration. M-TME infiltration by immune effector cells had little impact on patient survival. Accordingly, M-TME-infiltrating T cells were functionally impaired, but not upon checkpoint activation. Conversely, cytokine signaling in favor of M2-like TAMs activity appeared to underlie inhibited immunity in the M-TME and poor disease outcome. CONCLUSIONS: Immunosuppressive M2-like TAM infiltrating metastatic sites limit clinically relevant immune responses against HGSCs.

ShinySOM: graphical SOM-based analysis of single-cell cytometry data.

SUMMARY: ShinySOM offers a user-friendly interface for reproducible, high-throughput analysis of high-dimensional flow and mass cytometry data guided by self-organizing maps. The software implements a FlowSOM-style workflow, with improvements in performance, visualizations and data dissection possibilities. The outputs of the analysis include precise statistical information about the dissected samples, and R-compatible metadata useful for the batch processing of large sample volumes. AVAILABILITY AND IMPLEMENTATION: ShinySOM is free and open-source, available online at gitlab.com/exaexa/ShinySOM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

High-Content Immunophenotyping and Hierarchical Clustering Reveal Sources of Heterogeneity and New Surface Markers of Human Blood Monocyte Subsets.

OBJECTIVE: Blood monocyte subsets are emerging as biomarkers of cardiovascular inflammation. However, our understanding of human monocyte heterogeneity and their immunophenotypic features under healthy and inflammatory conditions is still evolving. RATIONALE: In this study, we sought to investigate the immunophenome of circulating human monocyte subsets. METHODS: Multiplexed, high-throughput flow cytometry screening arrays and computational data analysis were used to analyze the expression and hierarchical relationships of 242 specific surface markers on circulating classical (CD14++CD16-), intermediate (CD14++CD16+), and nonclassical (CD14+CD16++) monocytes in healthy adults. RESULTS: Using generalized linear models and hierarchical cluster analysis, we selected and clustered epitopes that most reliably differentiate between monocyte subsets. We validated existing transcriptional profiling data and revealed potential new surface markers that uniquely define the classical (e.g., BLTR1, CD35, CD38, CD49e, CD89, CD96), intermediate (e.g., CD39, CD275, CD305, CDw328), and nonclassical (e.g., CD29, CD132) subsets. In addition, our analysis revealed phenotypic cell clusters, identified by dendritic markers CMRF-44 and CMRF-56, independent of the traditional monocyte classification. CONCLUSION: These results reveal an advancement of the clinically applicable multiplexed screening arrays that may facilitate monocyte subset characterization and cytometry-based biomarker selection in various inflammatory disorders.

CD Maps-Dynamic Profiling of CD1-CD100 Surface Expression on Human Leukocyte and Lymphocyte Subsets.

CD molecules are surface molecules expressed on cells of the immune system that play key roles in immune cell-cell communication and sensing the microenvironment. These molecules are essential markers for the identification and isolation of leukocytes and lymphocyte subsets. Here, we present the results of the first phase of the CD Maps study, mapping the expression of CD1-CD100 (n = 110) on 47 immune cell subsets from blood, thymus, and tonsil using an eight-color standardized EuroFlow approach and quantification of expression. The resulting dataset included median antibody binding capacities (ABCs) and percentage of positivity for all markers on all subsets and was developed into an interactive CD Maps web resource. Using the resource, we examined differentially expressed proteins between granulocyte, monocyte, and dendritic cell subsets, and profiled dynamic expression of markers during thymocyte differentiation, T-cell maturation, and between functionally distinct B-cell subset clusters. The CD Maps resource will serve as a benchmark of antibody reactivities ensuring improved reproducibility of flow cytometry-based research. Moreover, it will provide a full picture of the surfaceome of human immune cells and serves as a useful platform to increase our understanding of leukocyte biology, as well as to facilitate the identification of new biomarkers and therapeutic targets of immunological and hematological diseases.

Toll-like receptor 2 expression on c-kit+ cells tracks the emergence of embryonic definitive hematopoietic progenitors.

Hematopoiesis in mammalian embryos proceeds through three successive waves of hematopoietic progenitors. Since their emergence spatially and temporally overlap and phenotypic markers are often shared, the specifics regarding their origin, development, lineage restriction and mutual relationships have not been fully determined. The identification of wave-specific markers would aid to resolve these uncertainties. Here, we show that toll-like receptors (TLRs) are expressed during early mouse embryogenesis. We provide phenotypic and functional evidence that the expression of TLR2 on E7.5 c-kit+ cells marks the emergence of precursors of erythro-myeloid progenitors (EMPs) and provides resolution for separate tracking of EMPs from primitive progenitors. Using in vivo fate mapping, we show that at E8.5 the Tlr2 locus is already active in emerging EMPs and in progenitors of adult hematopoietic stem cells (HSC). Together, this data demonstrates that the activation of the Tlr2 locus tracks the earliest events in the process of EMP and HSC specification.

Molecular Basis of Cisplatin Resistance in Testicular Germ Cell Tumors.

The emergence of cisplatin (CDDP) resistance is the main cause of treatment failure and death in patients with testicular germ cell tumors (TGCT), but its biologic background is poorly understood. To study the molecular basis of CDDP resistance in TGCT we prepared and sequenced CDDP-exposed TGCT cell lines as well as 31 primary patients' samples. Long-term exposure to CDDP increased the CDDP resistance 10 times in the NCCIT cell line, while no major resistance was achieved in Tera-2. Development of CDDP resistance was accompanied by changes in the cell cycle (increase in G1 and decrease in S-fraction), increased number of acquired mutations, of which 3 were present within ATRX gene, as well as changes in gene expression pattern. Copy number variation analysis showed, apart from obligatory gain of 12p, several other large-scale gains (chr 1, 17, 20, 21) and losses (chr X), with additional more CNVs found in CDDP-resistant cells (e.g., further losses on chr 1, 4, 18, and gain on chr 8). In the patients' samples, those who developed CDDP resistance and died of TGCT (2/31) showed high numbers of acquired aberrations, both SNPs and CNVs, and harbored mutations in genes potentially relevant to TGCT development (e.g., TRERF1, TFAP2C in one patient, MAP2K1 and NSD1 in another one). Among all primary tumor samples, the most commonly mutated gene was NSD1, affected in 9/31 patients. This gene encoding histone methyl transferase was also downregulated and identified among the 50 most differentially expressed genes in CDDP-resistant NCCIT cell line. Interestingly, 2/31 TGCT patients harbored mutations in the ATRX gene encoding a chromatin modifier that has been shown to have a critical function in sexual differentiation. Our research newly highlights its probable involvement also in testicular tumors. Both findings support the emerging role of altered epigenetic gene regulation in TGCT and CDDP resistance development.

Nuclear pore protein TPR associates with lamin B1 and affects nuclear lamina organization and nuclear pore distribution.

The organization of the nuclear periphery is crucial for many nuclear functions. Nuclear lamins form dense network at the nuclear periphery and play a substantial role in chromatin organization, transcription regulation and in organization of nuclear pore complexes (NPCs). Here, we show that TPR, the protein located preferentially within the nuclear baskets of NPCs, associates with lamin B1. The depletion of TPR affects the organization of lamin B1 but not lamin A/C within the nuclear lamina as shown by stimulated emission depletion microscopy. Finally, reduction of TPR affects the distribution of NPCs within the nuclear envelope and the effect can be reversed by simultaneous knock-down of lamin A/C or the overexpression of lamin B1. Our work suggests a novel role for the TPR at the nuclear periphery: the TPR contributes to the organization of the nuclear lamina and in cooperation with lamins guards the interphase assembly of nuclear pore complexes.

ERG deletions in childhood acute lymphoblastic leukemia with DUX4 rearrangements are mostly polyclonal, prognostically relevant and their detection rate strongly depends on screening method sensitivity.

ERG-deletions occur recurrently in acute lymphoblastic leukemia, especially in the DUX4-rearranged subtype. The ERG-deletion was shown to positively impact prognosis of patients with IKZF1-deletion and its presence precludes assignment into IKZF1 plus group, a novel high-risk category on AIEOP-BFM ALL trials. We analyzed the impact of different methods on ERG-deletion detection rate, evaluated ERG-deletion as a potential marker for DUX4-rearranged leukemia, studied its associations with molecular and clinical characteristics within this leukemia subtype, and analyzed its clonality. Using single-nucleotide-polymorphism array, genomic polymerase chain reaction (PCR) and amplicon-sequencing we found ERG-deletion in 34% (16 of 47), 66% (33 of 50) and 78% (39 of 50) of DUX4-rearranged leukemia, respectively. False negativity of ERG-deletion by single-nucleotide-polymorphism array caused IKZF1 plus misclassification in 5 patients. No ERG-deletion was found outside the DUX4-rearranged cases. Within DUX4-rearranged leukemia, the ERG-deletion was associated with higher total number of copy-number aberrations, and, importantly, the ERG-deletion positivity by PCR was associated with better outcome [5-year event-free survival (EFS), ERG-deletion-positive 93% vs. ERG-deletion-negative 68%, P=0.022; 5-year overall survival (OS), ERG-deletion-positive 97% vs. ERG-deletion-negative 75%, P=0.029]. Ultra-deep amplicon-sequencing revealed distinct co-existing ERG-deletions in 22 of 24 patients. In conclusion, our data demonstrate inadequate sensitivity of single-nucleotide-polymorphism array for ERG-deletion detection, unacceptable for proper IKZF1 plus classification. Even using more sensitive methods (PCR/amplicon-sequencing) for its detection, ERG-deletion is absent in 22-34% of DUX4-rearranged leukemia and does not represent an adequately sensitive marker of this leukemia subtype. Importantly, the ERG-deletion potentially stratifies the DUX4-rearranged leukemia into biologically/clinically distinct subsets. Frequent polyclonal pattern of ERG-deletions shows that late origin of this lesion is more common than has been previously described.

Genomic landscape of pediatric B-other acute lymphoblastic leukemia in a consecutive European cohort.

Novel biological subtypes and clinically important genetic aberrations (druggable lesions, prognostic factors) have been described in B-other acute lymphoblastic leukemia (ALL) during the last decade; however, due to a lack of studies on unselected cohorts, their population frequency and mutual associations still have to be established. We studied 110 consecutively diagnosed and uniformly treated childhood B-other patients using single nucleotide polymorphism arrays and whole exome/transcriptome sequencing. The frequency of DUX4-rearranged, BCR-ABL1-like, ZNF384-rearranged, ETV6-RUNX1-like, iAMP21 and MEF2D-rearranged subtypes was 27%, 15%, 5%, 5%, 4%, and 2%, respectively; 43% of cases were not classified into any of these subtypes (B-rest). We found worse early response to treatment in DUX4-rearranged leukemia and a strong association of ZNF384-rearranged leukemia with B-myeloid immunophenotype. Of the druggable lesions, JAK/STAT-class and RAS/RAF/MAPK-class aberrations were found in 21% and 43% of patients, respectively; an ABL-class aberration was found in one patient. A recently described negative prognostic factor, IKZF1plus , was found in 14% of patients and was enriched in (but not exclusive for) BCR-ABL1-like subtype. PAX5 fusions (including 4 novel), intragenic amplifications and P80R mutations were mutually exclusive and only occurred in the B-rest subset, altogether accounting for 20% of the B-other group. PAX5 P80R was associated with a specific gene expression signature, potentially defining a novel leukemia subtype. Our study shows unbiased European population-based frequencies of novel ALL subtypes, recurrent (cyto)genetic aberrations and their mutual associations. This study also strengthens and widens the current knowledge of B-other ALL and provides an objective basis for optimization of current genetic diagnostics.

Inhalation of ZnO Nanoparticles: Splice Junction Expression and Alternative Splicing in Mice.

Despite the wide application of nanomaterials, toxicity studies of nanoparticles (NP) are often limited to in vitro cell models, and the biological impact of NP exposure in mammals has not been thoroughly investigated. Zinc oxide (ZnO) NPs are commonly used in various consumer products. To evaluate the effects of the inhalation of ZnO NP in mice, we studied splice junction expression in the lungs as a proxy to gene expression changes analysis. Female ICR mice were treated with 6.46 × 104 and 1.93 × 106 NP/cm3 for 3 days and 3 months, respectively. An analysis of differential expression and alternative splicing events in 298 targets (splice junctions) of 68 genes involved in the processes relevant to the biological effects of ZnO NP was conducted using next-generation sequencing. Three days of exposure resulted in the upregulation of IL-6 and downregulation of BID, GSR, NF-kB2, PTGS2, SLC11A2, and TXNRD1 splice junction expression; 3 months of exposure increased the expression of splice junctions in ALDH3A1, APAF1, BID, CASP3, DHCR7, GCLC, GCLM, GSR, GSS, EHHADH, FAS, HMOX-1, IFNγ, NF-kB1, NQO-1, PTGS1, PTGS2, RAD51, RIPK2, SRXN1, TRAF6, and TXNRD1. Alternative splicing of TRAF6 and TXNRD1 was induced after 3 days of exposure to 1.93 × 106 NP/cm3. In summary, we observed changes of splice junction expression in genes involved in oxidative stress, apoptosis, immune response, inflammation, and DNA repair, as well as the induction of alternative splicing in genes associated with oxidative stress and inflammation. Our data indicate the potential negative biological effects of ZnO NP inhalation.

flowIO: Flow cytometry standard conformance testing, editing, and export tool.

The Flow Cytometry Standard (FCS) format is a widely accepted norm for storing Flow Cytometry (FCM) data. Its goal as a standard is to allow FCM data sharing and re-analysis. Over more than three decades of its existence FCS has evolved into a well-defined, flexible file format reflecting technical changes in the FCM field. Its flexibility as well as rising numbers of instrument vendors leads to suboptimal implementations of FCS in some cases. Such situations compromise the primary goal of the standard and hinder the ability to reproduce FCM analyses. It is further underlined by rapid rise of advanced FCM analyses, often carried out outside traditional software tools and heavily relying on standard data storage and presentation. We have developed flowIO, an R package which tests FCS file conformance with the standard as defined by International Society for Advancement of Cytometry (ISAC) normative. Along with the package we provide a web based application (also at http://bioinformin.cesnet.cz/flowIO/) allowing user friendly access to the conformance testing as well as FCS file editing and export for further analysis.

Low HOX gene expression in PML-RARα-positive leukemia results from suppressed histone demethylation.

Homeobox (HOX) genes are frequently dysregulated in leukemia. Previous studies have shown that aberrant HOX gene expression accompanies leukemogenesis and affects disease progression and leukemia patient survival. Patients with acute myeloid leukemia (AML) bearing PML-RARα fusion gene have distinct HOX gene signature in comparison to other subtypes of AML patients, although the mechanism of transcription regulation is not completely understood. We previously found an association between the mRNA levels of HOX genes and those of the histone demethylases JMJD3 and UTX in PML-RARα- positive leukemia patients. Here, we demonstrate that the release of the PML-RARα-mediated block in PML-RARα-positive myeloid leukemia cells increased both JMJD3 and HOX gene expression, while inhibition of JMJD3 using the specific inhibitor GSK-J4 reversed the effect. This effect was driven specifically through PML-RARα fusion protein since expression changes did not occur in cells with mutated RARα and was independent of differentiation. We confirmed that gene expression levels were inversely correlated with alterations in H3K27me3 histone marks localized at HOX gene promoters. Furthermore, data from chromatin immunoprecipitation followed by sequencing broaden a list of clustered HOX genes regulated by JMJD3 in PML-RARα-positive leukemic cells. Interestingly, the combination of GSK-J4 and all-trans retinoic acid (ATRA) significantly increased PML-RARα-positive cell apoptosis compared with ATRA treatment alone. This effect was also observed in ATRA-resistant NB4 clones, which may provide a new therapeutic opportunity for patients with acute promyelocytic leukemia (APL) resistant to current treatment. The results of our study reveal the mechanism of HOX gene expression regulation and contribute to our understanding of APL pathogenesis.

What matters in chronic Burkholderia cenocepacia infection in cystic fibrosis: Insights from comparative genomics.

Burkholderia cenocepacia causes severe pulmonary infections in cystic fibrosis (CF) patients. Since the bacterium is virtually untreatable by antibiotics, chronic infections persist for years and might develop into fatal septic pneumonia (cepacia syndrome, CS). To devise new strategies to combat chronic B. cenocepacia infections, it is essential to obtain comprehensive knowledge about their pathogenesis. We conducted a comparative genomic analysis of 32 Czech isolates of epidemic clone B. cenocepacia ST32 isolated from various stages of chronic infection in 8 CF patients. High numbers of large-scale deletions were found to occur during chronic infection, affecting preferentially genomic islands and nonessential replicons. Recombination between insertion sequences (IS) was inferred as the mechanism behind deletion formation; the most numerous IS group was specific for the ST32 clone and has undergone transposition burst since its divergence. Genes functionally related to transition metal metabolism were identified as hotspots for deletions and IS insertions. This functional category was also represented among genes where nonsynonymous point mutations and indels occurred parallelly among patients. Another category exhibiting parallel mutations was oxidative stress protection; mutations in catalase KatG resulted in impaired detoxification of hydrogen peroxide. Deep sequencing revealed substantial polymorphism in genes of both categories within the sputum B. cenocepacia ST32 populations, indicating extensive adaptive evolution. Neither oxidative stress response nor transition metal metabolism genes were previously reported to undergo parallel evolution during chronic CF infection. Mutations in katG and copper metabolism genes were overrepresented in patients where chronic infection developed into CS. Among professional phagocytes, macrophages use both hydrogen peroxide and copper for their bactericidal activity; our results thus tentatively point to macrophages as suspects in pathogenesis towards the fatal CS.