Semaphorin-plexin signalling genes associated with human breast tumourigenesis

Introduction

Gene expression profiling has enabled us to demonstrate the heterogeneity of breast cancers. The potential of a tumour to grow and metastasise is partly dependant on its ability to initiate angiogenesis or growth and remodelling of new blood vessels, usually from a pre-existing vascular network, to ensure delivery of oxygen, nutrients, and growth factors to rapidly dividing transformed cells along with access to the systemic circulation. Cell-cell signalling of semaphorin ligands through interaction with their plexin receptors is important for the homeostasis and morphogenesis of many tissues and has been widely studied for a role in neural connectivity, cancer, cell migration and immune responses. This study investigated the role of four semaphorin/plexin signalling genes in human breast cancers in vivo and in vitro.

Materials and methods

mRNA was extracted from formalin fixed paraffin embedded archival breast invasive ductal carcinoma tissue samples of progressive grades (grades I-III) and compared to tissue from benign tumours. Gene expression profiles were determined by microarray using the Affymetrix GeneChip® Human Genome U133 Plus 2.0 Arrays and validated by Q-PCR using a Corbett RotorGene 6000. Following validation, the gene expression profile of the identified targets was correlated with those of the human breast cancer cell lines MCF-7 and MDA-MD-231.

Results

The array data revealed that 888 genes were found to be significantly (p ≤ 0.05) differentially expressed between grades I and II tumours and 563 genes between grade III and benign tumours. From these genes, we identified four genes involved in semaphorin-plexin signalling including SEMA4D which has previously been identified as being involved in increased angiogenesis in breast cancers, and three other genes, SEMA4F, PLXNA2 and PLXNA3, which in the literature were associated with tumourigenesis, but not directly in breast tumourigenesis. The microarray analysis revealed that SEMA4D was significantly (P = 0.0347) down-regulated in the grade III tumours compared to benign tumours; SEMA4F, was significantly (P = 0.0159) down-regulated between grades I and II tumours; PLXNA2 was significantly (P = 0.036) down-regulated between grade III and benign tumours and PLXNA3 significantly (P = 0.042) up-regulated between grades I and II tumours. Gene expression of SEMA4D was validated using Q-PCR, demonstrating the same expression profile in both data sets. When the sample set was increased to incorporate more cases, SEMA4D continued to follow the same expression profile, including statistical significance for the differences observed and small standard deviations. In vitro the same pattern was present where expression for SEMA4D was significantly higher in MDA-MB-231 cells when compared to MCF-7 cells. The expression of SEMA4F, PLXNA2 and PLXNA3 could not be validated using Q-PCR, however in vitro analysis of these three genes revealed that both SEMA4F and PLXNA3 followed the microarray trend in expression, although they did not reach significance. In contrast, PLXNA2 demonstrated statistical significance and was in concordance with the literature.

Discussion

We, and others, have proposed SEMA4D to be a gene with a potentially protective effect in benign tumours that contributes to tumour growth and metastatic suppression. Previous data supports a role for SEMA4F as a tumour suppressor in the peripheral nervous system but our data seems to indicate that the gene is involved in tumour progression in breast cancer. Our in vitro analysis of PLXNA2 revealed that the gene has higher expression in more aggressive breast cancer cell types. Finally, our in vitro analysis on PLXNA3 also suggest that this gene may have some form of growth suppressive role in breast cancer, in addition to a similar role for the gene previously reported in ovarian cancer. From the data obtained in this study, SEMA4D may have a role in more aggressive and potentially metastatic breast tumours.

Conclusions

Semaphorins and their receptors, the plexins, have been implicated in numerous aspects of neural development, however their expression in many other epithelial tissues suggests that the semaphorin-plexin signalling system also contributes to blood vessel growth and development. These findings warrant further investigation of the role of semaphorins and plexins and their role in normal and tumour-induced angiogenesis in vivo and in vitro. This may represent a new front of attack in anti-angiogenic therapies of breast and other cancers.     

Identification of Gene-Expression Signatures and Protein Markers for Breast Cancer Grading and Staging

