Integration of genomic, transcriptomic and proteomic data identifies two biologically distinct subtypes of invasive lobular breast cancer.

Invasive lobular carcinoma (ILC) is the second most frequently occurring histological breast cancer subtype after invasive ductal carcinoma (IDC), accounting for around 10% of all breast cancers. The molecular processes that drive the development of ILC are still largely unknown. We have performed a comprehensive genomic, transcriptomic and proteomic analysis of a large ILC patient cohort and present here an integrated molecular portrait of ILC. Mutations in CDH1 and in the PI3K pathway are the most frequent molecular alterations in ILC. We identified two main subtypes of ILCs: (i) an immune related subtype with mRNA up-regulation of PD-L1, PD-1 and CTLA-4 and greater sensitivity to DNA-damaging agents in representative cell line models; (ii) a hormone related subtype, associated with Epithelial to Mesenchymal Transition (EMT), and gain of chromosomes 1q and 8q and loss of chromosome 11q. Using the somatic mutation rate and eIF4B protein level, we identified three groups with different clinical outcomes, including a group with extremely good prognosis. We provide a comprehensive overview of the molecular alterations driving ILC and have explored links with therapy response. This molecular characterization may help to tailor treatment of ILC through the application of specific targeted, chemo- and/or immune-therapies.

Integrated genome-wide genotyping and gene expression profiling reveals BCL11B as a putative oncogene in acute myeloid leukemia with 14q32 aberrations.

Acute myeloid leukemia is a neoplasm characterized by recurrent molecular aberrations traditionally demonstrated by cytogenetic analyses. We used high density genome-wide genotyping and gene expression profiling to reveal acquired cryptic abnormalities in acute myeloid leukemia. By genome-wide genotyping of 137 cases of primary acute myeloid leukemia, we disclosed a recurrent focal amplification on chromosome 14q32, which included the genes BCL11B, CCNK, C14orf177 and SETD3, in two cases. In the affected cases, the BCL11B gene showed consistently high mRNA expression, whereas the expression of the other genes was unperturbed. Fluorescence in situ hybridization on 40 cases of acute myeloid leukemia with high BCL11B mRNA expression [2.5-fold above median; 40 out of 530 cases (7.5%)] revealed 14q32 abnormalities in two additional cases. In the four BCL11B-rearranged cases the 14q32 locus was fused to different partner chromosomes. In fact, in two cases, we demonstrated that the focal 14q32 amplifications were integrated into transcriptionally active loci. The translocations involving BCL11B result in increased expression of full-length BCL11B protein. The BCL11B-rearranged acute myeloid leukemias expressed both myeloid and T-cell markers. These biphenotypic acute leukemias all carried FLT3 internal tandem duplications, a characteristic marker of acute myeloid leukemia. BCL11B mRNA expression in acute myeloid leukemia appeared to be strongly associated with expression of other T-cell-specific genes. Myeloid 32D(GCSF-R) cells ectopically expressing Bcl11b showed decreased proliferation rate and less maturation. In conclusion, by an integrated approach involving high-throughput genome-wide genotyping and gene expression profiling we identified BCL11B as a candidate oncogene in acute myeloid leukemia.

Loss of E-cadherin is not a necessity for epithelial to mesenchymal transition in human breast cancer.

