Prognostic value of microvessel density in stage II and III colon cancer patients:a retrospective cohort study

Background Microvessel density (MVD), as a derived marker for angiogenesis, has been associated with poor outcome in several types of cancer. This study aimed to evaluate the prognostic value of MVD in stage II and III colon cancer and its relation to tumour-stroma-percentage (TSP) and expression of HIF1A and VEGFA. Methods Formalin-fixed paraffin-embedded (FFPE) colon cancer tissues were collected from 53 stage II and 54 (5-fluorouracil-treated) stage III patients. MVD was scored by digital morphometric analysis of CD31-stained whole tumour sections. TSP was scored using haematoxylin-eosin stained slides. Protein expression of HIF1A and VEGFA was determined by immunohistochemical evaluation of tissue microarrays. Results Median MVD was higher in stage III compared to stage II colon cancers (11.1% versus 5.6% CD31-positive tissue area, p <0.001). High MVD in stage II patients tended to be associated with poor disease free survival (DFS) in univariate analysis (p = 0.056). In contrast, high MVD in 5FU-treated stage III patients was associated with better DFS (p = 0.006). Prognostic value for MVD was observed in multivariate analyses for both cancer stages. Conclusions MVD is an independent prognostic factor associated with poor DFS in stage II colon cancer patients, and with better DFS in stage III colon cancer patients treated with adjuvant chemotherapy

Fusion transcripts and their genomic breakpoints in polyadenylated and ribosomal RNA-minus RNA sequencing data

BACKGROUND: Fusion genes are typically identified by RNA sequencing (RNA-seq) without elucidating the causal genomic breakpoints. However, non–poly(A)-enriched RNA-seq contains large proportions of intronic reads that also span genomic breakpoints. RESULTS: We have developed an algorithm, Dr. Disco, that searches for fusion transcripts by taking an entire reference genome into account as search space. This includes exons but also introns, intergenic regions, and sequences that do not meet splice junction motifs. Using 1,275 RNA-seq samples, we investigated to what extent genomic breakpoints can be extracted from RNA-seq data and their implications regarding poly(A)-enriched and ribosomal RNA–minus RNA-seq data. Comparison with whole-genome sequencing data revealed that most genomic breakpoints are not, or minimally, transcribed while, in contrast, the genomic breakpoints of all 32 TMPRSS2-ERG–positive tumours were present at RNA level. We also revealed tumours in which the ERG breakpoint was located before ERG, which co-existed with additional deletions and messenger RNA that incorporated intergenic cryptic exons. In breast cancer we identified rearrangement hot spots near CCND1 and in glioma near CDK4 and MDM2 and could directly associate this with increased expression. Furthermore, in all datasets we find fusions to intergenic regions, often spanning multiple cryptic exons that potentially encode neo-antigens. Thus, fusion transcripts other than classical gene-to-gene fusions are prominently present and can be identified using RNA-seq. CONCLUSION: By using the full potential of non–poly(A)-enriched RNA-seq data, sophisticated analysis can reliably identify expressed genomic breakpoints and their transcriptional effects

Consensus molecular subtype classification of colorectal adenomas

Consensus molecular subtyping is an RNA expression-based classification system for colorectal cancer (CRC). Genomic alterations accumulate during CRC pathogenesis, including the premalignant adenoma stage, leading to changes in RNA expression. Only a minority of adenomas progress to malignancies, a transition that is associated with specific DNA copy number aberrations or microsatellite instability (MSI). We aimed to investigate whether colorectal adenomas can already be stratified into consensus molecular subtype (CMS) classes, and whether specific CMS classes are related to the presence of specific DNA copy number aberrations associated with progression to malignancy. RNA sequencing was performed on 62 adenomas and 59 CRCs. MSI status was determined with polymerase chain reaction-based methodology. DNA copy number was assessed by low-coverage DNA sequencing (n = 30) or array-comparative genomic hybridisation (n = 32). Adenomas were classified into CMS classes together with CRCs from the study cohort and from The Cancer Genome Atlas (n = 556), by use of the established CMS classifier. As a result, 54 of 62 (87%) adenomas were classified according to the CMS. The CMS3 ‘metabolic subtype’, which was least common among CRCs, was most prevalent among adenomas (n = 45; 73%). One of the two adenomas showing MSI was classified as CMS1 (2%), the ‘MSI immune’ subtype. Eight adenomas (13%) were classified as the ‘canonical’ CMS2. No adenomas were classified as the ‘mesenchymal’ CMS4, consistent with the fact that adenomas lack invasion-associated stroma. The distribution of the CMS classes among adenomas was confirmed in an independent series. CMS3 was enriched with adenomas at low risk of progressing to CRC, whereas relatively more high-risk adenomas were observed in CMS2. We conclude that adenomas can be stratified into the CMS classes. Considering that CMS1 and CMS2 expression signatures may mark adenomas at increased risk of progression, the distribution of the CMS classes among adenomas is consistent with the proportion of adenomas expected to progress to CRC