The grade of a cancer is a measure of the cancer''s malignancy level, and the stage of a cancer refers to the size and the extent that the cancer has spread. Here we present a computational method for prediction of gene signatures and blood/urine protein markers for breast cancer grades and stages based on RNA-seq data, which are retrieved from the TCGA breast cancer dataset and cover 111 pairs of disease and matching adjacent noncancerous tissues with pathologists-assigned stages and grades. By applying a differential expression and an SVM-based classification approach, we found that 324 and 227 genes in cancer have their expression levels consistently up-regulated vs. their matching controls in a grade- and stage-dependent manner, respectively. By using these genes, we predicted a 9-gene panel as a gene signature for distinguishing poorly differentiated from moderately and well differentiated breast cancers, and a 19-gene panel as a gene signature for discriminating between the moderately and well differentiated breast cancers. Similarly, a 30-gene panel and a 21-gene panel are predicted as gene signatures for distinguishing advanced stage (stages III-IV) from early stage (stages I-II) cancer samples and for distinguishing stage II from stage I samples, respectively. We expect these gene panels can be used as gene-expression signatures for cancer grade and stage classification. In addition, of the 324 grade-dependent genes, 188 and 66 encode proteins that are predicted to be blood-secretory and urine-excretory, respectively; and of the 227 stage-dependent genes, 123 and 51 encode proteins predicted to be blood-secretory and urine-excretory, respectively. We anticipate that some combinations of these blood and urine proteins could serve as markers for monitoring breast cancer at specific grades and stages through blood and urine tests.     

DNA content and p53 protein expression in ductal breast cancer

DNA content and p53 protein expression in ductal breast cancer
The DNA content of 85 ductal breast cancers of different histological grades was evaluated using static cytometry and correlated with immunocytochemical expression of p53 protein in tumour cells in cytological material. A statistically significant difference was observed between p53 protein expression and grade of malignancy ( P <0.001). The percentage of euploid tumours significantly decreased from grade I through grade II to grade III tumours ( P <0.001). Clonal DNA heterogeneity was observed in 26.6% of cases analysed and was correlated with p53 protein expression ( P <0.001). These changes probably reflect genomic alterations which may affect potential malignancy of breast cancer.     

Association of Single Nucleotide Polymorphisms in Wnt Signaling Pathway Genes with Breast Cancer in Saudi Patients

Breast cancer is a complex heterogeneous disease involving genetic and epigenetic alterations in genes encoding proteins that are components of various signaling pathways. Candidate gene approach have identified association of genetic variants in the Wnt signaling pathway genes and increased susceptibility to several diseases including breast cancer. Due to the rarity of somatic mutations in key genes of Wnt pathway, we investigated the association of genetic variants in these genes with predisposition to breast cancers. We performed a case-control study to identify risk variants by examining 15 SNPs located in 8 genes associated with Wnt signaling. Genotypic analysis of individual locus showed statistically significant association of five SNPs located in β-catenin, AXIN2, DKK3, SFRP3 and TCF7L2 with breast cancers. Increased risk was observed only with the SNP in β-catenin while the other four SNPs conferred protection against breast cancers. Majority of these associations persisted after stratification of the cases based on estrogen receptor status and age of on-set of breast cancer. The rs7775 SNP in exon 6 of SFRP3 gene that codes for either arginine or glycine exhibited very strong association with breast cancer, even after Bonferroni''s correction. Apart from these five variants, rs3923086 in AXIN2 and rs3763511 in DKK4 that did not show any association in the overall population were significantly associated with early on-set and estrogen receptor negative breast cancers, respectively. This is the first study to utilize pathway based approach to identify association of risk variants in the Wnt signaling pathway genes with breast cancers. Confirmation of our findings in larger populations of different ethnicities would provide evidence for the role of Wnt pathway as well as screening markers for early detection of breast carcinomas.     

Gene expression patterns in ovarian carcinomas

We used DNA microarrays to characterize the global gene expression patterns in surface epithelial cancers of the ovary. We identified groups of genes that distinguished the clear cell subtype from other ovarian carcinomas, grade I and II from grade III serous papillary carcinomas, and ovarian from breast carcinomas. Six clear cell carcinomas were distinguished from 36 other ovarian carcinomas (predominantly serous papillary) based on their gene expression patterns. The differences may yield insights into the worse prognosis and therapeutic resistance associated with clear cell carcinomas. A comparison of the gene expression patterns in the ovarian cancers to published data of gene expression in breast cancers revealed a large number of differentially expressed genes. We identified a group of 62 genes that correctly classified all 125 breast and ovarian cancer specimens. Among the best discriminators more highly expressed in the ovarian carcinomas were PAX8 (paired box gene 8), mesothelin, and ephrin-B1 (EFNB1). Although estrogen receptor was expressed in both the ovarian and breast cancers, genes that are coregulated with the estrogen receptor in breast cancers, including GATA-3, LIV-1, and X-box binding protein 1, did not show a similar pattern of coexpression in the ovarian cancers.     