Epithelial to mesenchymal transition (EMT) is typically defined by the acquisition of a spindle cell morphology in combination with loss of E-cadherin and upregulation of mesenchymal markers. However, by studying E-cadherin inactivation in 38 human breast cancer cell lines, we noted that not all cell lines that had undergone EMT had concomitantly lost E-cadherin expression. We further investigated this discrepancy functionally and in clinical breast cancer specimens. Interestingly, reconstitution of wild-type E-cadherin cDNA in a E-cadherin negative cell line that had undergone EMT (MDA-MB-231) did not revert the spindle morphology back to an epithelial morphology. Neither were changes observed in the expression of several markers known to be involved in the EMT process. Similarly, upregulation of E-cadherin via global DNA demethylation in eleven cell lines that had undergone EMT did not induce a change in cell morphology, nor did it alter the expression of EMT markers in these cells. Next, we extracted genes differentially expressed between cell lines that had undergone EMT versus cell lines that had not undergone EMT. Caveolin-1 was identified to be an excellent marker for EMT, irrespective of E-cadherin status (specificity and sensitivity of 100 %). Consistent with our observations in the breast cancer cell lines, expression of Caveolin-1 identified a subset of basal breast cancers, particularly of metaplastic pathology, and only 50 % of these lacked E-cadherin expression. The discrepancy between E-cadherin loss and EMT was thus reproduced in clinical samples. Together, these results indicate that in human breast cancer loss of E-cadherin is not causal for EMT and even not a necessity.

The human GFI136N variant induces epigenetic changes at the Hoxa9 locus and accelerates K-RAS driven myeloproliferative disorder in mice.

The coding single nucleotide polymorphism GFI136N in the human gene growth factor independence 1 (GFI1) is present in 3%-7% of whites and increases the risk for acute myeloid leukemia (AML) by 60%. We show here that GFI136N, in contrast to GFI136S, lacks the ability to bind to the Gfi1 target gene that encodes the leukemia-associated transcription factor Hoxa9 and fails to initiate histone modifications that regulate HoxA9 expression. Consistent with this, AML patients heterozygous for the GFI136N variant show increased HOXA9 expression compared with normal controls. Using ChipSeq, we demonstrate that GFI136N specific epigenetic changes are also present in other genes involved in the development of AML. Moreover, granulomonocytic progenitors, a bone marrow subset from which AML can arise in humans and mice, show a proliferative expansion in the presence of the GFI136N variant. In addition, granulomonocytic progenitors carrying the GFI136N variant allele have altered gene expression patterns and differ in their ability to grow after transplantation. Finally, GFI136N can accelerate a K-RAS driven fatal myeloproliferative disease in mice. Our data suggest that the presence of a GFI136N variant allele induces a preleukemic state in myeloid precursors by deregulating the expression of Hoxa9 and other AML-related genes.

Gene expression profiling assigns CHEK2 1100delC breast cancers to the luminal intrinsic subtypes.

CHEK2 1100delC is a moderate-risk cancer susceptibility allele that confers a high breast cancer risk in a polygenic setting. Gene expression profiling of CHEK2 1100delC breast cancers may reveal clues to the nature of the polygenic CHEK2 model and its genes involved. Here, we report global gene expression profiles of a cohort of 155 familial breast cancers, including 26 CHEK2 1100delC mutant tumors. In line with previous work, all CHEK2 1100delC mutant tumors clustered among the hormone receptor-positive breast cancers. In the hormone receptor-positive subset, a 40-gene CHEK2 signature was subsequently defined that significantly associated with CHEK2 1100delC breast cancers. The identification of a CHEK2 gene signature implies an unexpected biological homogeneity among the CHEK2 1100delC breast cancers. In addition, all 26 CHEK2 1100delC tumors classified as luminal intrinsic subtype breast cancers, with 8 luminal A and 18 luminal B tumors. This biological make-up of CHEK2 1100delC breast cancers suggests that a relatively limited number of additional susceptibility alleles are involved in the polygenic CHEK2 model. Identification of these as-yet-unknown susceptibility alleles should be aided by clues from the 40-gene CHEK2 signature.

Retroviral integration mutagenesis in mice and comparative analysis in human AML identify reduced PTP4A3 expression as a prognostic indicator.

