Recurrent amplification of RTEL1 and ABCA13 and its synergistic effect associated with clinicopathological data of gastric adenocarcinoma.

BACKGROUND: Despite progression in treatment of gastric cancer, prognosis of patients remains poor, in part due to the low rate of diagnosis during its early stages. This paradigm implies the necessity to identify molecular biomarkers for early gastric cancer diagnosis, as well as for disease monitoring, thus contributing to the development of new therapeutic approaches. In a previous study, performed by array-Comparative Genomic Hybridization, we described for the first time in literature recurrent amplification of RTEL1 and ABCA13 genes in gastric cancer. Thus, the aim of the present study was to validate recurrent amplification of RTEL1 and ABCA13 genes and associate CNV status with clinicopathological data.
FINDINGS: Results showed RTEL1 and ABCA13 amplification in 38 % of samples. Statistical analysis demonstrated that RTEL amplification is more common in older patients and more associated with intestinal type and ABCA13 amplification increases the risk of lymph node metastasis and is more common in men. Co-amplification of these genes showed a significant association with advanced staging.
CONCLUSIONS: aCGH is a very useful tool for investigating novel genes associated with carcinogenesis and RTEL1 amplification may be important for the development of gastric cancer in older patients, besides being a probable event contributing for chromosomal instability in intestinal gastric carcinogenesis. ABCA13 amplification may have age-specific function and could be considered a useful marker for predicting lymph node metastasis in resected gastric cancer patients in early stage. Lastly, RTEL1 and ABCA13 synergistic effect may be considered as a putative marker for advanced staging in gastric cancer patients.

Background

Gastric cancer is the fifth most frequent type of cancer and the third cause of cancer mortality worldwide [1].

The estimate for Brazil (2016 and 2017) indicates the occurrence of about 600,000 new cases of cancer. Except for non-melanoma skin cancer (approximately 180,000 new cases), there will be about 420,000 new cases of cancer. Among all different types of cancer that affect humans, gastric cancer, excluding non-melanoma skin cancer, ranks fourth as the most frequent tumor type in men and fifth in women [2]. In the Northern Brazil, excluding non-melanoma skin cancer, gastric cancer is the second most frequent cancer among men and the fourth among women [2].

The incidence rate of gastric cancer have decreased overall in recent years, however, it remains the leading cause of cancer-related mortality in developing countries [3]. Despite progression in treatment of advanced gastric cancer, the prognosis of gastric cancer patients remains poor, in part due to the low rate of diagnosis during its early stages [4]. This paradigm implies the necessity to search and identify molecular biomarkers for early gastric cancer diagnosis, as well as for disease monitoring, thus contributing to the development of new therapeutic approaches [5].

With regard to gastric cancer treatment, with the exception of trastuzumab (therapy based on the overexpression of HER2 protein and/or the amplification of its gene ERBB2), chemotherapy of localized and advanced gastric cancer still does not consider genotypic tumor characteristics. This implies that part of the patients, if not the majority, receives medical treatment with suboptimal or even lacking efficacy [6, 7].

Recently, numerous studies have investigated the molecular basis of gastric cancer involving the alteration of pathogenesis, including the mechanisms of invasion and metastasis. With the development of modern technologies, various novel biomarkers had been identified that appear to possess diagnostic and prognostic value [4]. Among the types of biomarkers, copy number variation (CNV) in key genes has been found in gastric cancer [8-10].

In a previous study, preformed by array-Comparative Genomic Hybridization (aCGH), we observed many genomic alterations in 22 patients and described for the first time in literature the recurrent amplification of RTEL1 (Regulator of Telomere Length 1) gene, located on 20q13.33 (OMIM: 608833), and ABCA13 (ATP-Binding Cassette, Sub-Family A, Member 13) gene, located on 7p12.3 (OMIM: 607807), in gastric adenocarcinoma, with a frequency of 50 and 23 %, respectively [11].

RTEL1 is an essential helicase that has been demonstrated to be required for the maintenance of telomere length and genomic stability. Thus, RTEL1 dysfunction is dramatically mutagenic and plays an important role in tumor initiation and progression [12, 13].

ABCA13 gene is a member of ABC (ATP-binding cassette) family of transporters that plays a crucial role in the development of resistance by the efflux of anticancer agents outside of cancer cells [14]. Recently, several studies have associated overexpression of ABCA13 with poor prognosis of cancer [15, 16].

Therefore, the aim of the present study was to validate the recurrent amplification of RTEL1 and ABCA13 genes observed previously by aCGH and associate CNV status with clinicopathological data of patients.

It is noteworthy to mention that we have chosen RTEL1 gene due to high frequency of its amplification in gastric adenocarcinoma samples observed by aCGH. On the other hand, ABCA13 was chosen due to its significant association with serosal extravasation. For both genes, we have also taken into account the agreement of results with data present in literature.

