Integrative expression analysis identifies a novel interplay between CFTR and linc-SUMF1-2 that involves CF-associated gene dysregulation.

Cystic fibrosis transmembrane regulator (CFTR) is a cyclic AMP-dependent Cl channel, and its dysfunction, due to CFTR gene mutations, causes the lethal inherited disorder cystic fibrosis (CF). To date, widespread dysregulation of certain coding genes in CF airway epithelial cells is well studied and considered as the driver of pulmonary abnormality. However, the involvement of non-coding genes, novel classes of functional RNAs with little or no protein-coding capacity, in the regulation of CF-associated gene dysregulation is poorly understood. Here, we utilized integrative analyses of human transcriptome array (HTA) and characterized 99 coding and 91 non-coding RNAs that are dysregulated in CFTR-defective CF bronchial epithelial cell line CFBE41o-. Among these genes, the expression level of linc-SUMF1-2, an intergenic non-coding RNA (lincRNA) whose function is unknown, was inversely correlated with that of WT-CFTR and consistently higher in primary human CF airway epithelial cells (DHBE-CF). Further integrative analyses under linc-SUMF1-knockdown condition determined MXRA5, SEMA5A, CXCL10, AK022877, CTGF, MYC, AREG and LAMB3 as both CFTR- and linc-SUMF1-2-dependent dysregulated gene sets in CF airway epithelial cells. Overall, our analyses reveal linc-SUMF1-2 as a dysregulated non-coding gene in CF as well as CFTR-linc-SUMF1-2 axis as a novel regulatory pathway involved in CF-associated gene dysregulation.