Next Article in Journal
Chitosan Functionalization: Covalent and Non-Covalent Interactions and Their Characterization
Next Article in Special Issue
Numerical Study of Customized Artificial Cornea Shape by Hydrogel Biomaterials on Imaging and Wavefront Aberration
Previous Article in Journal
Synthesis and Characterization of Hybrid Fiber-Reinforced Polymer by Adding Ceramic Nanoparticles for Aeronautical Structural Applications
Previous Article in Special Issue
A Review of the Effects of Collagen Treatment in Clinical Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Collagen–Glycosaminoglycan Mesh on Gene Expression as Determined by Using Principal Component Analysis-Based Unsupervised Feature Extraction

1
Department of Physics, Chuo University, Tokyo 112-8551, Japan
2
Department of Computer Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(23), 4117; https://doi.org/10.3390/polym13234117
Submission received: 20 September 2021 / Revised: 18 November 2021 / Accepted: 23 November 2021 / Published: 26 November 2021
(This article belongs to the Special Issue Biomaterials in Medical Applications)

Abstract

:
The development of the medical applications for substances or materials that contact cells is important. Hence, it is necessary to elucidate how substances that surround cells affect gene expression during incubation. In the current study, we compared the gene expression profiles of cell lines that were in contact with collagen–glycosaminoglycan mesh and control cells. Principal component analysis-based unsupervised feature extraction was applied to identify genes with altered expression during incubation in the treated cell lines but not in the controls. The identified genes were enriched in various biological terms. Our method also outperformed a conventional methodology, namely, gene selection based on linear regression with time course.

1. Introduction

Several factors are known to affect cell division; one such effective factor is contact with solid materials (or substance) [1]. Regulating the cell division process using biomaterials is the central theme of tissue engineering. The effect of tissue engineering scaffolds is especially important because tissue engineering cannot be conducted without equipment that can store cell lines. Collagen–glycosaminoglycan mesh is one such important biomaterial because it is used to aid wound healing [2]. Although Klappericha and Bertozzi [3] once investigated the effect of collagen–glycosaminoglycan mesh on cell division cycles using microarray analysis, the small number of samples studied prevented them from identifying genes whose expression significantly varied during development and whose expression profiles were distinct between controls and treated cells. Although they selected genes associated with p-values of less than 0.001, considering the number of genes as 10 4 , it is far below significant.
The recently proposed principal component analysis (PCA)-based unsupervised feature extraction (FE) [4] has the ability to identify genes with expression profiles that are significantly different using a small number of samples. In this study, we successfully applied PCA-based unsupervised FE to determine gene expression profiles during the cell division of cells in control conditions and in contact with collagen–glycosaminoglycan mesh. The identified genes were found to be associated with several enrichment terms with considerable biological significance.

2. Materials and Methods

2.1. Gene Expression Profiles

Gene expression profiles were downloaded from the Gene Expression Omnibus (GEO) database (GEO ID: GSE6432). The dataset in GSE6432_series_matrix.txt.gz is available in the Series Matrix File(s) section. It consists of 32 gene expression profiles of the IMR90 cell lines, and the relevant details are provided in Table 1.

2.2. PCA-Based Unsupervised FE

Gene expression profiles are formatted as matrices x i j R 22,283 × 19 for treated cells and x i j R 22,283 × 13 for control cells, where x i j denotes the gene expression of the ith probe at the jth sample. Before applying singular value decomposition (SVD), they were standardized as
i x i j = 0
i x i j 2 = 22,283
After applying SVD, we obtained the left-hand singular value vector u i , which corresponded to the principal component score attributed to the probes, and the right-hand singular value vector v j , which corresponded to the principal component loadings attributed to the samples, if we interpreted the application of SVD as PCA.
In order to see which v j is coincident with time points, we applied linear regression as
v j = a + b t j ,
where a and b are regression coefficients and t j is the time point (hours in Table 1) associated with the jth sample. We used the lm function in R [5], and the obtained p-values were corrected using the Benjamini–Hochberg criterion [4]. v 3 j for treated cell is associated with the adjusted p-values less than 0.05, whereas v j s for control cell is not associated with adjusted p-values less than 0.05. This result is appropriate because the simple cell division process may not be associated with any time development other than cell senescence [6], which might not be detected in only 24 h.
Probes are selected by assuming that u 3 i , associated with v 3 j , obeys the Gaussian distribution (null hypothesis) by assigning p-values to the probes as
P i = P χ 2 > u 3 i σ 3 2 ,
where P χ 2 [ > x ] is the cumulative χ 2 distribution, the argument is larger than x, and σ 3 is the standard deviation. Thus, 324 probes associated with the adjusted p-value less than 0.01 were selected for the treated cell lines.

