Next Article in Journal
Impact of Genomics on Clarifying the Evolutionary Relationships amongst Mycobacteria: Identification of Molecular Signatures Specific for the Tuberculosis-Complex of Bacteria with Potential Applications for Novel Diagnostics and Therapeutics
Previous Article in Journal
Current High-Throughput Approaches of Screening Modulatory Effects of Xenobiotics on Cytochrome P450 (CYP) Enzymes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
High-Throughput
Volume 7
Issue 4
10.3390/ht7040030
high-throughput-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
Please note that, as of 21 September 2020, High-Throughput has been renamed to BioTech and is now published here.
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Lesions of Insulator CTCF and Its Paralogue CTCFL (BORIS) in Cancer: An Analysis from Published Genomic Studies

by
Ioannis A. Voutsadakis
1,2
1
Algoma District Cancer Program, Sault Area Hospital, Sault Ste. Marie, ON P6B 0A8, Canada
2
Section of Internal Medicine, Division of Clinical Sciences, Northern Ontario School of Medicine, Sudbury, ON P3E 2C6, Canada
High-Throughput 2018, 7(4), 30; https://doi.org/10.3390/ht7040030
Submission received: 8 August 2018 / Revised: 10 September 2018 / Accepted: 26 September 2018 / Published: 1 October 2018
Browse Figures
Versions Notes

Abstract

:
CTCF (CCCTC-binding factor) is a transcription regulator with hundreds of binding sites in the human genome. It has a main function as an insulator protein, defining together with cohesins the boundaries of areas of the genome called topologically associating domains (TADs). TADs contain regulatory elements such as enhancers which function as regulators of the transcription of genes inside the boundaries of the TAD while they are restricted from regulating genes outside these boundaries. This paper will examine the most common genetic lesions of CTCF as well as its related protein CTCFL (CTCF-like also called BORIS) in cancer using publicly available data from published genomic studies. Cancer types where abnormalities in the two genes are more common will be examined for possible associations with underlying repair defects or other prevalent genetic lesions. The putative functional effects in CTCF and CTCFL lesions will also be explored.

1. Introduction

The three-dimensional organisation of DNA in interphase cell nuclei is important for the regulation of gene transcription [1]. DNA in human chromosomes is organised into higher-order domains called topologically associating domains (TADs) and these have subdomains called insulated neighbourhoods [2,3]. Protein CTCF (CCCTC-binding factor, also known as MRD21, Mental Retardation 21) defines borders of these domains in human genome.
The CTCF gene (Gene ID: 10664, Ensembl gene: ENSG00000102974) is located at human chromosome 16q22.1 region and has 13 exons spanning over 76,779 nucleotides according to GenBank. Its specific coordinates in the human genome assembly GRCh38 are chromosome 16: 67.562.407–67.639.185, forward strand. Alternative mRNA splicing produces different isoforms with 17 variants listed in Ensembl. Two promoters of the human gene are listed in the Eukaryotic Promoter database [4].
CTCF protein (UniProtKB-P49711) has a length of 727 amino acids and is comprised of 11 C2H2-type zinc fingers (ZFs) occupying the central portion of the protein from amino acid 266 to amino acid 577 and flanked by an aminoterminal and a carboxyterminal domain that are both unstructured (Figure 1) [5]. CTCF ZFs are 23 to 24 amino acids long and accommodate the zinc atom through two cysteine and two histidine residues. CTCF is expressed ubiquitously in all human adult tissues. CTCF functions as a specific DNA sequence binding protein and has a role in gene transcription regulation both as a suppressor and an activator and regulates gene imprinting. CTCF plays also a major additional role in gene regulation by acting as a DNA topologically associating domains (TADs) insulator. In this capacity, it associates with the cohesin complex and another CTCF factor bound at a distance creating DNA loops. In many occasions, transcriptional activities mediated by enhancers are restricted inside these loops, insulating these enhancers from acting in transcription of genes outside the loop. TADs may also possess functional subdomains termed insulated neighbourhoods also defined by CTCF borders. Mutations affecting TADs or neighbourhoods borders in cancer may have profound effects in gene regulation by creating new influences by enhancers that were originally outside the TAD or outside the neighbourhood or conversely restricting enhancers from exerting normal regulations. There are thousands of potential CTCF binding sites in the human genome (possibly in the range of 11,000 to 14,000 and up to 60,000 in some studies), although it is unclear what percentage of these identify potential insulated neighbourhood borders [6,7,8]. The function of CTCF in gene imprinting is commonly illustrated with the example of the IGF2/H19 locus where, under the influence of CTCF differential binding, the IGF2 gene is only expressed by the paternal allele while H19 is exclusively expressed by the maternal allele. This involves also the insulator function of CTCF. Binding of CTCF to the locus is dependent on the DNA methylation status of CTCF binding sites in the neighbourhood of the two alleles which modifies the length of the insulated neighbourhood that includes the two genes and includes or excludes different enhancers from regulating each of the two genes [9].
CTCF paralogue, CTCFL (CCCTC-binding factor-Like, also known as BORIS, Brother of the Regulator of Imprinted Sites or CT27: Cancer/Testis antigen 27) is transcribed from a gene in human chromosome 20q13.31 (Gene ID: 140690, Ensembl gene: ENSG00000124092). CTCFL specific coordinates in the human genome assembly GRCh38 are chromosome 20: 57.495.966–57.525.652, reverse strand. It has 16 exons and 22 splice variants according to Ensembl.
The CTCFL protein (UniProtKB-Q8NI51) consists of 663 amino acids and, in contrast to the ubiquitous expression of its paralogue CTCF, is normally expressed only during spermatogenesis [10]. Similarly to CTCF, CTCFL possesses 11 C2H2-type zinc fingers that are highly homologous to those of CTCF, are 23 to 24 amino acids long, and occupy the middle portion of the protein from amino acid 257 to amino acid 568. The DNA binding sequence of CTCFL and CTCF is very similar but their capacity to interact with partner proteins is not conserved due to significant divergence in their aminoterminal and carboxyterminal domains. Most CTCFL binding sites are shared with CTCF but the reverse is not true as CTCF binds to ten times more sites [11]. The two paralogues may have differences in their ability to bind methylated sites with CTCFL being able to bind methylated sites while CTCF preferring unmethylated sites [12]. However, mere loss of methylation in a target site is not sufficient for CTCF binding in most occasions [13]. Normal CTCFL expression is restricted to spermatogonia and preleptotene spermatocytes. CTCFL expression becomes silenced in late spermatogenesis through promoter methylation and remains absent or very low in most adult tissues. Expression is reactivated in some cancer cases through promoter hypomethylation. Thus, CTCFL belongs to the category of so called cancer testis antigens, alternatively termed cancer-germline antigens [14]. Although details of the interaction are not known, CTCFL has been reported to be part of the CTCF interactome [15]. Whether the two proteins interact directly or indirectly in cells where they are co-expressed remains to be confirmed experimentally. The two paralogues have also been reported to cooperate in binding on tandem sites and thus it is possible that CTCF binding facilitates preferential CTCFL binding in an adjacent site [16].
This paper will investigate molecular lesions of CTCF and CTCFL in various cancers from published sources. Underlying molecular defects putatively associated with development of CTCF and CTCFL anomalies as well as prognostic implications of CTCF and CTCFL mRNA dosage will be explored.

2. Methods