Toward the precision breast cancer survival prediction utilizing combined whole genome-wide expression and somatic mutation analysis

Background

Breast cancer is the most common type of invasive cancer in woman. It accounts for approximately 18% of all cancer deaths worldwide. It is well known that somatic mutation plays an essential role in cancer development. Hence, we propose that a prognostic prediction model that integrates somatic mutations with gene expression can improve survival prediction for cancer patients and also be able to reveal the genetic mutations associated with survival.

Method

Differential expression analysis was used to identify breast cancer related genes. Genetic algorithm (GA) and univariate Cox regression analysis were applied to filter out survival related genes. DAVID was used for enrichment analysis on somatic mutated gene set. The performance of survival predictors were assessed by Cox regression model and concordance index(C-index).

Results

We investigated the genome-wide gene expression profile and somatic mutations of 1091 breast invasive carcinoma cases from The Cancer Genome Atlas (TCGA). We identified 118 genes with high hazard ratios as breast cancer survival risk gene candidates (log rank p?<? 0.0001 and c-index?=?0.636). Multiple breast cancer survival related genes were found in this gene set, including FOXR2, FOXD1, MTNR1B and SDC1. Further genetic algorithm (GA) revealed an optimal gene set consisted of 88 genes with higher c-index (log rank p?<? 0.0001 and c-index?=?0.656). We validated this gene set on an independent breast cancer data set and achieved a similar performance (log rank p?<? 0.0001 and c-index?=?0.614). Moreover, we revealed 25 functional annotations, 15 gene ontology terms and 14 pathways that were significantly enriched in the genes that showed distinct mutation patterns in the different survival risk groups. These functional gene sets were used as new features for the survival prediction model. In particular, our results suggested that the Fanconi anemia pathway had an important role in breast cancer prognosis.

Conclusions

Our study indicated that the expression levels of the gene signatures remain the effective indicators for breast cancer survival prediction. Combining the gene expression information with other types of features derived from somatic mutations can further improve the performance of survival prediction. The pathways that were associated with survival risk suggested by our study can be further investigated for improving cancer patient survival.

     

Differentiation of tumours of ductal and lobular origin: II. Genomics of invasive ductal and lobular breast carcinomas

Breast cancer is considered to be a multifactorial disorder caused by both genetic and non-genetic factors. Different histological types of breast cancer differ in response to treatment and may have a divergent clinical course. Breast tissue is heterogeneous, with components of epithelial, mesenchymal, endothelial and lymphopoietic derivation. The genetic heterogeneity of invasive breast cancer is reflected by the wide spectrum of histological types and differentiation grades. Nevertheless, the influences of these cell types on the tumour's total pattern of gene expression can be estimated analytically. Microarrays permit total tissue analysis and provide a stable molecular portrait of tumours. Some investigations suggest differences in the gene expression profiling for ductal and lobular carcinomas. It has been reported that inactivating mutations of the E-cadherin gene are very frequent in infiltrating lobular breast carcinomas. Other than altered expression of E-cadherin, little is known about the underlying biology that distinguishes ductal and lobular tumour subtypes. However, about 8 genes have been identified differentially which are expressed in lobular and ductal cancers: E-CD, survivin, cathepsin B, TPI1, SPRY1, SCYA14, TFAP2B, and thrombospondin 4, osteopontin, HLA-G, and CHC1. Expression profiling of breast cancers can be used diagnostically to distinguish individual histologic subclassifications and may guide the selection of target therapeutics. However, future approaches will need to include methods for high throughput clinical validation and the ability to analyze microscopic samples.     

Novel markers of gene methylation and expression in breast cancer

