Tissue specific human fibroblast differential expression based on RNAsequencing analysis

Physical forces, such as mechanical stress, are essential for tissue homeostasis and influence gene expression of cells. In particular, the fibroblast has demonstrated sensitivity to extracellular matrices with assumed adaptation upon various mechanical loads. The purpose of this study was to compare the vocal fold fibroblast genotype, known for its unique mechanically stressful tissue environment, with cellular counterparts at various other anatomic locales to identify differences in functional gene expression profiles.

By using RNA-seq technology, we identified differentially expressed gene programs (DEseq2) among seven normal human fibroblast primary cell lines from healthy cadavers, which included: vocal fold, trachea, lung, abdomen, scalp, upper gingiva, and soft palate. Unsupervised gene expression analysis yielded 6216 genes differentially expressed across all anatomic sites. Hierarchical cluster analysis revealed grouping based on anatomic site origin rather than donor, suggesting global fibroblast phenotype heterogeneity. Sex and age-related effects were negligible. Functional enrichment analyses based on separate post-hoc 2-group comparisons revealed several functional themes within the vocal fold fibroblast related to transcription factors for signaling pathways regulating pluripotency of stem cells and extracellular matrix components such as cell signaling, migration, proliferation, and differentiation potential.

Human fibroblasts display a phenomenon of global topographic differentiation, which is maintained in isolation via in vitro assays. Epigenetic mechanical influences on vocal fold tissue may play a role in uniquely modelling and maintaining the local environmental cellular niche during homeostasis with vocal fold fibroblasts distinctly specialized related to their anatomic positional and developmental origins established during embryogenesis.