Exon array analysis of head and neck cancers identifies a hypoxia related splice variant of LAMA3 associated with a poor prognosis

Basement membrane-associated gene expression as a predictor of survival in oral cancer

BACKGROUND

This study sought to investigate the prognostic value of basement membrane (BM)-associated gene expressions in oral cancer.

METHODS

We harvested and integrated data on BM-associated genes (BMGs), the oral cancer transcriptome, and clinical information from public repositories. After identifying differentially expressed BMGs, we used Cox and Lasso regression analyses to create a BMG-based risk score for overall survival at various intervals. We then validated this score using the GSE42743 cohort as a validation set. The prognostic potential of the risk scores and their relations to clinical features were assessed. Further, we conducted functional pathway enrichment, immune cell infiltration, and immune checkpoint analyses to elucidate the immunological implications and therapeutic potential of the BMG-based risk score and constituent genes. To confirm the expression levels of the BMG LAMA3 in clinical samples of oral cancer tissue, we performed quantitative real-time PCR (qRT-PCR) and immunohistochemical staining.

RESULTS

The BMGs LAMA3, MMP14, and GPC2 demonstrated notable prognostic significance, facilitating the construction of a BMG-based risk score. A higher risk score derived from BMGs correlated with a poorer survival prognosis for oral cancer patients. Moreover, the risk-associated BMGs exhibited a significant relationship with immune function variability (P < 0.05), discrepancies in infiltrating immune cell fractions, and immune checkpoint expressions (P < 0.05). The upregulated expression levels of LAMA3 in oral cancer tissues were substantiated through qRT-PCR and immunohistochemical staining.

CONCLUSION

The BMG-based risk score emerged as a reliable prognostic tool for oral cancer, meriting further research for validation and potential clinical application.

Comprehensive Genomic Analysis of Puerarin in Inhibiting Bladder Urothelial Carcinoma Cell Proliferation and Migration

BACKGROUND

Bladder urothelial carcinoma (BUC) ranks second in the incidence of urogenital system tumors, and the treatment of BUC needs to be improved. Puerarin, a traditional Chinese medicine (TCM), has been shown to have various effects such as anti-cancer effects, the promotion of angiogenesis, and anti-inflammation. This study investigates the effects of puerarin on BUC and its molecular mechanisms.

METHODS

Through GeneChip experiments, we obtained differentially expressed genes (DEGs) and analyzed these DEGs using the Ingenuity® Pathway Analysis (IPA®), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway enrichment analyses. The Cell Counting Kit 8 (CCK8) assay was used to verify the inhibitory effect of puerarin on the proliferation of BUC T24 cells. String combined with Cytoscape® was used to create the Protein-Protein Interaction (PPI) network, and the MCC algorithm in cytoHubba plugin was used to screen key genes. Gene Set Enrichment Analysis (GSEA®) was used to verify the correlation between key genes and cell proliferation.

RESULTS

A total of 1617 DEGs were obtained by GeneChip. Based on the DEGs, the IPA® and pathway enrichment analysis showed they were mainly enriched in cancer cell proliferation and migration. CCK8 experiments proved that puerarin inhibited the proliferation of BUC T24 cells, and its IC50 at 48 hours was 218μmol/L. Through PPI and related algorithms, 7 key genes were obtained: ITGA1, LAMA3, LAMB3, LAMA4, PAK2, DMD, and UTRN. GSEA showed that these key genes were highly correlated with BUC cell proliferation. Survival curves showed that ITGA1 upregulation was associated with poor prognosis of BUC patients.

CONCLUSION

Our findings support the potential antitumor activity of puerarin in BUC. To the best of our knowledge, bioinformatics investigation suggests that puerarin demonstrates anticancer mechanisms via the upregulation of ITGA1, LAMA3 and 4, LAMB3, PAK2, DMD, and UTRN, all of which are involved in the proliferation and migration of bladder urothelial cancer cells.

Alternative Splicing Landscape of Head and Neck Squamous Cell Carcinoma

Head and neck malignancies are a significant global health concern, with head and neck squamous cell carcinoma (HNSCC) being the sixth most common cancer worldwide accounting for > 90% of cases. In recent years, there has been growing recognition of the potential role of alternative splicing (AS) in the etiology of cancer. Increasing evidence suggests that AS is associated with various aspects of cancer progression, including tumor occurrence, invasion, metastasis, and drug resistance. Additionally, AS is involved in shaping the tumor microenvironment, which plays a crucial role in tumor development and response to therapy. AS can influence the expression of factors involved in angiogenesis, immune response, and extracellular matrix remodeling, all of which contribute to the formation of a supportive microenvironment for tumor growth. Exploring the mechanism of AS events in HNSCC could provide insights into the development and progression of this cancer, as well as its interaction with the tumor microenvironment. Understanding how AS contributes to the molecular changes in HNSCC cells and influences the tumor microenvironment could lead to the identification of new therapeutic targets. Targeted chemotherapy and immunotherapy strategies tailored to the specific AS patterns in HNSCC could potentially improve treatment outcomes and reduce side effects. This review explores the concept, types, processes, and technological advancements of AS, focusing on its role in the initiation, progression, treatment, and prognosis of HNSCC.

Upregulated LAMA3 modulates proliferation, adhesion, migration and epithelial‑to‑mesenchymal transition of cholangiocarcinoma cells