Genomic studies of common cancers were interrogated in the cBioportal platform [17,18,19] for genetic lesions and mRNA dysregulation of the two paralogous genes of interest, CTCF and CTCFL. cBioportal contains several of the most extensive series of genomic studies performed by The Cancer Genome Atlas (TCGA) and other groups. The platform currently contains 226 genomic studies and allows for interrogation of each study for genetic lesions in any gene of interest. Genomic studies included in the cBioportal platform were examined for frequency and specific characteristics of cases with CTCF and CTCFL mutations and copy number alterations. Series with the higher absolute number of CTCF and CTCFL lesions were identified and examined in more detail to establish correlations with protein domain localisation of mutations and resulting total mutation burden. Studies selected for more detailed scrutiny included TCGA studies for endometrial, bladder, colon, and gastroesophageal carcinomas and the METABRIC breast cancer study [20,21,22,23,24]. These studies either contain the higher percentage of cases with CTCF and CTCFL lesions or, despite a lower percentage of lesions in these genes, the absolute numbers of defective cases were substantial in order to facilitate analysis. Studies of several other common types of cancers available in cBioportal were reviewed to determine frequency of CTCF and CTCFL defects [25,26,27,28,29,30,31,32].
Identified mutations were mapped in the different regions of each gene and assessed for their putative functional significance using the mutation assessor server, Computational Biology Center, Memorial Sloan Kettering Cancer Center, New York, NY, U.S.A. [33] which uses a multiple sequence alignment (msa) algorithm to assign a prediction score of functional significance to each mutation [34]. Additional investigations performed included identification of presence of MSI-related genes (MSH2, MSH6, PMS2, and MLH1) and polymerase δ and ε (POLD1 and POLE) defects in mutated samples and identification of the most commonly amplified region (amplicon) in samples with amplifications.
Survival of patients with high expression of CTCF and CTCFL mRNA versus those with low CTCF and CTCFL mRNA expression in examples of gastric, breast, and ovarian cancers was compared using the online tool Kaplan Meier Plotter [35,36].
Promoter methylation status was examined using published TCGA data and the online database for DNA methylation in cancer [37,38]. This database provides comparisons of methylation of promoter sequences of cancer samples of various cancers with corresponding sequences of respective normal tissues.
Categorical and continuous data were compared with the Fisher’s exact test and the t-test respectively. Correlations were explored with the Pearson correlation coefficient. All statistical comparisons were considered significant if p < 0.05.

3. Results

3.1. Molecular Lesions of CTCF in Cancer

An overview of lesions of CTCF gene in various cancers studied by TCGA and in the METABRIC breast cancer study shows that the most common type of genetic lesions is mutations (1.97% of all samples examined), while amplifications and deep deletions were very rare (0.18% and 0.49%, respectively) (Figure 2a and Table 1). Cancers presenting with the higher percentage of lesions in CTCF: Uterine endometrial carcinomas (37.25% of total samples have CTCF lesions), ovarian serous carcinomas (16.5% of total samples with CTCF lesions), bladder carcinomas (13.61% of total samples with CTCF lesions), colorectal carcinomas (11.24% of total samples with CTCF lesions), and prostate cancers (10.79% of total samples with CTCF lesions) (Table 1). However, even mutations are rare and observed in less than 2% of samples in most types of cancer with the exception of endometrial, gastroesophageal, colorectal, bladder, and breast cancers (Figure 2b). Importantly, the type of cancer that stands out as having the highest mutation rate of CTCF is endometrial cancer, where CTCF is mutated in more than one fourth of tumours (27.45%). Uterine carcinosarcomas and gastric cancers display a mutation rate of CTCF of approximately 5% followed by colorectal, bladder, and breast cancers which display mutation rates of CTCF in 4.87%, 3.22%, and 2.21% of cases examined, respectively (Figure 2b).
Most mutations (59%) of CTCF gene in endometrial carcinomas in the uterine TCGA PanCancer study were located in the area encoding the 11 ZFs of the protein (23 of 39 samples with CTCF mutations in cases that had complete mutations, copy number alterations, and mRNA expression analysis data available, Table 2) [23]. The rest of the samples, except one that had a mutation in the C-terminal domain, had mutations in the N-terminal domain. Twenty-six of 39 CTCF gene mutated samples (66.7%) had one or more mutations in one of the four microsatellite instability (MSI)-associated genes (MSH2, MSH6, PMS2, and MLH1). An alternative cause of hypermutation in cancer is mutations in polymerases epsilon (POLE) and delta 1 (POLD1) [39]. Among the 39 CTCF mutated samples 23 samples had a concomitant mutation in one of these polymerases (Table 2). Overall 29 samples of the 39 CTCF mutated samples (74.4%) had mutations in MSI or the two polymerases. The mean number of total mutations in samples with mutations in MSI genes or the two polymerases was over 5000 while the 10 samples without mutations in these genes had a mean number of 301 mutations. These data suggest that CTCF mutations are commonly but not exclusively seen in MSI-associated endometrial carcinomas. A recurrent CTCF mutation in the N-terminal domain of CTCF in endometrial carcinomas observed in six samples was a frameshift mutation at codon T204 producing truncation of the protein after 18 or 26 amino acids. This recurrent mutation was not always associated with MSI or polymerase mutations. In only three of these cases there were concomitant MSI-associated gene mutations and one had polymerase mutations (Table 2) suggesting that even the exact same mutation may be associated or caused by various underlying molecular defects. Interestingly, no association was observed with APOBEC3 mRNA upregulation, which is also a cause of mutation induction in cancer [40,41]. This is a gene encoding for a DNA cytosine deaminase physiologically involved in the innate immune system-mediated protection against retroviruses and retrotransposons. Its function promotes mutagenesis through deamination of cytidines to uracils [40].
In colorectal cancer CTCF alterations are observed overall in 11.24% of cases and mutations are less frequent (4.87%). In the 13 samples with CTCF mutations in the Colorectal TCGA PanCancer study cohort [21] most mutations (12 of 13, 92.3%) were located in the ZFs or the N-terminal domain (Table 3). Most samples (10 of 13, 76.9%) had also mutations in one of the MSI-associated genes or polymerases POLE and POLD1 or both. Seven samples had over a thousand mutations per sample and all of these seven samples had mutations in one or both polymerases. Interestingly, many of the CTCF mutated samples, including two of the three samples without MSI-associated/polymerases mutations, had mutations in APOBEC genes (Table 3).
Urothelial carcinomas from the Bladder TCGA PanCancer study [22] was analysed in more detail as an example of cancer not usually associated with MSI (Table 4). The total number of samples with MSI-associated mutations in this study was 35 (8.6%). Among the 13 samples with CTCF mutations only three (23.1%) had mutations in MSI-associated genes (one of those with concomitant POLE mutation). Five of the 13 samples with CTCF mutations had mRNA upregulation of one of the APOBEC genes or of AID (Activation Induced Deaminase), a deaminase of the same family.
The putative functional significance of CTCF mutations in endometrial cancer and in breast cancer (as an example of a cancer not associated with MSI-associated mutations) were evaluated using mutation accessor and OncoKB. Among 91 different CTCF mutations found in endometrial carcinoma in TCGA, 45 (49.45%) are listed as likely oncogenic (the rest has unknown oncogenic potential and some of them may prove to be oncogenic as data accumulate). The METABRIC breast cancer study identified 44 different CTCF mutations (2.1% of samples) of which 20 mutations (45.5%) are considered likely oncogenic. The OncoKB database of mutations maintained by Memorial Sloan Kettering Cancer Center (oncokb.org) lists six point mutations of CTCF (H284N/Y/P, R339W, R377H, and P378L) and truncating mutations and deletions as oncogenic because of likely protein loss of function (or switch of function in the case of R339W).
The prognostic significance for survival of mRNA levels of CTCF in gastric, breast, and ovarian cancer were checked using Kaplan Meier plotter [31,35]. Gastric cancer patients with a high CTCF mRNA expression had an improved OS compared with counterparts with low CTCF mRNA expression (HR = 0.76, 95% CI = 0.64–0.91, p = 0.0022, Figure 3a). In breast cancer, overall across subtypes, CTCF mRNA levels are not associated with differences in OS (Figure 3b). However, in HER2-positive disease patients with a high CTCF mRNA expression displayed a worse OS compared with counterparts with low CTCF mRNA expression (HR = 1.66, 95% CI = 1.07–2.59, p = 0.023, Figure 3c). Similarly, in ER-positive patients high CTCF mRNA expression is associated with worse OS (not shown). In stage I and II ovarian cancer patients with a high CTCF mRNA expression also displayed a worse OS compared with counterparts with low CTCF mRNA expression (HR = 1.96, 95% CI = 1.02–3.78, p = 0.039, not shown). In contrast, in stage III and IV ovarian cancers there is no difference in OS between high and low CTCF mRNA levels groups (p = 0.13). These data suggest that CTCF levels are not directly associated with prognosis or may have different prognostic implications depending on the particular tumour.

3.2. Lesions of CTCFL in Cancer