2.3. Gene Selection Using Linear Regression

As an alternative method to PCA-based unsupervised FE, we utilized linear regression-based FE. Linear regression is applied to x i j as
x i j = a i + b i t j ,
where a i and b i are regression coefficients and t j is the time point (hours in Table 1) associated with the jth sample. Subsequently, the adjusted p-values that were less than 0.01 were selected. The number of probes selected for treated cell lines was 813, and no probes were selected for the control cell lines.

2.4. Enrichment Analysis

The IDs of the selected probes were converted to gene symbol using the ID converter in DAVID [7]. Then, the gene symbols converted from the probe IDs were uploaded to Enrichr [8].

3. Results

As mentioned in the Materials and Methods Section, genes associated with the 318 probes for the treated cell lines (contact with collagen–glycosaminoglycan mesh) were uploaded to Enrichr (no probes were selected for control cell lines using this method). The full list of probes, genes, and enrichment analysis is provided in the Supplementary Materials (Data S1). Several enriched biological terms were determined.
The top ranked term in the GO biological process (BP) (Table 2) is “regulation of apoptotic process”. Na et al. reported [9] that collagen–glycosaminoglycan has an anti-apoptosis effect. Thus, the fact that this term is ranked first is reasonable.
“Focal adhesion” is the top ranked term in “GO Cellular Component 2021” (Table 3) and the nineth ranked in “KEGG 2021 Human” (Table 4); moreover, Murphy et al. [10] reported that the collagen–glycosaminoglycan scaffold plays critical roles in focal adhesion.
Other than these three categories, there are some additional categories that support the suitability of our analysis. For example, “ARCHS4 Cell-lines” lists IMR90, which is the cell line used in this study, as the top ranked cell line (Table 5).
Moreover, although it is not the top ranked term, “FETAL LUNG”, from which IMR90 cell lines were derived, is ranked within the top 10 ranked terms in “ARCHS4 Tissues” (Table 6).
Although we provide only a few examples, our results suggest that our analysis was robust.

4. Discussion

