Network analysis reveals a major role for 14q32 cluster miRNAs in determining transcriptional differences between IGHV-mutated and unmutated CLL

Abstract

Tumour cells from patients with chronic lymphocytic leukaemia (CLL) can express unmutated (U-CLL) or mutated (M-CLL) immunoglobulin heavy chain (IGHV) genes with differing clinical behaviours, variable B cell receptor (BCR) signalling capacity and distinct transcriptional profiles. As it remains unclear to what extent these differences reflect the tumour cells innate pre/post germinal centre origin or their BCR signalling competence, we applied RNA sequencing, small RNA sequencing and DNA methylation array analysis to 38 CLL cases categorised into three groups by IGHV mutational status and BCR signalling capacity. We identified 492 mRNAs and 38 miRNAs differentially expressed between U-CLL and M-CLL, but only 9 mRNAs and 0 miRNAs associated with BCR competence within M-CLL. A significant proportion of the IGHV-associated miRNAs derived from chr14q32 clusters (14/38 (37%)), where all miRNAs were co-expressed with the MEG3 lncRNA, as part of the DLK1-DIO3 genomic imprinted region, a locus of known importance in the pathogenesis of other human tumours. Integrative in silico analysis of miRNA/mRNA data revealed pronounced regulatory potential for the 14q32 miRNAs, potentially accounting for up to 25% of the IGHV-related transcriptome signature. GAB1, a positive regulator of BCR signalling, was predicted to be regulated by five 14q32 miRNAs and we confirmed that two of these (miR-409-3p and miR-411-3p) significantly repressed activity of the GAB1 3UTR. Our analysis demonstrates a potential key role of the 14q32 miRNA locus in the regulation of CLL-related gene regulation.