An optimized methylation-sensitive restriction fingerprinting technique was used to search for differentially methylated CpG islands in the tumor genome and detected seven genes subject to abnormal epigenetic regulation in breast cancer: SEMA6B, BIN1, VCPIP1, LAMC3, KCNH2, CACNG4, and PSMF1. For each gene, the rate of promoter methylation and changes in expression were estimated in tumor and morphologically intact paired specimens of breast tissue (N = 100). Significant methylation rates of 38, 18, and 8% were found for SEMA6B, BIN1, and LAMC3, respectively. The genes were not methylated in morphologically intact breast tissue. The expression of SEMA6B, BIN1, VCPIP1, LAMC3, KCNH2, CACNG4, and PSMF1 was decreased in 44–94% of tumor specimens by the real-time RT-PCR assay. The most profound changes in SEMA6B and LAMC3 suggest that these genes can be included in biomarker panels for breast cancer diagnosis. Fine methylation mapping of the most frequently methylated CpG islands (SEMA6B, BIN1, and LAMC3) provides a fundamental basis for developing efficient methylation tests for these genes.     

LIN28A Modulates Splicing and Gene Expression Programs in Breast Cancer Cells

LIN28 is an evolutionarily conserved RNA-binding protein with critical functions in developmental timing and cancer. However, the molecular mechanisms underlying LIN28''s oncogenic properties are yet to be described. RNA-protein immunoprecipitation coupled with genome-wide sequencing (RIP-Seq) analysis revealed significant LIN28 binding within 843 mRNAs in breast cancer cells. Many of the LIN28-bound mRNAs are implicated in the regulation of RNA and cell metabolism. We identify heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1), a protein with multiple roles in mRNA metabolism, as a LIN28-interacting partner. Subsequently, we used a custom computational method to identify differentially spliced gene isoforms in LIN28 and hnRNP A1 small interfering RNA (siRNA)-treated cells. The results reveal that these proteins regulate alternative splicing and steady-state mRNA expression of genes implicated in aspects of breast cancer biology. Notably, cells lacking LIN28 undergo significant isoform switching of the ENAH gene, resulting in a decrease in the expression of the ENAH exon 11a isoform. The expression of ENAH isoform 11a has been shown to be elevated in breast cancers that express HER2. Intriguingly, analysis of publicly available array data from the Cancer Genome Atlas (TCGA) reveals that LIN28 expression in the HER2 subtype is significantly different from that in other breast cancer subtypes. Collectively, our data suggest that LIN28 may regulate splicing and gene expression programs that drive breast cancer subtype phenotypes.     

<Emphasis Type="Italic">FAT4</Emphasis> hypermethylation and grade dependent downregulation in gastric adenocarcinoma

Gastric cancer is one of the major causes of death due to cancer in the world. It is a multi-factorial disease with epigenetic factors being also involved in its development. FAT4 is a tumor suppressor gene exerting an important role in cell adhesion. This study aimed at analyzing FAT4 expression and promoter methylation in gastric cancer. FAT4 expression was studied in 30 tumoral tissues and their non-tumoral counterparts using Taqman real time PCR method. Promoter methylation was assessed using bisulfite conversion method followed by sequencing. Tumor tissues showed reduced FAT4 expression (P = 0.04). FAT4 downregulation was associated with tumor grade, with higher repression at advanced grades. Significant increase of promoter methylation was observed in tumoral tissues. Reduced expression of FAT4 and increased methylation of its promoter may be one of the effective processes in turning a healthy stomach tissue into a tumor tissue.     

Expression of the CD15 antigen (Lewis x) in breast cancer

Summary The expression of the cell adhesion molecule CD15 (also known as Lewis x) by breast cancers and by adjacent normal and benign breast epithelium was investigated in a series of 98 tumours. Immunohistochemistry was performed on paraffin sections using the anti-CD15 monoclonal mouse IgM antibody Dako-M1. Some 35% of cancers expressed CD15, as did 45% of normal and 60% of hyperplasia. No association was observed between cancer cell staining, or any epithelial staining (cancer, benign and normal), and tumour size, histological grade, nodal status, age at diagnosis or the frequency of ‘events’ (recurrence or death). Chi-squared tests in each case were non-significant. The pattern of CD15 expression by breast cancer was frequently associated with the leading edge of invading tumour or with the outer edge of boli of carcinomain situ, possibly suggesting a potential role in invasiveness, and with cancer cells trapped intravascularly, possibly suggesting a role in metastasis.     

