Gene co-expression network analysis is an effective method for predicting gene functions and disease biomarkers. However, few studies have systematically identified co-expressed genes involved in the molecular origin and development of various types of tumors. In this study, we used a network mining algorithm to identify tightly connected gene co-expression networks that are frequently present in microarray datasets from 33 types of cancer which were derived from 16 organs/tissues. We compared the results with networks found in multiple normal tissue types and discovered 18 tightly connected frequent networks in cancers, with highly enriched functions on cancer-related activities. Most networks identified also formed physically interacting networks. In contrast, only 6 networks were found in normal tissues, which were highly enriched for housekeeping functions. The largest cancer network contained many genes with genome stability maintenance functions. We tested 13 selected genes from this network for their involvement in genome maintenance using two cell-based assays. Among them, 10 were shown to be involved in either homology-directed DNA repair or centrosome duplication control including the well- known cancer marker MKI67. Our results suggest that the commonly recognized characteristics of cancers are supported by highly coordinated transcriptomic activities. This study also demonstrated that the co-expression network directed approach provides a powerful tool for understanding cancer physiology, predicting new gene functions, as well as providing new target candidates for cancer therapeutics.

Ran Binding Protein 9 (RanBP9, also known as RanBPM) is an evolutionary conserved scaffold protein present both in the nucleus and the cytoplasm of cells whose biological functions remain elusive.We show that active ATM phosphorylates RanBP9 on at least two different residues (S181 and S603). In response to IR, RanBP9 rapidly accumulates into the nucleus of lung cancer cells, but this nuclear accumulation is prevented by ATM inhibition. RanBP9 stable silencing in three different lung cancer cell lines significantly affects the DNA Damage Response (DDR), resulting in delayed activation of key components of the cellular response to IR such as ATM itself, Chk2, γH2AX, and p53. Accordingly, abrogation of RanBP9 expression reduces homologous recombination-dependent DNA repair efficiency, causing an abnormal activation of IR-induced senescence and apoptosis.In summary, here we report that RanBP9 is a novel mediator of the cellular DDR, whose accumulation into the nucleus upon IR is dependent on ATM kinase activity. RanBP9 absence hampers the molecular mechanisms leading to efficient repair of damaged DNA, resulting in enhanced sensitivity to genotoxic stress. These findings suggest that targeting RanBP9 might enhance lung cancer cell sensitivity to genotoxic anti-neoplastic treatment.

Genes associated with hereditary breast and ovarian cancer (HBOC) are often sequenced in search of mutations that are predictive of susceptibility to these cancer types, but the sequence results are frequently ambiguous due to detection of missense substitutions for which the clinical impact is unknown. The BARD1 protein is the heterodimeric partner of BRCA1 and is included on clinical gene panels for testing for susceptibility to HBOC. Like BRCA1, it is required for homology-directed DNA repair (HDR). We measured the HDR function of 29 BARD1 missense variants, 27 culled from clinical test results and two synthetic variants. 23 of the assayed variants were functional for HDR; of these, four are known neutral variants. Three variants showed intermediate function, and three others were defective in HDR. When mapped to BARD1 domains, residues crucial for HDR were located in the N- and C- termini of BARD1. In the BARD1 RING domain, critical residues mapped to the zinc-coordinating amino acids and to the BRCA1-BARD1 binding interface, highlighting the importance of interaction between BRCA1 and BARD1 for HDR activity. Based on these results, we propose that the HDR assay is a useful complement to genetic analyses to classify BARD1 variants of unknown clinical significance.

The BARD1 protein, which heterodimerizes with BRCA1, is encoded by a known breast cancer susceptibility gene. While several BARD1 variants have been identified as pathogenic, many more missense variants exist that do not occur frequently enough to assign a clinical risk. In this paper, whole exome sequencing of over 10,000 cancer samples from 33 cancer types identified from somatic mutations and loss of heterozygosity in tumors 76 potentially cancer-associated BARD1 missense and truncation variants. These variants were tested in a functional assay for homology-directed repair (HDR), as HDR deficiencies have been shown to correlate with clinical pathogenicity for BRCA1 variants. From these 76 variants, 4 in the ankyrin repeat domain and 5 in the BRCT domain were found to be non-functional in HDR. Two known benign variants were found to be functional in HDR, and three known pathogenic variants were non-functional, supporting the notion that the HDR assay can be used to predict the clinical risk of BARD1 variants. The identification of HDR-deficient variants in the ankyrin repeat domain indicates there are DNA repair functions associated with this domain that have not been closely examined. In order to examine whether BARD1-associated loss of HDR function results in DNA damage sensitivity, cells expressing non-functional BARD1 variants were treated with ionizing radiation or cisplatin. These cells were found to be more sensitive to DNA damage, and variations in the residual HDR function of non-functional variants did not correlate with variations in sensitivity. These findings improve the understanding of BARD1 functional domains in DNA repair and support that this functional assay is useful for predicting the cancer association of BARD1 variants.