PCR-Free Shallow Whole Genome Sequencing for Chromosomal Copy Number Detection from Plasma of Cancer Patients Is an Efficient Alternative to the Conventional PCR-Based Approach

Somatic copy number alterations can be detected in cell-free DNA (cfDNA) by shallow whole genome sequencing (sWGS). PCR is typically included in library preparations, but a PCR-free method could serve as a high-throughput alternative. To evaluate a PCR-free method for research and diagnostics, archival peripheral blood or bone marrow plasma samples, collected in EDTA- or lithium-heparin-containing tubes, were collected from patients with non-small-cell lung cancer (n = 10 longitudinal samples; 4 patients), B-cell lymphoma (n = 31), and acute myeloid leukemia (n = 15), or from healthy donors (n = 14). sWGS was performed on PCR-free and PCR library preparations, and the mapping quality, percentage of unique reads, genome coverage, fragment lengths, and copy number profiles were compared. The percentage of unique reads was significantly higher for PCR-free method compared with PCR method, independent of the type of collection tube: EDTA PCR-free method, 96.4% (n = 35); EDTA PCR method, 85.1% (n = 32); heparin PCR-free method, 94.5% (n = 25); and heparin PCR method, 89.4% (n = 10). All other evaluated metrics were highly comparable for PCR-free and PCR library preparations. These results demonstrate the feasibility of somatic copy number alteration detection by PCR-free sWGS using cfDNA from plasma collected in EDTA- or lithium-heparin-containing tubes and pave the way for an automated cfDNA analysis workflow for samples from cancer patients.

Engineering large-scale chromosomal deletions by CRISPR-Cas9

Large-scale chromosomal deletions are a prevalent and defining feature of cancer. A high degree of tumor-type and subtype specific recurrencies suggest a selective oncogenic advantage. However, due to their large size it has been difficult to pinpoint the oncogenic drivers that confer this advantage. Suitable functional genomics approaches to study the oncogenic driving capacity of large-scale deletions are limited. Here, we present an effective technique to engineer large-scale deletions by CRISPR-Cas9 and create isogenic cell line models. We simultaneously induce double-strand breaks (DSBs) at two ends of a chromosomal arm and select the cells that have lost the intermittent region. Using this technique, we induced large-scale deletions on chromosome 11q (65 Mb) and chromosome 6q (53 Mb) in neuroblastoma cell lines. A high frequency of successful deletions (up to 30% of selected clones) and increased colony forming capacity in the 11q deleted lines suggest an oncogenic advantage of these deletions. Such isogenic models enable further research on the role of large-scale deletions in tumor development and growth, and their possible therapeutic potential.

Expression of let-7i and miR-192 is associated with resistance to cisplatin-based chemoradiotherapy in patients with larynx and hypopharynx cancer

OBJECTIVES

The majority of patients with locally advanced larynx or hypopharynx squamous cell carcinoma are treated with organ-preserving chemoradiotherapy (CRT). Clinical outcome following CRT varies greatly. We hypothesized that tumor microRNA (miRNA) expression is predictive for outcome following CRT.

METHODS

Next-generation sequencing (NGS) miRNA profiling was performed on 37 formalin-fixed paraffin-embedded (FFPE) tumor samples. Patients with a recurrence-free survival (RFS) of less than 2 years and patients with late/no recurrence within 2 years were compared by differential expression analysis. Tumor-specific miRNAs were selected based on normal mucosa miRNA expression data from The Cancer Genome Atlas database. A model was constructed to predict outcome using group-regularized penalized logistic ridge regression. Candidate miRNAs were validated by RT-qPCR in the initial sample set as well as in 46 additional samples.

RESULTS

Thirteen miRNAs were differentially expressed (p < 0.05, FDR < 0.1) according to outcome group. Initial class prediction in the NGS cohort (n = 37) resulted in a model combining five miRNAs and disease stage, able to predict CRT outcome with an area under the curve (AUC) of 0.82. In the RT-qPCR cohort (n = 83), 25 patients (30%) experienced early recurrence (median RFS 8 months; median follow-up 42 months). Class prediction resulted in a model combining let-7i-5p, miR-192-5p and disease stage, able to discriminate patients with good versus poor clinical outcome (AUC:0.80).

CONCLUSION

The combined miRNA expression and disease stage prediction model for CRT outcome is superior to using either factor alone. This study indicates NGS miRNA profiling using FFPE specimens is feasible, resulting in clinically relevant biomarkers.

Circulating Tumor DNA as a Preoperative Marker of Recurrence in Patients with Peritoneal Metastases of Colorectal Cancer: A Clinical Feasibility Study

Cytoreductive Surgery and Hyperthermic Intraperitoneal Chemotherapy (CRS-HIPEC) may be curative for colorectal cancer patients with peritoneal metastases (PMs) but it has a high rate of morbidity. Accurate preoperative patient selection is therefore imperative, but is constrained by the limitations of current imaging techniques. In this pilot study, we explored the feasibility of circulating tumor (ct) DNA analysis to select patients for CRS-HIPEC. Thirty patients eligible for CRS-HIPEC provided blood samples preoperatively and during follow-up if the procedure was completed. Targeted Next-Generation Sequencing (NGS) of DNA from PMs was used to identify bespoke mutations that were subsequently tested in corresponding plasma cell-free (cf) DNA samples using droplet digital (dd) PCR. CtDNA was detected preoperatively in cfDNA samples from 33% of patients and was associated with a reduced disease-free survival (DFS) after CRS-HIPEC (median 6.0 months vs median not reached, p = 0.016). This association could indicate the presence of undiagnosed systemic metastases or an increased metastatic potential of the tumors. We demonstrate the feasibility of ctDNA to serve as a preoperative marker of recurrence in patients with PMs of colorectal cancer using a highly sensitive technique. A more appropriate treatment for patients with preoperative ctDNA detection may be systemic chemotherapy in addition to, or instead of, CRS-HIPEC.

Clonality analysis of pulmonary tumors by genome-wide copy number profiling

Multiple tumors in patients are frequently diagnosed, either synchronous or metachronous. The distinction between a second primary and a metastasis is important for treatment. Chromosomal DNA copy number aberrations (CNA) patterns are highly unique to specific tumors. The aim of this study was to assess genome-wide CNA-patterns as method to identify clonally related tumors in a prospective cohort of patients with synchronous or metachronous tumors, with at least one intrapulmonary tumor. In total, 139 tumor pairs from 90 patients were examined: 35 synchronous and 104 metachronous pairs. Results of CNA were compared to histological type, clinicopathological methods (Martini-Melamed-classification (MM) and ACCP-2013-criteria), and, if available, EGFR- and KRAS-mutation analysis. CNA-results were clonal in 74 pairs (53%), non-clonal in 33 pairs (24%), and inconclusive in 32 pairs (23%). Histological similarity was found in 130 pairs (94%). Concordance between histology and conclusive CNA-results was 69% (74 of 107 pairs: 72 clonal and two non-clonal). In 31 of 103 pairs with similar histology, genetics revealed non-clonality. In two out of four pairs with non-matching histology, genetics revealed clonality. The subgroups of synchronous and metachronous pairs showed similar outcome for the comparison of histological versus CNA-results. MM-classification and ACCP-2013-criteria, applicable on 34 pairs, and CNA-results were concordant in 50% and 62% respectively. Concordance between mutation matching and conclusive CNA-results was 89% (8 of 9 pairs: six clonal and two non-clonal). Interestingly, in one patient both tumors had the same KRAS mutation, but the CNA result was non-clonal. In conclusion, although some concordance between histological comparison and CNA profiling is present, arguments exist to prefer extensive molecular testing to determine whether a second tumor is a metastasis or a second primary.