Acute myeloid leukemia (AML) results from multiple genetic and epigenetic aberrations, many of which remain unidentified. Frequent loss of large chromosomal regions marks haplo-insufficiency as one of the major mechanisms contributing to leukemogenesis. However, which haplo-insufficient genes (HIGs) are involved in leukemogenesis is largely unknown and powerful experimental strategies aimed at their identification are currently lacking. Here, we present a new approach to discover HIGs, using retroviral integration mutagenesis in mice in which methylated viral integration sites and neighbouring genes were identified. In total we mapped 6 genes which are flanked by methylated viral integration sites (mVIS). Three of these, i.e., Lrmp, Hcls1 and Prkrir, were up regulated and one, i.e., Ptp4a3, was down regulated in the affected tumor. Next, we investigated the role of PTP4A3 in human AML and we show that PTP4A3 expression is a negative prognostic indicator, independent of other prognostic parameters. In conclusion, our novel strategy has identified PTP4A3 to potentially have a role in AML, on one hand as a candidate HIG contributing to leukemogenesis in mice and on the other hand as a prognostic indicator in human AML.

Characterization of CEBPA mutations and promoter hypermethylation in pediatric acute myeloid leukemia.

BACKGROUND: Dysfunctioning of CCAAT/enhancer binding protein α (C/EBPα) in acute myeloid leukemia can be caused, amongst others, by mutations in the encoding gene (CEBPA) and by promoter hypermethylation. CEBPA-mutated acute myeloid leukemia is associated with a favorable outcome, but this may be restricted to the case of double mutations in CEBPA in adult acute myeloid leukemia. In pediatric acute myeloid leukemia, data on the impact of these mutations are limited to one series, and data on promoter hypermethylation are lacking. Our objective was to investigate the characteristics, gene expression profiles and prognostic impact of the different CEBPA aberrations in pediatric acute myeloid leukemia. DESIGN AND METHODS: We screened a large pediatric cohort (n=252) for CEBPA single and double mutations by direct sequencing, and for promoter hypermethylation by methylation-specific polymerase chain reaction. Furthermore, we determined the gene-expression profiles (Affymetrix HGU133 plus 2.0 arrays) of this cohort (n=237). RESULTS: Thirty-four mutations were identified in 20 out of the 252 cases (7.9%), including 14 double-mutant and 6 single-mutant cases. CEBPA double mutations conferred a significantly better 5-year overall survival compared with single mutations (79% versus 25%, respectively; P=0.04), and compared with CEBPA wild-type acute myeloid leukemia excluding core-binding factor cases (47%; P=0.07). Multivariate analysis confirmed that the double mutations were an independent favorable prognostic factor for survival (hazard ratio 0.23, P=0.04). The combination of screening for promoter hypermethylation and gene expression profiling identified five patients with silenced CEBPA, of whom four cases relapsed. All cases characteristically expressed T-lymphoid markers. Moreover, unsupervised clustering of gene expression profiles showed a clustering of CEBPA double-mutant and silenced cases, pointing towards a common hallmark of abrogated C/EBPα-functioning in these acute myeloid leukemias. CONCLUSIONS: We showed the independent favorable outcome of patients with CEBPA double-mutant acute myeloid leukemia in a large pediatric series. This molecular marker may, therefore, improve risk-group stratification in pediatric acute myeloid leukemia. For the first time, CEBPA-silenced cases are suggested to confer a poor outcome in pediatric acute myeloid leukemia, indicating that further investigation of this aberration is needed. Furthermore, clustering of gene expression profiles provided insight into the biological similarities and diversities of the different aberrations in CEBPA in pediatric acute myeloid leukemia.

Leukemic IDH1 and IDH2 mutations result in a hypermethylation phenotype, disrupt TET2 function, and impair hematopoietic differentiation.

Cancer-associated IDH mutations are characterized by neomorphic enzyme activity and resultant 2-hydroxyglutarate (2HG) production. Mutational and epigenetic profiling of a large acute myeloid leukemia (AML) patient cohort revealed that IDH1/2-mutant AMLs display global DNA hypermethylation and a specific hypermethylation signature. Furthermore, expression of 2HG-producing IDH alleles in cells induced global DNA hypermethylation. In the AML cohort, IDH1/2 mutations were mutually exclusive with mutations in the α-ketoglutarate-dependent enzyme TET2, and TET2 loss-of-function mutations were associated with similar epigenetic defects as IDH1/2 mutants. Consistent with these genetic and epigenetic data, expression of IDH mutants impaired TET2 catalytic function in cells. Finally, either expression of mutant IDH1/2 or Tet2 depletion impaired hematopoietic differentiation and increased stem/progenitor cell marker expression, suggesting a shared proleukemogenic effect.