Although we successfully applied our methodology to the dataset, one might wonder whether more conventional methods can achieve similar performance. Since this dataset was generated using archaic technology, namely, microarray, more modernized methodologies adapted to high-throughput sequencing technology (e.g., edgeR [11] or DESeq2 [12]) cannot be employed. Moreover, the archaic technologies adapted to microarray (e.g., SAM [13] and limma [14]) cannot be employed, because they can only deal with categorical classification, whereas we need to identify genes whose expressions are altered as a numerical variable (hours). Thus, we decided to employ more conventional methodology than SAM or limma, namely, gene selection using linear regression.
As described in the Materials and Methods Section, we identified 813 probes using linear regression-based FE and uploaded the gene symbols associated with the identified probes to Enrichr. When considering only the number of probes selected, it performed better than the PCA-based unsupervised FE, which could only identify 324 probes. Selecting no probes for the control cell lines is the same as PCA-based unsupervised FE. Thus, it seems that the application of PCA-based unsupervised FE, instead of linear regression, was not productive.
Nevertheless, if we consider the performance of the enrichment analysis more carefully, this impression is reversed. A full list of the probes, genes, and the results of enrichment analysis are provided in the Supplementary Materials (Data S2). First, for “GO BP 2021”, in which PCA-based unsupervised FE ranked apoptosis first (Table 2 and Table 7), although the top ranked term “regulation of apoptotic process” in Table 2 is associated with the adjusted p-value as small as 2.52 × 10 13 , the top ranked term in Table 7 is associated with adjusted p-value as large as 4.56 × 10 2 , which is much less significant. Even the tenth ranked term in Table 2 is more significant than the top ranked term in Table 7. Generally, more genes uploaded have more opportunities to be associated with more significant enrichment. Nevertheless, genes associated with 813 probes, which were greater than the 324 probes identified using PCA-based unsupervised FE, could be associated with the less significant terms. This clearly suggests the inferiority of linear regression as compared to PCA-based unsupervised FE.
Regarding the comparison of the “GO Cellular Component 2021” in Table 3 and Table 8, we have a similar impression.
Although “focal adhesion” is ranked first in both Tables, its significance is very distinct. It is associated with an adjusted P-value as small as 4.77 × 10 27 in Table 3, whereas it is associated with that as large as 3.39 × 10 7 in Table 8. The number of overlapping genes is only 39 in Table 8, whereas it is higher (43) in Table 3, despite the fact that a higher total number of genes was uploaded to Enrichr, as shown in Table 8. Thus, the performance of linear regression is again poorer than that of PCA-based unsupervised FE.
For KEGG, not only are the generally adjusted p-values larger (i.e., less significant) in Table 9 than those in Table 4, but also “Glycolysis/Gluconeogenesis” and “Focal adhesion”, which are ranked within the top 10 in Table 4, are not even listed in Table 9, and no other terms seemingly related to the experiments are mentioned. Thus, the performance of linear regression is again poorer than that of PCA-based unsupervised FE.
For “ARCHS4 Cell-lines” and “ARCHS4 Tissue”, the results are similar. In Table 10, not only are the adjusted p-values generally larger (i.e., less significant) than those in Table 5, but the adjusted p-values attributed to IMR90 in Table 10 ( 1.06 × 10 5 ) are also much larger (i.e., less significant) than those in Table 5. The number of overlapping genes for IMR90 is only 128 in Table 5, whereas that in Table 10 is 89, despite the fact that more than twice the total number of genes were uploaded to Enrichr, as shown in Table 5. However, the number of overlapping genes for HUVEC, which is the wrong one, is as large as 113 in Table 10, whereas that in Table 5 is only 64. Thus, the increased number of genes selected using linear regression substantially contributes to the increase in overlapping genes assigned to the wrong answer. Moreover, lower ranked terms failed to demonstrate an association with significant p-values (e.g., less than 0.015). These finding suggest the inferiority of linear regression as compared to PCA-based unsupervised FE.
Although “FETAL LUNG” is fourth ranked in Table 11, its adjusted p-value is 1.05 × 10 3 , which is much less significant than that in Table 6 ( 9.58 × 10 9 ). Thus, overall, PCA-based unsupervised FE performed better than linear regression.
Finally, we attempted to conduct a time-series analysis, which is more widely used than linear regression for time course data. To this end, we used the fsMTS [15] package implemented in R [5] that included multiple methods, such as correlation-based, lasso-based, mutual information-based, and random forest-based methods. Nevertheless, none of the fsMTS methods could be performed. This was because time-series analysis requires auto/cross-correlations that require the memory size proportional to the square of the number of features. Since the number of features in this analysis was as high as 10 4 , it was computationally infeasible to execute the methods in fsMTS. Thus, our strategy, PCA-based unsupervised FE, was the only one applicable to the present data set.
The limitation of our methodology is that because of its unsupervised nature, when it fails to select biologically reasonable genes, there are no ways to improve it, although it occasionally worked effectively in the present study.

5. Conclusions

In the current study, we applied PCA-based unsupervised FE to gene expression profiles for IMR90 cell lines incubated in collagen–glycosaminoglycan mesh. Whereas no genes whose expressions vary over time were detected in control cell lines, the expression profiles of several genes were altered during the cell division process. These genes are associated with several enriched biological terms. One conventional method, linear regression, was employed for comparison. Although it could select several hundred genes whose expressions vary over time, their enrichment was inferior to that seen using PCA-based unsupervised FE. Thus, not only can PCA-based unsupervised FE achieve a good performance, but it can also outperform a conventional method. We demonstrated that collagen–glycosaminoglycan is an effective medium that could be used for cell culture.

Supplementary Materials