Abstract 1801: MicroRNA expression profiling predicts clinical outcome in patients with locally advanced larynx and hypopharynx cancer treated with chemoradiotherapy

Background

For patients with locally advanced larynx and hypopharynx cancer, organ preservation chemoradiotherapy (CRT) protocols are frequently applied. However, some tumors are resistant to CRT or will recur within a short time frame. For these patients upfront surgical resection may result in an improved cure rate. Currently, there are no molecular markers that can predict CRT response. In this study, tumor microRNA (miRNA) profiling was performed to predict clinical outcome in patients with locally advanced larynx or hypopharynx cancer treated with CRT.

Methods

First, matched fresh frozen (FF) and formalin-fixed paraffin-embedded (FFPE) tumor biopsies from 12 patients were profiled by Next Generation Sequencing (NGS) to ensure valid miRNA results from FFPE. Next, miRNA profiles of FFPE tumor biopsies from 37 patients with larynx or hypopharynx cancer, treated with primary CRT, were determined. Clinicopathological data, including response to CRT, tumor stage and nodal stage, was collected of all patients. Differential expression analysis was performed to compare patients with poor clinical outcome, i.e residual tumor after treatment or progression free survival (PFS) &lt; 2 years versus patients with good clinical outcome, i.e. complete response and PFS &gt; 2 years. The median follow up of this cohort was 60 months (range 6 - 118). The Cancer Genome Atlas (TCGA) head and neck miRNA expression data of larynx and hypopharynx tumors and normal squamous cell tissues were used to select tumor specific miRNAs. By group-regularized penalized logistic ridge regression analysis, a model to predict clinical outcome was generated using tumor specific miRNAs combined with clinical covariates.

Results

In total, 324 novel candidate and 2,262 mature miRNAs were detected. A comparative analysis of miRNA profiles generated using matched FF and FFPE tumor biopsies showed a high correlation (mean r2=0.75). Using only FFPE samples, 13 miRNAs were significantly differentially expressed between patients with good versus poor clinical outcome (p&lt;0.05, FDR&lt;0.1). Class prediction resulted in a model of three miRNAs combined with AJCC stage, able to predict recurrence with an area under the curve of 0.83. Further validation of the model is ongoing in an additional cohort of patients with larynx or hypopharynx cancer treated with CRT.

Conclusion

NGS miRNA profiling using FFPE tumor biopsies from patients with head and neck cancer is technically feasible, resulting in high quality miRNA expression data comparable to fresh frozen tissue. Based on miRNA expression profiling, candidate miRNAs were identified predictive for clinical outcome following treatment with CRT. A validation study is ongoing.

Citation Format: Dennis Poel, Francois Rustenburg, Daoud Sie, Hendrik F. van Essen, Paul P. Eijk, Elisabeth Bloemena, Teresita E. Benites, Bauke Ylstra, Brakenhoff H. Ruud, René C. Leemans, Tineke E. Buffart, Henk M. Verheul, Jens Voortman. MicroRNA expression profiling predicts clinical outcome in patients with locally advanced larynx and hypopharynx cancer treated with chemoradiotherapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1801.

Combination of a six microRNA expression profile with four clinicopathological factors for response prediction of systemic treatment in patients with advanced colorectal cancer

Background

First line chemotherapy is effective in 75 to 80% of patients with metastatic colorectal cancer (mCRC). We studied whether microRNA (miR) expression profiles can predict treatment outcome for first line fluoropyrimidine containing systemic therapy in patients with mCRC.

Methods

MiR expression levels were determined by next generation sequencing from snap frozen tumor samples of 88 patients with mCRC. Predictive miRs were selected with penalized logistic regression and posterior forward selection. The prediction co-efficients of the miRs were re-estimated and validated by real-time quantitative PCR in an independent cohort of 81 patients with mCRC.

Results

Expression levels of miR-17-5p, miR-20a-5p, miR-30a-5p, miR-92a-3p, miR-92b-3p and miR-98-5p in combination with age, tumor differentiation, adjuvant therapy and type of systemic treatment, were predictive for clinical benefit in the training cohort with an AUC of 0.78. In the validation cohort the addition of the six miR signature to the four clinicopathological factors demonstrated a significant increased AUC for predicting treatment response versus those with stable disease (SD) from 0.79 to 0.90. The increase for predicting treatment response versus progressive disease (PD) and for patients with SD versus those with PD was not significant. in the validation cohort. MiR-17-5p, miR-20a-5p and miR-92a-3p were significantly upregulated in patients with treatment response in both the training and validation cohorts.

Conclusion

A six miR expression signature was identified that predicted treatment response to fluoropyrimidine containing first line systemic treatment in patients with mCRC when combined with four clinicopathological factors. Independent validation demonstrated added predictive value of this miR-signature for predicting treatment response versus SD. However, added predicted value for separating patients with PD could not be validated. The clinical relevance of the identified miRs for predicting treatment response has to be further explored.

Loss of Chromosome 18q11.2-q12.1 Is Predictive for Survival in Patients With Metastatic Colorectal Cancer Treated With Bevacizumab

Purpose Patients with metastatic colorectal cancer (mCRC) have limited benefit from the addition of bevacizumab to standard chemotherapy. However, a subset probably benefits substantially, highlighting an unmet clinical need for a biomarker of response to bevacizumab. Previously, we demonstrated that losses of chromosomes 5q34, 17q12, and 18q11.2-q12.1 had a significant correlation with progression-free survival (PFS) in patients with mCRC treated with bevacizumab in the CAIRO2 clinical trial but not in patients who did not receive bevacizumab in the CAIRO trial. This study was designed to validate these findings. Materials and Methods Primary mCRC samples were analyzed from two cohorts of patients who received bevacizumab as first-line treatment; 96 samples from the European multicenter study Angiopredict (APD) and 81 samples from the Italian multicenter study, MOMA. A third cohort of 90 samples from patients with mCRC who did not receive bevacizumab was analyzed. Copy number aberrations of tumor biopsy specimens were measured by shallow whole-genome sequencing and were correlated with PFS, overall survival (OS), and response. Results Loss of chromosome 18q11.2-q12.1 was associated with prolonged PFS most significantly in both the cohorts that received bevacizumab (APD: hazard ratio, 0.54; P = .01; PFS difference, 65 days; MOMA: hazard ratio, 0.55; P = .019; PFS difference, 49 days). A similar association was found for OS and overall response rate in these two cohorts, which became significant when combined with the CAIRO2 cohort. Median PFS in the cohort of patients with mCRC who did not receive bevacizumab and in the CAIRO cohort was similar to that of the APD, MOMA, and CAIRO2 patients without an 18q11.2-q12.1 loss. Conclusion We conclude that the loss of chromosome 18q11.2-q12.1 is consistently predictive for prolonged PFS in patients receiving bevacizumab. The predictive value of this loss is substantiated by a significant gain in OS and overall response rate.

Abstract 4928: MiR expression profiles can predict response to systemic treatment in patients with advanced colorectal cancer

Background and aim: Patients with advanced colorectal cancer (mCRC) are commonly treated with systemic treatment consisting of fluoropyrimidine-based regimens being ineffective in 20-25% of the patients. Currently, selection criteria for patients to predict who will respond to this treatment is lacking. The aim of this study is to identify which patients will respond to first line fluoropyrimidine-based treatment based using microRNA (miR) expression profiles in order to avoid ineffective treatment.

Material and methods: Total RNA was isolated from 88 fresh frozen colorectal cancer tissue samples consisting of ≥ 70% tumor cells, collected from patients with mCRC. MiR expression profiles were generated by next generation sequencing using the Illumina High Seq 2000 platform. Of all patients clinical and pathological data, including treatment response based on RECIST criteria, were collected. Class prediction and miR selection were performed using the GRidge package in R. Penalized selection and internal cross validation were used to select miRs predictive for treatment response.