HAT: hypergeometric analysis of tiling-arrays with application to promoter-GeneChip data.

BACKGROUND: Tiling-arrays are applicable to multiple types of biological research questions. Due to its advantages (high sensitivity, resolution, unbiased), the technology is often employed in genome-wide investigations. A major challenge in the analysis of tiling-array data is to define regions-of-interest, i.e., contiguous probes with increased signal intensity (as a result of hybridization of labeled DNA) in a region. Currently, no standard criteria are available to define these regions-of-interest as there is no single probe intensity cut-off level, different regions-of-interest can contain various numbers of probes, and can vary in genomic width. Furthermore, the chromosomal distance between neighboring probes can vary across the genome among different arrays. RESULTS: We have developed Hypergeometric Analysis of Tiling-arrays (HAT), and first evaluated its performance for tiling-array datasets from a Chromatin Immunoprecipitation study on chip (ChIP-on-chip) for the identification of genome-wide DNA binding profiles of transcription factor Cebpa (used for method comparison). Using this assay, we can refine the detection of regions-of-interest by illustrating that regions detected by HAT are more highly enriched for expected motifs in comparison with an alternative detection method (MAT). Subsequently, data from a retroviral insertional mutagenesis screen were used to examine the performance of HAT among different applications of tiling-array datasets. In both studies, detected regions-of-interest have been validated with (q)PCR. CONCLUSIONS: We demonstrate that HAT has increased specificity for analysis of tiling-array data in comparison with the alternative method, and that it accurately detects regions-of-interest in two different applications of tiling-arrays. HAT has several advantages over previous methods: i) as there is no single cut-off level for probe-intensity, HAT can detect regions-of-interest at various thresholds, ii) it can detect regions-of-interest of any size, iii) it is independent of probe-resolution across the genome, and across tiling-array platforms and iv) it employs a single user defined parameter: the significance level. Regions-of-interest are detected by computing the hypergeometric-probability, while controlling the Family Wise Error. Furthermore, the method does not require experimental replicates, common regions-of-interest are indicated, a sequence-of-interest can be examined for every detected region-of-interest, and flanking genes can be reported.

A variant allele of Growth Factor Independence 1 (GFI1) is associated with acute myeloid leukemia.

The GFI1 gene encodes a transcriptional repressor, which regulates myeloid differentiation. In the mouse, Gfi1 deficiency causes neutropenia and an accumulation of granulomonocytic precursor cells that is reminiscent of a myelodysplastic syndrome. We report here that a variant allele of GFI1 (GFI1(36N)) is associated with acute myeloid leukemia (AML) in white subjects with an odds ratio of 1.6 (P < 8 x 10(-5)). The GFI1(36N) variant occurred in 1806 AML patients with an allele frequency of 0.055 compared with 0.035 in 1691 healthy control patients in 2 independent cohorts. We observed that both GFI1 variants maintain the same activity as transcriptional repressors but differ in their regulation by the AML1/ETO (RUNX1/RUNX1T1) fusion protein produced in AML patients with a t(8;21) translocation. AML1/ETO interacts and colocalizes with the more common GFI1(36S) form in the nucleus and inhibits its repressor activity. However, the variant GFI1(36N) protein has a different subnuclear localization than GFI1(36S). As a consequence, AML1/ETO does not colocalize with GFI1(36N) and is unable to inhibit its repressor activity. We conclude that both variants of GFI1 differ in their ability to be regulated by interacting proteins and that the GFI1(36N) variant form exhibits distinct biochemical features that may confer a predisposition to AML.

Distinct gene mutation profiles among luminal-type and basal-type breast cancer cell lines.

Breast cancer has for long been recognized as a highly diverse tumor group, but the underlying genetic basis has been elusive. Here, we report an extensive molecular characterization of a collection of 41 human breast cancer cell lines. Protein and gene expression analyses indicated that the collection of breast cancer cell lines has retained most, if not all, molecular characteristics that are typical for clinical breast cancers. Gene mutation analyses identified 146 oncogenic mutations among 27 well-known cancer genes, amounting to an average of 3.6 mutations per cell line. Mutations in genes from the p53, RB and PI3K tumor suppressor pathways were widespread among all breast cancer cell lines. Most important, we have identified two gene mutation profiles that are specifically associated with luminal-type and basal-type breast cancer cell lines. The luminal mutation profile involved E-cadherin and MAP2K4 gene mutations and amplifications of Cyclin D1, ERBB2 and HDM2, whereas the basal mutation profile involved BRCA1, RB1, RAS and BRAF gene mutations and deletions of p16 and p14ARF. These subtype-specific gene mutation profiles constitute a genetic basis for the heterogeneity observed among human breast cancers, providing clues for their underlying biology and providing guidance for targeted pharmacogenetic intervention in breast cancer patients.

Micro-RNA-15a and micro-RNA-16 expression and chromosome 13 deletions in multiple myeloma.

We have used copy number variation (CNV) analysis with SNP mapping arrays for miRNA-15a and miRNA-16-1 expression analysis in patients with multiple myeloma (MM) with or without deletion of chromosome 13q14. MiRNA-15a and miRNA-16 display a range of expression patterns in MM patients, independent of the chromosome 13 status. These findings suggest that genes other than miR-15a and miR-16 may explain the prognostic significance of 13q14 deletions.

Unusually high prevalence of panton-valentine leukocidin genes among methicillin-sensitive Staphylococcus aureus strains carried in the Indonesian population.

Few data on the molecular characteristics and epidemiology of Staphylococcus aureus from Indonesia are available. The purpose of the present study was to define S. aureus reservoirs in both the Indonesian community and hospital using a collection of 329 nasal carriage isolates obtained during a survey of 3,995 healthy individuals and patients from Java, Indonesia. Only one strain (0.3%) was identified as methicillin-resistant S. aureus by mecA gene PCR. The Panton-Valentine leukocidin (PVL) genes were detected in 35 methicillin-sensitive S. aureus strains (10.6%). Molecular typing by pulsed-field gel electrophoresis of the 329 isolates showed extensive genetic diversity among both PVL-positive and PVL-negative strains. In Surabaya, Indonesia, however, a cluster was identified that was strongly associated with the presence of the PVL locus (P < 0.0001). As determined by high-throughput amplified fragment length polymorphism, PVL-positive strains occurred throughout all major AFLP clusters (I to IV). Multilocus sequence typing of a subset of isolates showed that most PVL-positive strains belonged to sequence type (ST) 188, while most PVL-negative isolates belonged to ST45. The high prevalence of PVL-positive S. aureus strains in certain regions of Indonesia is of concern since these strains may cause severe infections in the community and in hospitals.

Methicillin-resistant and -susceptible Staphylococcus aureus sequence type 398 in pigs and humans.

Methicillin-resistant Staphylococcus aureus sequence type 398 (ST398 MRSA) was identified in Dutch pigs and pig farmers. ST398 methicillin-susceptible S. aureus circulates among humans at low frequency (0.2%) but was isolated in 3 human cases of bacteremia (2.1%; p = 0.026). Although its natural host is probably porcine, ST398 MRSA likely causes infections in humans.

Host-microbe interplay in persistent Staphylococcus aureus nasal carriage in HIV patients.

It has been shown that persistent Staphylococcus aureus nasal carriage results in increased bacterial dispersal and a higher risk of infection compared to non-or-intermittent S. aureus carriage. Although many studies investigated S. aureus nasal carriage in HIV patients, none compared persistent carriage to non-persistent carriage nor were studies performed in the HAART era. We investigated the host-microbe interplay of persistent S. aureus nasal carriage in HIV-infected patients by studying host determinants of persistent carriage as well as the genetic structure of S. aureus strains isolated. We compared this genetic structure with the previously determined population structure of S. aureus isolates obtained from healthy individuals. Between February 2004 and June 2005 all HIV patients visiting the outpatient department of Erasmus MC (Rotterdam, The Netherlands) were asked to participate in this study. Participants were interviewed and screened for persistent S. aureus carriage using two semi-quantitative nasal swab cultures. For 443 patients two cultures were available, 131 (29.6%) were persistent carriers, which is significantly higher as compared to healthy individuals from the same geographic region (17.6%; P<0.0001). Male sex (odds ratio [OR], 2.22; 95% confidence interval [CI], 1.32-3.73), current smoking (OR, 0.58; 95% CI, 0.38-0.90), Pneumocystis jiroveci pneumonia (PCP) prophylaxis (OR, 0.39; 95% CI, 0.16-0.97) and antiretroviral therapy (OR, 0.61; 95% CI, 0.38-0.98) were independent determinants of persistent carriage. Only two strains were mecA positive (1.2%) and no PVL positive strains were detected. The population structure of S. aureus strains isolated from HIV patients appeared to be strongly overlapping with that of S. aureus isolates from healthy individuals.

Overlapping population structures of nasal isolates of Staphylococcus aureus from healthy Dutch and American individuals.

To understand Staphylococcus aureus nasal carriage and its relationship with subsequent disease, insight into the natural (nonclinical) bacterial population structure is essential. This study investigated whether the distributions of S. aureus genotypes that cause colonization differ by geographic locales. High-throughput amplified fragment length polymorphism (AFLP) analysis was performed on nasal isolates of S. aureus from healthy American (n = 391) and Dutch (n = 829) volunteers. In total, 164,970 binary outcomes, covering 135 different markers per isolate, were scored. Methicillin resistance was defined for all strains; pulsed-field gel electrophoresis typing was performed for the American isolates. The overall population structures of the American and Dutch S. aureus isolates were comparable. The same four major AFLP clusters (I to IV) and subclusters were identified for both collections. However, the Dutch methicillin-susceptible S. aureus (MSSA) isolates were overrepresented in AFLP cluster III (P = 0.0016). Furthermore, the majority of the American methicillin-resistant S. aureus isolates (90.5%) were located in AFLP cluster I (P < 0.0001). This result identifies differences in the local prevalence of certain S. aureus genotypes. AFLP clusters II and III, which represent multilocus sequence typing clonal complexes 30 and 45, respectively, account for 46.4% of all MSSA isolates in the study, suggesting that these two lineages have evolved as extremely successful pandemic colonizers of humans. In conclusion, the overall population structures of American and Dutch nasal carriage isolates of S. aureus are surprisingly similar, despite subtle geographic differences in the prevalence of certain S. aureus genotypes.

High-resolution genome profiling differentiated Staphylococcus epidermidis isolated from patients with ocular infections and normal individuals.

PURPOSE: To investigate the potential phenotypic and genetic differences among the Staphylococcus epidermidis isolates obtained from control subjects (lower conjunctival sac; n = 14) with those from patients with keratitis (corneal scrapings; n = 18) or endophthalmitis (vitreous; n = 24). METHODS: Biofilm-forming capacity was detected by PCR for the icaAB gene and phenotyping by microtiter plate assay and congo red agar plate. Genotyping was performed by using fluorescence-amplified fragment length polymorphism (FAFLP) and in silico analysis of the FAFLP profiles. RESULTS: Biofilm phenotyping (congo red agar/microtiter plate) differentiated disease-causing strains from control subjects. PCR assays (mecA, icaAB) were not useful in differentiating disease-causing strains from that of control subjects. The biofilm-forming capability appeared more critical in the pathogenesis of keratitis than in that of endophthalmitis. Cluster analysis of FAFLP data generated 11 clusters comprising 4 major clusters (I, II, III, and V) and 7 minor ones. FAFLP analysis clearly showed clustering of most of the commensal isolates in cluster I, separate from keratitis and endophthalmitis isolates. In silico analysis mapped signature bands to genes such as ebh, tagD, ptsI, and sepA, which might have a significant role in transforming less virulent populations of S. epidermidis to more virulent ones. CONCLUSIONS: The population dynamics of S. epidermidis revealed that there are significant genetic variations that can be detected through FAFLP between ocular disease causing isolates and the commensal population.

Comparison of multilocus sequence typing (MLST), pulsed-field gel electrophoresis (PFGE), and amplified fragment length polymorphism (AFLP) for genetic typing of Staphylococcus aureus.

We compared multilocus sequence typing (MLST), pulsed-field gel electrophoresis (PFGE), and amplified fragment length polymorphism (AFLP) for typing of Staphylococcus aureus and show that the methods yield similar results, although with differences in resolving power and reproducibility. Epidemiological conditions should determine which is the optimal typing method to be used.

Nasopharyngeal co-colonization with Staphylococcus aureus and Streptococcus pneumoniae in children is bacterial genotype independent.

Bacterial interference between Staphylococcus aureus and Streptococcus pneumoniae in the nasopharynx has been observed during colonization, which might have important clinical implications for the widespread use of pneumococcal conjugate vaccine in young children. This study aimed to determine whether the capacity of Staph. aureus to compete with Strep. pneumoniae is dependent on bacterial genotype. Demographic and microbiological determinants of carriage of specific genotypes of Staph. aureus in children were also studied. Children (n=3198) were sampled in the nasopharynx to detect carriage of Staph. aureus, Strep. pneumoniae and Neisseria meningitidis. Staph. aureus genotypes and pneumococcal sero- and genotypes were determined. Age, gender, zip code, active smoking and co-colonization with N. meningitidis or Strep. pneumoniae, both vaccine- and non-vaccine types, were not associated with colonization by specific Staph. aureus genotypes. Based on the whole-genome typing data obtained, there was no obvious correlation between staphylococcal and pneumococcal genotypes during co-colonization. Passive smoking showed a significant association (P=0.003) with carriage of a specific Staph. aureus cluster. This study suggests that there are no major differences between Staph. aureus clones (with different disease-invoking potential) in their capacity to compete with Strep. pneumoniae subtypes. Further studies should demonstrate whether differences in bacterial interference are due to more subtle genetic changes.

Exon expression arrays as a tool to identify new cancer genes.

BACKGROUND: Identification of genes that are causally implicated in oncogenesis is a major goal in cancer research. An estimated 10-20% of cancer-related gene mutations result in skipping of one or more exons in the encoded transcripts. Here we report on a strategy to screen in a global fashion for such exon-skipping events using PAttern based Correlation (PAC). The PAC algorithm has been used previously to identify differentially expressed splice variants between two predefined subgroups. As genetic changes in cancer are sample specific, we tested the ability of PAC to identify aberrantly expressed exons in single samples. PRINCIPAL FINDINGS: As a proof-of-principle, we tested the PAC strategy on human cancer samples of which the complete coding sequence of eight cancer genes had been screened for mutations. PAC detected all seven exon-skipping mutants among 12 cancer cell lines. PAC also identified exon-skipping mutants in clinical cancer specimens although detection was compromised due to heterogeneous (wild-type) transcript expression. PAC reduced the number of candidate genes/exons for subsequent mutational analysis by two to three orders of magnitude and had a substantial true positive rate. Importantly, of 112 randomly selected outlier exons, sequence analysis identified two novel exon skipping events, two novel base changes and 21 previously reported base changes (SNPs). CONCLUSIONS: The ability of PAC to enrich for mutated transcripts and to identify known and novel genetic changes confirms its suitability as a strategy to identify candidate cancer genes.