A poor outcome for cholangiocarcinoma (CCA) patients is still a clinical challenge. CCA is typically recognized by the desmoplastic nature, which accounts for its malignancy. Among various extracellular matrix proteins, laminin is the most potent inducer for CCA migration. Herein, we accessed the expression profiles of laminin gene family and explored the significance of the key laminin subunit on CCA aggressiveness. Of all 11 laminin genes, LAMA3, LAMA5, LAMB3 and LAMC2 were concordantly upregulated based on the analysis of multiple public transcriptomic datasets and also overexpressed in Thai CCA cell lines and patient tissues in which LAMA3A upregulated in the highest frequency (97%) of the cases. Differential expression genes (DEGs) analysis of low and high laminin signature groups revealed LAMA3 as the sole common DEG in all investigated datasets. Restratifying CCA samples according to LAMA3 expression indicated the association of LAMA3 in the focal adhesion pathway. Silencing LAMA3 revealed that it plays important roles in CCA cell proliferation, adhesion, migration and epithelial-to-mesenchymal transition. Taken together, this research signifies the roles of dysregulated ECM homeostasis in CCA malignancy and highlights, for the first time, the potential usage of LAMA3 as the diagnostic biomarker and the therapeutic target to tackle the CCA stromal.

Laminin N-terminus (LaNt) proteins, laminins and basement membrane regulation

Basement membranes (BMs) are structured regions of the extracellular matrix that provide multiple functions including physical support and acting as a barrier, as a repository for nutrients and growth factors, and as biophysical signalling hubs. At the core of all BMs is the laminin (LM) family of proteins. These large heterotrimeric glycoproteins are essential for tissue integrity, and differences between LM family members represent a key nexus in dictating context and tissue-specific functions. These variations reflect genetic diversity within the family, which allows for multiple structurally and functionally distinct heterotrimers to be produced, each with different architectures and affinities for other matrix proteins and cell surface receptors. The ratios of these LM isoforms also influence the biophysical properties of a BM owing to differences in their relative ability to form polymers or networks. Intriguingly, the LM superfamily is further diversified through the related netrin family of proteins and through alternative splicing leading to the generation of non-LM short proteins known as the laminin N-terminus (LaNt) domain proteins. Both the netrins and LaNt proteins contain structural domains involved in LM-to-LM interaction and network assembly. Emerging findings indicate that one netrin and at least one LaNt protein can potently influence the structure and function of BMs, disrupting the networks, changing physical properties, and thereby influencing tissue function. These findings are altering the way that we think about LM polymerisation and, in the case of the LaNt proteins, suggest a hitherto unappreciated form of LM self-regulation.

Chemerin effect on transcriptome of the porcine endometrium during implantation determined by RNA-sequencing†

It is well known that the body's metabolism and reproduction are closely related. Chemerin is one of many biologically active proteins secreted by the adipose tissue involved in the regulation of the energy homeostasis of the organism. In the present study, RNA-Sequencing (RNA-Seq) was performed to investigate the differentially expressed genes (DEGs), long non-coding RNAs (lncRNAs) and alternatively spliced (AS) transcripts in the cultured in vitro porcine endometrium exposed to chemerin for 24 hours (CHEM; 400 ng/ml) collected during the implantation period (15 to 16 days of gestation). High-throughput sequencing of transcriptomes was performed on the Illumina NovaSeq 6000 platform (Illumina, USA). In the current study, among all 130 DEGs, 58 were up-regulated and were 72 down-regulated in the CHEM-treated group. DEGs were assigned to 73 functional annotations. Twelve identified lncRNAs indicated a difference in the expression profile after CHEM administration. Additionally, we detected 386 differentially AS events encompassed 274 protein-coding genes and 2 lncRNAs. All AS events were divided into 5 alternative splicing types: alternative 3' splice site (A3SS), 5' splice site (A5SS), mutually exclusive exons (MXE), retention intron (RI), and skipping exon (SE). Within all AS events, we identified 42 A3SS, 43 A5SS, 53 MXE, 9 RI, and 239 SE. In summary, CHEM affects the transcriptomic profile of the porcine endometrium, controlling the expression of numerous genes, including those involved in the cell migration and adhesion, angiogenesis, inflammation, and steroidogenesis. It can be assumed that CHEM may be an important factor for a proper course of gestation and embryo development.

LAMA3 DNA methylation and transcriptome changes associated with chemotherapy resistance in ovarian cancer

Objective

LAMA3 is a widely studied methylated gene in multiple tumors, but the relationship between chemotherapy resistance in ovarian cancer is unclear. In this study, LAMA3 methylation was predicted by bioinformatics, and the ability of LAMA3 methylation to predict the chemotherapy resistance and prognosis of ovarian cancer was confirmed in experiments.

Methods

Multiple databases have performed the bioinformatics analysis of methylation and transcription factor binding site (TFBS) on the promoter region of LAMA3 gene. Pyrosequencing detected the methylation of LAMA3. QRT-PCR and immunohistochemistry detected the expression of LAMA3. Real Time Cell Analyzer (RTCA) detects changes in cell proliferation, migration and invasion ability. Flow cytometry was used to detect apoptosis.

Results

CPG islands of 176 bp, 134 bp, 125 bp and 531 bp were predicted in the promoter region of LAMA3 gene. The 4 prediction results are basically overlapped. 7 transcription factor binding sites were predicted, and the one with the highest score was on the predicted CpG island located in the proximal promoter region. LAMA3 hypermethylation and low expression are both associated with chemotherapy resistance and poor prognosis in ovarian cancer. LAMA3 methylation was negatively correlated with expression. After upregulation of LAMA3, the proliferation ability of chemoresistant ovarian cancer cell decreased, while the ability of apoptosis, invasion and migration increased.

Conclusion

LAMA3 hypermethylation is associated with chemotherapy resistance and poor prognosis. As a typical CpG island gene, LAMA3(cg20937934) and LAMA3(cg13270625) hypermethylation is negatively correlated with low expression. LAMA3 promotes the invasion, migration and apoptosis of SKOV3DDP. In the future, the mechanism of LAMA3 methylation in ovarian cancer will need to be further studied.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13048-021-00807-y.

Putative N-glycoprotein markers of MCF-7/ADR cancer stem cells from N-glycoproteomics characterization of the whole cell lysate