RESULTS: Next generation sequencing resulted in a mean of 10.087.107 (range 6.114.932 to 74.313.067) reads per sample corresponding to 2567 unique mature miR sequences, including 457 novel candidate and 2110 known miRs sequences (miRbase version 19). Penalized regression analysis on tumor specific miRs identified an expression profile which was predictive for clinical benefit (defined as response and stable disease) from first line treatment.

Conclusion: With miR profiling of CRC tissue samples response prediction to first line fluoropyrimidine-based treatment in patients with mCRC is possible. We foresee that selection of treatment using miR expression profiling will avoid unnecessary treatment related toxicity and improve outcome for patients with mCRC

Citation Format: Dennis Poel, Maarten Neerincx, Daoud L.S. Sie, Nicole C.T. van Grieken, R Shankaraiah, F. S.W. van der Wolf, J. H.T.M van Waesberghe, J. D. Burggraaf, Paul P. Eijk, Bauke Ylstra, Cees Verhoef, Mark A. van de Wiel, Henk M.W. Verheul, Tineke E. Buffart. MiR expression profiles can predict response to systemic treatment in patients with advanced colorectal cancer. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 4928.

Abstract 3426: A novel approach to copy number assessment by whole genome sequencing reveals extensive spatial heterogeneity in diffuse low-grade glioma

Introduction: The disease course of patients with diffuse low-grade gliomas (LGGs) is notoriously unpredictable. Detailed analysis of the genetic make-up is therefore key to understand this wide variety in overall survival of patients diagnosed with LGGs. Reports on the prognostic value of CNAs other than 1p/19q co-deletion are conflicting, which may be explained by spatial copy number heterogeneity.

Aim: Prognostic implications of spatial and temporal copy number aberrations in LGGs by whole genome sequencing.

Methods: Approximately 25% of LGG patients have a life expectation of more than 20 years following diagnosis, which necessitates collection of samples with long clinical follow up information. Therefore archival material was selected for this project and WGS was developed to access these samples without the requirement for a normal reference. We collected 156 formalin-fixed and paraffin-embedded (FFPE) samples in a discovery cohort of 98 LGG patients with extensive clinical follow-up data, which included recurrent tumours and spatially distinct regions. An independent cohort of 126 samples was obtained from a recently published study (Alentorn et al, Neuro-oncology 2014). Technical challenges for genome-wide inference of copy number aberrations include repetitive and common DNA sequences in the genome and sequence variation across the general population, as well as the compromised and variable quality of DNA obtained from archived tissues. We developed a robust and cost effective method that infers copy number aberrations from WGS data of approximately 0.1x coverage, without the need for a reference signal. This method implements (1) a combined LOESS correction for mappability and GC content, and improves on previous methods by (2) comprehensive filtering based on public genome project data, 1000-Genomes project and ENCODE blacklists.

Results: Both prognostic value, temporal evolution and spatial heterogeneity of CNA were assessed by WGS. We confirmed prognostic favourable value of 1p/19q co-deletion, and demonstrated loss of 10q to be an unfavourable marker. In paired recurrences 10q loss was invariably maintained and surfaced in 4 additional recurrences of the discovery cohort. In spatial regions of LGGs we recognized extensive copy number heterogeneity; 15 of 17 LGGs show spatial variability of CNAs. 1p/19q co-deletion is homogeneous, while loss of 10q is heterogeneously present.

Conclusions: we present clinically relevant CNAs, but also demonstrate extensive spatial copy number heterogeneity in diffuse LGGs that might complicate unequivocal biomarker discovery.

Citation Format: Ilari Scheinin, Hinke F. van Thuijl, Daoud Sie, Hendrik F. van Essen, Paul P. Eijk, Francois Rustenburg, Ahmed Idbaih, Agusti Alentorn, Gerrit A. Meijer, Mark A. van der Wiel, Henrik Bengtsson, Adam Olshen, Eleonora Aronica, Jan J. Heimans, Jaap C. Reijneveld, Pieter Wesseling, Donna G. Albertson, Dan Pinkel, Bauke Ylstra. A novel approach to copy number assessment by whole genome sequencing reveals extensive spatial heterogeneity in diffuse low-grade glioma. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3426. doi:10.1158/1538-7445.AM2014-3426

Performance of amplicon-based next generation DNA sequencing for diagnostic gene mutation profiling in oncopathology

PURPOSE

Next generation DNA sequencing (NGS) holds promise for diagnostic applications, yet implementation in routine molecular pathology practice requires performance evaluation on DNA derived from routine formalin-fixed paraffin-embedded (FFPE) tissue specimens. The current study presents a comprehensive analysis of TruSeq Amplicon Cancer Panel-based NGS using a MiSeq Personal sequencer (TSACP-MiSeq-NGS) for somatic mutation profiling.

METHODS

TSACP-MiSeq-NGS (testing 212 hotspot mutation amplicons of 48 genes) and a data analysis pipeline were evaluated in a retrospective learning/test set approach (n = 58/n = 45 FFPE-tumor DNA samples) against 'gold standard' high-resolution-melting (HRM)-sequencing for the genes KRAS, EGFR, BRAF and PIK3CA. Next, the performance of the validated test algorithm was assessed in an independent, prospective cohort of FFPE-tumor DNA samples (n = 75).

RESULTS

In the learning set, a number of minimum parameter settings was defined to decide whether a FFPE-DNA sample is qualified for TSACP-MiSeq-NGS and for calling mutations. The resulting test algorithm revealed 82% (37/45) compliance to the quality criteria and 95% (35/37) concordant assay findings for KRAS, EGFR, BRAF and PIK3CA with HRM-sequencing (kappa = 0.92; 95% CI = 0.81-1.03) in the test set. Subsequent application of the validated test algorithm to the prospective cohort yielded a success rate of 84% (63/75), and a high concordance with HRM-sequencing (95% (60/63); kappa = 0.92; 95% CI = 0.84-1.01). TSACP-MiSeq-NGS detected 77 mutations in 29 additional genes.

CONCLUSION

TSACP-MiSeq-NGS is suitable for diagnostic gene mutation profiling in oncopathology.

Targeted sequencing by proximity ligation for comprehensive variant detection and local haplotyping

Despite developments in targeted gene sequencing and whole-genome analysis techniques, the robust detection of all genetic variation, including structural variants, in and around genes of interest and in an allele-specific manner remains a challenge. Here we present targeted locus amplification (TLA), a strategy to selectively amplify and sequence entire genes on the basis of the crosslinking of physically proximal sequences. We show that, unlike other targeted re-sequencing methods, TLA works without detailed prior locus information, as one or a few primer pairs are sufficient for sequencing tens to hundreds of kilobases of surrounding DNA. This enables robust detection of single nucleotide variants, structural variants and gene fusions in clinically relevant genes, including BRCA1 and BRCA2, and enables haplotyping. We show that TLA can also be used to uncover insertion sites and sequences of integrated transgenes and viruses. TLA therefore promises to be a useful method in genetic research and diagnostics when comprehensive or allele-specific genetic information is needed.

An exploration of pathways involved in lung carcinoid progression using gene expression profiling