Molecular lesions of CTCFL gene were observed in less than 10% of common tumours examined (Figure 4). However, several types of cancer had lesions in more than 6% of cases (Figure 4a). In addition to mutations, amplifications were common in CTCFL. Amplifications were the almost exclusive CTCFL lesion in breast cancer and constituted a significant percentage of CTCFL molecular lesions in ovarian, colon, and gastric carcinomas as well as uterine carcinosarcomas. On the other hand, melanomas and endometrial cancers showed mutations as the dominant type of lesion (Figure 4).
In the 10 samples with CTCFL mutations among samples that had complete mutations and copy number alterations information in the Colorectal PanCancer Atlas cohort [15] (Table 5), all had concomitant mutations in one of the four MSI associated genes or the POLE or POLD1 genes, while only one of the 25 samples with CTCFL amplifications (used as a control) had such lesions in MSI-associated genes or POLE/POLD1 genes (in POLD1) (Fisher’s two-tailed exact test p < 0.0001). In the respective endometrial carcinoma TCGA PanCancer Atlas study, which included 509 cases with complete mutations and copy number alterations information [23], among the 27 cases with CTCFL mutations 22 had mutations in one of the four MSI associated genes or the POLE or POLD1 genes. In the nine samples with CTCFL amplifications two had mutations in those genes (one of the two had a concomitant CTCFL mutation) (Fisher’s two-tailed exact test p = 0.0025). In the cutaneous melanoma PanCancer Atlas study (provisional), 11 of 27 samples (40.7%) with CTCFL mutations had lesions in one of the four MSI associated genes or the POLE or POLD1 genes. Thus, it appears that CTCFL mutations are often produced by underlying MSI or POLE or POLD1 defects in both cancers commonly associated with these defects (colorectal and endometrial) and in other cancers less commonly associated with them (melanoma).
Amplifications of CTCFL do not always correlate with increased CTCFL mRNA expression. For example, mRNA expression in breast and ovarian cancers is not significantly increased in amplified cases (Figure 5a,b). In contrast in colon cancer CTCFL amplified samples display a higher CTCFL mRNA expression (Figure 5c). The mean normalised CTCFL mRNA expression of diploid colon cases was 19.46 (SD: 7.97) and the mean normalised CTCFL mRNA expression of amplified cases was 29.79 (SD: 10.84, t = 5.4, p = 0.001). However, there was no correlation of mean normalised CTCFL mRNA expression values with the Log2 copy number values in either breast or ovarian or colon cancers (Pearson correlation p 0.06, 0.07, and 0.26, respectively). As a comparison, amplifications of the ERBB2 gene (encoding for the HER2 protein) in breast cancer result in increased mRNA expression compared with nonamplified tumours for ERBB2 (Figure 5d). Moreover, CTCFL protein is rarely expressed in cancers, in contrast to CTCF that is ubiquitously expressed. The Human Protein Atlas [42] records an absence of CTCFL expression in all cancers examined, including cancers with comparatively high rate of CTCFL amplification such as gastroesophageal, breast, and colon (Figure 4b).
The amplicon of CTCFL at chromosome 20q13.31–32 contains, in addition to CTCFL, genes: PCK1, PMEPA1 (also known as STAG1), ZBP1, BMP7, MIR4325, MTRNR2L3, RAE1, RBM38, SPO11, MIR4532, C20ORF85, and ANKRD60. All genes in the amplicon are amplified in a similar percentage of cases in different series, albeit in variable levels across cancers. In breast cancer, for example, all amplicon genes are amplified in approximately 6% to 7% in the TCGA PanCancer Atlas study and in 7% to 8% of cases in the METABRIC study (Figure S1) [19]. Thus, there is no clear indication of whether there is a driver gene among the amplicon genes that favours the amplification by promoting cancer cell fitness. Interestingly, some genes of the amplicon for which data are available in the human protein atlas, such as PCK1, PMEPA1, and RBM38, are expressed in the protein level at low to moderate levels in several cancers. Of additional interest is that the locus of 20q13 chromosomal region is commonly amplified in cancers but various subregions in this locus may be part of different amplicons. As an example, from the METABRIC study, Figure 6 shows that occurrence of CTCFL amplifications and amplifications of zinc finger transcription factor ZNF217 (another zinc finger transcription factor located in a neighbouring locus at 20q13.2 and proposed to be an oncogene) are only partially overlapping in samples of breast cancer despite both being present in approximately 8% of cases. These amplifications are also partially overlapping with amplifications of ERBB2 encoding for HER2 (Figure 6).
Using the Kaplan Meier plotter [35] the prognostic significance of mRNA levels of CTCFL in gastric, breast, and ovarian cancer were interrogated, similarly to the respective levels of CTCF. In gastric cancer high CTCFL mRNA expression levels are associated with worse OS compared with gastric cancers having low CTCFL mRNA expression (HR = 1.7, 95% CI = 1.35–2.13, p = 0.000, Figure 7a). In breast cancer, independently of subtypes, survival of patients with high CTCFL mRNA levels is not different from OS of patients with low CTCFL mRNA levels (Figure 7b). HER2-positive breast cancers with a high CTCFL mRNA expression have a trend towards worse OS compared with counterparts with low CTCFL mRNA expression (HR = 1.67, 95% CI = 0.96–2.89, p = 0.065, Figure 7c). Patients with ER-positive breast cancer and high CTCFL mRNA expression have no difference in survival compared with low CTCFL mRNA expression counterparts (not shown). In stage I and II ovarian cancer patients with a high CTCFL mRNA expression displayed a worse OS compared with counterparts with low CTCF mRNA expression (HR = 2.66, 95% CI = 1.3–5.46, p = 0.0055, not shown). Similarly, stage III and IV ovarian cancers suffer worse survival when the CTCFL mRNA level of their tumours is high compared with patients with low levels (HR = 1.26, 95% CI = 1–1.59, p = 0.05). These data suggest that high CTCFL mRNA levels are commonly associated with adverse prognosis in various tumours, although there are exceptions such as ER-positive breast cancer.
The most common mechanism causing re-expression of CTCFL in cancer is promoter hypomethylation. Data from the online platform MethHC comparing the methylation status of CTCFL promoters in various cancers with the status of these promoters in corresponding normal tissues disclose that compared with normal corresponding tissues, several common cancers such as carcinomas of the bladder, clear cell kidney, and squamous carcinomas of the lung, head, and neck display hypomethylation of CTCFL promoter (Figure 8a–d). However other common cancers such as colon and lung adenocarcinomas show no difference in their CTCFL promoter methylation status compared to normal colon and lung tissues (Figure 8f,g), while others, such as breast cancer and melanoma, even have promoter hypermethylation (Figure 8e,h).

4. Discussion