Breast cancer is one of the most malignant diseases among females. N-glycoproteomics studies have shown that N-glycosylation alteration of tumor cells is the key player of cancer progression, multidrug resistance (MDR) and high mortality. Cancer stem cells (CSCs) have the remarkable potential of self-renewing and differentiation which leads to drug resistance and metastasis. To investigate the differentially expressed N-glycosylation in adriamycin-resistant breast cancer stem cells MCF-7/ADR CSCs (relative to MCF-7 CSCs) and find the putative biomarkers, 1:1 paired ZIC-HILIC-enriched and stable isotopic diethyl labelled (SIDE) intact N-glycopeptides from MCF-7/ADR CSCs and MCF-7 CSCs were analyzed with C18-RPLC-ESI-MS/MS (HCD with stepped NCE); differentially expressed intact N-glycopeptides (DEGPs) were identified and quantified via search engine GPSeeker. With control of spectrum-level FDR≤1%, 5515 intact N-glycopeptides were identified (1737 N-glycosites, 1705 peptide backbones and 1516 intact N-glycoproteins; 181 putative N-glycan linkages and 68 monosaccharide compositions). Among 5515 intact N-glycopeptide IDs, 3864 were identified with glycoform score≥1, i.e., one or more structure-diagnostic fragment ions were observed to distinguish sequence isomers. With the three technical replicates and the criteria of fold change≥1.5 and p value<0.05, 380 DEGPs (corresponding to 153 intact N-glycoproteins) were found along with 293 down-regulated and 87 up-regulated. For these 153 intact N-glycoproteins, the molecular functions and biological processes of were comprehensively discussed, and side-to-side comparison of differential expression results with other method were also made.

Hypoxia-induced alternative splicing: the 11th Hallmark of Cancer

Hypoxia-induced alternative splicing is a potent driving force in tumour pathogenesis and progression. In this review, we update currents concepts of hypoxia-induced alternative splicing and how it influences tumour biology. Following brief descriptions of tumour-associated hypoxia and the pre-mRNA splicing process, we review the many ways hypoxia regulates alternative splicing and how hypoxia-induced alternative splicing impacts each individual hallmark of cancer. Hypoxia-induced alternative splicing integrates chemical and cellular tumour microenvironments, underpins continuous adaptation of the tumour cellular microenvironment responsible for metastatic progression and plays clear roles in oncogene activation and autonomous tumour growth, tumor suppressor inactivation, tumour cell immortalization, angiogenesis, tumour cell evasion of programmed cell death and the anti-tumour immune response, a tumour-promoting inflammatory response, adaptive metabolic re-programming, epithelial to mesenchymal transition, invasion and genetic instability, all of which combine to promote metastatic disease. The impressive number of hypoxia-induced alternative spliced protein isoforms that characterize tumour progression, classifies hypoxia-induced alternative splicing as the 11th hallmark of cancer, and offers a fertile source of potential diagnostic/prognostic markers and therapeutic targets.

The clinical significance and potential molecular mechanism of integrin subunit beta 4 in laryngeal squamous cell carcinoma

The relationship between integrin beta 4 (ITGB4) expression and laryngeal squamous cell carcinoma (LSCC) remains unclarified. The object of the present study was to explore the clinical significance and potential molecular mechanism of ITGB4 in LSCC. The protein level of ITGB4 was significantly higher in 46 LSCC patients than in 26 non-LSCC tissues detected by in-house immunohistochemistry. Consistently, ITGB4 mRNA level was also greatly upregulated based on microarray and RNA-seq data (standard mean difference, SMD = 1.62, 95 % CI: 1.23-2.00). And the area under curves (AUC) of summary receiver operator characteristic (SROC) was 0.87 (95 % CI: 0.84-0.90) based on 172 cases of LSCC and 59 cases of non-cancerous controls. Ninety genes were intersected by the ITGB4 related genes and LSCC differential expressed genes (DEGs) from all available microarray and RNA-seq datasets. Based on Gene Ontology (GO) analysis, the top terms of biological process (BP), cellular component (CC) and molecular function (MF) for the 90 ITGB4 related DEGs were extracellular matrix organization, basement membrane and extracellular matrix structural constituent, respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that ITGB4 related DEGs mainly participated in the pathways of ECM-receptor interaction, Focal adhesion and Small cell lung cancer. Moreover, the Protein-Protein Interaction (PPI) network indicated that ITGA3, ITGA5, ITGB4, MET, LAMA3, and COL4A1 might be the core genes of LSCC development related to ITGB4. In conclusion, high ITGB4 expression may lead to the occurrence and development of LSCC via various signaling pathways.

Epigenetic changes in fibroblasts drive cancer metabolism and differentiation

Genomic changes that drive cancer initiation and progression contribute to the co-evolution of the adjacent stroma. The nature of the stromal reprogramming involves differential DNA methylation patterns and levels that change in response to the tumor and systemic therapeutic intervention. Epigenetic reprogramming in carcinoma-associated fibroblasts are robust biomarkers for cancer progression and have a transcriptional impact that support cancer epithelial progression in a paracrine manner. For prostate cancer, promoter hypermethylation and silencing of the RasGAP, RASAL3 that resulted in the activation of Ras signaling in carcinoma-associated fibroblasts. Stromal Ras activity initiated a process of macropinocytosis that provided prostate cancer epithelia with abundant glutamine for metabolic conversion to fuel its proliferation and a signal to transdifferentiate into a neuroendocrine phenotype. This epigenetic oncogenic metabolic/signaling axis seemed to be further potentiated by androgen receptor signaling antagonists and contributed to therapeutic resistance. Intervention of stromal signaling may complement conventional therapies targeting the cancer cell.