Pulmonary carcinoids comprise a well-differentiated subset of neuroendocrine tumors usually associated with a favorable prognosis, but mechanisms underlying disease progression are poorly understood. In an explorative approach to identify pathways associated with progression, we compared gene expression profiles of tumors from five patients with a favorable and five with a poor disease outcome. Differentially expressed genes were validated using quantitative real-time PCR on 65 carcinoid tumors, in combination with survival analysis. One of the identified pathways was further examined using immunohistochemistry. As compared with other chromosomal locations, a significantly higher number of genes downregulated in carcinoids with a poor prognosis were located at chromosome 11q (P = 0.00017), a region known to be frequently lost in carcinoids. In addition, a number of upregulated genes were found involved in the mitotic spindle checkpoint, the chromosomal passenger complex (CPC), mitotic kinase CDC2 activity and the BRCA-Fanconi anemia pathway. At the individual gene level, BIRC5 (survivin), BUB1, CD44, IL20RA, KLK12 and OTP were independent predictors of patient outcome. For survivin, the number of positive nuclei was also related to poor prognosis within the group of carcinoids. Aurora B kinase and survivin, major components of the CPC, were particularly upregulated in high-grade carcinomas and may therefore comprise therapeutic targets for these tumors. To our knowledge, this is the first expression profiling study focusing specifically on pulmonary carcinoids and progression. We have identified novel pathways underlying malignant progression and validated several genes as being strong prognostic indicators, some of which could serve as putative therapeutic targets.

Detection limits of DNA copy number alterations in heterogeneous cell populations

BACKGROUND

Array Comparative Genomic Hybridization (aCGH) is a widely used technique to assess chromosomal copy number alterations. Chromosomal content, however, is often not uniform throughout cell populations. Here we evaluated to what extent aCGH can detect DNA copy number alterations in heterogeneous cell populations. A systematic evaluation is currently lacking, despite its importance in diagnostics and research. The detection limits reported are a compound of analytical software and laboratory techniques and do not account for the number of probes in relation to sample homogeneity.

METHODS

Detection limits were explored with DNA isolated from a patient with intellectual disability (ID) and from tumor cell line BT474. Both were diluted with increasing amounts of normal DNA to simulate different levels of cellularity. Samples were hybridized on microarrays containing 180,880 oligonucleotides evenly distributed over the genome (spacing ~17 kb).

RESULTS

Single copy number alterations, represented by down to 249 probes (4 Mb) and present in 10 % of a cell population, could be detected. Alterations encompassing as few as 14 probes (~238 Kb) could also be detected, but for this a 35 % mosaic level was required.

CONCLUSIONS

DNA copy number alterations can be detected in cell populations containing 10 % abnormal cells. Detection of sub-megabase alterations requires a higher percentage of abnormal cells or microarrays with a higher probe density.

CGH arrays compared for DNA isolated from formalin-fixed, paraffin-embedded material

Formalin-fixed, paraffin-embedded (FFPE) archival tissue is an important source of DNA material. The most commonly used technique to identify copy number aberrations from chromosomal DNA in tumorigenesis is array comparative genomic hybridization (aCGH). Although copy number analysis using DNA from FFPE archival tissue is challenging, several research groups have reported high quality and reproducible DNA copy number results using aCGH. Aim of this study is to compare the commercially available aCGH platforms suitable for high-resolution copy number analysis using FFPE-derived DNA. Two dual channel aCGH platforms (Agilent and NimbleGen) and a single channel SNP-based platform (Affymetrix) were evaluated using seven FFPE colon cancer samples, and median absolute deviation (MAD), deflection, signal-to-noise ratio, and DNA input requirements were used as quality criteria. Large differences were observed between platforms; Agilent and NimbleGen showed better MAD values (0.13 for both) compared with Affymetrix (0.22). On the contrary, Affymetrix showed a better deflection of 0.94, followed by 0.71 for Agilent and 0.51 for NimbleGen. This resulted in signal-to-nose ratios that were comparable between the three commercially available platforms. Interestingly, DNA input amounts from FFPE material lower than recommended still yielded high quality profiles on all platforms. Copy number analysis using DNA derived from FFPE archival material is feasible using all three high-resolution copy number platforms and shows reproducible results, also with DNA input amounts lower than recommended.

Integrative analysis of array-comparative genomic hybridisation and matched gene expression profiling data reveals novel genes with prognostic significance in oesophageal adenocarcinoma

BACKGROUND AND AIMS

The incidence of oesophageal adenocarcinoma (OAC) has been increasing rapidly with a dismal survival rate of less than 20%. Understanding the genomic aberrations and biology of this cancer may enhance disease interventions. This study aimed to use genome-wide genomic and expression data to enhance the understanding of OAC pathogenesis and identify groups with differential outcomes.

METHODS

Array-comparative genomic hybridisation (aCGH) analysis was carried out on 56 fresh frozen OAC resection samples with long-term clinical follow-up data. Samples with aberrations were further analysed with whole-genome single-nucleotide polymorphism arrays to confirm aCGH findings. Matched gene expression microarray data were used to identify genes with high copy number-expression correlations. Nested-multiplex PCR on DNA from microdissected specimens and fluorescence in situ hybridisation assays were used for target validation. Immunohistochemistry on the same cohort and independent samples (n=371) was used for subsequent validation. Kaplan-Meier survival analyses were performed based on aCGH data after unsupervised K-means clustering (K=5, 50 iterations) and immunohistochemistry data.

RESULTS

aCGH identified 17 common regions (>5% samples) of gains and 11 common regions of losses, including novel regions in OAC (loci 11p13 and 21q21.2). Integration of aCGH data with matched gene expression microarray data highlighted genes with high copy number-expression correlations: two deletions (p16/CDKN2A, MBNL1) and four gains (EGFR, WT1, NEIL2, MTMR9). Immunohistochemistry demonstrated protein over-expression of targets with gains: EGFR (10%), WT1 (20%), NEIL2 (14%) and MTMR9 (25%). These targets individually (p<0.060) and in combination had prognostic significance (p=0.008). On the genomic level, K-means clustering identified a cluster (32% of cohort) with differential log(2) ratios of 16 CGH probes (p<4×10(-7)) and a worse prognosis (median survival=1.37 years; p=0.015).

CONCLUSIONS

Integration of aCGH and gene expression data identified copy number aberrations and novel genes with prognostic potential in OAC.

Small bowel adenocarcinoma copy number profiles are more closely related to colorectal than to gastric cancers

BACKGROUND

Small bowel adenocarcinoma (SBA) is a rare cancer and consequently, the options for clinical trials are limited. As they are treated according to either a colorectal or a gastric cancer regimen and the molecular biology of a tumor is a pivotal determinant for therapy response, chromosomal copy number aberrations were compared with the colorectal and gastric adenocarcinomas.

MATERIALS AND METHODS

A total of 85 microsatellite stable (MSS) adenocarcinomas from the stomach, colorectum and small bowel were selected from existing array comparative genomic hybridization (aCGH) datasets. We compared the aCGH profiles of the three tumor sites by supervised analysis and hierarchical clustering.

RESULTS

Hierarchical clustering revealed substantial overlap of 27 SBA copy number profiles with matched colorectal adenocarcinomas but less overlap with profiles of gastric adenocarcinomas. DNA copy number aberrations located at chromosomes 1p36.3-p34.3, 4p15.3-q35.2, 9p24.3-p11.1, 13q13.2-q31.3 and 17p13.3-p13.2 were the strongest features discriminating SBAs and colorectal adenocarcinomas from gastric adenocarcinomas.

CONCLUSIONS

We show that MSS SBAs are more similar to colorectal than to gastric cancer, based on the 27 genome-wide DNA copy number profiles that are currently available. These molecular similarities provide added support for treatment of MSS small bowel cancers according to colorectal cancer regimens.

Trifluorothymidine resistance is associated with decreased thymidine kinase and equilibrative nucleoside transporter expression or increased secretory phospholipase A2

Trifluorothymidine (TFT) is part of the novel oral formulation TAS-102, which is currently evaluated in phase II studies. Drug resistance is an important limitation of cancer therapy. The aim of the present study was to induce resistance to TFT in H630 colon cancer cells using two different schedules and to analyze the resistance mechanism. Cells were exposed either continuously or intermittently to TFT, resulting in H630-cTFT and H630-4TFT, respectively. Cells were analyzed for cross-resistance, cell cycle, protein expression, and activity of thymidine phosphorylase (TP), thymidine kinase (TK), thymidylate synthase (TS), equilibrative nucleoside transporter (hENT), gene expression (microarray), and genomic alterations. Both cell lines were cross-resistant to 2'-deoxy-5-fluorouridine (>170-fold). Exposure to IC(75)-TFT increased the S/G(2)-M phase of H630 cells, whereas in the resistant variants, no change was observed. The two main target enzymes TS and TP remained unchanged in both TFT-resistant variants. In H630-4TFT cells, TK protein expression and activity were decreased, resulting in less activated TFT and was most likely the mechanism of TFT resistance. In H630-cTFT cells, hENT mRNA expression was decreased 2- to 3-fold, resulting in a 5- to 10-fold decreased TFT-nucleotide accumulation. Surprisingly, microarray-mRNA analysis revealed a strong increase of secretory phospholipase-A2 (sPLA2; 47-fold), which was also found by reverse transcription-PCR (RT-PCR; 211-fold). sPLA2 inhibition reversed TFT resistance partially. H630-cTFT had many chromosomal aberrations, but the exact role of sPLA2 in TFT resistance remains unclear. Altogether, resistance induction to TFT can lead to different mechanisms of resistance, including decreased TK protein expression and enzyme activity, decreased hENT expression, as well as (phospho)lipid metabolism. Mol Cancer Ther; 9(4); 1047-57. (c)2010 AACR.

Anti-proliferative action of vitamin D in MCF7 is still active after siRNA-VDR knock-down

Background

The active form of Vitamin D, 1,25-dihydroxyvitamin D₃ (1,25D), has strong anti-proliferative effects, yet the molecular mechanisms underneath this effect remain unclear. In contrast, the molecular mechanism of 1,25D for the regulation of calcium homeostasis has principally been resolved, demonstrating a pivotal role for the vitamin D receptor (VDR).

Results

We first addressed the question whether the anti-proliferative effects of 1,25D are influenced by VDR. Knockdown of VDR by siRNA did not affect the anti-proliferative effects of 1,25D in MCF7 breast cancer cells. This unanticipated finding led us to take an alternative approach using genome wide screens to study the molecular mechanisms of 1,25D in proliferation. For that purpose, four independently developed and stable 1,25D resistant MCF7 cell lines were analyzed. Array CGH identified a copy number alteration in a region of 13.5 Mb at chromosome 11q13.4-14.1 common to all four 1,25D resistant cell lines. Expression arrays revealed that no single gene was differentially expressed between the sensitive and resistant cells, but multiple membrane receptor signaling pathways were altered in the 1,25D resistant cell lines. Importantly, in the genome wide experiments neither VDR, CYP24A1 nor other known vitamin D signaling pathway genes were associated with 1,25D resistance.

Conclusion

In conclusion, siRNA and genome wide studies both suggest that the anti-proliferative effects of 1,25D in MCF7 breast tumor cell lines do not rely on classical Vitamin D pathway per se.

Genomic profiles associated with early micrometastasis in lung cancer: relevance of 4q deletion

PURPOSE

Bone marrow is a common homing organ for early disseminated tumor cells (DTC) and their presence can predict the subsequent occurrence of overt metastasis and survival in lung cancer. It is still unclear whether the shedding of DTC from the primary tumor is a random process or a selective release driven by a specific genomic pattern.

EXPERIMENTAL DESIGN

DTCs were identified in bone marrow from lung cancer patients by an immunocytochemical cytokeratin assay. Genomic aberrations and expression profiles of the respective primary tumors were assessed by microarrays and fluorescence in situ hybridization analyses. The most significant results were validated on an independent set of primary lung tumors and brain metastases.

RESULTS

Combination of DNA copy number profiles (array comparative genomic hybridization) with gene expression profiles identified five chromosomal regions differentiating bone marrow-negative from bone marrow-positive patients (4q12-q32, 10p12-p11, 10q21-q22, 17q21, and 20q11-q13). Copy number changes of 4q12-q32 were the most prominent finding, containing the highest number of differentially expressed genes irrespective of chromosomal size (P=0.018). Fluorescence in situ hybridization analyses on further primary lung tumor samples confirmed the association between loss of 4q and bone marrow-positive status. In bone marrow-positive patients, 4q was frequently lost (37% versus 7%), whereas gains could be commonly found among bone marrow-negative patients (7% versus 17%). The same loss was also found to be common in brain metastases from both small and non-small cell lung cancer patients (39%).

CONCLUSIONS

Thus, our data indicate, for the first time, that early hematogenous dissemination of tumor cells might be driven by a specific pattern of genomic changes.

DNA copy number profiles of primary tumors as predictors of response to chemotherapy in advanced colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is biologically a heterogeneous disease, which may affect response to drug therapy. We investigated the correlation of genome-wide DNA copy number profiles of primary tumors with response to systemic chemotherapy in advanced CRC.

PATIENTS AND METHODS

DNA was isolated from formaldehyde-fixed paraffin-embedded primary tumors of 32 patients with advanced CRC, which were selected based on either a good response (n = 16) or a poor response (n = 16) to first-line combination therapy with capecitabine and irinotecan. High-resolution DNA copy number profiles were obtained by means of 30 K oligonucleotide-based array comparative genomic hybridization (aCGH).

RESULTS

Unsupervised hierarchical cluster analysis of the aCGH data revealed two clusters of 19 and 13 tumors, respectively, and cluster membership showed a significant correlation with response status (P < 0.03). The nonresponders had fewer chromosomal alterations compared with the responders, in particular less losses were found (P < 0.03). Most prominent differences between the two groups were losses of regions 18p11.32-q11.2 (P < 0.02) and 18q12.1-q23 (P < 0.03), which were more frequently observed in responders.

CONCLUSIONS

Differences in DNA copy number profiles of primary CRCs are associated with response to systemic combination chemotherapy with capecitabine and irinotecan. Responders overall had more chromosomal alterations, especially loss of chromosome 18.

Micro-array analysis of resistance for gemcitabine results in increased expression of ribonucleotide reductase subunits

To study in detail the relation between gene expression and resistance against gemcitabine, a cell line was isolated from a tumor for which gemcitabine resistance was induced in vivo. Similar to the in vivo tumor, resistance in this cell line, C 26-G, was not related to deficiency of deoxycytidine kinase (dCK). Micro-array analysis showed increased expression of ribonucleotide reductase (RR) subunits M1 and M2 as confirmed by real time PCR analysis (28- and 2.7-fold, respectively). In cell culture, moderate cross-resistance (about 2-fold) was observed to 1-ss-D-arabinofuranosylcytosine (ara-C), 2-chloro-2'deoxyadenosine (CdA), LY231514 (ALIMTA), and cisplatin (CDDP), and pronounced cross-resistance (>23-fold) to 2',2'-difluorodeoxyuridine (dFdU) and 2',2'-difluorodeoxyguanosine (dFdG). Culture in the absence of gemcitabine reduced resistance as well as RRM1 RNA expression, demonstrating a direct relationship of RRM1 RNA expression with acquired resistance to gemcitabine.

Profiling of apoptosis genes allows for clinical stratification of primary nodal diffuse large B-cell lymphomas