CNV analysis demonstrated a significant association between gene amplification and co-amplification with clinicopathological data of patients with gastric adenocarcinoma.

Results

RTEL1 gene amplification was observed in 38 % of samples (Fig. 1). Statistical analysis showed that this amplification is 2.6 times more common in patients older than 50 years (p = 0.045; 95 % CI = 1.011–5.128) and 2.3 times more associated with intestinal type (p = 0.034; 95 % CI = 1.057–6.3) (Table 1).

Fig. 1

Copy number bar plot of RTEL1 gene showing high frequency of amplification in gastric adenocarcinoma samples

Table 1

Clinicopathological data of patients with and without RTEL1 and ABCA13 amplification

	RTEL1 CNV		p-value	ABCA13 CNV		p-value	RTEL1 and ABCA13 amplification		p-value
	≥3 copies	Others		≥3 copies	Others		Present	Absent
Gender
Male	32	53	0.987	38	42	0.003*	17	61	0.031*
Female	15	25		7	30		2	34
Age
≤50 years	8	27	0.034*	10	21	0.408	2	29	0.074
>50 years	39	51		35	51		17	66
Histopathology
Intestinal	36	46	0.045*	34	43	0.079	15	61	0.214
Diffuse	11	32		11	29		4	34
Localization
Cardia	12	12	0.197	9	15	0.901	4	19	0.969
Non-cardia	35	63		35	55		15	73
Stage
I–II	17	37	0.218	16	34	0.215	4	43	0.05*
III–IV	30	41		29	38		15	52
pN
N0	9	19	0.395	5	19	0.033*	2	23	0.157
N1 or more	37	53		39	48		17	67
pT
T1–T3	34	62	0.359	36	53	0.431	13	74	0.375
T4	13	16		9	19		6	21
pM
M0	10	15	0.746	10	13	0.741	4	17	0.360
M1	12	15		13	14		8	18

M male, F female, pN lymph node metastasis status, N0 without lymph node metastasis, N1 or more metastasis in one ore more lymph nodes, pT extent of the primary tumor, T1–T3 without serosal extravasation, T4 with serosal extravasation, pM, distant metastasis status, M1 with distant metastasis, M0 without distant metastasis. *Significant difference between groups with and without amplification, p ≤ 0.05, Chi-square test

ABCA13 gene amplification was also observed in 38 % of samples (Fig. 2). Statistical analysis showed that ABCA13 amplification increases 3 times the risk of lymph node metastasis (p = 0.033; 95 % CI = 1.057–9.018). Additionally, this amplification is 4 times more common in men (p = 0.003; 95 % CI = 1.526–9.851) and demonstrated an inconclusive association with intestinal type (p = 0.079) (Table 1).

Fig. 2

Copy number bar plot of ABCA13 gene showing high frequency of amplification in gastric adenocarcinoma samples

Synergistic effect of these two amplifications was also evaluated. The results showed that RTEL1 and ABCA13 amplification in the same patient was 3 times more associated with advanced stage (III and IV) (p = 0.05; 95 % CI = 0.958–10.037) (Table 1).

Subcategorization of samples into intestinal and diffuse types did not result in any significant clinicopathological association.

It is important to note that 51 % (23/45) and 21 % (10/47) of ABCA13 and RTEL1 amplifications, respectively, were high-level amplifications (≥5 copies), but subdivision of samples into high-level amplifications and other amplifications (3 and 4 copies) or high-level amplifications and all other (1, 2, 3 and 4 copies) did not result in any significant clinicopathological association.

Regarding oral squamous cell carcinoma, we observed 30 and 25 % of RTEL1 and ABCA13 amplification, respectively, but the presence of amplification was not statistically associated with clinicopathological data of patients.

Discussion

High amplification frequencies observed in the present study corroborate with the previous study performed by aCGH, indicating that this technique is useful to investigate CNV and find novel genes associated with diseases, even with a low number of cases.

The statistical association of ABCA13 with serosal extravasation [11] was not confirmed by real time PCR and new significant associations were observed, including for RTEL1 gene. This discrepancy is probably due to differences in sample size of each investigation.

The region 20q, where RTEL1 gene is located, is amplified in several types of cancer [17-19], but we were the first group to describe RTEL1 amplification in gastric cancer [11].

Rtel1 is an essential helicase for telomere maintenance and the regulation of homologous recombination (HR) [12]. HR is one of the major pathways to maintain genomic stability and is involved in the repair of complex DNA damage, DSBs, interstrand crosslinks and DNA gaps. Thus, upregulated Rtel1 function might prevent HR when it is needed as a legitimate means for repair, leading to malfunction of repair system [20].