The CTCF transcription regulator is a main organizer of the human genome functioning as an insulator defining TADs borders. These act as physical barriers preventing function of remote enhancers from acting on genes outside the limits of the specific TAD. The insulating function of CTCF takes place through binding of the protein to specific DNA sequences that are ubiquitous throughout the human genome and recruitment of additional partner proteins interacting with the aminoterminal or the carboxyterminal domain of CTCF. Genetic lesions affecting either DNA binding of CTCF or interaction with partners could have severe implications for the function of CTCF as an insulator and lead to profound changes in the regulation of multiple genes through alterations in enhancer regulation, preventing enhancers from acting on normal target genes or creating new influences. These effects could be widespread throughout the genome. Upregulation of the expression of an oncogene or downregulation of a tumour suppressor under abnormal enhancer influences may promote cancer [43,44]. CTCF hemizygous mice displayed dysregulation of hundreds of cancer-related genes [45]. In this model of quantitative reduction of CTCF protein most affected were CTCF binding sites with weaker affinity for the protein. In another model of CTCF haploinsufficiency, using shRNA, decreased CTCF dose promotes cell survival and affects cell polarity, a hallmark of normal polarised epithelia [46]. CTCF mutations in endometrial carcinomas lead to nonsense-mediated decay of the transcripts or loss of function of the protein with missense mutations.
In this paper, CTCF DNA lesions were explored using published publicly available genomic studies and open platforms such as cBioportal available online. Several conclusions can be drawn from this investigation. First, CTCF lesions are rare across cancers, but mutations are much more common in certain cancers such as endometrial cancers than others. Second, underlying MSI-associated or polymerase mutations are common concomitant defects in CTCF mutant cases. This observation agrees with previous publications [47]. These authors have also observed that, similarly to the extensive published series included in the current report, MSI, although common, is not always present in CTCF-mutant cases. In some cases, CTCF mutations are associated with concomitant POLE or POLD1 mutations. The two polymerases are responsible for the synthesis of the leading and lagging strand respectively during DNA replication and mutations in them lead to a hypermutator phenotype [48]. Mutations in POLE or POLD1 lead to a polyposis syndrome called polymerase proofreading-associated polyposis (PPAP) [49]. Some cancers with lower MSI incidence may have other underlying defects that promote CTCF mutations, such as APOBEC deaminases abnormalities. Lastly, despite the fact that a significant proportion of CTCF mutations are considered oncogenic, the association of mRNA dose with tumour aggressiveness and prognosis is variable, suggesting that the protein is not a tumour suppressor in all contexts. This could be expected given the extensive role of CTCF in tertiary DNA organisation which leads to multiple gene dysregulations when defective. In addition, dosage of mRNA does not capture mutations which may have deleterious effects in protein function of an otherwise well-expressed protein.
An additional function of CTCF consists of its involvement in double strand DNA repair [50,51]. Although the mechanism is not entirely clear, and whether involvement of PARP1, BRCA2, and RAD51 as repair partners required in this function is debated, CTCF appears to promote homologous recombination over nonhomologous end-joining (NHEJ) as the mechanism of double strand repair, thus favouring error-free DNA repair. As a result, point mutations in the protein or haplo-insufficiency due to nonsense mutations may have deleterious influence in double strand DNA repair which would be forced to proceed through the error-prone NHEJ mechanism. This could pose an additional burden for cancer cells with other repair defects such as MSI or promote errors creation even in cells that are microsatellite stable.
CTCF function may be more commonly affected through mutations in its ubiquitous DNA binding sites instead of mutations or other DNA alterations affecting the locus of the protein itself. Mutations in specific binding sites of CTCF may have effects in specific TADs but would be expected to have much less widespread influence on genome regulations than DNA lesions affecting the CTCF protein. On the other hand, such binding site mutations may have very specific oncogenic effects which may, for example, lead to expression of an oncogene under the influence of a new enhancer after TAD reshuffling. In agreement with this discussion, CTCF binding sites have been reported to be mutated at high frequency (25% and 19%, respectively) in gastric and colorectal cancers [52]. Moreover, in these cases, CTCF binding site mutations are commonly seen concomitantly with MSI.
CTCF paralogue CTCFL is normally not expressed in adult tissues, besides specific stages of spermatogenesis, due to promoter methylation of its gene, but is re-expressed in some cancer cases. In addition to epigenetic promoter hypomethylation that could promote CTCFL re-expression in cancer, genetic lesions may contribute to CTCFL de-repression. Abnormal re-expression of CTCFL in these cancer cases could have functional implications by interfering with binding of CTCF in a subset of its sites or by binding to methylated sites where CTCF may be less apt or cannot bind. Amplifications of CTCFL is the most common genetic lesion overall and could lead to overexpression of the protein. Data from the published studies presented in the current paper show that amplification is not always associated with mRNA upregulation (Figure 5) and thus the implications of such amplifications remain unclear. In addition, no other gene in the 20q13 amplicon appears to be more commonly amplified, a fact that would suggest cancer cell survival benefit leading to clonal dominance. Breast cancers with the common ERBB2 amplification defining the HER2-positive subset have only partially overlapping amplifications with CTCFL suggesting that the two amplicons may be created by different underlying mechanisms and not by a common mechanism affecting the two chromosomes, 17 and 20.
CTCFL promoter hypomethylation was evident in some types of cancers compared with corresponding normal tissues but not in other cancers (Figure 8). Whether there is a correlation of CTCFL promoter hypomethylation with protein expression in vivo in cancer patients remains untested. Squamous cell carcinomas of lung and head and neck which show CTCFL promoter hypomethylation (Figure 8b,d) have been observed to express CTCFL transcripts but no evidence exist for the corresponding protein expression [53]. Conversely, breast cancer displays hypermethylation in the CTCFL promoter (Figure 8e) and data suggest that CTCFL protein may not be expressed in human breast cancers, although this is controversial [54,55].
Some cancers present more commonly with CTCFL mutations rather than amplifications (Figure 4), and these could be of significance if the defective protein is expressed. In most occasions this is not the case and re-expression of CTCFL remains a rare occurrence in cancers. Thus, CTCFL mutations seem to be in many cases a marker of underlying genetic instability without pathologic implications of the specific mutation per se. Genetic instability such as MSI, POLE, and POLD1 defects lead to an increased mutation burden which arises as a leading biomarker of response to immune blockade inhibitors, new drugs that have improved outcomes in several cancer types through immune stimulation. Besides responding better to these novel drugs, cancers with increased mutation burden tend also to have an improved prognosis. Given the significant effects that TAD border reshuffling may have for gene expression regulation, it would be of interest to further investigate the effect of CTCF and CTCFL lesions or their binding sites in TAD borders for refining the prognostic implications of mutation burden as a prognostic and immune checkpoint inhibitors predictive biomarker.

Supplementary Materials