Intrinsic resistance of lymphoma cells to apoptosis is a probable mechanism causing chemotherapy resistance and eventual fatal outcome in patients with diffuse large B cell lymphomas (DLBCL). We investigated whether microarray expression profiling of apoptosis related genes predicts clinical outcome in 46 patients with primary nodal DLBCL. Unsupervised cluster analysis using genes involved in apoptosis (n = 246) resulted in three separate DLBCL groups partly overlapping with germinal centre B-lymphocytes versus activated B-cells like phenotype. One group with poor clinical outcome was characterised by high expression levels of pro-and anti-apoptotic genes involved in the intrinsic apoptosis pathway. A second group, also with poor clinical outcome, was characterised by high levels of apoptosis inducing cytotoxic effector genes, possibly reflecting a cellular cytotoxic immune response. The third group showing a favourable outcome was characterised by low expression levels of genes characteristic for both other groups. Our results suggest that chemotherapy refractory DLBCL are characterised either by an intense cellular cytotoxic immune response or by constitutive activation of the intrinsic mediated apoptosis pathway with concomitant downstream inhibition of this apoptosis pathway. Consequently, strategies neutralising the function of apoptosis-inhibiting proteins might be effective as alternative treatment modality in part of chemotherapy refractory DLBCL.

Cross-platform array comparative genomic hybridization meta-analysis separates hematopoietic and mesenchymal from epithelial tumors

A series of studies have been published that evaluate the chromosomal copy number changes of different tumor classes using array comparative genomic hybridization (array CGH); however, the chromosomal aberrations that distinguish the different tumor classes have not been fully characterized. Therefore, we performed a meta-analysis of different array CGH data sets in an attempt to classify samples tested across different platforms. As opposed to RNA expression, a common reference is used in dual channel CGH arrays: normal human DNA, theoretically facilitating cross-platform analysis. To this aim, cell line and primary cancer data sets from three different dual channel array CGH platforms obtained by four different institutes were integrated. The cell line data were used to develop preprocessing methods, which performed noise reduction and transformed samples into a common format. The transformed array CGH profiles allowed perfect clustering by cell line, but importantly not by platform or institute. The same preprocessing procedures used for the cell line data were applied to data from 373 primary tumors profiled by array CGH, including controls. Results indicated that there is no apparent feature related to the institute or platform and that array CGH allows for unambiguous cross-platform meta-analysis. Major clusters with common tissue origin were identified. Interestingly, tumors of hematopoietic and mesenchymal origins cluster separately from tumors of epithelial origin. Therefore, it can be concluded that chromosomal aberrations of tumors from hematopoietic and mesenchymal origin versus tumors of epithelial origin are distinct, and these differences can be picked up by meta-analysis of array CGH data. This suggests the possibility of prospectively using combined analysis of diverse copy number data sets for cancer subtype classification.

Gastric cancers in young and elderly patients show different genomic profiles

Although most gastric cancers occur in elderly patients, a substantial number of cases of this common disease occur in young patients. Gastric cancer is a heterogeneous disease at the genomic level and different patterns of DNA copy number alterations are associated with different clinical behaviour. The aim of the present study was to explore differences in DNA copy number alterations in relation to age of onset of gastric cancer. DNA isolated from 46 paraffin-embedded gastric cancer tissue samples from 17 patients less than 50 years of age [median 43 (21-49) years] and 29 patients greater than or equal to 70 years of age [median 75 (70-83) years] was analysed by genome-wide microarray comparative genomic hybridization (array CGH) using an array of 5000 BAC clones. Patterns of DNA copy number aberrations were analysed by hierarchical cluster analysis of the mode-normalized and smoothed log(2) ratios of tumour to normal reference fluorescence signal intensities using TMEV software, after which cluster membership was correlated with age group. In addition, supervised analysis was performed using CGH Multi-array. Hierarchical cluster analysis of the array CGH data revealed three clusters with different genomic profiles that correlated significantly with age (p = 0.006). Cluster 1 mainly contained young patients, while elderly patients were divided over clusters 2 and 3. Chromosome regions 11q23.3 and 19p13.3 contributed most to age-related differences in tumour profiles. Gastric cancers of young and old patients belong to groups with different genomic profiles, which likely reflect different pathogenic mechanisms of the disease.

In vivo induction of resistance to gemcitabine results in increased expression of ribonucleotide reductase subunit M1 as the major determinant

Gemcitabine is a deoxycytidine (dCyd) analogue with activity against several solid cancers. Gemcitabine is activated by dCyd kinase (dCK) and interferes, as its triphosphate dFdCTP, with tumor growth through incorporation into DNA. Alternatively, the metabolite gemcitabine diphosphate (dFdCDP) can interfere with DNA synthesis and thus tumor growth through inhibition of ribonucleotide reductase. Gemcitabine can be inactivated by the enzyme dCyd deaminase (dCDA). In most in vitro models, resistance to gemcitabine was associated with a decreased dCK activity. In all these models, resistance was established using continuous exposure to gemcitabine with increasing concentrations; however, these in vitro models have limited clinical relevance. To develop in vivo resistance to gemcitabine, we treated mice bearing a moderately sensitive tumor Colon 26-A (T/C = 0.25) with a clinically relevant schedule (120 mg/kg every 3 days). By repeated transplant of the most resistant tumor and continuation of gemcitabine treatment for >1 year, the completely resistant tumor Colon 26-G (T/C = 0.96) was created. Initial studies focused on resistance mechanisms known from in vitro studies. In Colon 26-G, dCK activity was 1.7-fold decreased; dCDA and DNA polymerase were not changed; and Colon 26-G accumulated 1.5-fold less dFdCTP, 6 hours after a gemcitabine injection, than the parental tumor. Based on in vitro studies, these relative minor changes were considered insufficient to explain the completely resistant phenotype. Therefore, an expression microarray was done with Colon 26-A versus Colon 26-G. Using independently grown nonresistant and resistant tumors, a striking increase in expression of the RRM1 subunit gene was found in Colon 26-G. The expression of RRM1 mRNA was 25-fold increased in the resistant tumor, as measured by real-time PCR, which was confirmed by Western blotting. In contrast, RRM2 mRNA was 2-fold decreased. However, ribonucleotide reductase enzyme activity was only moderately increased in Colon 26-G. In conclusion, this is the first model with in vivo induced resistance to gemcitabine. In contrast to most in vitro studies, dCK activity was not the most important determinant of gemcitabine resistance. Expression microarray identified RRM1 as the gene with the highest increase in expression in the Colon 26-G, which might clarify its complete gemcitabine-resistant phenotype. This study is the first in vivo evidence for a key role for RRM1 in acquired gemcitabine resistance.

Microarray analysis reveals pivotal divergent mRNA expression profiles early in the development of either compensated ventricular hypertrophy or heart failure

Myocardial right ventricular (RV) hypertrophy due to pulmonary hypertension is aimed at normalizing ventricular wall stress. Depending on the degree of pressure overload, RV hypertrophy may progress to a state of impaired contractile function and heart failure, but this cannot be discerned during the early stages of ventricular remodeling. We tested whether critical differences in gene expression profiles exist between ventricles before the ultimate development of either a compensated or decompensated hypertrophic phenotype. Both phenotypes were selectively induced in Wistar rats by a single subcutaneous injection of either a low or a high dose of the pyrrolizidine alkaloid monocrotaline (MCT). Spotted oligonucleotide microarrays were used to investigate pressure-dependent cardiac gene expression profiles at 2 wk after the MCT injections, between control rats and rats that would ultimately develop either compensated or decompensated hypertrophy. Clustering of significantly regulated genes revealed specific expression profiles for each group, although the degree of hypertrophy was still similar in both. The ventricles destined to progress to failure showed activation of pro-apoptotic pathways, particularly related to mitochondria, whereas the group developing compensated hypertrophy showed blocked pro-death effector signaling via p38-MAPK, through upregulation of MAPK phosphatase-1. In summary, we show that, already at an early time point, pivotal differences in gene expression exist between ventricles that will ultimately develop either a compensated or a decompensated phenotype, depending on the degree of pressure overload. These data reveal genes that may provide markers for the early prediction of clinical outcome as well as potential targets for early intervention.