Long Non-coding RNA LINC00628 Interacts Epigenetically with the LAMA3 Promoter and Contributes to Lung Adenocarcinoma

Long non-coding RNAs (lncRNAs) have emerged as key regulators of cellular progress in lung adenocarcinoma. In this study, to identify cancer-related lncRNAs and genes, we screened for those lncRNAs that were differentially expressed in lung adenocarcinoma, which revealed LINC00628 overexpression and low expression of laminin subunit alpha 3 (LAMA3). This was further validated in the cancerous tissues from patients diagnosed with lung adenocarcinoma. Thereafter, we explored the functional relevance of LINC00628 and LAMA3 in lung adenocarcinoma by analyzing the recruitment of DNA methyltransferase (DNMT) and the cellular processes of lung adenocarcinoma cells following treatments that induced LINC00628 overexpression or LINC00628 silencing or with 5-azacytidine (5-Aza, a DNMT inhibitor). The results showed that LINC00628 silencing decreased cell proliferation, migration, and invasion as well as the drug resistance of lung adenocarcinoma cells to vincristine (VCR). The results were opposite in the cells with LAMA3 demethylation induced by 5-Aza treatment. Further research indicated that LINC00628 recruited DNMT1, DNMT3A, and DNMT3B to promote the methylation of LAMA3 promoter, thereby decreasing its expression. Moreover, an in vivo experiment was performed in nude mice to assess the tumor growth ability and drug resistance of human lung adenocarcinoma cells. It was observed that LINC00628 silencing or 5-Aza treatment inhibited the in vivo tumor growth ability of the human lung adenocarcinoma cells and reduced their resistance to VCR. Altogether, our results provide evidence of a mechanism by which LINC00628 silencing exerts an inhibitory role in lung adenocarcinoma by modulating the DNA methylation of LAMA3, indicative of a novel molecular target for treatment of lung adenocarcinoma patients showing resistance to VCR.

General regulatory effects of hypoxia on human cartilage endplate‑derived stem cells: A genome‑wide analysis of differential gene expression and alternative splicing events

Intervertebral disc (IVD) degeneration of is considered to be initiated by the degeneration of the cartilage endplate (CEP). CEP‑derived stem cells (CESCs) with the capacity for osteochondrogenic differentiation may be responsible for CEP cartilage restoration. As CEP is avascular and hypoxic, and hypoxia can greatly influence biological activities of stem cells, physiological hypoxia may serve important roles in regulating the physiological functions of CESCs. The aim of the present study was to investigate the mechanisms of hypoxia‑regulated CESCs fate by using the Human Transcriptome Array 2.0 system to identify differentially expressed genes (DEGs) and alternatively spliced genes (ASGs) in CESCs cultured under hypoxic and normoxic conditions. The high‑throughput analysis of both DEGs and ASGs were notably enriched in the immune response signal, which so far has not been investigated in IVD cells, due to their avascular nature and low immunogenicity. The present results provided a referential study direction of the mechanisms of hypoxia‑regulated CESC fate at the level of gene expression and alternative splicing, which may aid in our understanding of the processes of CEP degeneration.

A genome-wide analysis of the gene expression profiles and alternative splicing events during the hypoxia-regulated osteogenic differentiation of human cartilage endplate-derived stem cells

It has been hypothesized that intervertebral disc degeneration is initiated by degeneration of the cartilage endplate (CEP), which is characterized by cartilage ossification. CEP‑derived stem cells (CESCs), with the potential for chondro‑osteogenic differentiation, may be responsible for the balance between chondrification and ossification in the CEP. The CEP remains in an avascular and hypoxic microenvironment; the present study observed that hypoxia was able to markedly inhibit the osteogenic differentiation of CESCs. This tissue‑specific CESC differentiation in response to a hypoxic microenvironment was physiologically important for the prevention of ossification in the CEP. In order to study the hypoxia‑regulated mechanisms underlying osteogenic differentiation of CESCs, a Human Transcriptome Array 2.0 was used to detect differentially expressed genes (DEGs) and alternatively spliced genes (ASGs) during the osteogenic differentiation of CESCs under hypoxia, compared with those induced under normoxia. High‑throughput analysis of DEGs and ASGs demonstrated that genes in the complement pathway were enriched, which may be a potential mechanism underlying hypoxia inhibition of CESCs osteogenesis. The results of the present study may provide a basis for future mechanistic studies regarding gene expression levels and alternative splicing events during the hypoxia‑regulated inhibition of osteogenesis, which may be helpful in identifying targets for CEP degeneration therapy.

A saga of cancer epigenetics: linking epigenetics to alternative splicing

The discovery of an increasing number of alternative splicing events in the human genome highlighted that ∼94% of genes generate alternatively spliced transcripts that may produce different protein isoforms with diverse functions. It is now well known that several diseases are a direct and indirect consequence of aberrant splicing events in humans. In addition to the conventional mode of alternative splicing regulation by 'cis' RNA-binding sites and 'trans' RNA-binding proteins, recent literature provides enormous evidence for epigenetic regulation of alternative splicing. The epigenetic modifications may regulate alternative splicing by either influencing the transcription elongation rate of RNA polymerase II or by recruiting a specific splicing regulator via different chromatin adaptors. The epigenetic alterations and aberrant alternative splicing are known to be associated with various diseases individually, but this review discusses/highlights the latest literature on the role of epigenetic alterations in the regulation of alternative splicing and thereby cancer progression. This review also points out the need for further studies to understand the interplay between epigenetic modifications and aberrant alternative splicing in cancer progression.

Characterization of visceral and subcutaneous adipose tissue transcriptome in pregnant women with and without spontaneous labor at term: implication of alternative splicing in the metabolic adaptations of adipose tissue to parturition

OBJECTIVE

The aim of this study was to determine gene expression and splicing changes associated with parturition and regions (visceral vs. subcutaneous) of the adipose tissue of pregnant women.

