Journal of Oncology
Volume 2018 (2018), Article ID 2439527, 13 pages
https://doi.org/10.1155/2018/2439527

Research Article

RNA-seq Reveals the Overexpression of IGSF9 in Endometrial Cancer

Zonggao Shi,1,2 Chunyan Li,1,3 Laura Tarwater,¹ Jun Li,⁴ Yang Li,² William Kaliney,1,5 Darshan S. Chandrashekar,⁶ and M. Sharon Stack1,2

¹Harper Cancer Research Institute, University of Notre Dame, South Bend, IN 46617, USA
²Department of Chemistry and Biochemistry, University of Notre Dame, South Bend, IN 46617, USA
³Department of Gynecology and Obstetrics, Shandong Provincial Hospital Affiliated to Shandong University, Jinan, Shandong Province, China
⁴Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, South Bend, IN 46617, USA
⁵Department of Pathology, The Medical Foundation, South Bend, IN 46601, USA
⁶Department of Pathology, University of Alabama School of Medicine, Birmingham, AL 35294, USA

Correspondence should be addressed to Zonggao Shi; zshi2@nd.edu

Received 25 November 2017; Revised 3 January 2018; Accepted 17 January 2018; Published 14 February 2018

Academic Editor: Akira Hara

Copyright © 2018 Zonggao Shi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We performed RNA-seq on an Illumina platform for 7 patients with endometrioid endometrial carcinoma for which both tumor tissue and adjacent noncancer tissue were available. A total of 66 genes were differentially expressed with significance level at adjusted value < 0.01. Using the gene functional classification tool in the NIH DAVID bioinformatics resource, 5 genes were found to be the only enriched group out of that list of genes. The gene IGSF9 was chosen for further characterization with immunohistochemical staining of a larger cohort of human endometrioid carcinoma tissues. The expression level of IGSF9 in cancer cells was significantly higher than that in control glandular cells in paired tissue samples from the same patients () or in overall comparison between cancer and the control (). IGSF9 expression is higher in patients with myometrium invasion relative to those without invasion (). Reanalysis of RNA-seq dataset from The Cancer Genome Atlas shows higher expression of IGSF9 in endometrial cancer versus normal control and expression was associated with poor prognosis. These results suggest IGSF9 as a new biomarker in endometrial cancer and warrant further studies on its function, mechanism of action, and potential clinical utility.

1. Introduction

Precision medicine calls for better characterization of diseases at molecular level. As the cost of next-generation sequencing (NGS) technology becomes increasingly affordable, more and more efforts have been invested on genome wide observation of genetic changes and gene expression profiles in cancer. NGS-based RNA-sequencing (RNA-seq) is an established and effective screening method in identifying new biomarkers and better understanding cancer biology [1]. Nationwide and international consortia like The Cancer Genome Atlas (TCGA) and the International Cancer Genomics Consortium (ICGC) have generated tremendous amounts of new data and new insights regarding the molecular landscape of many major types of cancer using various high throughput profiling technologies including RNA-seq [2–4].