Assessment of Chlamydia trachomatis infection of semen specimens by ligase chain reaction

Diagnostic potential of the Chlamydia trachomatis ligase chain reaction system (LCx) to assess the presence of C. trachomatis in urine and semen specimens was evaluated. Paired urine and semen specimens from 153 asymptomatic male partners of subfertile couples attending our Center for Reproductive Medicine were examined by LCx. As controls, 19 semen samples from four donors who were participating in the programme for artificial insemination were used. Of these, 12 samples had previously been shown to be C. trachomatis-positive by an in-house PCR. C. trachomatis was detected by LCx in seven of 153 (5 %) urine samples. None of the 153 semen samples tested positive by LCx. Also, none of the 12 C. trachomatis-containing control semen samples were positive by LCx. By in-house PCR, seven urine specimens and two of 153 (1 %) semen samples tested positive. The corresponding urine samples of these male partners were also C. trachomatis-positive, as well as the 12 C. trachomatis-containing samples from donors. In conclusion, LCx is not sensitive enough to assess the presence of C. trachomatis in semen specimens; therefore, this method is not recommended to routinely screen semen specimens from donors who participate in programmes for artificial insemination or male partners of subfertile couples for C. trachomatis.

Normal IncA expression and fusogenicity of inclusions in Chlamydia trachomatis isolates with the incA I47T mutation

To investigate the correlation between the incA I47T mutation in Chlamydia trachomatis and the nonfusogenic phenotype, the incA genes of 25 isolates were sequenced. Four major sequence types were identified. Seven isolates (28%) had the I47T mutation. Isolates representing the four sequence types expressed IncA in the membrane of one large single inclusion. In conclusion, the incA I47T mutation is not associated with the nonfusogenic phenotype.

PCR assessment of Chlamydia trachomatis infection of semen specimens processed for artificial insemination

In order to ascertain the microbiological quality of stored semen specimens processed for artificial insemination by a donor (AID), we developed a PCR assay targeting the chlamydial plasmid to detect Chlamydia trachomatis in semen. The lower limit of detection of this assay corresponded to 2.5 to 5 elementary bodies per microl of semen. A total of 669 cryopreserved ejaculates from 97 asymptomatic donors were tested for C. trachomatis infection. Twelve ejaculates, originating from four donors, were found to be positive, indicating a 4% prevalence of C. trachomatis infection among the donor population studied. Cross-contamination between the cryopreserved specimens in the storage container was studied by typing using sequence analysis of PCR-amplified omp1 genes of the strains. Two donors were infected with serovar E, one was infected with serovar F, and one was infected with serovar K. For two donors, the duration of C. trachomatis positivity could be assessed. One donor donated C. trachomatis-positive semen for at least 4 successive months, and the other did so for at least 16 months. The occurrence of C. trachomatis infection in cryopreserved donor semen indicates that ejaculates from donors not tested for a C. trachomatis infection just prior to donation should be tested for infection by a direct test such as the PCR described here. Direct testing of semen specimens will detect not only donors with an active infection but also C. trachomatis-infected ejaculates already stored and will thus improve the microbiological quality of AID, since discrepancies in the presence of C. trachomatis in urine and semen specimens have been reported.

Stimulation of the adherence of Haemophilus influenzae to human lung epithelial cells by antimicrobial neutrophil defensins

Patients with chronic obstructive pulmonary disease (COPD) frequently have recurrent lower respiratory tract infections with nonencapsulated Haemophilus influenzae. The infected mucosa of these patients is infiltrated with neutrophils, which upon activation may release antimicrobial peptides, including defensins. It was shown that defensins isolated from neutrophils or from sputum samples of COPD patients did not kill H. influenzae from these patients, but they did stimulate its adherence to human bronchial epithelial cells in a time- and dose-dependent manner. Maximal stimulation was observed after 3 h in the presence of > or = 10 micrograms/mL defensins, resulting in 65 +/- 36 cfu/cell (61-fold increase). The enhanced adherence was not solely due to charge effects and was specifically blocked by alpha 1-proteinase inhibitor. Because adherence is the first step in the onset of respiratory tract infections, our findings indicate that neutrophil defensins likely contribute to the pathogenesis of H. influenzae infection in the lower respiratory tract.

The fimbria gene cluster of nonencapsulated Haemophilus influenzae

The occurrence of fimbria gene clusters in nonencapsulated Haemophilus influenzae strains from chronic bronchitis patients (n = 58), patients with acute otitis media (n = 13), and healthy carriers (n = 12) was determined by DNA hybridization and PCR, based on sequences of fimbriate H. influenzae type b. Although an average of 18% of all nonencapsulated strains had a fimbria gene cluster consisting of hifA to hifE inserted in the chromosome between purE and pepN, differences in the frequency of fimbria cluster-positive strains were observed, depending on the source of isolates. The compositions of the fimbria gene clusters of seven strains from chronic bronchitis patients and one strain from an otitis media patient were analyzed in more detail. After enrichment for fimbria expression, the promoter of the gene cluster contained 10 TA repeats (n = 2), leading to optimal positioning between the -10 and -35 promoter regions. The promoter regions of five fimbria-negative strains were sequenced; four were found to have nine TA repeats, and one had only four TA repeats. The protein sequence of three ganglioside GM1-specific HifA adhesins consisted of conserved regions intermingled with regions of sequence diversity. hifA appeared to be flanked by intergenic regions that varied between strains and contained both direct and inverted DNA repeats. Since noncoding DNA between hifA and purE has not been found in H. influenzae type b, these DNA sequences are probably not essential for fimbria expression. An analysis of strains lacking the gene cluster revealed the presence of similar sequences in 13 of 15 strains from chronic bronchitis patients, 5 of 5 strains from otitis media patients, and 3 of 5 strains from healthy carriers. The lengths of these intergenic regions were the same for multiple isolates of strains obtained during persistent infections. The presence or absence and the composition of the fimbria gene cluster and other sequences between the flanking genes purE and pepN suggest that the fimbria gene cluster was originally contained on a mobile element.

Molecular variation in the major outer membrane protein P5 gene of nonencapsulated Haemophilus influenzae during chronic infections

During the course of persistent infections by nonencapsulated Haemophilus influenzae in patients with chronic bronchitis, the major outer membrane protein (MOMP) P5 varies in molecular weight. The nature of this variability was determined by DNA sequence analysis of the P5 gene from five different H. influenzae strains and their seven MOMP P5 variants which were isolated from patients with chronic infections of the lower respiratory tract. Analysis of the P5 sequence data from the different strains revealed four well-defined, heterogeneous regions. These regions of variable sequence appeared to correspond to the regions of the gene encoding the putative surface-exposed loops of MOMP P5. The MOMP P5 variants with alterations in MOMP P5 were shown to result from DNA point mutations and codon deletions. In addition, in three variants derived sequentially from one H. influenzae strain, a frameshift mutation resulted in the formation of a stop codon in the region encoding the signal sequence of the MOMP P5 gene. Strikingly, all nucleotide substitutions in the MOMP P5 loop regions of variants were nonsynonymous, suggesting that variants with alterated amino acid compositions of the surface-exposed parts of MOMP P5 obtained a selective advantage during persistence of the infection by nonencapsulated H. influenzae in chronic bronchitis patients.

Haemophilus influenzae infections in patients with chronic obstructive pulmonary disease despite specific antibodies in serum and sputum