STUDY DESIGN

The transcriptome of visceral and abdominal subcutaneous adipose tissue from pregnant women at term with (n=15) and without (n=25) spontaneous labor was profiled with the Affymetrix GeneChip Human Exon 1.0 ST array. Overall gene expression changes and the differential exon usage rate were compared between patient groups (unpaired analyses) and adipose tissue regions (paired analyses). Selected genes were tested by quantitative reverse transcription-polymerase chain reaction.

RESULTS

Four hundred and eighty-two genes were differentially expressed between visceral and subcutaneous fat of pregnant women with spontaneous labor at term (q-value <0.1; fold change >1.5). Biological processes enriched in this comparison included tissue and vasculature development as well as inflammatory and metabolic pathways. Differential splicing was found for 42 genes [q-value <0.1; differences in Finding Isoforms using Robust Multichip Analysis scores >2] between adipose tissue regions of women not in labor. Differential exon usage associated with parturition was found for three genes (LIMS1, HSPA5, and GSTK1) in subcutaneous tissues.

CONCLUSION

We show for the first time evidence of implication of mRNA splicing and processing machinery in the subcutaneous adipose tissue of women in labor compared to those without labor.

Global profiling of the gene expression and alternative splicing events during hypoxia-regulated chondrogenic differentiation in human cartilage endplate-derived stem cells

The intervertebral disc (IVD) degeneration is initiated by cartilage endplate (CEP) degeneration and is characterised by reduced chondrification. Cartilage endplate-derived stem cells (CESCs) with chondrogenic differentiation abilities are responsible for the restoration of cartilage. CEP remains in an avascular and hypoxic microenvironment. In this study, we observed that the physiological hypoxia greatly promotes the chondrogenic differentiation of CESCs. This tissue specificity of the differentiation fate of CESCs in response to the hypoxic microenvironment was physiologically significant for the CEP to maintain the chondrification status. To investigate the mechanisms underlying the hypoxia-regulated chondrogenic differentiation of CESCs, we adopted a high-throughput scanning technology to detect the global profiling of gene expression and alternative splicing (AS) event changes during chondrogenic differentiation under hypoxia in CESCs compared to those induced under normoxia. An Affymetrix Human Transcriptome Array 2.0 was used to identify the differentially expressed genes (DEGs) and alternatively spliced genes (ASGs). After RT-PCR validation, GO and KEGG pathway analyses of both the DEGs and ASGs were performed. The enrichment of the GO functional terms and signalling pathways provided referential direction of the mechanism to study the gene expression and AS in the hypoxia-regulated chondrogenesis promotion, which could be helpful in understanding this physiological phenomenon, and it could also be instrumental in finding targets for CEP degeneration therapy.

Characterization of Visceral and Subcutaneous Adipose Tissue Transcriptome and Biological Pathways in Pregnant and Non-Pregnant Women: Evidence for Pregnancy-Related Regional-Specific Differences in Adipose Tissue

Objective

The purpose of this study was to compare the transcriptome of visceral and subcutaneous adipose tissues between pregnant and non-pregnant women.

Study Design

The transcriptome of paired visceral and abdominal subcutaneous adipose tissues from pregnant women at term and matched non-pregnant women (n = 11) was profiled with the Affymetrix Human Exon 1.0 ST array. Differential expression of selected genes was validated with the use of quantitative reverse transcription–polymerase chain reaction.

Results

Six hundred forty-four transcripts from 633 known genes were differentially expressed (false discovery rate (FDR) <0.1; fold-change >1.5), while 42 exons from 36 genes showed differential usage (difference in FIRMA scores >2 and FDR<0.1) between the visceral and subcutaneous fat of pregnant women. Fifty-six known genes were differentially expressed between pregnant and non-pregnant subcutaneous fat and three genes in the visceral fat. Enriched biological processes in the subcutaneous adipose tissue of pregnant women were mostly related to inflammation.

Conclusion

The transcriptome of visceral and subcutaneous fat depots reveals pregnancy-related gene expression and splicing differences in both visceral and subcutaneous adipose tissue. Furthermore, for the first time, alternative splicing in adipose tissue has been associated with regional differences and human parturition.

The translational significance of epithelial-mesenchymal transition in head and neck cancer

Positive markers of epithelial-mesenchymal transition (EMT) in head and neck cancers complicate clinical management and are associated with reduced survival. We discuss recent translational discoveries in EMT and suggest additional actionable molecular pathways, biomarkers, and clinical agents.

Severe hypoxia exerts parallel and cell-specific regulation of gene expression and alternative splicing in human mesenchymal stem cells

BACKGROUND

The endosteum of the bone marrow provides a specialized hypoxic niche that may serve to preserve the integrity, pluripotency, longevity and stemness of resident mesenchymal stem cells (MSCs). To explore the molecular genetic consequences of such a niche we subjected human (h) MSCs to a pO2 of 4 mmHg and analyzed global gene expression and alternative splicing (AS) by genome-exon microarray and RT-qPCR, and phenotype by western blot and immunostaining.

RESULTS

Out of 446 genes differentially regulated by >2.5-fold, down-regulated genes outnumbered up-regulated genes by 243:203. Exon analyses revealed 60 hypoxia-regulated AS events with splice indices (SI) >1.0 from 53 genes and a correlation between high SI and degree of transcript regulation. Parallel analyses of a publicly available AS study on human umbilical vein endothelial cells (HUVECs) showed that there was a strong cell-specific component with only 11 genes commonly regulated in hMSCs and HUVECs and 17 common differentially spliced genes. Only 3 genes were differentially responsive to hypoxia at the gene (>2.0) and AS levels in both cell types. Functional assignments revealed unique profiles of gene expression with complex regulation of differentiation, extracellular matrix, intermediate filament and metabolic marker genes. Antioxidant genes, striated muscle genes and insulin/IGF-1 signaling intermediates were down-regulated. There was a coordinate induction of 9 out of 12 acidic keratins that along with other epithelial and cell adhesion markers implies a partial mesenchymal to epithelial transition.

CONCLUSIONS

We conclude that severe hypoxia confers a quiescent phenotype in hMSCs that is reflected by both the transcriptome profile and gene-specific changes of splicosome actions. The results reveal that severe hypoxia imposes markedly different patterns of gene regulation of MSCs compared with more moderate hypoxia. This is the first study to report hypoxia-regulation of AS in stem/progenitor cells and the first molecular genetic characterization of MSC in a hypoxia-induced quiescent immobile state.

Expression microarray analysis reveals alternative splicing of LAMA3 and DST genes in head and neck squamous cell carcinoma

Purpose

Prior studies have demonstrated tumor-specific alternative splicing events in various solid tumor types. The role of alternative splicing in the development and progression of head and neck squamous cell carcinoma (HNSCC) is unclear. Our study queried exon-level expression to implicate splice variants in HNSCC tumors.

Experimental Design

We performed a comparative genome-wide analysis of 44 HNSCC tumors and 25 uvulopalatopharyngoplasty (UPPP) tissue samples at an exon expression level. In our comparison we ranked genes based upon a novel score—the Maximum-Minimum Exon Score (MMES) – designed to predict the likelihood of an alternative splicing event occurring. We validated predicted alternative splicing events using quantitative RT-PCR on an independent cohort.

Results

After MMES scoring of 17,422 genes, the top 900 genes with the highest scores underwent additional manual inspection of expression patterns in a graphical analysis. The genes LAMA3, DST, VEGFC, SDHA, RASIP1, and TP63 were selected for further validation studies because of a high frequency of alternative splicing suggested in our graphical analysis, and literature review showing their biological relevance and known splicing patterns. We confirmed TP63 as having dominant expression of the short DeltaNp63 isoform in HNSCC tumor samples, consistent with prior reports. Two of the six genes (LAMA3 and DST) validated by quantitative RT-PCR for tumor-specific alternative splicing events (Student's t test, P<0.001).

Conclusion

Alternative splicing events of oncologically relevant proteins occur in HNSCC. The number of genes expressing tumor-specific splice variants needs further elucidation, as does the functional significance of selective isoform expression.

Epidermal growth-factor-induced transcript isoform variation drives mammary cell migration

Signal-induced transcript isoform variation (TIV) includes alternative promoter usage as well as alternative splicing and alternative polyadenylation of mRNA. To assess the phenotypic relevance of signal-induced TIV, we employed exon arrays and breast epithelial cells, which migrate in response to the epidermal growth factor (EGF). We show that EGF rapidly – within one hour – induces widespread TIV in a significant fraction of the transcriptome. Importantly, TIV characterizes many genes that display no differential expression upon stimulus. In addition, similar EGF-dependent changes are shared by a panel of mammary cell lines. A functional screen, which utilized isoform-specific siRNA oligonucleotides, indicated that several isoforms play essential, non-redundant roles in EGF-induced mammary cell migration. Taken together, our findings highlight the importance of TIV in the rapid evolvement of a phenotypic response to extracellular signals.

Oncogenic alternative splicing switches: role in cancer progression and prospects for therapy

Alterations in the abundance or activities of alternative splicing regulators generate alternatively spliced variants that contribute to multiple aspects of tumor establishment, progression and resistance to therapeutic treatments. Notably, many cancer-associated genes are regulated through alternative splicing suggesting a significant role of this post-transcriptional regulatory mechanism in the production of oncogenes and tumor suppressors. Thus, the study of alternative splicing in cancer might provide a better understanding of the malignant transformation and identify novel pathways that are uniquely relevant to tumorigenesis. Understanding the molecular underpinnings of cancer-associated alternative splicing isoforms will not only help to explain many fundamental hallmarks of cancer, but will also offer unprecedented opportunities to improve the efficacy of anti-cancer treatments.

AKT activation drives the nuclear localization of CSE1L and a pro-oncogenic transcriptional activation in ovarian cancer cells

The human homolog of the yeast cse1 gene (CSE1L) is over-expressed in ovarian cancer. CSE1L forms complex with Ran and importin-α and has roles in nucleocytoplasmic traffic and gene expression. CSE1L accumulated in the nucleus of ovarian cancer cell lines, while it was localized also in the cytoplasm of other cancer cell lines. Nuclear localization depended on AKT, which was constitutively active in ovarian cancer cells, as the CSE1L protein translocated to the cytoplasm when AKT was inactivated. Moreover, the expression of a constitutively active AKT forced the translocation of CSE1L from the cytoplasm to the nucleus in other cancer cells. Nuclear accrual of CSE1L was associated to the nuclear accumulation of the phosphorylated Ran Binding protein 3 (RanBP3), which depended on AKT as well. Also in samples of human ovarian cancer, AKT activation was associated to nuclear accumulation of CSE1L and phosphorylation of RanBP3. Expression profiling of ovarian cancer cells after CSE1L silencing showed that CSE1L was required for the expression of genes promoting invasion and metastasis. In agreement, CSE1L silencing impaired motility and invasiveness of ovarian cancer cells. Altogether these data show that in ovarian cancer cells activated AKT by affecting RanBP3 phosphorylation determines the nuclear accumulation of CSE1L and likely the nuclear concentration of transcription factors conveying pro-oncogenic signals.

Isoform level expression profiles provide better cancer signatures than gene level expression profiles

BACKGROUND

The majority of mammalian genes generate multiple transcript variants and protein isoforms through alternative transcription and/or alternative splicing, and the dynamic changes at the transcript/isoform level between non-oncogenic and cancer cells remain largely unexplored. We hypothesized that isoform level expression profiles would be better than gene level expression profiles at discriminating between non-oncogenic and cancer cellsgene level.

METHODS

We analyzed 160 Affymetrix exon-array datasets, comprising cell lines of non-oncogenic or oncogenic tissue origins. We obtained the transcript-level and gene level expression estimates, and used unsupervised and supervised clustering algorithms to study the profile similarity between the samples at both gene and isoform levels.

RESULTS

Hierarchical clustering, based on isoform level expressions, effectively grouped the non-oncogenic and oncogenic cell lines with a virtually perfect homogeneity-grouping rate (97.5%), regardless of the tissue origin of the cell lines. However, gene levelthis rate was much lower, being 75% at best based on the gene level expressions. Statistical analyses of the difference between cancer and non-oncogenic samples identified the existence of numerous genes with differentially expressed isoforms, which otherwise were not significant at the gene level. We also found that canonical pathways of protein ubiquitination, purine metabolism, and breast-cancer regulation by stathmin1 were significantly enriched among genes thatshow differential expression at isoform level but not at gene level.

CONCLUSIONS

In summary, cancer cell lines, regardless of their tissue of origin, can be effectively discriminated from non-cancer cell lines at isoform level, but not at gene level. This study suggests the existence of an isoform signature, rather than a gene signature, which could be used to distinguish cancer cells from normal cells.

Post-transcriptional regulation in cancer progression : Microenvironmental control of alternative splicing and translation

The microenvironment acts as a conduit for cellular communication, delivering signals that direct development and sustain tissue homeostasis. In pathologies such as cancer, this integral function of the microenvironment is hijacked to support tumor growth and progression. Cells sense the microenvironment via signal transduction pathways culminating in altered gene expression. In addition to induced transcriptional changes, the microenvironment exerts its effect on the cell through regulation of post-transcriptional processes including alternative splicing and translational control. Here we describe how alternative splicing and protein translation are controlled by microenvironmental parameters such as oxygen availability. We also emphasize how these pathways can be utilized to support processes that are hallmarks of cancer such as angiogenesis, proliferation, and cell migration. We stress that cancer cells respond to their microenvironment through an integrated regulation of gene expression at multiple levels that collectively contribute to disease progression.

Alternative transcription and alternative splicing in cancer

In recent years, the notion of "one gene makes one protein that functions in one signaling pathway" in mammalian cells has been shown to be overly simplistic. Recent genome-wide studies suggest that at least half of the human genes, including many therapeutic target genes, produce multiple protein isoforms through alternative splicing and alternative usage of transcription initiation and/or termination. For example, alternative splicing of the vascular endothelial growth factor gene (VEGFA) produces multiple protein isoforms, which display either pro-angiogenic or anti-angiogenic activities. Similarly, for the majority of human genes, the inclusion or exclusion of exonic sequences enhances the generation of transcript variants and/or protein isoforms that can vary in structure and functional properties. Many of the isoforms produced in this manner are tightly regulated during normal development but are misregulated in cancer cells. Altered expression of transcript variants and protein isoforms for numerous genes is linked with disease and its prognosis, and cancer cells manipulate regulatory mechanisms to express specific isoforms that confer drug resistance and survival advantages. Emerging insights indicate that modulating the expression of transcript and protein isoforms of a gene may hold the key to impeding tumor growth and act as a model for efficient targeting of disease-associated genes at the isoform level. This review highlights the role and regulation of alternative transcription and splicing mechanisms in generating the transcriptome, and the misuse and diagnostic/prognostic potential of alternative transcription and splicing in cancer.

Hypoxia-induced alternative splicing in endothelial cells

Background

Adaptation to low oxygen by changing gene expression is vitally important for cell survival and tissue development. The sprouting of new blood vessels, initiated from endothelial cells, restores the oxygen supply of ischemic tissues. In contrast to the transcriptional response induced by hypoxia, which is mainly mediated by members of the HIF family, there are only few studies investigating alternative splicing events. Therefore, we performed an exon array for the genome-wide analysis of hypoxia-related changes of alternative splicing in endothelial cells.

Methodology/Principal findings

Human umbilical vein endothelial cells (HUVECs) were incubated under hypoxic conditions (1% O₂) for 48 h. Genome-wide transcript and exon expression levels were assessed using the Affymetrix GeneChip Human Exon 1.0 ST Array. We found altered expression of 294 genes after hypoxia treatment. Upregulated genes are highly enriched in glucose metabolism and angiogenesis related processes, whereas downregulated genes are mainly connected to cell cycle and DNA repair. Thus, gene expression patterns recapitulate known adaptations to low oxygen supply. Alternative splicing events, until now not related to hypoxia, are shown for nine genes: six which are implicated in angiogenesis-mediated cytoskeleton remodeling (cask, itsn1, larp6, sptan1, tpm1 and robo1); one, which is involved in the synthesis of membrane-anchors (pign) and two universal regulators of gene expression (cugbp1 and max).

Conclusions/Significance

For the first time, this study investigates changes in splicing in the physiological response to hypoxia on a genome-wide scale. Nine alternative splicing events, until now not related to hypoxia, are reported, considerably expanding the information on splicing changes due to low oxygen supply. Therefore, this study provides further knowledge on hypoxia induced gene expression changes and presents new starting points to study the hypoxia adaptation of endothelial cells.

Using an exon microarray to identify a global profile of gene expression and alternative splicing in K562 cells exposed to sodium valproate

To investigate the effect of valproate treatment on the K562 cell line, a model for chronic myelogenous leukaemia, the growth and survival of the K562 cell line were investigated using the Annexin-V/PI dual staining method, and global profiles of gene expression and alternative splicing in K562 cells were assessed using exon microarrays. A significant increase in cell apoptosis was observed in valproate-exposed K562 cells using flow cytometry. A total of 628 transcripts were identified as being significantly differentially expressed. The number of genes demonstrating increased expression levels was greater than the number of genes demonstrating decreased expression levels (445 genes vs. 183 genes, respectively). The significant enrichment analysis of GO terms for the differentially expressed genes revealed that these genes are involved in many important biological processes such as apoptosis. Six of the genes observed to be differentially expressed that might be involved in apoptosis were selected to undergo qRT-PCR validation. In total, 198 candidates of alternative splicing variants were identified. Among them, three alternative splicing events were selected for validation, and CBLC and TBX1 were confirmed to be alternatively spliced by semi-nested PCR. In conclusion, valproate exposure facilitated cell apoptosis, altered mRNA expression and alternative splicing events in the K562 cell line.

Exon array data analysis using Affymetrix power tools and R statistical software

The use of microarray technology to measure gene expression on a genome-wide scale has been well established for more than a decade. Methods to process and analyse the vast quantity of expression data generated by a typical microarray experiment are similarly well-established. The Affymetrix Exon 1.0 ST array is a relatively new type of array, which has the capability to assess expression at the individual exon level. This allows a more comprehensive analysis of the transcriptome, and in particular enables the study of alternative splicing, a gene regulation mechanism important in both normal conditions and in diseases. Some aspects of exon array data analysis are shared with those for standard gene expression data but others present new challenges that have required development of novel tools. Here, I will introduce the exon array and present a detailed example tutorial for analysis of data generated using this platform.

Comprehensive exon array data processing method for quantitative analysis of alternative spliced variants

Alternative splicing of pre-mRNA generates protein diversity. Dysfunction of splicing machinery and expression of specific transcripts has been linked to cancer progression and drug response. Exon microarray technology enables genome-wide quantification of expression levels of the majority of exons and facilitates the discovery of alternative splicing events. Analysis of exon array data is more challenging than the analysis of gene expression data and there is a need for reliable quantification of exons and alternatively spliced variants. We introduce a novel, computationally efficient methodology, Multiple Exon Array Preprocessing (MEAP), for exon array data pre-processing, analysis and visualization. We compared MEAP with existing pre-processing methods, and validation of six exons and two alternatively spliced variants with qPCR corroborated MEAP expression estimates. Analysis of exon array data from head and neck squamous cell carcinoma (HNSCC) cell lines revealed several transcripts associated with 11q13 amplification, which is related with decreased survival and metastasis in HNSCC patients. Our results demonstrate that MEAP produces reliable expression values at exon, alternatively spliced variant and gene levels, which allows generating novel experimentally testable predictions.

Comprehensive analysis of alternative splicing and functionality in neuronal differentiation of P19 cells

BACKGROUND

Alternative splicing, which produces multiple mRNAs from a single gene, occurs in most human genes and contributes to protein diversity. Many alternative isoforms are expressed in a spatio-temporal manner, and function in diverse processes, including in the neural system.

METHODOLOGY/PRINCIPAL FINDINGS

The purpose of the present study was to comprehensively investigate neural-splicing using P19 cells. GeneChip Exon Array analysis was performed using total RNAs purified from cells during neuronal cell differentiation. To efficiently and readily extract the alternative exon candidates, 9 filtering conditions were prepared, yielding 262 candidate exons (236 genes). Semiquantitative RT-PCR results in 30 randomly selected candidates suggested that 87% of the candidates were differentially alternatively spliced in neuronal cells compared to undifferentiated cells. Gene ontology and pathway analyses suggested that many of the candidate genes were associated with neural events. Together with 66 genes whose functions in neural cells or organs were reported previously, 47 candidate genes were found to be linked to 189 events in the gene-level profile of neural differentiation. By text-mining for the alternative isoform, distinct functions of the isoforms of 9 candidate genes indicated by the result of Exon Array were confirmed.

CONCLUSIONS/SIGNIFICANCE

Alternative exons were successfully extracted. Results from the informatics analyses suggested that neural events were primarily governed by genes whose expression was increased and whose transcripts were differentially alternatively spliced in the neuronal cells. In addition to known functions in neural cells or organs, the uninvestigated alternative splicing events of 11 genes among 47 candidate genes suggested that cell cycle events are also potentially important. These genes may help researchers to differentiate the roles of alternative splicing in cell differentiation and cell proliferation.

Exon array analysis using re-defined probe sets results in reliable identification of alternatively spliced genes in non-small cell lung cancer

Background

Treatment of non-small cell lung cancer with novel targeted therapies is a major unmet clinical need. Alternative splicing is a mechanism which generates diverse protein products and is of functional relevance in cancer.

Results

In this study, a genome-wide analysis of the alteration of splicing patterns between lung cancer and normal lung tissue was performed. We generated an exon array data set derived from matched pairs of lung cancer and normal lung tissue including both the adenocarcinoma and the squamous cell carcinoma subtypes. An enhanced workflow was developed to reliably detect differential splicing in an exon array data set. In total, 330 genes were found to be differentially spliced in non-small cell lung cancer compared to normal lung tissue. Microarray findings were validated with independent laboratory methods for CLSTN1, FN1, KIAA1217, MYO18A, NCOR2, NUMB, SLK, SYNE2, TPM1, (in total, 10 events) and ADD3, which was analysed in depth. We achieved a high validation rate of 69%. Evidence was found that the activity of FOX2, the splicing factor shown to cause cancer-specific splicing patterns in breast and ovarian cancer, is not altered at the transcript level in several cancer types including lung cancer.

Conclusions

This study demonstrates how alternatively spliced genes can reliably be identified in a cancer data set. Our findings underline that key processes of cancer progression in NSCLC are affected by alternative splicing, which can be exploited in the search for novel targeted therapies.