Elevation of SIPL1 (SHARPIN) Increases Breast Cancer Risk

SIPL1 (Sharpin) or Sharpin plays a role in tumorigenesis. However, its involvement in breast cancer tumorigenesis remains largely unknown. To investigate this issue, we have systemically analyzed SIPL1 gene amplification and expression data available from Oncomine datasets, which were derived from 17 studies and contained approximately 20,000 genes, 3438 breast cancer cases, and 228 normal individuals. We found a SIPL1 gene amplification in invasive ductal breast cancers compared to normal breast tissues and a significant elevation of SIPL1 mRNA in breast cancers in comparison to non-tumor breast tissues. These results collectively reveal that increases in SIPL1 expression occur during breast cancer tumorigenesis. To further investigate this association, we observed increases in the SIPL1 gene and mRNA in the breast cancer subtypes of estrogen receptor (ER)+, progesterone receptor (PR)+, HER2+, or triple negative. Additionally, a gain of the SIPL1 gene correlated with breast cancer grade and the levels of SIPL1 mRNA associated with both breast cancer stages and grades. Elevation of SIPL1 gene copy and mRNA is linked to a decrease in patient survival, especially for those with PR+, ER+, or HER2- breast cancers. These results are supported by our analysis of SIPL1 protein expression using a tissue microarray containing 224 breast cancer cases, in which higher levels of SIPL1 relate to ER+ and PR+ tumors and AKT activation. Furthermore, we were able to show that progesterone significantly reduced SIPL1 mRNA and protein expression in MCF7 cells. As progesterone enhances breast cancer tumorigenesis in a context dependent manner, inhibition of SIPL1 expression may contribute to progesterone''s non-tumorigenic function which might be countered by SIPL1 upregulation. Taken together, we demonstrate a positive correlation of SIPL1 with BC tumorigenesis.     

Grade-specific diagnostic and prognostic biomarkers in breast cancer

An integrative approach is presented to identify grade-specific biomarkers for breast cancer. Grade-specific molecular interaction networks were constructed with differentially expressed genes (DEGs) of cancer grade 1, 2, and 3. We observed that the molecular network of grade3 is predominantly associated with cancer-specific processes. Among the top ten connected DEGs in the grade3, the increase in the expression of UBE2C and CCNB2 genes was statistically significant across different grades. Along with UBE2C and CCNB2 genes, the CDK1, KIF2C, NDC80, and CCNB2 genes are also profoundly expressed in different grades and reduce the patient's survival. Gene set enrichment analysis of these six genes reconfirms their role in metastatic phenotype. Moreover, the coexpression network shows a strong association of these six genes promotes cancer specific biological processes and possibly drives cancer from lower to a higher grade. Collectively the identified genes can act as potential biomarkers for breast cancer diagnosis and prognosis.     

Assessment of HER-2 gene overexpression in Isfahan province breast cancer patients using Real Time RT-PCR and immunohistochemistry

Purpose

Overexpression of proto-oncogene HER-2 is one of the main molecular markers of breast cancer involved in prognosis and diagnosis and also in trastuzumab therapy. Thus, a request for the evaluation of HER-2 status in breast cancer has been increasing. The aim of our study was assessment of HER-2 overexpression in malignant and benign breast cancer specimens by Real Time RT-PCR technique and comparison of its results with IHC outcomes.

Methods

Twenty benign and sixty malignant breast cancers in addition to fifteen normal breast tissue specimens were analyzed by Real Time RT-PCR method. Fresh tissue samples were disrupted by mortar and pestle. A syringe and a needle were used for complete homogenization of the tissues. The RNA was then isolated from the samples and converted to cDNA. A standard curve was initially plotted using BioEasy SYBR Green I and then all 95 specimens were studied by Real Time RT-PCR using 2^− ΔΔCt method.

Results

23.3% of 60 malignant specimens showed HER-2 overexpression, while all of the benign samples represented the normal expression level of HER-2 gene. The concordance rate between the results of Real Time RT-PCR and IHC was 86.6%.

Conclusion

Real Time RT-PCR method is an almost reliable technique and at least can be used as a complementary method for confirming IHC results. This is emanated from relatively high rate of concordance between outcomes of IHC test, as a routine method of detecting the HER-2 gene expression status, and Real Time RT-PCR technique.