The titer and specificity of antibodies to the infecting Haemophilus influenzae was determined in sera and sputa from 27 patients with chronic obstructive pulmonary disease (COPD) to analyze the specific immune response. COPD patients had significantly higher serum IgG and IgA antibody titers than 13 healthy controls (mean IgG titers 12,302 and 5,623, respectively; mean IgA titers, 2,398 and 912; p less than 0.001). The mean IgM titers were comparable: 501 and 447, respectively. Specific IgA antibodies were also detectable in the sputum of the COPD patients (mean IgA antibody titer, 776). The local antibody production was determined by calculating the relative coefficient of excretion (RCE) to albumin. The mean RCE of 89.1 for IgA indicated statistically significant local production (p less than 0.02), in contrast to a nonsignificant increase for IgG (mean RCE of 3.6). Specific IgM was below the detection level. Immunoblotting experiments showed that the antibodies in sera from COPD patients and controls were directed against most of the outer membrane proteins of H. influenzae, with individual differences between IgG, IgA, and IgM. The IgA and IgG antibodies in serum had a similar specificity as those in sputum. The appearance or persistence of H. influenzae coincided with minor changes in antibody titer and specificity. From these results we conclude that COPD patients are infected with H. influenzae despite the presence of at least as many antibodies in sputum and serum as in controls and that these antibodies are directed against a variety of antigenic determinants of the infecting strain.

Endogenous and exogenous reinfections by Haemophilus influenzae in patients with chronic obstructive pulmonary disease: the effect of antibiotic treatment on persistence

To analyze whether exacerbations in chronic obstructive pulmonary disease (COPD) coincide with reinfection by Haemophilus influenzae, 16 COPD patients were studied longitudinally for 3 years. Exacerbations coincided with reinfection by H. influenzae, either endogenous, by a strain with a DNA fingerprint indistinguishable from the strain previously present but with another major outer membrane protein (MOMP) pattern (2 patients), or exogenous, by a strain with a different DNA fingerprint and MOMP pattern (3 patients). The other patients, remaining in an infectious state without clear exacerbations for longer periods, were persistently infected by a particular H. influenzae strain (median persistence time, 5.5 months; range, 2-23 months). Of 8 antibiotic-treated patients, 7 remained infected by H. influenzae with the same DNA fingerprint, although all strains were sensitive to the antibiotics prescribed. Results of the study suggested that exacerbations in COPD patients coincide with endogenous or exogenous reinfection by H. influenzae, persistently infected patients keep the same H. influenzae strain for longer periods, and antibiotic treatment was not effective in eradicating H. influenzae.

Antigenic drift of Haemophilus influenzae in patients with chronic obstructive pulmonary disease

Differences in the major outer membrane protein b,c (molecular weight, 39,000 to 41,000) of related Haemophilus influenzae strains isolated from the sputum of patients with chronic obstructive pulmonary disease were analyzed biochemically and immunologically. Protein b,c was isolated from a total of six related H. influenzae strains from two chronic obstructive pulmonary disease patients. After CNBr digestion of the proteins, the differences in size appeared in the larger of the two fragments. Trypsin and chymotrypsin digests of proteins from related H. influenzae strains showed that proteins differed by only a few peptides or were very similar, in contrast to the peptide maps of proteins from nonrelated strains. Peptide analysis of b,c proteins from related H. influenzae strains by high-performance liquid chromatography after Staphylococcus aureus V8 protease digestion and amino acid analysis of corresponding fractions revealed highly comparable patterns, indicating only minor differences in the amino acid sequences of these proteins. Immunization of rabbits with intact bacteria of four related H. influenzae strains resulted in a strong anti-protein b,c antibody response in all animals. In three of four rabbits, antibodies specific for the b,c protein of the strain used for immunization were elicited, indicating that the changed proteins contained specific immunodominant epitopes. Anti-protein b,c antibodies promoted strain-specific, complement-dependent, bactericidal activity. From these results, we conclude that H. influenzae shows antigenic drift in immunodominant epitopes, caused by small changes in amino acid composition of the b,c protein. Antibodies to these epitopes promote complement-dependent bactericidal activity.

Nonculture detection of Haemophilus influenzae in sputum with monoclonal antibodies specific for outer membrane lipoprotein P6

Isolation of Haemophilus influenzae from sputum is hampered by overgrowth by other microorganisms or by antibiotic treatment of the patient. To overcome this problem in the detection of H. influenzae, an in situ immunoperoxidase staining technique was developed with monoclonal antibody (MAb) 8BD9, immunoglobulin subclass G2a. MAb 8BD9 appeared to be directed to an epitope on the outer membrane lipoprotein P6 of H. influenzae. The species specificity of MAb 8BD9 was analyzed by staining isolates from different bacterial species. MAb 8BD9 reacted with all 300 H. influenzae strains tested and with H. aegyptius and H. haemolyticus. Twenty-six of 30 H. parainfluenzae strains, other Haemophilus species, and other bacterial species often isolated from sputum were not stained. The staining technique was compared with culture of 845 routinely obtained sputum samples from patients with respiratory tract infections. In 829 sputa (98.1%), the results of both techniques were in agreement; 173 were positive for H. influenzae, and 656 were negative. With 14 sputum samples, the staining method gave a positive result for H. influenzae, but the bacterium was not cultured. This could be ascribed to antibiotic treatment of the patient (n = 7), the presence of other MAb 8BD9-positive Haemophilus species in the sputum (n = 5), and overgrowth by swarming Proteus mirabilis or by Branhamella catarrhalis. In the immunoperoxidase- and Gram-stained smears of two sputum samples, no bacteria were seen, although some H. influenzae was cultured. On the basis of these results, we conclude that immunoperoxidase staining with MAb 8BD9 is a fast and reliable extension of the available detection techniques for H. influenzae.

Changes in outer membrane proteins of nontypable Haemophilus influenzae in patients with chronic obstructive pulmonary disease

Five individual colonies of Haemophilus influenzae were isolated from each of one to three cultures of sputum collected from 18 patients with chronic obstructive pulmonary disease (COPD). The isolates were studied to investigate whether the major outer membrane proteins (MOMPs) changed during persistence. The relationship between isolates was analyzed by fingerprinting their chromosomal DNA. The fingerprints of eight strains (isolated from eight patients) with various MOMP compositions were different, whereas fingerprints of isolates with identical MOMP compositions were indistinguishable. In 12 patients, two or more strains with different MOMP compositions were found; one strain was isolated from the sputum samples of each of the six remaining patients. In seven of the 12 patients, strains with different MOMPs but with indistinguishable fingerprints were found. The differences were found in proteins b,c (five patients) and d (five patients). In patients with COPD, the MOMPs of H. influenzae are subject to changes that may enable this bacterium to escape immunological defense mechanisms.

Comparison of two nonculture techniques for detection of Hemophilus influenzae in sputum. In situ hybridization and immunoperoxidase staining with monoclonal antibodies

Two nonculture methods, in situ hybridization and immunoperoxidase staining with monoclonal antibodies, were compared for the detection of Hemophilus influenzae in 184 sputa. For in situ hybridization, a biotin-labeled probe of total genomic DNA of H influenzae type b was prepared that hybridizes specifically with H influenzae, H parainfluenzae, H hemolyticus, and H parahemolyticus DNA. Immunoperoxidase staining was done with monoclonal antibody 8BD9 directed against outer membrane protein P6 of H influenzae. Both techniques detected Hemophilus in sputum equally well and were superior to culture: all 30 sputum samples culture-positive for H influenzae were positive on both nonculture tests, and 13 additional positive sputum samples were detected from which Hemophilus was not cultured. The higher sensitivity of the nonculture tests was mainly attributed to culture failure because of overgrowth of H influenzae by other bacteria, especially in patients with cystic fibrosis. The immunoperoxidase staining technique appeared slightly easier and quicker to perform than the in situ hybridization test. For the in situ DNA hybridization probe, DNA can be prepared from any strain of H influenzae. The immunoperoxidase test requires monoclonal antibody 8BD9 but has a higher specificity than the hybridization technique. Both techniques can be reliably applied, especially for the detection of Hemophilus in sputum of patients with cystic fibrosis.