The following are available online at https://www.mdpi.com/2571-5135/7/4/30/s1, Figure S1: The CTCFL amplicon in the METABRIC breast cancer study.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Kaiser, V.B.; Semple, C.A. When TADs go bad: Chromatin structure and nuclear organisation in human disease. F1000 Res. 2017, 6, 314. [Google Scholar] [CrossRef] [PubMed]
  2. Hnisz, D.; Day, D.S.; Young, R.A. Insulated neighborhoods: Structural and functional units of mammalian gene control. Cell 2016, 167, 1188–1200. [Google Scholar] [CrossRef] [PubMed]
  3. Dixon, J.R.; Gorkin, D.U.; Ren, B. Chromatin domains: The unit of chromosome organisation. Mol. Cell 2016, 62, 668–680. [Google Scholar] [CrossRef] [PubMed]
  4. Eukaryotic Promoter Database. Available online: https://epd.vital-it.ch (accessed on 7 July 2018).
  5. Wang, D.C.; Wang, W.; Zhang, L.; Wang, X. A tour of 3D genome with a focus on CTCF. Semin Cell Dev. Biol. 2018. [Google Scholar] [CrossRef] [PubMed]
  6. Chen, H.; Tian, Y.; Shu, W.; Bo, X.; Wang, S. Comprehensive identification and annotation of cell type-specific and ubiquitous CTCF-binding sites in the human genome. PLoS ONE 2012, 7, e41374. [Google Scholar] [CrossRef] [PubMed]
  7. Barski, A.; Cuddapah, S.; Cui, K.; Roh, T.Y.; Schones, D.E.; Wang, Z.; Wei, G.; Chepelev, I.; Zhao, K. High-resolution profiling of histone methylations in the human genome. Cell 2007, 129, 823–837. [Google Scholar] [CrossRef] [PubMed]
  8. Kim, T.H.; Abdullaev, Z.K.; Smith, A.D.; Ching, K.A.; Loukinov, D.I.; Green, R.D.; Zhang, M.Q.; Lobanenkov, V.V.; Ren, B. Analysis of the vertebrate insulator protein CTCF-binding sites in the human genome. Cell 2007, 128, 1231–1245. [Google Scholar] [CrossRef] [PubMed]
  9. Holwerda, S.J.B.; de Laat, W. CTCF: The protein, the binding partners, the binding sites and their chromatin loops. Phil. Trans. R. Soc. B 2013, 368, 20120369. [Google Scholar] [CrossRef] [PubMed]
  10. Martin-Kleiner, I. BORIS in human cancers—A review. Eur. J. Cancer 2012, 48, 929–935. [Google Scholar] [CrossRef] [PubMed]
  11. Marshall, A.D.; Bailey, C.G.; Rasko, J.E.J. CTCF and BORIS in genome regulation and cancer. Curr. Opin. Genet. Dev. 2014, 24, 8–15. [Google Scholar] [CrossRef] [PubMed]
  12. Nguyen, P.; Cui, H.; Bisht, K.S.; Sun, L.; Patel, K.; Lee, R.S.; Kugoh, H.; Oshimura, M.; Feinberg, A.P.; Gius, D. CTCFL/BORIS is a methylation-independent DNA-binding protein that preferentially binds to the paternal H19 differentially methylated region. Cancer Res. 2008, 68, 5546–5551. [Google Scholar] [CrossRef] [PubMed]
  13. Maurano, M.T.; Wang, H.; John, S.; Shafer, A.; Canfield, T.; Lee, K.; Stamatoyannopoulos, J.A. Role of DNA methylation in modulating transcription factor occupancy. Cell Rep. 2015, 12, 1184–1195. [Google Scholar] [CrossRef] [PubMed]
  14. Van Tongelen, A.; Loriot, A.; De Smet, C. Oncogenic roles of DNA hypomethylation through the activation of cancer-germline genes. Cancer Lett. 2017, 396, 130–137. [Google Scholar] [CrossRef] [PubMed]
  15. Jabbari, K.; Heger, P.; Sharma, R.; Wiehe, T. The diverging routes of BORIS and CTCF: An interactomic and phylogenomic analysis. Life 2018, 8, 4. [Google Scholar] [CrossRef] [PubMed]
  16. Pugacheva, E.M.; Rivero-Hinojosa, S.; Espinoza, C.A.; Méndez-Catalá, F.; Kang, S.; Suzuki, T.; Kosaka-Suzuki, N.; Robinson, S.; Nagarajan, V.; Ye, Z.; et al. Comparative analyses of CTCF and BORIS occupancies uncover two distinct classes of CTCF binding genomic regions. Genome. Biol. 2015, 16, 161. [Google Scholar] [CrossRef] [PubMed]
  17. Bioportal for Cancer Genomics. Available online: http://www.cbioportal.org (accessed on 7 July 2018).
  18. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio Cancer Genomics Portal: An open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci. Signal. 2013, 6, 269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Robertson, A.G.; Kim, J.; Al-Ahmadie, H.; Bellmunt, J.; Guo, G.; Cherniack, A.D.; Hinoue, T.; Laird, P.W.; Hoadley, K.A.; Akbani, R.; et al. Comprehensive molecular characterization of muscle-invasive bladder cancer. Cell 2017, 171, 540–556. [Google Scholar] [CrossRef] [PubMed]
  21. Cancer Genome Atlas Network. Comprehensive molecular characterization of human colon and rectal cancer. Nature 2012, 487, 330–337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Cancer Genome Atlas Research Network. Comprehensive molecular characterization of gastric adenocarcinoma. Nature 2014, 513, 202–209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Cancer Genome Atlas Network; Kandoth, C.; Schultz, N.; Cherniack, A.D.; Akbani, R.; Liu, Y.; Shen, H.; Robertson, A.G.; Pashtan, I.; Shen, R.; et al. Integrated genomic characterization of endometrial carcinoma. Nature 2013, 497, 67–73. [Google Scholar] [PubMed] [Green Version]
  24. Curtis, C.; Shah, S.P.; Chin, S.-F.; Turashvili, G.; Rueda, O.M.; Dunning, M.J.; Speed, D.; Lynch, A.G.; Samarajiwa, S.; Yuan, Y.; et al. The genomic and transcriptomic architecture of 2000 breast tumours reveals novel subgroups. Nature 2012, 486, 346–352. [Google Scholar] [CrossRef] [PubMed]
  25. Barretina, J.; Caponigro, G.; Stransky, N.; Venkatesan, K.; Margolin, A.A.; Kim, S.; Wilson, C.J.; Lehár, J.; Kryukov, G.V.; Sonkin, D.; et al. The Cancer Cell Line Encyclopedia enables predictive modeling of anticancer drug sensitivity. Nature 2012, 483, 603–607. [Google Scholar] [CrossRef] [PubMed]
  26. Ciriello, G.; Gatza, M.L.; Beck, A.H.; Wilkerson, M.D.; Rhie, S.K.; Pastore, A.; Zhang, H.; McLellan, M.; Yau, C.; Kandoth, C.; et al. Comprehensive molecular portraits of invasive lobular breast cancer. Cell 2015, 163, 506–519. [Google Scholar] [CrossRef] [PubMed]
  27. Brennan, C.W.; Verhaak, R.G.; McKenna, A.; Campos, B.; Noushmehr, H.; Salama, S.R.; Zheng, S.; Chakravarty, D.; Sanborn, J.Z.; Berman, S.H.; et al. The somatic genomic landscape of glioblastoma. Cell 2013, 155, 462–477. [Google Scholar] [CrossRef] [PubMed]
  28. Cancer Genome Atlas Network. Comprehensive molecular profiling of lung adenocarcinoma. Nature 2014, 511, 543–550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Cancer Genome Atlas Network. Comprehensive genomic characterization of squamous cell lung cancers. Nature 2012, 489, 519–525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Cancer Genome Atlas Research Network. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature 2013, 499, 43–49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Cancer Genome Atlas Research Network. Integrated genomic analyses of ovarian carcinoma. Nature 2011, 474, 609–615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Cancer Genome Atlas Research Network. The molecular taxonomy of primary prostate cancer. Cell 2015, 163, 1011–1025. [Google Scholar] [CrossRef] [PubMed]
  33. Mutation Assessor. Available online: www.mutationassessor.org (accessed on 29 July 2018).
  34. Reva, B.; Antipin, Y.; Sander, C. Predicting the functional impact of protein mutations: Application to cancer genomics. Nucleic Acids Res. 2011, 39, e118. [Google Scholar] [CrossRef] [PubMed]
  35. Kaplan Meier Plotter. Available online: www.kmplot.com (accessed on 14 July 2018).
  36. Szász, A.M.; Lánczky, A.; Nagy, Á.; Föster, S.; Hark, K.; Green, J.E.; Boussioutas, A.; Busuttil, R.; Szabó, A.; Győrffy, B. Cross-validation of survival associated biomarkers in gastric cancer using transcriptomic data of 1065 patients. Oncotarget 2016, 7, 49322–49333. [Google Scholar] [CrossRef] [PubMed]
  37. MethHC A database of DNA Methylation and Gene Expression in Human Cancer. Available online: www.methhc.mbc.nctu.edu.tw (accessed on 3 July 2018).
  38. Huang, W.Y.; Hsu, S.D.; Huang, H.Y.; Sun, Y.M.; Chou, C.H.; Weng, S.L.; Huang, H.D. MethHC: A database of DNA methylation and gene expression in human cancer. Nucleic Acids Res. 2015, 43, D856–D861. [Google Scholar] [CrossRef] [PubMed]
  39. Konstantinopoulos, P.A.; Matulonis, U.A. POLE mutations as an alternative pathway for Microsatellite Instability in endometrial cancer: Implications for Lynch syndrome testing. Cancer 2015, 121, 331–334. [Google Scholar] [CrossRef] [PubMed]
  40. Rebhandl, S.; Huemer, M.; Greil, R.; Geisberger, R. AID/APOBEC deaminases and cancer. Oncoscience 2015, 2, 320–333. [Google Scholar] [CrossRef] [PubMed]
  41. Venkatesan, S.; Rosenthal, R.; Kanu, N.; McGranahan, N.; Bartek, J.; Quezada, S.A.; Hare, J.; Harris, R.S.; Swanton, C. Perspective: APOBEC mutagenesis in drug resistance and immune escape in HIV and cancer evolution. Ann. Oncol. 2018, 29, 563–572. [Google Scholar] [CrossRef] [PubMed]
  42. Human Protein Atlas. Available online: www.proteinatlas.org (accessed on 10 July 2018).
  43. Lai, A.Y.; Fatemi, M.; Dhasarathy, A.; Malone, C.; Sobol, S.E.; Geigerman, C.; Jaye, D.L.; May, D.; Shah, R.; Li, L.; et al. DNA methylation prevents CTCF-mediated silencing of the oncogene BCL6 in B cell lymphomas. J. Exp. Med. 2010, 207, 1939–1950. [Google Scholar] [CrossRef] [PubMed]
  44. Witcher, M.; Emerson, B.M. Epigenetic silencing of the p16(INK4a) tumour suppressor is associated with loss of CTCF binding and a chromatin boundary. Mol. Cell 2009, 34, 271–284. [Google Scholar] [CrossRef] [PubMed]
  45. Aitken, S.J.; Ibarra-Soria, X.; Kentepozidou, E.; Flicek, P.; Feig, C.; Marioni, J.C.; Odom, D.T. CTCF maintains regulatory homeostasis of cancer pathways. Genome. Biol. 2018, 19, 106. [Google Scholar] [CrossRef] [PubMed]
  46. Marshall, A.D.; Bailey, C.G.; Champ, K.; Vellozzi, M.; O’Young, P.; Metierre, C.; Feng, Y.; Thoeng, A.; Richards, M.; Schmitz, U.; et al. CTCF genetic alterations in endometrial carcinoma are pro-tumourigenic. Oncogene 2017, 36, 4100–4110. [Google Scholar] [CrossRef] [PubMed]
  47. Zighelboim, I.; Mutch, D.G.; Knapp, A.; Ding, L.; Xie, M.; Cohn, D.E.; Goodfellow, P.J. High frequency strand slippage mutations in CTCF in MSI-positive endometrial cancers. Human Mutat. 2014, 35, 63–65. [Google Scholar] [CrossRef] [PubMed]
  48. Voutsadakis, I.A. Polymerase epsilon mutations and concomitant β2-microglobulin mutations in cancer. Gene 2018, 647, 31–38. [Google Scholar] [CrossRef] [PubMed]
  49. Palles, C.; Cazier, J.B.; Howarth, K.M.; Domingo, E.; Jones, A.M.; Broderick, P.; Kemp, Z.; Spain, S.L.; Guarino, E.; Salguero, I.; et al. Germline mutations affecting the proofreading domains of POLE and POLD1 predispose to colorectal adenomas and carcinomas. Nature Genet. 2012, 45, 136–144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Lang, F.; Li, X.; Zheng, W.; Li, Z.; Lu, D.; Chen, G.; Gong, D.; Yang, L.; Fu, J.; Shi, P.; et al. CTCF prevents genomic instability by promoting homologous recombination-directed DNA double-strand break repair. Proc. Natl. Acad. Sci. USA 2017, 114, 10912–10917. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Hilmi, K.; Jangal, M.; Marques, M.; Zhao, T.; Saad, A.; Zhang, C.; Luo, V.M.; Syme, A.; Rejon, C.; Yu, Z.; et al. CTCF facilitates DNA double-strand break repair by enhancing homologous recombination repair. Sci. Adv. 2017, 3, e1601898. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Guo, Y.A.; Chang, M.M.; Huang, W.; Ooi, W.F.; Xing, M.; Tan, P.; Skanderup, A.J. Mutation hotspots at CTCF binding sites coupled to chromosomal instability in gastrointestinal cancers. Nat. Commun. 2018, 9, 1520. [Google Scholar] [CrossRef] [PubMed]
  53. Smith, I.M.; Glazer, C.A.; Mithani, S.K.; Ochs, M.F.; Sun, W.; Bhan, S.; Vostrov, A.; Abdullaev, Z.; Lobanenkov, V.; Gray, A.; et al. Coordinated activation of candidate proto-oncogenes and cancer testes antigens via promoter demethylation in head and neck cancer and lung cancer. PLoS ONE 2009, 4, e4961. [Google Scholar] [CrossRef] [PubMed]
  54. Hines, W.C.; Bazarov, A.V.; Mukhopadhyay, R.; Yaswen, R. BORIS (CTCFL) is not expressed in most human breast cell lines and high grade breast carcinomas. PLoS ONE 2010, 5, e9738. [Google Scholar] [CrossRef] [PubMed]
  55. D’Arcy, V.; Pore, N.; Docquier, F.; Abdullaev, Z.K.; Chernukhin, I.; Kita, G.X.; Rai, S.; Smart, M.; Farrar, D.; Pack, S.; et al. BORIS, a paralogue of the transcription factor, CTCF, is aberrantly expressed in breast tumours. Br. J. Cancer 2008, 98, 571–579. [Google Scholar] [CrossRef] [PubMed] [Green Version]