Another hypothesis is that the excessive activity of Rtel1 would increase Rtel1 helicase activity, leading to disengaging of T-loop structure and, consequently, to telomere deprotection, which may result in end-to-end fusions and exonucleolytic attack [13, 21].

Several observations suggested a possible role for RTEL1 during DNA replication. Mouse cells deficient for RTEL1 exhibit reduced proliferative capacity, and worms and mammalian cells lacking RTEL1 are particularly sensitive to DNA damaging agents that hinder DNA replication, such as inter-strand crosslinking agents [13, 20, 22].

Wu et al. [12] showed that increased expression of Rtel1 in mouse hepatocytes induced the development of liver tumors. This finding is consistent with human genetic data that showed that amplification of RTEL1 genomic locus is not only a common genetic alteration in human hepatocellular carcinoma, but also closely associated with its malignancy and progression [23-27].

In this study, we found for the first time recurrent RTEL1 amplification statistically associated with advanced age and intestinal type of gastric adenocarcinoma. Thus, we suggest that amplification of RTEL1 may have age-specific function and an important role in adenocarcinoma of intestinal type, which corroborates with the hypothesis that these two histological types have different genetic pathways [28, 29].

Consistent with these results, El-Rifai et al. [30] and Kokkola et al. [31] found a significant association between 20q amplification and intestinal type of gastric cancer. Interestingly, a molecular classification of gastric cancer, proposed by Cancer Genome Atlas Research Network [32], categorized intestinal gastric cancer as correlated with chromosomal instability. Thus, RTEL1 may be a key gene of 20q region, since its upregulation triggers chromosomal instability [13, 21].

Noteworthy, we observed that the frequency of RTEL1 amplification is almost equal in all stages of cancer (I-41 %, II-30 %, III-44 % and IV-44 %), corroborating with the hypothesis proposed by Tabach et al. [33] that amplification of the 20q chromosomal arm occurs early in tumorigenic transformation and may initiate cancer.

The human ABC transporters are encoded by a large transporter gene superfamily, which is composed of 49 members grouped into seven subfamilies (A–G) according to the sequence homology. ABC proteins facilitate translocation of heterogeneous substrates including metabolic products, lipids and sterols, peptides and proteins, saccharides, amino acids and drugs across the cell membrane. To transport these substrates across extracellular and intracellular membranes against a concentration gradient, ABCs use energy acquired by the hydrolysis of ATP [34].

ABCA13 is a member of ABC gene subfamily A (ABCA) that plays a crucial role in the development of resistance by the efflux of anticancer agents outside of cancer cells [14] and overexpression of one or more membrane-bound ATP-binding cassette (ABC) transportes has been associated with such mechanism of drug resistance [35].

There are few studies in literature regarding the role of ABCA13 in cancer, but they demonstrate a positive association between ABCA13 upregulation and unfavorable outcomes.

Upregulation of ABCA12, ABCA13, ABCB6, ABCC1, ABCC2 and ABCE1 genes was found by Hlavata et al. [14] in colorectal cancer samples when compared to normal tissues and Nymoen et al. [15] observed that ABCA13 mRNA overexpression was significantly related to shorter overall survival in metastatic ovarian serous carcinoma.

Importantly, Hlaváč et al. [16] stated that ABCA13, ABCB3 and ABCC1 levels were significantly higher in patients with grade 3 than in patients with grade 1 or 2 of breast carcinoma, suggesting that overexpression of these genes may be associated with poor prognosis.

In the present study, we found for the first time recurrent amplification of ABCA13 statistically associated with lymph node metastasis in gastric carcinogenesis. Therefore, we suggest that amplification of ABCA13 gene has an important role in development of lymph node metastasis, which is associated with poor outcomes [36].

Moreover, we observed a significant association between ABCA13 amplification and male gender. In this regard, it is important to note that ABC transporters can move substrates in (influx) or out (efflux) of cells, such as inorganic anions, metal ions and a large number of hydrophobic compounds and metabolites across the plasma membrane [37], which may include substances that could be harmful to health. Since men are more exposed to carcinogens than woman [38-43], the influx of carcinogenic substances by ABCA13 channel could explain the high frequency of amplification of this gene in male gender related to gastric cancer.

Conclusion

aCGH is a very useful tool for investigating novel genes associated with carcinogenesis. Through this technique, we were able to identify recurrent amplification of RTEL1 and ABCA13 and this observation was validated by real-time PCR for copy number analysis on a larger number of samples and in other type of cancer, demonstrating that these genes may have important roles in the carcinogenesis process. Thus, RTEL1 amplification may be important for the development of gastric cancer in older patients, besides being a probable event contributing for chromosomal instability in intestinal gastric carcinogenesis. Moreover, ABCA13 amplification may have age-specific function and could be considered a useful marker for predicting lymph node metastasis in resected gastric cancer patients in early stage. Lastly, RTEL1 and ABCA13 synergistic effect may be considered as a putative marker for advanced staging in gastric cancer patients.