Endometrial cancer of the uterus is responsible for about 74,000 deaths of women each year worldwide [5, 6]. It is the most common gynecological malignancy in America and other western countries. Among the histologic subtypes, endometrioid carcinoma is the most prevalent, accounting for 75%–85% of all patients. Other types, like serous, mucinous, clear cell and squamous cell carcinomas, are much less common. According to the data from SEER (Surveillance, Epidemiology, and End Results) program at the National Cancer Institute (https://seer.cancer.gov/), although the current survival rate of women with endometrial cancer is as high as 81%, annual deaths from endometrial cancer still numbered over 10,000 patients in the United States alone. Furthermore, for over three decades, the age-adjusted death rate of endometrial cancer in the general population has not decreased. Thus, it is imperative to better characterize the molecular details of endometrial cancer using the latest high throughput profiling technology to identify new biomarkers that may enhance our understanding of disease progression and may impact clinical management of this cancer.

The largest effort in interrogating the transcriptome of endometrial cancer is from TCGA [7]; however both the initial comprehensive report from TCGA on endometrial cancer and the follow-up report by others who focused on the RNA-seq data of endometrial cancer in TCGA did not provide transcriptome-wide information on the differentially expressed genes in the scenario of endometrial cancer versus normal endometrial tissue [7, 8]. Another effort by Xiong et al. examined the transcriptome profiles of endometrial cancer versus adjacent noncancer tissue, but with 3 patients only [9]. In the current study, we screened paired human endometrioid carcinoma tissues and the adjacent nontumor tissues from 7 patients using RNA-seq technology and compared their gene expression profiles. Further validation was performed by quantitative PCR and immunohistochemistry. Focusing on the overexpressed gene IGSF9, we explored its associations with clinicopathologic characteristics in endometrioid endometrial carcinoma.

2. Materials and Methods

2.1. Human Tissue Specimens

Two different approaches were used in human tissue specimen procurement. For RNA-seq, freshly harvested, RNAlater-preserved cancer and adjacent noncancer tissues were obtained from 7 patients through frozen section rooms at The Medical Foundation (South Bend, IN). For the follow-up immunohistochemistry study, formalin-fixed and paraffin-embedded (FFPE) blocks from a total of 56 patients were obtained from The Medical Foundation surgical archive. All specimens were collected during the period from 2013 to 2015 via the Harper Cancer Research Institute Tissue Biorepository project. All collection of human tissues from surgically removed uterine endometrial cancer specimens were conducted by The Medical Foundation pathologists. Pathologic diagnoses and classification were initially made by the pathologist on-duty and further reviewed by board-certified pathologists. Informed consent was obtained from all patients. Protocols were approved by the University of Notre Dame Institutional Review Board (IRB) and local hospital IRBs.

2.2. RNA-Sequencing and Data Processing

For RNA preparation, ~30 mg tissue from each specimen was homogenized using an Omni-TH tissue homogenizer (OMNI International, Kennesaw, GA) in prechilled lysis buffer from the AllPrep RNA/DNA/protein mini kit (Qiagen, Germantown, MD) according to the manufacturer’s instructions. RNA elution was then quantified with a Nanodrop spectrophotometer (Thermo Fisher) and further analyzed by BioAnalyzer (Agilent) for quality control. Samples with were used in RNA-seq library preparation by the Notre Dame Genomics Core Facility. RNA-seq libraries were prepared with nonstranded TruSeq RNA kit from Illumina (San Diego, CA) and sequenced on a HiSeq 4000 sequencer (BGI Hong Kong). The reads generated were paired-end and of length 100-nt. Raw reads were cleaned up by BGI with SOAPnuke to remove adapters and low-quality reads. All other command line software tools were installed and used on the high performance computer cluster at University of Notre Dame Center for Research Computing. FastQC (version 0.11.2) was used to examine the read quality. No further read trimming or filtering was performed. Splicing aware mapping tool, STAR (version 2.4.2a) was used for alignment in 2-pass mode. The resulting BAM files were examined by Qualimap (version 2.2) and then used to count reads per gene with the htseq-count function from HTSeq software (version 0.6.1).

R/Bioconductor packages including DESeq2 were used for gene expression analysis [10–12]. Dealing with paired samples, we used a multifactor design which includes the sample information by putting the condition of interest (cancer or noncancer) at the end of the design. When the data was treated as from two independent groups, only single factor design was used. For comparison, we also downloaded and reanalyzed the raw RNA-seq data published by Xiong et al. [9] from SRA (NIH Sequence Read Archive) with our own data processing protocol as mentioned above. To identify the most overrepresented biological terms and genes, we made use of the gene classification tool from DAVID Bioinformatics Resources 6.8 [13] (https://david.ncifcrf.gov). High classification stringency was used.

To access tumor subgroup gene expression and survival data in The Cancer Genome Atlas (TCGA) project, we used the web interface of UALCAN [14] and part of the processed data from UALCAN. UALCAN uses TCGA level 3 RNA-seq and clinical data from 33 cancer types. The expression level was normalized as transcripts per million reads (TMP) for comparison across groups and individuals.

2.3. Quantitative PCR

Taqman® gene expression assays were used in quantitative PCR. The primers and probe set for the IGSF9 gene and reference control gene ACTB were bought from Life Technologies (Carlsbad, CA). RNA sample preparation and cDNA synthesis used TRIzol reagent and SuperScript II reagent, respectively, from Life Technologies. The q-PCR reaction was run on an ABI StepOnePlus real time thermal cycler (Life Technologies). Quantification was calculated with the ΔΔCt method using control samples and the expression of the reference gene ACTB for normalization [15].

2.4. Immunohistochemistry

The immunohistochemistry procedure is as previously described [16]. FFPE tissue sections (5 μm in thickness) were dewaxed in xylene and then rehydrated. Heat induced epitope retrieval was performed with citrate buffer (pH 6.0). Primary antibodies for IGSF9 (catalog number: NBP1-93676, 1 : 200 diluted) and Ki67 (catalog number: ab15580, 1 : 500 diluted) were bought from NOVUS Biologicals (Littleton, CO) and Abcam (Cambridge, MA), respectively. Polymer-based chromogen visualization system used the ImmPRESS™ HRP anti-rabbit IgG product from Vector Laboratories (Burlingame, CA).

2.5. Image Analysis and Quantification

All immunohistochemistry slides were scanned and analyzed with the digital pathology system, including Aperio ScanScope CS, eSlide manager, and image analysis tools, made by Leica Biosystem (Chicago, Illinois). Quantification of IGSF9 used a macro built from the pixel counting algorithm (version 1), which reports the percentage of total number of positively stained pixels over total number of pixels analyzed. Cellular quantification nuclear stain algorithm (version 9) was used to make a macro to count positive cell percentage of total number of tumor cells analyzed in Ki67 immunostains.

2.6. Statistical Analysis

Rstudio and related R packages were used to perform all the statistical analyses involving RNA-seq data downstream processing, clinical pathological data, and IHC results. Significance level was set at value < 0.05, while in screening the long gene lists from RNA-seq, adjusted values were described in the results. Principle component analysis was used in analyzing the variance in overall RNA-seq counts. Wald test was used when comparing differential gene expression as implemented in DESeq2 package. When comparing means of two groups in clinical and IHC, nonparametric method Wilcoxon rank sum test (also known as Mann–Whitney test) was used. Correlation tests used Pearson’s method as implemented in base R and the psych package. Survival analysis used Kaplan-Meier method and log-rank test.

3. Results

3.1. Gene Expression Profiles of Endometrioid Endometrial Carcinoma and Benign Adjacent Tissues via RNA-seq

RNA-seq was performed for a total of 7 patients with both cancer tissues and adjacent benign endometrial tissues available. Patient ages ranged from 48 to 78. All specimens were from the same major histological type, endometrioid endometrial carcinoma; two patients were at histological grade I and five at grade II, all at clinical stage I (FIGO, 2009). On average, 17.2 million reads were produced from each of the 14 libraries. Sequencing read quality was examined by FastQC, all with a Phred score above 20. Alignment or mapping quality metrics as produced by Qualimap over the output BAM files from alignment with STAR are listed in Suppl. Table 1.

While exploring the gene count data, plotting the top two components from principle component analysis did not reveal any clustering trend in the 7 cancer samples versus the 7 controls (Figure 1(a)). Using DESeq2 as the main tool, we sought to find the differential expression of genes in cancer versus the control with two different approaches: (1) using the pairing information in the samples, (2) without considering the pairing information.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1: RNA-seq identifies differentially expressed genes in endometrioid carcinoma. (a) Principle component analysis showing the distance of variance among all the 14 samples from 7 patients, no obvious clustering for cancer versus noncancerous control; (b) volcano plot depicting the number of differentially expressed genes based on their values and log2 fold change in the analysis of 7 pairs of samples; (c), (d), (e), and (f) normalized counts from the top 4 genes obtained by comparing the control and cancer tissues as two independent groups.

With the first approach, when comparing cancer tissues with noncancer controls, a total of 665 genes (357 up and 308 down) were identified as differentially expressed at the statistical significance level of adjusted (for multiple comparison) value < 0.05. The entire list of differentially expressed genes is given in Suppl. Table 2. On the top of gene list ranked by ascending adjusted value there are EGR1, FOSB, PTPRT, SELE, FOS, ENPP2, TFPI, ZFP36, and others. Those with adjusted value < 0.01 and log2 fold change (L2FC) >2 were labeled in the volcano plot in Figure 1(b). TUBB3 is the most increased (L2FC = 2.177). With a cutoff at adjusted value < 0.01, we obtained a list of 66 genes (Figure 2), which was used for analysis with NIH DAVID as described below.

Figure 2: Heat map of the top 66 differentially expressed genes. Genes shown were selected by the adjusted value < 0.01 in ascending order in the output of DESeq2 analysis on RNA-seq data from 7 pairs of samples (cancer or control). Colored legend indicates patient code, condition (cancer or control), and the relative expression level.

Using the second method, the 14 samples were treated as two independent groups, cancer versus control, with pairing information ignored, then only 4 genes were found to be differentially expressed at adjusted value < 0.05, including IGSF9, c10orf35, and ZNF710 genes overexpressed and the HHATL gene downregulated in cancer versus control. The normalized counts of these 4 genes were plotted in Figures 1(c)–1(f) and the full list of genes is given in Suppl. Table 3. Reanalysis of the RNA-seq study by Xiong et al. on paired cancer and normal tissues from 3 endometrial cancer patients produced another list of differentially expressed genes (Suppl. Table 4) [9].

With the aforementioned list of top 66 genes that were taken from the differentially expressed genes cutoff at adjusted value < 0.01, we performed the analysis with the gene function classification tool from NIH DAVID [17], which has the capability of addressing the enriched and redundant relationships among many-genes-to-many-terms and can reduce the complexity of gene list. Out of the 66 genes, 5 genes, that is, VSIG2, PTPRT, PRTG, IGSF9, and PTPRF (Figure 3), were found to share functional similarity and form the only enriched cluster, with an enrichment sore of 0.96. Among them, IGSF9 is the one that increased most in terms of log2 folds. This result from the gene function classification assay with DAVID plus the short list of differentially expressed genes we obtained when treating the samples as 2 independent groups prompted us to further evaluate the expression of IGSF9 with immunohistochemistry on human endometrial cancer tissues.

Figure 3: Gene functional classification. With either medium or high stringency, the gene functional classification tool from NIH DAVID bioinformatics resource produced only one cluster consisting of 5 genes when fed with 66 genes at adjust value < 0.01 as shown in the Supplementary Table 2. Columns are genes names and rows show the related functional categories shared by the genes indicated by the blue dots.

3.2. IGSF9 Expression in Endometrial Tissue and Endometrioid Carcinoma

The first study on IGSF9 expression was on the developing nervous system in mouse and human [18]; however to date there are no published data detailing IGSF9 expression in endometrial tissue. We first performed quantitative PCR to confirm the differential expression of IGSF9 in cancer versus noncancer endometrial tissues. With 3 paired RNA samples from the RNA-seq project, the relative expression in cancer tissue was found to be over 6 folds compared to adjacent tissue (Figure 4(a)). IGSF9 is also one of differentially expressed genes identified by our reanalysis of the RNA-seq data from Xiong et al. For side-by-side comparison, we plotted their normalized counts (Figure 4(b)).

(a)

(b)

(c)

Figure 4: Expression of IGSF9 in normal endometrium and endometrial carcinoma. (a) Q-PCR measurement of IGSF9 from RNA samples of three paired cancer and normal tissues; the Wilcoxon signed rank test gives value = 0.25, while the paired -test gives value = 0.038. (b) Reanalysis of published RNA-seq data from Xiong et al. [9]; the Wald test gives value = and adjusted value = . (c) IHC staining of IGSF9, the right half is noncancerous endometrial glands with weak brown stain, while the left side cancer cells stained strongly brown indicating higher level of expression.

As further characterization of IGSF9 expression in relation to cancer microscopic features requires immunohistochemistry, we evaluated IGSF9 expression at protein level in FFPE materials from 56 patients. Normal endometrial stroma and glands usually stain negative or weak while endometrioid carcinoma stained stronger (Figure 4(c)). The staining pattern is both membranous and cytoplasmic, but not the entire membrane, with uneven distribution inside the cells, often concentrated on one end of the stained cells. Some specimens show a punctate and particulate perinuclear staining pattern.

3.3. Correlation of IGSF9 Expression with Clinical Pathologic Features

The association of IGSF9 expression with clinicopathologic characteristics of these patients was systematically examined in our collection of 56 patients with endometrial cancer, which were aged from 39 to 94 years, with a median of 64. All patients were diagnosed as endometrioid carcinoma, 28 at grade I, and 28 at grade II in histologic grading (Table 1, Figures 5(a)-5(b)) [19]. As for clinical stages, 46 patients were at stage I, three at stage II, and seven at stage III. Tumor size, the depth of myometrial invasion, tumor cell labeling of proliferation marker Ki67 (Figures 5(c)-5(d)), and the presence or absence of lymphovascular invasion were recorded as part of clinical and pathologic data.

Table 1: Clinicopathological characteristics of patients in IHC study.

Parameters
Category/measurement

Age, mean ± SD, range (year)
(39–94)

Tumor size, mean ± SD, range (cm)
(0.3–12)

Clinical stage
I
46

II
3

III
7

Histologic type
Endometrioid
56

Serous
0

Histologic grade
G1
28

G2
28

G3
0

Myometrium invasion
Absent
3

Present
53

Lymphovascular invasion
Absent
49

Present
7

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5: Histopathological and immunohistochemical evaluation of endometrial carcinoma. (a) and (b) H&E stain of endometrioid endometrial carcinoma, histological grades I and II, respectively; (c) and (d) Ki67 stain on grade I and II sections, respectively; (e) and (f) IGSF9 stain on grade I and II sections, respectively. Both immunostains were quantified with the image analysis tools from Aperio digital pathology system. Numeric data were used in statistical analysis.

The positivity of IGSF9 expression as determined by immunohistochemistry (Figures 5(e)-5(f)) was presented in decimal format. Among the cases evaluated were eight patients for which both cancer and control (normal/noncancerous) tissues were available. Paired comparison of relative IGSF9 expression shows significantly enhanced staining in cancer tissues relative to controls with values of 0.448 versus 0.200, respectively ( value = 0.008 by the Wilcoxon test) (Figure 6(a)). Comparison of the all the cancer tissues with the 8 control tissues showed a similar trend (Figure 6(b)) with statistical significance ( value = 0.002 by the Wilcoxon test).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6: IGSF9 expression in association with categorical clinicopathologic features. (a) Comparison of IGSF9 positivity in 8 pairs endometrioid carcinoma and noncancerous endometrial epithelia ( value = 0.008). (b) Comparison between noncancerous tissues and all the endometrioid carcinoma tissues from 56 patients ( value = 0.002). (c) Histological grade I (G1) versus grade II (G2). (d) Clinical stages. (e) Lymphovascular invasion. (f) Myometrial invasion ( value = 0.013).

A trend was observed in which IGSF9 positivity is higher in histological grade II (G2) relative to grade I (G1) (Figure 6(c)), but the results are not statistically significant. No association was found in relation to clinical staging and lymphovascular invasion (Figures 6(d) and 6(e)). Interestingly, however, statistically significant difference was obtained when comparing the IGSF9 positivity between patients without any degree of myometrial invasion and those with myometrial invasion (0.151 versus 0.395; value = 0.013 by the Wilcoxon test) (Figure 6(f)).

The correlation of IGSF9 positivity with other continuous variables, such as age, tumor size in cm, depth of myometrium invasion as percentage of total myometrium thickness, and Ki67 labeling index were plotted in Figure 7. IGSF9 expression has somewhat of a correlation with Ki67 labeling index, but not statistically significant. However the positive correlation of tumor size with depth of myometrium invasion (coefficient = 0.412, value = 0.002) and that of Ki67 with depth of myometrium invasion (coefficient = 0.360, value = 0.006) were statistically significant.

Figure 7: IGSF9 expression in association with continuous clinicopathologic features. Scatter plots showing the pairwise correlations of age (years), tumor size (diameter in cm), depth of myometrium invasion (percentage of the entire myometrium thickness), Ki67 label index (positive percentage of tumor cells), and IGSF9 positivity (positive percentage of analyzed pixels) with each other in the off-diagonals. Loess smoothed (red) and linear (green) fit lines were superimposed on these plots. Each parameter was named in the principal diagonal overlaid with density and rug plots of that parameter. Depth of myometrial invasion correlated with tumor size (coefficient = 0.412, value = 0.020) and Ki67 labeling index (coefficient = 0.360, value = 0.006).

3.4. IGSF9 Expression Status in RNA-seq Data from TCGA

As a way of in silico validation, we surveyed the expression of IGSF9 in TCGA RNA-seq dataset via the web portal of UALCAN [14]. As of November 2017, normalized RNA expression data from 33 types of cancer are available, but not all have proper amount of normal controls for comparison. IGSF9 had higher expression in uterine corpus endometrial cancer (UCEC) tissues and 7 other types of cancers compared to their corresponding normal controls (Table 2). In contrast, IGSF9 expression is lower in cancer versus normal control in rectum adenocarcinoma and colon adenocarcinoma (Table 2).

Table 2: IGSF9 expression and survival correlation in TCGA RNA-seq dataset.

Abbrev.
Cancer type
Sample size⁽¹⁾
Expression level⁽²⁾
Overall survival⁽³⁾

BLCA
Bladder urothelial carcinoma
Normal ()
Tumor ()
Up
-

BRCA
Breast invasive carcinoma
Normal ()
Tumor ()
Up
-

CESC
Cervical squamous cell carcinoma and endocervical adenocarcinoma
Normal ()
Tumor ()
Up
-

CHOL
Cholangiocarcinoma
Normal ()
Tumor ()
Up
-

COAD
Colon adenocarcinoma
Normal ()
Tumor ()
Down
-

HNSC
Head and neck squamous cell carcinoma
Normal ()
Tumor ()
Up
-

LGG
Brain lower grade glioma
Grade II ()
Grade III ()
-
Yes

LUSC
Lung squamous cell carcinoma
Normal ()
Tumor ()
Up
-

READ
Rectum adenocarcinoma
Normal ()
Tumor ()
Down
-

SKCM
Skin cutaneous melanoma
Normal ()
Tumor ()
-
Yes

THCA
Thyroid carcinoma
Normal ()
Tumor ()
Up
-

THYM
Thymoma
Normal ()
Tumor ()
-
Yes

UCEC
Uterine corpus endometrial cancer
Normal ()
Tumor ()
Up
Yes

Notes. ⁽¹⁾The number of normal versus tumor samples, except for LGG. ⁽²⁾“Up” or “Down” indicates significant expression level changes, tumor versus normal; dash indicates no significant change. ⁽³⁾Kaplan-Meier analysis; “Yes” indicates that higher expression of IGSF9 is related to poor overall survival.

UALCAN contains reprocessed RNA-seq data from 546 UCEC tissues and 35 normal endometrial tissues. Out of that, 23 patients provided both cancer and normal tissues. IGSF9 mRNA expression was found to be higher in both the comparison of all cancer versus normal control ( by Wilcoxon test) and comparison of paired cancer versus normal from the same patients ( by Wilcoxon test) (Figures 8(a) and 8(b)).

(a)

(b)

(c)

Figure 8: IGSF9 expression in TCGA endometrial cancer project. (a) Normalized RNA-seq data for IGSF9, transcripts per million (TMP), clinical data of 546 patients cancer specimen versus 35 normal tissues ( value < 0.001). (b) comparison of 23 patients with RNA-seq data from both cancer tissue and normal control ( value < 0.001). (c) Kaplan-Meier survival analysis of 543 patients with IGSF9 expression and outcome data. Log-rank test, value = 0.017.

A total of 543 endometrial cancer patients with RNA-seq and overall outcome data were used for Kaplan-Meier survival analysis. Overexpression of IGSF9 (top 25%) was associated with poor survival (Figure 8(c), log-rank test, ). A statistically significant association of IGSF9 overexpression and poor overall survival was also observed in thymoma, skin cutaneous melanoma, and brain lower grade glioma (Table 2).

4. Discussion

To find differentially expressed genes between groups of specimens, RNA-seq is a very powerful tool as it is high throughput and reliable in quantification [10, 20]. For gene expression profiling with RNA-seq, TCGA has data from the largest number of patients [2]. The reports and datasets from TCGA are a rich resource for evaluation of genetic variations or molecular markers of endometrial cancer. Based on their comprehensive findings on somatic copy number alterations and exome sequencing, a new molecular classification system was proposed for endometrial cancer [7]. However, limited evaluation of the mRNA expression data was provided [7, 8]. By design, the raw RNA-seq data and clinicopathologic characteristics in all the TCGA projects are open to researchers who are interested in further analyses. Noticeably, many secondary online resources based on TCGA datasets (e.g., the UALCAN website) have recently appeared and greatly eased the use of these data [14, 21]. The report by Dellinger and colleagues utilized TCGA data to show that L1CAM is an independent predictor of poor survival in endometrial cancer [22].

The lack of information from normal tissues or benign lesions in RNA-seq data from the initial TCGA endometrial cancer publication piqued our interest and others as well. Xiong and colleagues used RNA-seq to evaluate both mRNA and microRNA expression profiles in three patients with both cancer and noncancer tissues [9]. In addition to mRNA and microRNA expression profiles, their study also performed genetic variation calling from RNA-seq data, which is a feasible but less established approach [23].

In the current study, our RNA-seq with paired cancer and control tissues from 7 endometrioid carcinoma patients discovered a unique set of differentially expressed genes. For example, our reanalysis of Xiong’s data identified WISP2 (CCN5) (log2 fold change −7.12; Suppl. Table 4) as a highly differentially expressed gene. It is reported to be a tumor suppressor; the deficiency of WISP2 promotes breast cancer growth [24] and WISP2 RNA expression was decreased in 79% of human colon cancers [25], but no studies of WISP2 in endometrial cancer have been reported. The TUBB3 gene is the most increased gene on our list of genes (Figure 1(b)) from our own patient materials. TUBB3 encodes class III β-tubulin and its overexpression has been linked to resistance to paclitaxel and correlated with poor survival in ovarian, breast, gastric, and non-small-cell lung cancers and unknown primary tumors [26, 27]. This finding was not confirmed in endometrial cancer and class III β-tubulin expression was also not correlated with clinicopathologic characteristics [27, 28].

Facing the long list of differentially expressed genes, tools for functional profiling are often used to extract biological meaning from such data. RNA-seq specific analysis packages like GOseq and SeqGSEA are available [29]. Our use of the gene functional classification tool in the NIH DAVID package and our top 66 genes provided only one cluster of genes, and treating our paired samples as two independent groups further narrowed the gene list. As IGSF9 is the only common gene from both analyses, we focused on IGSF9 and explored the association of IGSF9 protein expression with clinicopathological characteristics in endometrioid endometrial carcinoma.

IGSF9 was first cloned and characterized as a member of the immunoglobulin superfamily expressed in a wide variety of tissues in human and mouse [18]. Structurally, it contains five Ig domains, two fibronectin III domains, a single transmembrane domain, and an intracellular C-terminal and it was found to be related to dendrite arborization in the brain [30]. To date, there are no published reports on IGSF9 in any type of cancer tissue. Existing studies of IGSF9 are very limited (10 papers in PubMed as of January 2018). Our study is the first to investigate the association of its expression with cancer, but its expression and mutation profiles could be assessed via searching online databases, such as the Human Protein Atlas project (http://www.proteinatlas.org), GTEx (https://GTExPortal.org), and the Catalog of Somatic Mutations in Cancer (COSMIC, http://cancer.sanger.ac.uk/cosmic).

The expression pattern as revealed by our immunohistochemical analyses in endometrioid carcinoma confirms the localization of IGSF9 on the cellular membrane and in the cytoplasm. We confirmed the differential expression of IGSF9 in cancer versus noncancer tissues with materials from a cohort of 56 patients. The expression of IGSF9 is higher in cancer than that in normal or benign endometrial tissue. The positivity of IGSF9 is also higher in cancer with myometrium invasion than those without, likely an indicator of cancer aggressiveness. However, as these observations were based on the small and retrospective cohort of patients available, further studies on larger and prospective cohorts were needed. Analysis of the expression data from TCGA as we surveyed with UALCAN provided the necessary support. These data confirmed that IGSF9 is overexpressed in cancerous versus normal endometrial tissue as well as in 7 additional major cancer types as well. Interestingly we observed the opposite relationship in colorectal cancer, wherein IGSF9 was expressed at lower levels in cancer versus normal, suggesting a topic to be further studied.

A limitation of the current study is our lack of access to survival data for our local cohort of 56 patients. However, analysis of publicly available TCGA RNA-seq data via UALCAN provided further support for our initial finding and enabled survival analyses. Overexpression of IGSF9 is an indicator of poor prognosis in endometrial cancer. Moreover, IGSF9 overexpression is also associated with poor prognosis in several other cancer types, including brain lower grade glioma, skin cutaneous melanoma, and thymoma, suggesting that the significance of IGSF9 in cancer pathology is likely not limited to endometrial cancer.

5. Conclusions

We identified via RNA-seq a list of differentially expressed genes in endometrioid endometrial carcinomas versus noncancer controls. Focusing on IGSF9, an adhesion molecule that has not been characterized in connection with cancer pathology, we confirmed the overexpression of IGSF9 in endometrial cancer and found an association with myometrial invasion and poor outcome. These studies provide preliminary support for consideration of IGSF9 as a new biomarker in endometrial cancer with potential utility in diagnostics, prognostics, and therapeutics. Further studies on the role of IGSF9 in cancer pathology are warranted.

Abbreviations

IGSF9: | Immunoglobulin Super Family 9
NGS: | Next-generation sequencing
RNA-seq: | RNA-sequencing
SEER: | Surveillance, Epidemiology, and End Results
ICGC: | International Cancer Genomics Consortium
TCGA: | The Cancer Genome Atlas
FFPE: | Formalin-fixed and paraffin-embedded
IHC: | Immunohistochemistry
Q-PCR: | Quantitative Polymerase Chain Reaction
COSMIC: | Catalog of Somatic Mutations in Cancer.

Ethical Approval

Protocols were approved by the University of Notre Dame Institutional Review Board (IRB).

Consent

Informed consent was obtained from all patients.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Zonggao Shi, Jun Li, and M. Sharon Stack conceived, designed, and managed the project. Zonggao Shi, Chunyan Li, and Laura Tarwater performed sample collection, nucleic acid extraction, q-PCR, and immunohistochemistry. William Kaliney, Chunyan Li, and Zonggao Shi did pathology review. Zonggao Shi, Chunyan Li, Yang Li, M. Sharon Stack, Jun Li, and Darshan S. Chandrashekar did data collection and analysis. Zonggao Shi and M. Sharon Stack wrote the manuscript. All authors participated in the revision of this manuscript and approved the submission.

Acknowledgments

The authors are grateful to following individuals and affiliated institutions for their generous help: Maggie Li for initial data collection and analysis; Genomics Core Facility of University of Notre Dame for RNA-seq library preparation and sequencing contract coordination; Diane Enos and Dr. Qing Li from TMF for FFPE block retrieval. The study was supported by the Walther Cancer Foundation (M. Sharon Stack).

Supplementary Materials

Supplementary 1. Supplementary Table : QC and alignment results.

Supplementary 2. Supplementary Table : differentially expressed genes as paired tissues.

Supplementary 3. Supplementary Table : differentially expressed genes as independent groups.

Supplementary 4. Supplementary Table : differentially expressed genes from reanalyzing data of Xiong et al.

References

J. C. Marioni, C. E. Mason, S. M. Mane, M. Stephens, and Y. Gilad, “RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays,” Genome Research, vol. 18, no. 9, pp. 1509–1517, 2008.

K. Tomczak, P. Czerwińska, and M. Wiznerowicz, “The Cancer Genome Atlas (TCGA): An immeasurable source of knowledge,” Wspolczesna Onkologia, vol. 1A, pp. A68–A77, 2015.

L. Chin, W. C. Hahn, G. Getz, and M. Meyerson, “Making sense of cancer genomic data,” Genes & Development, vol. 25, no. 6, pp. 534–555, 2011.

C. International Cancer, T. J. Genome, and W. Anderson, “International network of cancer genome projects, Nature,” in Proceedings of the International network of cancer genome projects, Nature, vol. 464, p. 993, 2010.

M. Le Gallo, A. J. O'Hara, M. L. Rudd et al., “Exome sequencing of serous endometrial tumors identifies recurrent somatic mutations in chromatin-remodeling and ubiquitin ligase complex genes,” Nature Genetics, vol. 44, no. 12, pp. 1310–1315, 2012.

L. A. Torre, F. Bray, R. L. Siegel, J. Ferlay, and J. Lortet-Tieulent, “Global cancer statistics, 2012,” CA: A Cancer Journal for Clinicians, vol. 65, no. 2, pp. 87–108, 2015.

C. Kandoth, N. Schultz, A. D. Cherniack et al., “Integrated genomic characterization of endometrial carcinoma,” Nature, vol. 497, no. 7447, pp. 67–73, 2013.

H. Sun, L. Yan, R. Tu et al., “Expression Profiles of Endometrial Carcinoma by Integrative Analysis of TCGA Data,” Gynecologic and Obstetric Investigation, vol. 82, no. 1, pp. 30–38, 2017.

H. Xiong, Q. Li, S. Liu et al., “Integrated microRNA and mRNA transcriptome sequencing reveals the potential roles of miRNAs in stage I endometrioid endometrial carcinoma,” PLoS ONE, vol. 9, no. 10, Article ID e110163, 2014.

S. Anders and W. Huber, “Differential expression analysis for sequence count data,” Genome Biology, vol. 11, no. 10, article R106, 2010.

S. Anders, D. J. McCarthy, Y. Chen et al., “Count-based differential expression analysis of RNA sequencing data using R and Bioconductor,” Nature Protocols, vol. 8, no. 9, pp. 1765–1786, 2013.

G. Yu and Q.-Y. He, “ReactomePA: An R/Bioconductor package for reactome pathway analysis and visualization,” Molecular BioSystems, vol. 12, no. 2, pp. 477–479, 2016.

D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.

D. S. Chandrashekar, B. Bashel, S. A. H. Balasubramanya et al., “UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses,” Neoplasia (United States), vol. 19, no. 8, pp. 649–658, 2017.

Z. Shi, Y. Liu, J. J. Johnson, and M. S. Stack, “Urinary-type plasminogen activator receptor (uPAR) modulates oral cancer cell behavior with alteration in p130cas,” Molecular and Cellular Biochemistry, vol. 357, no. 1-2, pp. 151–161, 2011.

Z. Shi and M. Stack, “An Update on Immunohistochemistry in Translational Cancer Research,” Cancer Translational Medicine, vol. 1, no. 4, p. 115, 2015.

D. W. Huang, B. T. Sherman, Q. Tan et al., “The DAVID gene functional classification tool: a novel biological module-centric algorithm to functionally analyze large gene lists,” Genome Biology, vol. 8, no. 9, article R183, 2007.

K. Doudney, J. N. Murdoch, C. Braybrook et al., “Cloning and characterization of Igsf9 in mouse and human: A new member of the immunoglobulin superfamily expressed in the developing nervous system,” Genomics, vol. 79, no. 5, pp. 663–670, 2002.

R. A. Reichert, Diagnostic gynecologic and obstetric pathology: an atlas and text, LWW, Philadelphia, 2012.

A. Conesa, P. Madrigal, S. Tarazona et al., “A survey of best practices for RNA-seq data analysis,” Genome Biology, vol. 17, no. 1, article no. 181, p. 13, 2016.

Y. Zhu, P. Qiu, and Y. Ji, “TCGA-assembler: Open-source software for retrieving and processing TCGA data,” Nature Methods, vol. 11, no. 6, pp. 599-600, 2014.

T. H. Dellinger, D. D. Smith, C. Ouyang, C. D. Warden, J. C. Williams, and E. S. Han, “L1CAM is an independent predictor of poor survival in endometrial cancer - An analysis of The Cancer Genome Atlas (TCGA),” Gynecologic Oncology, vol. 141, no. 2, pp. 336–340, 2016.

Z. Sun, A. Bhagwate, N. Prodduturi, P. Yang, and J. A. Kocher, “Indel detection from RNA-seq data: tool evaluation and strategies for accurate detection of actionable mutations,” Briefings in Bioinformatics, p. bbw069.

A. Das, K. Dhar, G. Maity et al., “Deficiency of CCN5/WISP-2-Driven Program in breast cancer Promotes Cancer Epithelial cells to mesenchymal stem cells and Breast Cancer growth,” Scientific Reports, vol. 7, no. 1, article no. 1220, 2017.

D. Pennica, T. A. Swanson, J. W. Welsh et al., “WISP genes are members of the connective tissue growth factor family that are up-regulated in wnt-1-transformed cells and aberrantly expressed in human colon tumors,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 95, no. 25, pp. 14717–14722, 1998.

D. P. English, D. M. Roque, and A. D. Santin, “Class III b-tubulin overexpression in gynecologic tumors: Implications for the choice of microtubule targeted agents?” Expert Review of Anticancer Therapy, vol. 13, no. 1, pp. 63–74, 2013.

I. Vandenput, A. Capoen, L. Coenegrachts et al., “Expression of ERCC1, p53, and Class III β-tubulin do not reveal chemoresistance in endometrial cancer results from an immunohistochemical study,” International Journal of Gynecological Cancer, vol. 21, no. 6, pp. 1071–1077, 2011.

C. Zhu, J. Luo, and H. Shi, “Expression of tubulin, p53, ki67, receptors for estrogen, and progesterone in endometrial cancer,” Eur J Gynaecol Oncol, vol. 30, no. 5, p. 514, 2009.

M. D. Young, M. J. Wakefield, G. K. Smyth, and A. Oshlack, “Gene ontology analysis for RNA-seq: accounting for selection bias,” Genome Biology, vol. 11, no. 2, article R14, 2010.

S.-H. Shi, D. M. Cox, D. Wang, L. Y. Jan, and Y.-N. Jan, “Control of dendrite arborization by an Ig family member, dendrite arborization and synapse maturation 1 (Dasm1),” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 101, no. 36, pp. 13341–13345, 2004.