High throughput 07 00030 g001 550
Figure 1. Schematic representation of CTCF protein with domain organisation and examples of the most frequent amino acid locations of recurrent mutations. Numbers below the protein schema represent margins of the protein domains. Domains are as follows: amino acids 1–266: amino-terminal domain, amino acids 267–577: Zinc finger domain, amino acids 578–727: carboxy-terminal domain. The amino acid locations of recurrent mutations are represented above the schematic. N: aminoterminal, C: carboxy-terminal.
Figure 1. Schematic representation of CTCF protein with domain organisation and examples of the most frequent amino acid locations of recurrent mutations. Numbers below the protein schema represent margins of the protein domains. Domains are as follows: amino acids 1–266: amino-terminal domain, amino acids 267–577: Zinc finger domain, amino acids 578–727: carboxy-terminal domain. The amino acid locations of recurrent mutations are represented above the schematic. N: aminoterminal, C: carboxy-terminal.
High throughput 07 00030 g001
High throughput 07 00030 g002 550
Figure 2. Percentage of CTCF mutations and amplifications combined (a) and mutations (b) in 12 common cancers. The total number of samples examined was 6043.
Figure 2. Percentage of CTCF mutations and amplifications combined (a) and mutations (b) in 12 common cancers. The total number of samples examined was 6043.
High throughput 07 00030 g002
High throughput 07 00030 g003a 550High throughput 07 00030 g003b 550
Figure 3. Overall survival of patients with high CTCF mRNA expression compared with patients with low CTCF mRNA expression. (a) Gastric cancer. (b) Breast cancer across all subtypes. (c) HER2+ Breast cancer. Low (black lines) and high (red lines) curves represent cases with mRNA levels below and above the mean of the groups.
Figure 3. Overall survival of patients with high CTCF mRNA expression compared with patients with low CTCF mRNA expression. (a) Gastric cancer. (b) Breast cancer across all subtypes. (c) HER2+ Breast cancer. Low (black lines) and high (red lines) curves represent cases with mRNA levels below and above the mean of the groups.
High throughput 07 00030 g003aHigh throughput 07 00030 g003b
High throughput 07 00030 g004 550
Figure 4. Percentage of CTCFL mutations and amplifications combined (a) and amplifications (b) in 12 common cancers. Please note that studies and types of cancer presented in parts (a) and (b) of the figure are only partially overlapping as the two parts include studies and cancer types with the higher prevalence of total CTCFL lesions and CTCFL amplifications which are different in each occasion.
Figure 4. Percentage of CTCFL mutations and amplifications combined (a) and amplifications (b) in 12 common cancers. Please note that studies and types of cancer presented in parts (a) and (b) of the figure are only partially overlapping as the two parts include studies and cancer types with the higher prevalence of total CTCFL lesions and CTCFL amplifications which are different in each occasion.
High throughput 07 00030 g004
High throughput 07 00030 g005a 550High throughput 07 00030 g005b 550High throughput 07 00030 g005c 550
Figure 5. CTCFL gene amplification versus mRNA expression in various cancers. (a) Breast cancer. (b) Ovarian serous Cystadenocarcinoma. (c) Colon adenocarcinoma. (d) Gene amplification versus mRNA expression of ERBB2 in breast cancer. In (a,b) the mean mRNA levels in diploid and amplified cases are similar at 4 and 7, respectively. In (c,d) amplified cases present a higher mean mRNA level than diploid cases. GISTIC: Genomic Identification of Significant Targets in Cancer. The copy number analysis algorithm according to GISTIC defines a copy number below −2 as deep deletion (possible homozygous deletion), copy number between −2 and −1 as swallow deletion (possible heterozygous deletion), copy number between −1 and 1 as diploid, copy number between 1 and 2 as low-level gain, and copy number above 2 as amplification.
Figure 5. CTCFL gene amplification versus mRNA expression in various cancers. (a) Breast cancer. (b) Ovarian serous Cystadenocarcinoma. (c) Colon adenocarcinoma. (d) Gene amplification versus mRNA expression of ERBB2 in breast cancer. In (a,b) the mean mRNA levels in diploid and amplified cases are similar at 4 and 7, respectively. In (c,d) amplified cases present a higher mean mRNA level than diploid cases. GISTIC: Genomic Identification of Significant Targets in Cancer. The copy number analysis algorithm according to GISTIC defines a copy number below −2 as deep deletion (possible homozygous deletion), copy number between −2 and −1 as swallow deletion (possible heterozygous deletion), copy number between −1 and 1 as diploid, copy number between 1 and 2 as low-level gain, and copy number above 2 as amplification.
High throughput 07 00030 g005aHigh throughput 07 00030 g005bHigh throughput 07 00030 g005c
High throughput 07 00030 g006 550
Figure 6. Amplifications of CTCFL, ZNF217, and ERBB2 in the METABRIC breast cancer study samples. Percentages represent the percentage of cases of each gene amplified in the study.
Figure 6. Amplifications of CTCFL, ZNF217, and ERBB2 in the METABRIC breast cancer study samples. Percentages represent the percentage of cases of each gene amplified in the study.
High throughput 07 00030 g006
High throughput 07 00030 g007a 550High throughput 07 00030 g007b 550
Figure 7. Overall survival of patients with high CTCFL mRNA expression compared with patients with low CTCFL mRNA expression. (a) Gastric cancer. (b) Breast cancer across all subtypes. (c) HER2+ Breast cancer. Low (black lines) and high (red lines) curves represent cases with mRNA levels below and above the mean of the groups.
Figure 7. Overall survival of patients with high CTCFL mRNA expression compared with patients with low CTCFL mRNA expression. (a) Gastric cancer. (b) Breast cancer across all subtypes. (c) HER2+ Breast cancer. Low (black lines) and high (red lines) curves represent cases with mRNA levels below and above the mean of the groups.
High throughput 07 00030 g007aHigh throughput 07 00030 g007b
High throughput 07 00030 g008 550
Figure 8. Methylation status of the CTCFL promoter in various cancers compared with corresponding normal tissues. Red symbols represent tumour samples and green normal samples. Numbers on the left of each comparison represent average beta values. ** p < 0.005. (a) Bladder carcinoma (BLCA). (b) Squamous carcinoma of the lung (LUSQ). (c) Clear cell Renal cell carcinoma (KIRC). (d) Head and neck cancer (HNSC). (e) Breast cancer (BRCA). (f) Colon cancer (COAD). (g) Adenocarcinoma of the lung (LUAD). (h) Melanoma (SKCM).
Figure 8. Methylation status of the CTCFL promoter in various cancers compared with corresponding normal tissues. Red symbols represent tumour samples and green normal samples. Numbers on the left of each comparison represent average beta values. ** p < 0.005. (a) Bladder carcinoma (BLCA). (b) Squamous carcinoma of the lung (LUSQ). (c) Clear cell Renal cell carcinoma (KIRC). (d) Head and neck cancer (HNSC). (e) Breast cancer (BRCA). (f) Colon cancer (COAD). (g) Adenocarcinoma of the lung (LUAD). (h) Melanoma (SKCM).
High throughput 07 00030 g008
Table
Table 1. Summary of studies with CTCF lesions in various cancer studies. In the second column the number of samples with lesions is in the numerator and the total number of samples in the series is in the denominator. GBM: Glioblastoma multiforme, RCC: Renal cell carcinoma, HCC: Hepatocellular carcinoma, TGCTs: Testicular germ cell tumours, NA: Not available.
Table 1. Summary of studies with CTCF lesions in various cancer studies. In the second column the number of samples with lesions is in the numerator and the total number of samples in the series is in the denominator. GBM: Glioblastoma multiforme, RCC: Renal cell carcinoma, HCC: Hepatocellular carcinoma, TGCTs: Testicular germ cell tumours, NA: Not available.
Type of Cancer (Reference)All LesionsAmplificationsDeep DeletionsmRNA UpregulationmRNA DownregulationMutationsMultiple Lesions
Cell line encyclopedia [23]57/877 (6.5%)13 (1.48%)6 (0.68%)24 (2.74%) 14 (1.6%) --
Uterine endometrial TCGA [23]38/102 (37.25%)-1 (0.98%)2 (1.96%)3 (2.94%)28 (27.45%)4 (3.92%)
Ovarian serous TCGA [31]33/200 (16.5%)-2 (1%)3 (1.5%)23 (11.5%)1 (0.5%)4 (2%), two fusions
Bladder TCGA [20]55/404 (13.61%)3 (0.74%)1 (0.25%)25 (6.19%)8 (1.98%)13 (3.22%)5 (1.24%)
Colorectal TCGA [21]30/267 (11.24%)--15 (5.62%)2 (0.75%)13 (4.87%)-
Prostate TCGA [32]53/491 (10.79%)-14 (2.85%)12 (2.44%)22 (4.48%)3 (0.61%)2 (0.41%)
Melanoma TCGA PanCancer (Provis.)36/363 (9.92%)--14 (3.86)11 (3.03%)7 (1.93%)4 (1.1%)
Lung adenocarcinoma TCGA [28]40/503 (7.95%)-1 (0.2%)21 (4.17%)12 (2.39%)5 (0.99%)1 (0.2%)
Lung squamous TCGA [29]34/466 (7.3)--19 (4.08%)9 (1.93%)4 (0.86%)2 (0.42%), one fusion
Pancreatic TCGA PanCancer (Provis.)12/168 (7.14%)--6 (3.57%)5 (2.98%)1 (0.6%)-
TGCTs TCGA PanCancer (Provis.)10/144 (6.94%)--6 (4.17%)3 (2.08%)-1 (0.69%)
HCC TCGA PanCancer (Provis.)20/345 (5.8%)1 (0.29%)-6 (1.74%)9 (2.61%)4 (1.16%)-
RCC TCGA [30]20/352 (5.68%)--6 (1.7%)13 (3.69%)1 (0.28%)-
Uterine serous TCGA [23]3/53 (5.66%)--2 (3.77%)1 (1.89%)--
GBM TCGA [27]5/142 (3.52%)--3 (2.11%)2 (1.41%)--
Breast TCGA [26]26/816 (3.19%)3 (0.37%)5 (0.61%)--18 (2.21%)-
Gastric adenocarcinoma TCGA [22]4/188 (2.13%)1 (0.53%)1 (0.53%)NANA2 (1.06%)-
Oesophageal adenocarcinoma TCGA [22]1/77 (1.3%)1 (1.3%)-NANA--
All (Not cell lines)420/5081 (8.27%)9/5081 (0.18%)25/5081 (0.49%)140/4816 (2.9%)123/4816 (2.55%)100/5081 (1.97%)22/5081 (0.43%)
Table
Table 2. Endometrial cancer CTCF mutations in the TCGA PanCa study. All CTCF mutated samples were diploid. ZF: Zinc finger, N: Aminoterminal, C: Carboxyterminal, FS: Frameshift. No: No mutation, *: Stop codon.
Table 2. Endometrial cancer CTCF mutations in the TCGA PanCa study. All CTCF mutated samples were diploid. ZF: Zinc finger, N: Aminoterminal, C: Carboxyterminal, FS: Frameshift. No: No mutation, *: Stop codon.
Number of SampleProtein ChangeDomainMutation TypeAllele FrequencyNumber of Mutations in SampleMutations in MSI-Associated GenesMutations in POLE/POLD1 Genes
1G48Vfs*14NFS del0.37774NoNo
2G48Vfs*14NFS del0.28716MSH6No
3A88VNMissense0.175737MSH2, MSH6, MLH1, PMS2POLE, POLD1
4G173 *NNonsense0.257390MSH2, MSH6, PMS2POLE
5Q180 *NNonsense0.4594MSH6No
6E182Gfs*9NFS ins0.339662MSH2, PMS2POLE
7T204Qfs*18NFS del0.19597MLH1No
8T204Nfs*26NFS ins0.26538MSH6, PMS2No
9T204Nfs*26NFS ins0.35323NoNo
10T204Nfs*26NFS ins0.5864PMS2POLE
11T204Nfs*26NFS ins0.31451NoNo
12T204Nfs*26NFS ins0.2556NoNo
13R213CNMissense0.3312218MSH6, MLH1, PMS2POLE, POLD1
14D222Mfs*28NFS del0.32218NoNo
15G261CNMissense0.29451NoNo
16Q267PZFMissense0.3213840MSH2, PMS2, PMS1POLE, POLD1
17H312NZFMissense0.393190MSH6, PMS2POLE
18T317Rfs*91ZFFS del0.5669PMS2No
19G318Qfs*16ZFFS ins0.22562NoNo
20C324 *ZFNonsense0.32611NoPOLE
21R341HZFMissense0.459662MSH2, PMS2POLE
22R342CZFMissense0.49440MSH2, MSH6, PMS1POLE
23Y343CZFMissense0.347644MSH6, MLH1POLE
24H369RZFMissense0.171326NoPOLE
25R377CZFMissense0.341307MSH6POLE, POLD1
26P378LZFMissense0.213840MSH2, PMS2POLE, POLD1
27L394delZFIF del0.09716MSH6No
28G420CZFMissense0.261307MSH6POLE, POLD1
29E432Gfs*10ZFFS del0.3841NoNo
30A447Vfs*61ZFFS del0.454346PMS1, MLH1POLE, POLD1
31R448 *ZFNonsense0.4681NoNo
32K449TZFMissense0.348511MSH2, MLH1, PMS1POLE
33S450Kfs*2ZFFS ins0.264096MSH2, MSH6POLE
34X453_spliceZFSplice0.2858NoNo
35R457 *ZFNonsense0.3710061MSH6, MSH2, PMS2POLE, POLD1
36Q499 *ZFNonsense0.353925MSH6POLE
37R533HZFMissense0.47611NoPOLE
38R566HZFMissense0.2884MSH2, MSH6POLD1
39E631 *CNonsense0.413387MSH6, MLH1POLE
Notes: ZF: Zinc finger, N: Aminoterminal, C: Carboxyterminal, FS: Frameshift. No: No mutation, *: Stop codon.
Table
Table 3. Colorectal cancer CTCF mutations. All samples were diploid except one which had CTCF DNA gain.
Table 3. Colorectal cancer CTCF mutations. All samples were diploid except one which had CTCF DNA gain.
NumberProtein ChangeDomainMutation TypeAllele Frequency Number of Mutations in SampleMutations MSIPOLE/POLD1 MutationsAPOBEC Mutations
1R29WNMissense0.44170NoNoNo
2E182Gfs*9NFS ins0.291858NoPOLD1APOBEC2
3D194Rfs*36NFS ins0.25814NoNoAPOBEC4
4T204Qfs*18NFS del0.1530MLH1, PMS1NoAPOBEC4
5T204Nfs*26, D194_A201delinsTQTISNFS ins, Missense0.11, 0.171917PMS1POLEAPOBEC1
6K260 *NNonsense0.25181NoNoAPOBEC3B
7R278CZFMissense0.241002MSH2POLE, POLD1APOBEC3C
8R368CZFMissense0.58171NoNoNo
9R377CZFMissense0.412139MSH6, PMS2POLE, POLD1APOBEC4, APOBEC3C
10E464KZFMissense0.281107NoPOLD1APOBEC3A
11G519RZFMissense0.18662NoPOLD1No
12A524TZFMissense0.294195MSH6POLE, POLD1APOBEC3G
13E691Sfs*30CFS del0.331057NoPOLD1APOBEC3C
Notes: ZF: Zinc finger, N: Aminoterminal, C: Carboxyterminal, No: No mutation, *: Stop codon.
Table
Table 4. CTCF mutations in urothelial cancer. Most mutations were in Zinc Fingers (ZF) domain of the protein. Seven samples were diploid and three each had gains and shallow deletions.
Table 4. CTCF mutations in urothelial cancer. Most mutations were in Zinc Fingers (ZF) domain of the protein. Seven samples were diploid and three each had gains and shallow deletions.
NumberProtein ChangeDomainMutation TypeNumber of Mutations in SampleMutations in MSIMutations in POLE/POLD1
1E104 *, E145Q, K264NNNonsense580NoNo
2Q117 *NNonsense58NoNo
3D290NZFMissense291NoNo
4H322YZFMissense588NoNo
5T346NZFMissense881NoNo
6F351LZFMissense508NoNo
7S354FZFMissense182NoNo
8S354YZFMissense124NoNo
9G375AZFMissense857MSH2No
10E376 *ZFNonsense133NoNo
11S388NZFMissense283NoNo
12E631 *CNonsense3545MSH2, MLH1POLE
13E687 *CNonsense766MSH6No
Notes: N: Aminoterminal, C: Carboxyterminal, *: Stop codon.
Table
Table 5. Summary of CTCFL lesions in various cancer studies. GBM: Glioblastoma multiforme, RCC: Renal cell carcinoma, HCC: Hepatocellular carcinoma, TGCTs: Testicular germ cell tumours, NA: Not available.
Table 5. Summary of CTCFL lesions in various cancer studies. GBM: Glioblastoma multiforme, RCC: Renal cell carcinoma, HCC: Hepatocellular carcinoma, TGCTs: Testicular germ cell tumours, NA: Not available.
Type of Cancer (Reference)All LesionsAmplificationsDeep DeletionsmRNA UpregulationmRNA DownregulationMutationsMultiple Lesions
Cell line encyclopedia [25]100/877 (11.4%)79 (9.01%)-19 (2.17%) --2 (0.23%)
Uterine serous TCGA [23]17/53 (32.08%)4 (7.55%)-13 (24.53%)---
Ovarian serous TCGA [31]37/200 (18.5%)12 (6%)-22 (11%)-1 (0.5%)2 (1%) 1 fusion
Gastric adenocarcinoma TCGA [22]26/188 (13.83%)24 (12.77%)-NANA2 (1.06%)-
Melanoma TCGA PanCancer (Provisional)49/363 (13.5%)4 (1.1%)-17 (4.68%)-27 (7.44%)1 (0.28%)
Colorectal TCGA [21]33/267 (12.36%)9 (3.37%)-13 (4.87%)-5 (1.87%)6 (2.25%)
Uterine endometrial [23]12/102 (11.76%)--3 (2.94%)-7 (6.86%)2 (1.96%)
Lung adenocarcinoma TCGA [28]54/503 (10.74%)19 (3.78%)1 (0.2%)13 (2.58%)-19 (3.78%)2 (0.4%)
Esophageal adenocarcinoma TCGA [22]7/77 (9.09%)7/77 (9.09%)-NANA--
Breast TCGA [26]62/816 (7.6%)56 (6.86%)2 (0.25%)--3 (0.37%)-
Lung squamous TCGA [29]33/466 (7.08%)1 (0.21%)1 (0.2%)22 (4.72%)-8 (1.72%)1 (0.2%)
Prostate TCGA [32]32/491 (6.52%)4 (0.81%)1 (0.2%)13 (2.65%)12 (2.44%)2 (0.41%)-
TGCTs TCGA PanCancer [Provisional]9/144 (6.25%)--9 (6.25%)---
Pancreatic TCGA PanCancer (Provisional)10/168 (5.95%)1 (0.6%)-6 (3.57%)1 (0.6%)-2 (1.19%)
HCC TCGA PanCancer (Provisional)19/345 (5.51%)5 (1.45%)-12 (3.48%)-2 (0.58%)-
GBM TCGA [27]6/142 (4.23%)--6/142 (4.23%)---
RCC TCGA [30]10/352 (2.84%)--8 (2.27%)1 (0.28%)1 (0.28%)-
Bladder TCGA [20]10/404 (2.48%)4 (0.99%)-3 (0.74%)-3 (0.74%)-
All (Not lines)426/5081 (8.38%)150/5081 (2.95%)5/5081 (0.1%)160/4816 (3.3%)14/4816 (0.29%)80/5081 (1.57%)16/5081 (0.31%)

Share and Cite

MDPI and ACS Style

Voutsadakis, I.A. Molecular Lesions of Insulator CTCF and Its Paralogue CTCFL (BORIS) in Cancer: An Analysis from Published Genomic Studies. High-Throughput 2018, 7, 30. https://doi.org/10.3390/ht7040030

AMA Style

Voutsadakis IA. Molecular Lesions of Insulator CTCF and Its Paralogue CTCFL (BORIS) in Cancer: An Analysis from Published Genomic Studies. High-Throughput. 2018; 7(4):30. https://doi.org/10.3390/ht7040030

Chicago/Turabian Style

Voutsadakis, Ioannis A. 2018. "Molecular Lesions of Insulator CTCF and Its Paralogue CTCFL (BORIS) in Cancer: An Analysis from Published Genomic Studies" High-Throughput 7, no. 4: 30. https://doi.org/10.3390/ht7040030

APA Style

Voutsadakis, I. A. (2018). Molecular Lesions of Insulator CTCF and Its Paralogue CTCFL (BORIS) in Cancer: An Analysis from Published Genomic Studies. High-Throughput, 7(4), 30. https://doi.org/10.3390/ht7040030

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop