Combined deletion and DNA methylation result in silencing of FAM107A gene in laryngeal tumors

Larynx squamous cell carcinoma (LSCC) is characterized by complex genotypes, with numerous abnormalities in various genes. Despite the progress in diagnosis and treatment of this disease, 5-year survival rates remain unsatisfactory. Therefore, the extended studies are conducted, with the aim to find genes, potentially implicated in this cancer. In this study, we focus on the FAM107A (3p14.3) gene, since we found its significantly reduced expression in LSCC by microarray profiling (Affymetrix U133 Plus 2.0 array). By RT-PCR we have confirmed complete FAM107A downregulation in laryngeal cancer cell lines (15/15) and primary tumors (21/21) and this finding was further supported by FAM107A protein immunohistochemistry (15/15). We further demonstrate that a combined two hit mechanism including loss of 3p and hypermethylation of FAM107A promoter region (in 9/15 cell lines (p < 0.0001) and in 15/21 primary tumors (p < 0.0001)) prevails in the gene transcriptional loss. As a proof of principle, we show that Decitabine - a hypomethylating agent – restores FAM107A expression (5 to 6 fold increase) in the UT-SCC-29 cell line, characterized by high DNA methylation. Therefore, we report the recurrent inactivation of FAM107A in LSCC, what may suggest that the gene is a promising tumor suppressor candidate involved in LSCC development.

about genetic background of laryngeal cancers, but also may help to estimate the risk of disease burden or the response to the applied therapy.
The application of novel high throughput technologies based on microarrays and next generation sequencing has accelerated tumor genetic studies in recent years 10,11 . In this study we have used expression profiles from Affymetrix GeneChip Human Genome U133 Plus 2.0 arrays, performed previously by our group 12 . To identify novel, potential tumor suppressor genes, we screened this data for genes showing reduced or lack of the mRNA expression. This enabled us to indicate the FAM107A gene as a promising tumor suppressor gene candidate involved in larynx cancer development and to identify the main mechanisms of its inactivation in this type of cancer.

Results
The microarray based expression analysis revealed downregulation of FAM107A in laryngeal tumor cell lines and primary samples. Using the expression profiles established previously with the application of Affymetrix GeneChip Human Genome U133 Plus 2.0 array we have searched for genes differentially expressed in all laryngeal cancer samples (12 cell lines and 5 primary tumors; n = 17) in comparison to 3 non-tumor controls. For this purpose, we have screened the microarray data for tags carrying the "absent" call in each laryngeal cancer sample and the "present" call in each non-tumor control. This filtering resulted in 11 out of 54 675 tags that fulfilled these criteria and included 209074_s_at, that corresponded to FAM107A gene localized in 3p14.3 chromosomal region (chr3:58,549,845-58,563,491; UCSC Genome Browser GRCh37/hg19). This gene was selected for further study because in its case, the difference in mean expression level between tumor and non-tumor samples was the highest. Another tag, namely 207547_s_at also corresponded to this gene. Figure 1 presents the relative expression of FAM107A as indicated by both tags in the analyzed samples, showing their statistically significant downregulation in the laryngeal cancer samples (both primary tumors and cell lines) in comparison to non-tumor controls (Mann-Whitney test). Using the 209074_s_at tag we have identified a 48.6 fold (p = 0.004) and 13.5 fold (p = 0.036) decrease in expression in cell lines and primary tumor samples, respectively as compared to non-tumor controls. Likewise, 10 and 8.5 fold decrease of expression was observed for 207547_s_at tag (p = 0.004 and p = 0.036, respectively). Intrigued by this finding, we have further evaluated the status of FAM107A copy number on Agilent Human Genome CGH 244A and 44K Microarrays (aCGH), performed previously on 13 cell lines. Homozygous deletions targeting this gene were excluded earlier 12,13 , however basing on the copy number plots presenting chromosome 3 in 13 cell lines we have observed deletions of the short arm of this chromosome, resulting in loss of one copy of the FAM107A gene (Supplementary Figure S1).

RT-PCR confirms the complete loss of FAM107A expression in laryngeal cancer samples.
To confirm the downregulation of FAM107A indicated by the expression microarray we have amplified the whole coding region of this gene by RT-PCR in the presence of GAPDH gene as the internal control. Because the expression microarray was performed on both, cell lines and primary tumors we found it to be reasonable to include in this experiment both types of samples. The results, confirmed the complete lack of FAM107A expression in laryngeal cancer cell lines (15/15; 100%; Fig. 2b,c) and primary tumors (21/21; 100%; Fig. 2c,d). On the contrary, two bands, specific for both, FAM107A and GAPDH genes were visible in all applied non-tumor controls (9/9; Fig. 2a). This finding encouraged us to search for the potential mechanism inactivating the remaining copy of FAM107A.
Sequencing of coding exons and intron/exon junctions of FAM107A. To find the potential loss of function mutations in FAM107A we have performed Sanger sequencing in 15 LSCC cell lines. The coding sequences of FAM107A (NM_001076778.2), covering exon 2, 3 and 4, together with 5′ and 3′ splicing sites (up to 8 intronic bp) were analyzed. Four different single nucleotide polymorphisms (SNP) were detected in four different cell lines (Table 1 and Supplementary Figure S2). All identified variants were homozygous, further corroborating the deletion of one copy observed by aCGH. Three of them were missense variants (rs1043942, rs11539086 and rs141609445) -resulted in amino acid change and one was synonymous change (rs1139701). The rs1043942 and rs141609445 missense variants were indicated during in silico analysis with the application of PolyPhen-2 software as probably damaging for the protein structure. Additionally, the analysis of population data for the latter revealed, that in European-American population only the reference allele was observed (http:// browser.1000genomes.org).
While no clear novel inactivating mutations were detected in this analysis, we have discontinued searching for molecular variants in primary laryngeal tumor samples.

DNA methylation analysis in the FAM107A promoter region by bisulfite pyrosequencing.
To study further the potential mechanisms underlying FAM107A downregulation we have analyzed the DNA methylation level of this gene in 15 cell lines, 21 primary tumors and 8 non-tumor controls (Fig. 3). The data obtained for non-tumor controls were used to calculate the cut-off value (41%), above which the samples were assigned as hypermethylated. The dilution series of methylated and unmethylated templates were used to determine the assay sensitivity and PCR bias. It showed lowered selectivity of the primers towards methylated template and thus the cut-off point assigned at 41% of methylation corresponded to the mix of standards containing 70% of methylated template. Therefore, all methylation results above the cut-off point were considered as hypermethylated (Supplementary Figure S3).
The methylation level ranged from 15 to 28% for non-tumor controls, 27% to 87% for cell lines and 14% to 69% for primary tumors. In detail, the pyrosequencing analysis revealed significant differences in DNA methylation level between tumor samples (both laryngeal cancer cell lines and primary tumors) and non-tumor controls (Fig. 3). The hypermethylation was detected in 9/15 cell lines (60%; p < 0.0001) and in 15/21 primary tumors (71%; p < 0.0001).

Demethylation of UT-SCC-29 cell line with Decitabine. The pyrosequencing analysis revealed that
UT-SCC-29 cell line showed the highest DNA methylation level for FAM107A in our study group. Therefore, we decided to demethylate this cell line with the use of Decitabine (DAC; 5-aza-2′-deoxycytidine; Sigma-Aldrich)a known hypomethylating agent. Further, we have analyzed DNA methylation level in the promoter region of FAM107A gene, revealing its significant reduction in samples treated with DAC (Fig. 4a). Differences between mean methylation (MM) values of untreated UT-SCC-29 cell line (MM = 87%) and those treated with 0.1 μM DAC (MM = 31%) and 0.3 μM DAC (MM = 38%) were 56% and 49%, respectively. On the contrary, there was no noticeable difference between mean methylation (MM) values of untreated UT-SCC-29 cell line (MM = 87%) and those treated with 0.1 μM and 0.3 μM acetic acid (MM = 88% and 86%, respectively) -a solvent of DAC, having no influence on methylation or demethylation process (Fig. 4a).
Next, cDNA from the same set of samples was applied for RT-qPCR to analyze FAM107A expression after the demethylation. The obtained data revealed, that in UT-SCC-29 cell line treated with either 0.1 μM DAC or 0.3 μM DAC, FAM107A expression is restored (Fig. 4b). The relative expression values increased 5 to 6 fold in comparison to UT-SCC-29 samples untreated ("mock") or to treated with acetic acid, where the expression is barely observed (Fig. 4b), suggesting an important role of FAM107A epigenetic silencing in LSCC.
Evaluation of FAM107A protein expression by immunohistochemistry. The nuclear/cytoplasmic expression of FAM107A protein was found in glandular cells of fallopian tube (Fig. 5a), which served as the positive control for antibody binding. The FAM107A immunostaining revealed lack of the protein expression in all analyzed primary LSCC samples (  The relationship between DNA methylation of FAM107A and tumor stage. We have assessed the relationship between DNA methylation of FAM107A and: tumor size (T), nodal involvement (N) and differentiation degree (G), however, there was no significant association between any of the analyzed parameters (Table 2). We were not able to analyze the impact of FAM107A hypermethylation on patients' survival due to the lack of relevant data.

Discussion
The solid tumor formation is complex and associated with acquisition of multiple genetic alterations, which influence and deregulate the different cellular processes 14 . Within this, tumor suppressor genes (TSGs) are frequently inactivated at the early stages of disease development. Their loss of anticancer functions is preceded by two events leading to inactivation of each allele, as described by Knudson 15 . 20 years ago Califano et al. have shown that one of the earliest steps of head and neck tumor development, namely the transition from hyperplastic mucosa to dysplasia is related to the loss of heterozygosity (LOH) at 3p chromosomal arm 4 . In line with the observation by Califano et. al, several genes localized in 3p arm, known to be engaged in head and neck cancer appearance, e.g. FHIT or RASSF1A have been identified [16][17][18] .
In this study we have focused on a novel TSG candidate localized in 3p14.3 region, named FAM107A (Homo sapiens family with sequence similarity 107, member A; aliases TU3A or DRR1), potentially involved in laryngeal carcinogenesis. We have shown recurrent transcriptional loss of FAM107A in LSCC using microarray data and RT-PCR (Figs 1 and 2). This is in concordance with the results obtained by other authors, who showed that FAM107A (previously known as TU3A or DRR1) is downregulated in various tumors and cell lines [19][20][21][22] . However, we are the first who show that FAM107A is also downregulated in laryngeal tumors. Moving forward, we have evaluated the DNA copy number of FAM107A using the array CGH data generated previously 12,13 . Although the homozygous deletions including this gene were previously excluded 12 , it is well known, that LOH affecting 3p chromosomal arm is a frequent event in the early step of HNSCC development 4,16,18,23 . The analysis of copy number profiles of chromosome 3 in the analyzed cell lines supports this assumption (Supplementary Figure S1).
Given the above-mentioned arrangements we were interested in finding the potential mechanism underlying FAM107A inactivation in laryngeal tumors, considering that other mechanisms, like point mutations or DNA methylation may affect the remaining gene copy. At first we have investigated the mutational status of this gene in 15 laryngeal cancer cell lines. We found four different sequence variants in four different cell lines -three missense and one synonymous (Table 1 and Supplementary Figure S2). One of them, namely rs141609445 resulting  As the mutation analyses did not reveal a potential mechanism behind FAM107A silencing, in the next step we have analyzed the methylation level of FAM107A promoter region. With the application of bisulfite pyrosequencing we have shown that 60% of cell lines (9/15) and 71% of primary tumors (15/21) are hypermethylated and that the differences in mean methylation value between analyzed samples and non-tumor controls are  statistically significant (p < 0.001; Fig. 3). The hypermethylation was detected among others in cell lines UT-SCC 42B and UT-SCC 116, in which, by Sanger sequencing we have revealed two, probably damaging sequence variants (Table 1). On the other hand, cell lines UT-SCC 19B and UT-SCC 23, harboring synonymous and missense (benign by PolyPhen-2) variants respectively, were not hypermethylated (Table 1 and Fig. 3). Both findings further supports that FAM107A is not inactivated by mutation or SNP. The high prevalence of recurrent FAM107A hypermethylation in the analyzed samples indicates with high degree of certainty that it is the predominant mechanism responsible for the transcriptional silencing of the remaining copy of the gene. Nevertheless, the available data concerning FAM107A methylation in tumors are inconsistent. Awakura et al. showed that more than 40% of renal cell carcinomas presents methylated FAM107A promoter 21 . This gene was also indicated as hypermethylated in hepatocellular carcinoma 26 . On the contrary, in the group of 60 lung cancer patients with decreased FAM107A expression, DNA methylation was observed in minor amount of samples and thus was excluded by the authors as the mechanism of its inactivation 27 . However, in the latter case, the use of the less sensitive MSP method (methylation-specific PCR) to primary tumors might introduce the bias in detecting methylation.
Intrigued by the fact, that the UT-SCC-29 cell line was almost completely methylated in our study (Fig. 3), we took an attempt to demethylate this cell line using well-known hypomethylating agent: 5-aza-2′-deoxycytidine (DAC) and subsequently to verify FAM107A methylation and expression. The used reagent -Decitabine -is currently applied in the epigenetic therapy of various hematologic malignancies. It acts through the whole genome (not selectively) as an inhibitor of DNA methyltransferase 1 (DNMT1), leading to the reactivation of silenced genes 28 . In our study, treatment with DAC resulted in significant decrease of FAM107A methylation, while it was preserved at the high level in control samples, untreated with DAC (Fig. 4a). This shows, that FAM107A methylation is reversible. Moreover, the subsequent RT-qPCR revealed that FAM107A expression was restored in samples treated with DAC (Fig. 4b). This result confirms that DNA methylation is responsible for FAM107A transcriptional inactivation and raises direct question whether DNA methylation is associated with cancer stage. In the study of Awakura et al. the authors found that DNA hypermethylation of TU3A (former name of FAM107A) corresponded with tumor stage (>T2) in primary renal cell carcinoma, however they did not observe such relationship for primary bladder and testicular cancers 21 . Moreover, the tumor differentiation did not correlate with DNA methylation status in any of the analyzed tumors. Our study confirmed lack of relationship between tumor extension, nodal involvement and tumor differentiation with FAM107A hypermethylation (Table 2). But it should be emphasized, that in our study group, the primary samples originated predominantly from patients with advanced tumor stage, thus if FAM107A silencing took place at early stage of carcinogenesis, we were not able to observe such association. Thus, it might be fruitful to assess FAM107A expression or methylation status in premalignant lesions of the larynx.
The shown recurrent inactivation of the gene in LSCC is especially interesting in relation to the potential function of FAM107A. Although it is not completely understood so far, Wang et al. have shown that the protein encoded by the gene contains a nuclear localization signal and a coiled domain, suggesting that it may play a role in regulation of gene transcription and signal transduction 19 . In our study we confirmed the nuclear and cytoplasmic expression of this protein in non-tumor control mucosa (Fig. 5e-h) and importantly showed the lack of FAM107A expression in laryngeal tumor samples (Fig. 5b-d). Functional studies of Wang et al. and Liu et al. revealed that the re-expression of DRR1 gene (former name of FAM107A) in renal and lung cancer cell line respectively, resulted in growth suppression and apoptosis 19,29 . The same effect, i.e. tumor growth inhibition and elevation of apoptosis was observed in the xenograft model after injection of nude mice bearing A549 tumor cells with vector containing DRR1 cDNA 29 .
In summary, our study indicates, that in the course of laryngeal squamous cell carcinoma, FAM107A is inactivated by combined deletions and DNA methylation events. Considering its downregulation in multiple tumors we postulate, that its inactivation is not specific to a given cancer (like laryngeal cancer) but is in contrast, a widespread phenomenon.  Eagle's medium supplemented with 10% of fetal bovine serum at 37 °C under 5% CO 2 atmosphere. For DNA and RNA isolation cell lines were cultured to reach 80% of confluence and subsequently harvested with 0.1% trypsin and 0.2% EDTA. The detailed characteristic of cell lines is presented in Supplementary Table S1.

Primary tumors. The group of 21 primary laryngeal tumor samples derived from patients who underwent surgery in the Department of Otolaryngology and Laryngological Oncology, Poznan University of Medical
Sciences was used in this study. All samples collected during the surgery were cut into three parts: (I) for histopathological assessment and IHC analysis (into the probe filled with 10% of buffered formalin), (II) for DNA isolation (required immediate freezing at −80 °C in an empty tube) and (III) for RNA analysis (stored in RNAlater solution (Sigma-Aldrich, Saint Louis, USA) until the isolation step).
Samples containing more than 60% of tumor cells were included into the study. All experiments were performed in accordance with relevant guidelines and regulations approved by the Ethics Review Board of Poznan University of Medical Sciences (decision no. 164/10 and 502/15) and informed consent was obtained from all donors. The clinical data of patient and tumor characteristics are presented in Supplementary Table S2.
Non-tumor controls. Depending on the performed analysis, various non-tumor samples were applied as controls in this study: (I) for microarray expression analysis three non-tumor controls were applied, i.e. human larynx total RNA (Stratagene, Agilent Technologies, Santa Clara, USA) and RNA derived from bronchial airway epithelia reconstituted in vitro (two donors) (EC, Epithelix Sarl, Geneve, Switzerland) -both commercially available, as well as LX10 -histologically normal mucosa from surgical margin, (II) for RT-PCR analysis nine non-tumor controls were applied, i.e. eight samples obtained during the uvulopalatoplasty procedure (nonmalignant ailments connected with snoring) as well as commercially available total RNA collected from different adult human tissues (qPCR Human Reference Total RNA, Clontech Laboratories, Mountain View, USA), thus providing the broadest coverage of the expressed genes, (III) for bisulfite pyrosequencing analysis eight non-tumor controls collected during the uvulopalatoplasty procedure were applied. Aditionally, the fully methylated standard (MET, Millipore, Hilden, Germany) and unmethylated DNA (UMET), i.e. the whole genome amplified DNA from pooled peripheral blood lymphocytes were used in each run, (IV) for immunohistochemistry (IHC) five samples of Reinke's edema (nonmalignant ailments of the larynx) and one sample of fallopian tube were applied.
DNA and RNA isolation. The nucleic acids from all cell lines, primary tumor samples and non-tumor controls were isolated with the application of standard methods: DNA was obtained using phenol/chloroform extraction and ethanol precipitation and RNA was isolated according to Chomczynski's method with application of Trizol reagent 32 . The whole genome amplified DNA from pooled peripheral blood lymphocytes was prepared with GenomePlex ® Whole Genome Amplification Kit (Sigma-Aldrich), according to manufacturer's procedure.
The concentration and purity of nucleic acids were assessed using NanoDrop-1000 spectrophotometer and RNA integrity was estimated with the use of RNA 6000 Nano Kit on Agilent 2100 BioAnalyzer (Agilent). DNA was stored at −20 °C and RNA at −80 °C until the laboratory use.

Microarray-based gene expression and copy number analysis.
The microarray-based analyses were performed by our group previously and the detailed description of the respective procedures is published elsewhere 12,13 . The expression analysis was performed using Affymetrix GeneChip Human Genome U133 Plus 2.0 array and the genome-wide analysis of DNA copy number changes was performed with the use of Agilent Human Genome CGH 244 A and 44 K Microarray Kits (Agilent Technologies, Waldbronn, Germany).
Reverse transcription and PCR. 1 µg of total RNA from 15 cell lines, 21 primary tumors and 9 non-tumor controls was reverse transcribed to cDNA with the application of Maxima First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, Massachusetts, USA), according to the manufacturer's procedure. In the next step, the level of FAM107A expression was analyzed by multiplex RT-PCR. Two primer pairs were designed: one, covering the whole coding sequence of FAM107A gene and second, covering exons 5-8 of GAPDH gene (Table 3). GAPDH fragment constituted an internal control of RNA integrity and appropriate PCR reaction. In both cases, the primers were designed with the use of Primer3   Fig. 6), localized in the promoter region of FAM107A gene were analyzed for methylation level. The detailed sequencing protocol was described elsewhere 33 . Mean methylation level was calculated from the analyzed CG dinucleotides and was treated as FAM107A methylation in a given sample. Due to the technical impediments we quantified mean methylation from only two CG nucleotide repeats in the primary tumor samples, while all three CG repeats were analyzed in the remaining samples.
To evaluate which samples are hypermethylated the cut off value was determined, based on the results collected for non-tumor controls (three times the standard deviation of methylation for control samples + the highest value of DNA methylation observed in controls). Additionally, dilution series of commercially available methylated DNA template in unmethylated DNA template (WGA -from 0% to 100%, every 10%; Supplementary Figure S3) were used to estimate the extent of PCR bias and the sensitivity of the assay to measure DNA methylation levels in the analyzed samples.  Decitabine -induced DNA demethylation and validation by pyrosequencing and quantitative real time PCR. The cell line UT-SCC-29 was seeded on the 6-well plate and cultured to reach 10-20% confluence under the same conditions as described in section 1.1. The cell line was incubated with 0.1 μM and 0.3 μM concentrations of 5-Aza-2′-deoxycytidine (Decitabine; DAC) in the culture medium. The medium, supplemented with freshly prepared DAC solution was replaced every 24 hours. For the control culture the same conditions, except the DAC solution replaced by equal volume of 50% acetic acid (DAC solvent) were used. Additionally, the "mock" control, i.e. the UT-SCC-29 cell line cultured only with DMEM medium, without any additional treatment was included. After reaching 80% of confluence, cells were harvested and DNA and RNA were isolated (as described above) and the effect of Decitabine application was assessed by bisulfite pyrosequencing and real-time qPCR.
DNA was isolated from the UT-SCC-29 cell line culture, treated either with DAC or acetic acid or untreated ("mock") and bisulfite pyrosequencing of FAM107A promoter was performed using the same conditions as earlier in this study.
Likewise, RNA isolated from the respective cultures was reverse transcribed using Maxima First Strand cDNA Synthesis Kit (Thermo Fisher Scientific), according to manufacturer's protocol. Primers for RT-qPCR were designed with the application of Beacon Designer TM 7.5 software (PRIMER Biosoft International) in a way to include the intron localized between two coding exons of FAM107A. The primer sequences were verified with primer BLAST database (http://blast.ncbi.nlm.nih.gov/Blast.cgi) to confirm their specificity. As the reference, GAPDH gene was used. Primer sequences are presented in Table 3. The RT-qPCR reaction was performed with the use of HOT FIREPol ® EvaGreen ® qPCR Mix Plus (no ROX) (Solis BioDyne, Estonia), according to manufacturer's protocol. No template control (NTC) was included to indicate lack of PCR contaminations. RT-qPCR was performed with the use of CFX96 Real-Time System (BIO-RAD) and the reaction conditions are described in Supplementary Table S3. The melting curve was generated to verify the product specificity and the PCR efficiency was appointed to estimate the relative expression value of FAM107A, as described previously 34 . All calculations were performed using Gene Expression MacroTM 1.10 software.

Analysis of protein expression by immunohistochemistry. The samples designated for FAM107A
protein analysis by IHC consisted of 15 primary LSCC and 5 non-tumor controls. All samples were processed in the Department of Clinical Pathology, Collegium Medicum, Nicolaus Copernicus University in Bydgoszcz, Poland and the previously described procedures of immunohistochemistry staining were applied 35,36 . The primary rabbit polyclonal anti-FAM107A antibody (Thermo Fisher Scientific; cat. No: PA5-50409) was applied to estimate FAM107A protein expression.
The protocol was standardized using a series of positive and negative control reactions. The positive control reaction was performed on a model tissue selected according to The Human Protein Atlas (http://www.proteinatlas.org) and the antibody datasheet 37 . Therefore, consecutive 3 µm tissue sections of non-tumor fallopian tube sample were cut and used subsequently for IHC staining. The presence of the analyzed antigen was evaluated in glandular cells of fallopian tube. The nuclear/cytoplasmic/membranous expression was considered as positive FAM107A protein expression, according to reference sources. All negative control reactions were performed on additionally analyzed tissue sections, by substituting the primary antibody with a solution of 1% BSA (bovine serum albumin) diluted in PBS (phosphate buffered saline).
The paraffin blocks were cut on the manual rotary microtome (AccuCut, Sakura, Torrance, USA). 3 µm paraffin sections were prepared and mounted onto the extra adhesive slides (SuperFrostPlus, MenzelGlasser, Braunschweig, Germany). Deparaffinization, rehydratation and antigen retrieval were performed by heating sections in Epitope Retrieval Solution high-pH at 95-98 °C for 20 min. (Dako, Agilent Technologies, USA) in PT-Link (Dako). Subsequently, endogenous peroxidase activity was blocked with the use of 3% H 2 O 2 solution for 15 minutes in room temperature (RT) and the non-specific binding was blocked using 5% solution of BSA for 15 minutes in RT. Incubation with the primary rabbit polyclonal anti-FAM107A antibody (dilution 1:200) was performed overnight at 4 °C. The antibody complex was detected using EnVisionFlex Anti-Mouse/Rabbit HRP-Labeled Polymer (Dako, Agilent Technologies) and localized using 3-3′diaminobenzidine (DAB) as chromogen. Finally, tissue sections were counterstained in hematoxylin, and subsequently dehydrated, cleared in series of xylenes, and coverslipped using mounting medium (Dako, Agilent Technologies).
The evaluation of protein expression was performed at 20x original objective magnification, in the light microscope ECLIPSE E400 (Nikon Instruments Europe, Amsterdam, Netherlands). For evaluation of FAM107A expression, immunohistochemical reactions were scored on a two-point qualitative scale: 0 -absence of FAM107A staining; 1 -presence of nuclear/cytoplasmic/membranous staining of FAM107A protein. Statistical analysis. The statistical calculations were performed using non-parametric Mann-Whitney test (GraphPad Prism 7 demo version) and chi-square test (http://quantpsy.org) 38 . Significance level was defined as p value less than 0.05 (p < 0.05). The box plot charts (Fig. 1) were prepared with the use of GraphPad Prism 7 demo version.