Killer immunoglobulin-like Receptors (KIR) haplogroups A and B track with Natural Killer Cells and Cytokine Profile in Aged Subjects: Observations from Octo/Nonagenarians in the Belfast Elderly Longitudinal Free-living Aging STudy (BELFAST)

BACKGROUND: Natural Killer Cells (NK) play an important role in detection and elimination of virus-infected, damaged or cancer cells. NK cell function is guided by expression of Killer Immunoglobulin-like Receptors (KIRs) and contributed to by the cytokine milieu. KIR molecules are grouped on NK cells into stimulatory and inhibitory KIR haplotypes A and B, through which NKs sense and tolerate HLA self-antigens or up-regulate the NK-cytotoxic response to cells with altered HLA self-antigens, damaged by viruses or tumours. We have previously described increased numbers of NK and NK-related subsets in association with sIL-2R cytokine serum levels in BELFAST octo/nonagenarians. We hypothesised that changes in KIR A and B haplotype gene frequencies could explain the increased cytokine profiles and NK compartments previously described in Belfast Elderly Longitudinal Free-living Aging STudy (BELFAST) octo/nonagenarians, who show evidence of ageing well. RESULTS: In the BELFAST study, 24% of octo/nonagenarians carried the KIR A haplotype and 76% KIR B haplotype with no differences for KIR A haplogroup frequency between male or female subjects (23% v 24%; p=0.88) or for KIR B haplogroup (77% v 76%; p=0.99). Octo/nonagenarian KIR A haplotype carriers showed increased NK numbers and percentage compared to Group B KIR subjects (p=0.003; p=0.016 respectively). There were no KIR A/ B haplogroup-associated changes for related CD57+CD8 ((high or low)) subsets. Using logistic regression, KIR B carriers were predicted to have higher IL-12 cytokine levels compared to KIR A carriers by about 3% (OR 1.03, confidence limits CI 0.99–1.09; p=0.027) and 14% higher levels for TGF-β (active), a cytokine with an anti-inflammatory role, (OR 1.14, confidence limits CI 0.99–1.09; p=0.002). CONCLUSION: In this observational study, BELFAST octo/nonagenarians carrying KIR A haplotype showed higher NK cell numbers and percentage compared to KIR B carriers. Conversely, KIR B haplotype carriers, with genes encoding for activating KIRs, showed a tendency for higher serum pro-inflammatory cytokines compared to KIR A carriers. While the findings in this study should be considered exploratory they may serve to stimulate debate about the immune signatures of those who appear to age slowly and who represent a model for good quality survivor-hood

GeneBreak: detection of recurrent DNA copy number aberration-associated chromosomal breakpoints within genes

Development of cancer is driven by somatic alterations, including numerical and structural chromosomal aberrations. Currently, several computational methods are available and are widely applied to detect numerical copy number aberrations (CNAs) of chromosomal segments in tumor genomes. However, there is lack of computational methods that systematically detect structural chromosomal aberrations by virtue of the genomic location of CNA-associated chromosomal breaks and identify genes that appear non-randomly affected by chromosomal breakpoints across (large) series of tumor samples. ‘GeneBreak’ is developed to systematically identify genes recurrently affected by the genomic location of chromosomal CNA-associated breaks by a genome-wide approach, which can be applied to DNA copy number data obtained by array-Comparative Genomic Hybridization (CGH) or by (low-pass) whole genome sequencing (WGS). First, ‘GeneBreak’ collects the genomic locations of chromosomal CNA-associated breaks that were previously pinpointed by the segmentation algorithm that was applied to obtain CNA profiles. Next, a tailored annotation approach for breakpoint-to-gene mapping is implemented. Finally, dedicated cohort-based statistics is incorporated with correction for covariates that influence the probability to be a breakpoint gene. In addition, multiple testing correction is integrated to reveal recurrent breakpoint events. This easy-to-use algorithm, ‘GeneBreak’, is implemented in R (www.cran.r-project.org) and is available from Bioconductor (www.bioconductor.org/packages/release/bioc/html/GeneBreak.html)