The following are available at https://www.mdpi.com/article/10.3390/polym13234117/s1, Data S1: The full list of probes, genes, and enrichment analysis is provided in the supplementary material for PCA based unsupervised FE, Data S2: The full list of probes, genes, and enrichment analysis is provided in the supplementary material for linear regression.

Author Contributions

Y.-h.T. planned the research and performed analyses. Y.-h.T. and T.T. evaluated the results, discussions, outcomes and wrote and reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by KAKENHI (grant numbers 19H05270, 20H04848, and 20K12067) to Y.-h.T.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used in this study are available in GEO ID GSE6432.

Conflicts of Interest

The authors declare no conflict of interest. The founders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. McMurray, R.J.; Dalby, M.J.; Tsimbouri, P.M. Using biomaterials to study stem cell mechanotransduction, growth and differentiation. J. Tissue Eng. Regen. Med. 2015, 9, 528–539. [Google Scholar] [CrossRef] [PubMed]
  2. Ansari, M.; Kordestani, S.S.; Nazralizadeh, S.; Eslami, H. Biodegradable Cell-Seeded Collagen Based Polymer Scaffolds for Wound Healing and Skin Reconstruction. J. Macromol. Sci. Part B 2018, 57, 100–109. [Google Scholar] [CrossRef]
  3. Klapperich, C.M.; Bertozzi, C.R. Global gene expression of cells attached to a tissue engineering scaffold. Biomaterials 2004, 25, 5631–5641. [Google Scholar] [CrossRef] [PubMed]
  4. Taguchi, Y.-h. Unsupervised Feature Extraction Applied to Bioinformatics; Springer International Publishing: Berlin/Heidelberg, Germany, 2020. [Google Scholar] [CrossRef]
  5. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  6. Micco, R.D.; Krizhanovsky, V.; Baker, D.; d’Adda di Fagagna, F. Cellular senescence in ageing: From mechanisms to therapeutic opportunities. Nat. Rev. Mol. Cell Biol. 2020, 22, 75–95. [Google Scholar] [CrossRef] [PubMed]
  7. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2008, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  8. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Na, J.; Bak, D.; Im, S.I.; Choi, H.; Hwang, J.H.; Kong, S.Y.; No, Y.A.; Lee, Y.; Kim, B.J. Anti-apoptotic effects of glycosaminoglycans via inhibition of ERK/AP-1 signaling in TNF-α-stimulated human dermal fibroblasts. Int. J. Mol. Med. 2018. [Google Scholar] [CrossRef] [PubMed]
  10. Murphy, C.M.; Matsiko, A.; Haugh, M.G.; Gleeson, J.P.; O’Brien, F.J. Mesenchymal stem cell fate is regulated by the composition and mechanical properties of collagen-glycosaminoglycan scaffolds. J. Mech. Behav. Biomed. Mater. 2012, 11, 53–62. [Google Scholar] [CrossRef] [PubMed]
  11. McCarthy, D.J.; Chen, Y.; Smyth, G.K. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012, 40, 4288–4297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 1–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Tusher, V.G.; Tibshirani, R.; Chu, G. Significance analysis of microarrays applied to the ionizing radiation response. Proc. Natl. Acad. Sci. USA 2001, 98, 5116–5121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  15. Pavlyuk, D. fsMTS: Feature Selection for Multivariate Time Series; R package Version 0.1.5. 2020. Available online: https://cran.r-project.org/web/packages/fsMTS/fsMTS.pdf (accessed on 6 April 2020).