Methods

Samples

We analyzed gastric adenocarcinoma samples obtained from primary gastric tumors of patients from João de Barros Barreto University Hospital (HUJBB), located in Pará State, Brazil.

All samples were obtained before administration of chemical treatments or radiotherapy. This study was approved by HUJBB ethics committee (CAAE: 42999115.7.0000.5634) and all individuals signed a Consent Form allowing the use of biological samples and clinical data.

For RTEL1 copy number investigation we used 125 fresh frozen samples (Table 1), 68 % from male and 32 % from female patients, with a mean age of 59 years (±13). Regarding tumor site, 80 % were obtained from tumors located outside from cardia region and 20 % of tumors located in the cardia region. Of the total of samples, 57 % were collected from patients with advanced stage (III and IV) and 43 % from patients with early stage (I and II), 66 % of tumors belonging to the intestinal type and 34 % to the diffuse type of Laurén. Also with respect to staging, it was observed that 23 % of patients presented serosal extravasation (T4), 76 % presented lymph node metastasis and 52 % presented distant metastasis. For ABCA13 copy number investigation we used 117 fresh frozen samples (Table 1), because eight reactions were unsuccessful.

Additionally, we evaluated CNV status of these genes in 47 oral squamous cell carcinoma samples, in an attempt to investigate if the frequency of RTEL1 and ABCA13 amplification is also high in other type of neoplasia. Sample composition was 64 % male and 36 % female patients, with a mean age of 61 years (±13). Regarding tumor site, 80 % were obtained from tumors located in tongue or floor of mouth. Of the total of samples, 56 % were collected from patients with advanced stage (III and IV) and 44 % from patients with early stage (I and II).

Histopathology

Tumor samples were included only if at least 80 % of the sample consisted of cancer cells. Histopathological data, such as histological subtype, degree of differentiation, depth of invasion, lymph node involvement and distant metastasis were taken from pathology reports of the Department of Pathology of HUJBB. The histopathological analysis of tumor fragments was performed according to Laurén’s classification [44].

Quantitative analysis of copy number variants based on real-time PCR

Genomic DNA extraction was carried out using Gentra Puregene Kit (Qiagen, Germantown, MD, USA), according to the manufacturer’s instructions. After extraction, we evaluated quantity and quality of each sample using Nanodrop 1000 Spectophotometer (NanoDrop Technologies, Houston, TX, USA). The volume for each sample was adjusted accordingly to achieve 10 ng/μl using nuclease free water.

We performed TaqMan Copy Number Assay (Applied Biosystems, Foster City, CA). Briefly, 1 μl of 10 ng DNA was added to 5 μL of TaqMan Genotyping Master Mix (Applied Biosystems, Foster City, CA), with 0.5 μl of RTEL1 or ABCA13 probe and 3 μL of water. We measured copy number gain using the following profile: denaturation at 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. We determined relative quantification using the 7500 Rreal-time PCR system (Applied Biosystems, Foster City, CA) in quadruplicate. RNaseP (Applied Biosystems, Foster City, CA) was used as a control. After amplification, we imported the experiment results containing threshold-cycle values for the copy number and reference assay into the CopyCaller Software (Applied Biosystems, Foster City, CA) for post-PCR data analysis as previously described by Graziano et al. [45].

Statistics

Statistical analysis for comparisons of categorical variables between groups were done by means of Chi-square test and were performed using PASW Statistics program. Odds Ratio (OR) and Confidence Interval (CI = 95 %) were also calculated. A two-tailed probability value p ≤ 0.05 was considered to be statistically significant.

Abbreviations

ABC, ATP-binding cassette; ABCA12, ATP-Binding Cassette, Sub-Family A, Member 12; ABCA13, ATP-Binding Cassette, Sub-Family A, Member 13; ABCB3, ATP-Binding Cassette, Sub-Family B, Member 3; ABCB6, ATP-Binding Cassette, Sub-Family B, Member 6; ABCC1, ATP-Binding Cassette, Sub-Family C, Member 1; ABCC2, ATP-Binding Cassette, Sub-Family C, Member 2; ABCE1, ATP-Binding Cassette, Sub-Family E, Member 1; aCGH, array Comparative Genomic Hybridization; CNV, copy number variation; DSBs, Double Strand Breaks; ERBB2, Erb-B2 Receptor Tyrosine Kinase 2; HER2, human epidermal growth factor receptor 2; HR, homologous recombination; HUJBB, João de Barros Barreto University Hospital; PCR, polymerase chain reaction; RTEL1, Regulator of Telomere Length 1