Table 1. The number of samples with the gene expression profiles. Treated means contact with the collagen–glycosaminoglycan mesh.
Table 1. The number of samples with the gene expression profiles. Treated means contact with the collagen–glycosaminoglycan mesh.
Conditions1 h2 h4 h8 h12 h24 h48 hTotal
treated323432219
control123231113
Table 2. The top 10 enriched terms in “GO Biological Process 2021” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
Table 2. The top 10 enriched terms in “GO Biological Process 2021” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
TermOverlapp-ValueAdjusted p-Value
regulation of apoptotic process (GO:0042981)41/742 1.18 × 10 16 2.52 × 10 13
SRP-dependent cotranslational protein targeting to membrane (GO:0006614)17/90 4.06 × 10 16 4.33 × 10 13
cotranslational protein targeting to membrane (GO:0006613)17/94 8.78 × 10 16 6.24 × 10 13
protein targeting to ER (GO:0045047)17/103 4.36 × 10 15 2.33 × 10 12
cytoplasmic translation (GO:0002181)16/93 1.44 × 10 14 6.15 × 10 12
nuclear-transcribed mRNA catabolic process, nonsense-mediated decay (GO:0000184)17/113 2.16 × 10 14 7.70 × 10 12
cellular protein metabolic process (GO:0044267)27/417 7.60 × 10 13 2.10 × 10 10
peptide biosynthetic process (GO:0043043)18/162 7.89 × 10 13 2.10 × 10 10
negative regulation of programmed cell death (GO:0043069)25/381 4.24 × 10 12 1.00 × 10 9
nuclear-transcribed mRNA catabolic process (GO:0000956)17/171 2.11 × 10 11 4.51 × 10 9
Table 3. The top 10 enriched terms in “GO Cellular Component 2021” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
Table 3. The top 10 enriched terms in “GO Cellular Component 2021” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
TermOverlapp-ValueAdjusted p-Value
focal adhesion (GO:0005925)43/387 2.22 × 10 29 4.77 × 10 27
cell-substrate junction (GO:0030055)43/394 4.69 × 10 29 5.04 × 10 27
intracellular organelle lumen (GO:0070013)40/848 5.63 × 10 14 4.03 × 10 12
collagen-containing extracellular matrix (GO:0062023)25/380 4.00 × 10 12 2.15 × 10 10
endoplasmic reticulum lumen (GO:0005788)21/285 2.78 × 10 11 1.19 × 10 9
cytosolic large ribosomal subunit (GO:0022625)10/55 7.95 × 10 10 2.85 × 10 8
large ribosomal subunit (GO:0015934)10/59 1.64 × 10 9 5.03 × 10 8
ribosome (GO:0005840)10/62 2.72 × 10 9 7.30 × 10 8
secretory granule lumen (GO:0034774)19/316 7.48 × 10 9 1.79 × 10 7
ficolin-1-rich granule lumen (GO:1904813)12/123 2.50 × 10 8 5.37 × 10 7
Table 4. The top 10 enriched terms in “KEGG 2021 Human” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
Table 4. The top 10 enriched terms in “KEGG 2021 Human” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
TermOverlapp-ValueAdjusted p-Value
Coronavirus disease21/232 5.36 × 10 13 1.22 × 10 10
Ribosome17/158 5.88 × 10 12 6.67 × 10 10
Legionellosis9/57 2.09 × 10 8 1.58 × 10 6
Salmonella infection16/249 4.75 × 10 8 2.70 × 10 6
IL-17 signaling pathway10/94 1.64 × 10 7 7.46 × 10 6
Glycolysis/Gluconeogenesis8/67 1.20 × 10 6 4.52 × 10 5
Lipid and atherosclerosis13/215 1.75 × 10 6 5.69 × 10 5
Protein digestion and absorption9/103 3.65 × 10 6 1.04 × 10 4
Focal adhesion12/201 5.03 × 10 6 1.17 × 10 4
PI3K-Akt signaling pathway16/354 5.14 × 10 6 1.17 × 10 4
Table 5. The top 10 enriched terms in “ARCHS4 Cell-lines” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
Table 5. The top 10 enriched terms in “ARCHS4 Cell-lines” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
TermOverlapp-ValueAdjusted p-Value
IMR9089/2395 1.56 × 10 24 1.95 × 10 22
NHDF79/2395 2.70 × 10 18 1.69 × 10 16
BJ CELL72/2395 2.02 × 10 14 8.41 × 10 1 3
HUVEC64/2395 1.61 × 10 10 5.04 × 10 9
T2462/2395 1.23 × 10 9 3.09 × 10 8
T98G59/2395 2.22 × 10 8 4.63 × 10 7
BT54956/2395 3.28 × 10 7 5.12 × 10 6
DU14556/2395 3.28 × 10 7 5.12 × 10 6
CAKI155/2395 7.68 × 10 7 9.60 × 10 6
U8755/2395 7.68 × 10 7 9.60 × 10 6
Table 6. The top 10 enriched terms in “ARCHS4 Tissues” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
Table 6. The top 10 enriched terms in “ARCHS4 Tissues” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using PCA-based unsupervised FE.
TermOverlapp-ValueAdjusted p-Value
MYOBLAST87/2316 3.16 × 10 24 3.38 × 10 22
FIBROBLAST86/2316 1.45 × 10 23 7.76 × 10 22
FORESKIN FIBROBLAST72/2316 3.58 × 10 15 1.28 × 10 13
BLOOD DENDRITIC CELLS66/2316 4.31 × 10 12 1.15 × 10 10
DENDRITIC CELL63/2316 1.13 × 10 10 2.02 × 10 9
OSTEOBLAST63/2316 1.13 × 10 10 2.02 × 10 9
KUPFFER CELL62/2316 3.22 × 10 10 4.30 × 10 9
VASCULAR SMOOTH MUSCLE62/2316 3.22 × 10 10 4.30 × 10 9
FETAL LUNG61/2316 8.95 × 10 10 9.58 × 10 9
LIVER (BULK TISSUE)61/2316 8.95 × 10 10 9.58 × 10 9
Table 7. The top 10 enriched terms in “GO Biological Process 2021” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
Table 7. The top 10 enriched terms in “GO Biological Process 2021” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
TermOverlapp-ValueAdjusted p-Value
actin polymerization or depolymerization (GO:0008154)9/50 4.05 × 10 5 4.56 × 10 2
rRNA-containing ribonucleoprotein complex export from nucleus (GO:0071428)4/7 4.22 × 10 5 4.56 × 10 2
positive regulation of protein modification process (GO:0031401)20/214 4.32 × 10 5 4.56 × 10 2
transmembrane receptor protein tyrosine kinase signaling pathway (GO:0007169)30/404 5.44 × 10 5 4.56 × 10 2
protein stabilization (GO:0050821)17/179 1.31 × 10 4 7.08 × 10 2
regulation of lipid biosynthetic process (GO:0046890)7/35 1.46 × 10 4 7.08 × 10 2
regulation of cellular metabolic process (GO:0031323)8/47 1.63 × 10 4 7.08 × 10 2
positive regulation of cellular protein metabolic process (GO:0032270)12/102 1.74 × 10 4 7.08 × 10 2
regulation of stress fiber assembly (GO:0051492)10/74 1.90 × 10 4 7.08 × 10 2
regulation of mRNA catabolic process (GO:0061013)13/122 2.58 × 10 4 7.65 × 10 2
Table 8. The top 10 enriched terms in “GO Cellular Component 2021” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
Table 8. The top 10 enriched terms in “GO Cellular Component 2021” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
TermOverlapp-ValueAdjusted p-Value
focal adhesion (GO:0005925)39/387 1.41 × 10 9 3.39 × 10 7
cell-substrate junction (GO:0030055)39/394 2.34 × 10 9 3.39 × 10 7
actin cytoskeleton (GO:0015629)25/316 8.27 × 10 5 7.97 × 10 3
intracellular organelle lumen (GO:0070013)49/848 2.02 × 10 4 1.46 × 10 2
nucleus (GO:0005634)186/4484 1.06 × 10 3 5.51 × 10 2
collagen-containing extracellular matrix (GO:0062023)25/380 1.30 × 10 3 5.51 × 10 2
cytoplasmic stress granule (GO:0010494)8/65 1.54 × 10 3 5.51 × 10 2
intracellular membrane-bounded organelle (GO:0043231)210/5192 1.66 × 10 3 5.51 × 10 2
endoplasmic reticulum lumen (GO:0005788)20/285 1.81 × 10 3 5.51 × 10 2
intracellular non-membrane-bounded organelle (GO:0043232)58/1158 1.91 × 10 3 5.51 × 10 2
Table 9. The top 10 enriched terms in “KEGG 2021 Human” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
Table 9. The top 10 enriched terms in “KEGG 2021 Human” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
TermOverlapp-ValueAdjusted p-Value
PI3K-Akt signaling pathway28/354 3.17 × 10 5 9.30 × 10 3
Sphingolipid signaling pathway13/119 2.01 × 10 4 1.66 × 10 2
Arrhythmogenic right ventricular cardiomyopathy10/77 2.65 × 10 4 1.66 × 10 2
Antigen processing and presentation10/78 2.95 × 10 4 1.66 × 10 2
Hepatitis C15/157 2.99 × 10 4 1.66 × 10 2
Salmonella infection20/249 3.39 × 10 4 1.66 × 10 2
Hippo signaling pathway15/163 4.47 × 10 4 1.87 × 10 2
Tight junction15/169 6.54 × 10 4 2.29 × 10 2
Protein processing in endoplasmic reticulum15/171 7.39 × 10 4 2.29 × 10 2
AMPK signaling pathway12/120 7.81 × 10 4 2.29 × 10 2
Table 10. The top 10 enriched terms in “ARCHS4 Cell-lines” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
Table 10. The top 10 enriched terms in “ARCHS4 Cell-lines” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
TermOverlapp-ValueAdjusted p-Value
IMR90128/2395 8.49 × 10 8 1.06 × 10 5
HUVEC113/2395 1.54 × 10 4 9.65 × 10 3
NHDF112/2395 2.35 × 10 4 9.78 × 10 3
BT549103/2395 6.32 × 10 3 1.62 × 10 1
BJ CELL101/2395 1.17 × 10 2 1.62 × 10 1
HNSCC101/2395 1.17 × 10 2 1.62 × 10 1
KNS42101/2395 1.17 × 10 2 1.62 × 10 1
NHBE101/2395 1.17 × 10 2 1.62 × 10 1
U87101/2395 1.17 × 10 2 1.62 × 10 1
DAOY99/2395 2.07 × 10 2 2.35 × 10 1
Table 11. The top 10 enriched terms in “ARCHS4 Tissues” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
Table 11. The top 10 enriched terms in “ARCHS4 Tissues” using Enrichr. Overlap is the number of common genes between the genes uploaded and the genes in the category divided by the number of genes in the category. Probes, whose associated genes were uploaded to Enrichr, were identified using linear regression.
TermOverlapp-ValueAdjusted p-Value
FIBROBLAST116/2316 9.27 × 10 6 5.01 × 10 4
VENTRICLE116/2316 9.27 × 10 6 5.01 × 10 4
ADIPOSE (BULK TISSUE)113/2316 3.89 × 10 5 1.05 × 10 3
FETAL LUNG113/2316 3.89 × 10 5 1.05 × 10 3
RESPIRATORY SMOOTH MUSCLE112/2316 6.14 × 10 5 1.33 × 10 3
OMENTUM111/2316 9.58 × 10 5 1.73 × 10 3
SUBCUTANEOUS ADIPOSE TISSUE106/2316 7.60 × 10 4 1.17 × 10 2
MYOBLAST104/2316 1.61 × 10 3 2.18 × 10 2
NEURONAL EPITHELIUM103/2316 2.31 × 10 3 2.77 × 10 2
ASTROCYTE99/2316 8.75 × 10 3 8.54 × 10 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Taguchi, Y.-h.; Turki, T. Effects of Collagen–Glycosaminoglycan Mesh on Gene Expression as Determined by Using Principal Component Analysis-Based Unsupervised Feature Extraction. Polymers 2021, 13, 4117. https://doi.org/10.3390/polym13234117

AMA Style

Taguchi Y-h, Turki T. Effects of Collagen–Glycosaminoglycan Mesh on Gene Expression as Determined by Using Principal Component Analysis-Based Unsupervised Feature Extraction. Polymers. 2021; 13(23):4117. https://doi.org/10.3390/polym13234117

Chicago/Turabian Style

Taguchi, Y-h., and Turki Turki. 2021. "Effects of Collagen–Glycosaminoglycan Mesh on Gene Expression as Determined by Using Principal Component Analysis-Based Unsupervised Feature Extraction" Polymers 13, no. 23: 4117. https://doi.org/10.3390/polym13234117

APA Style

Taguchi, Y. -h., & Turki, T. (2021). Effects of Collagen–Glycosaminoglycan Mesh on Gene Expression as Determined by Using Principal Component Analysis-Based Unsupervised Feature Extraction. Polymers, 13(23), 4117. https://doi.org/10.3390/polym13234117

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop