MXRA5 is a TGF-β1-regulated human protein with anti-inflammatory and anti-fibrotic properties

Identification of HTRA1, DPT and MXRA5 as potential biomarkers associated with osteoarthritis progression and immune infiltration

BACKGROUND

Our study aimed to identify potential specific biomarkers for osteoarthritis (OA) and assess their relationship with immune infiltration.

METHODS

We utilized data from GSE117999, GSE51588, and GSE57218 as training sets, while GSE114007 served as a validation set, all obtained from the GEO database. First, weighted gene co-expression network analysis (WGCNA) and functional enrichment analysis were performed to identify hub modules and potential functions of genes. We subsequently screened for potential OA biomarkers within the differentially expressed genes (DEGs) of the hub module using machine learning methods. The diagnostic accuracy of the candidate genes was validated. Additionally, single gene analysis and ssGSEA was performed. Then, we explored the relationship between biomarkers and immune cells. Lastly, we employed RT-PCR to validate our results.

RESULTS

WGCNA results suggested that the blue module was the most associated with OA and was functionally associated with extracellular matrix (ECM)-related terms. Our analysis identified ALB, HTRA1, DPT, MXRA5, CILP, MPO, and PLAT as potential biomarkers. Notably, HTRA1, DPT, and MXRA5 consistently exhibited increased expression in OA across both training and validation cohorts, demonstrating robust diagnostic potential. The ssGSEA results revealed that abnormal infiltration of DCs, NK cells, Tfh, Th2, and Treg cells might contribute to OA progression. HTRA1, DPT, and MXRA5 showed significant correlation with immune cell infiltration. The RT-PCR results also confirmed these findings.

CONCLUSIONS

HTRA1, DPT, and MXRA5 are promising biomarkers for OA. Their overexpression strongly correlates with OA progression and immune cell infiltration.

MXRA7 is involved in monocyte-to-macrophage differentiation

Macrophages are critical in mediating immune and inflammatory responses, while monocyte-to-macrophage differentiation is one of the main macrophage resources that involves various matrix proteins. Matrix remodeling associated 7 (MXRA7) was recently discovered to affect a variety of physiological and pathological processes related to matrix biology. In the present study, we investigated the role of MXRA7 in monocyte-to-macrophage differentiation in vitro. We found that knockdown of MXRA7 inhibited the proliferation of THP-1 human monocytic cells. Knockdown of MXRA7 increased the adhesion ability of THP-1 cells through upregulation the expression of adhesion molecules VCAM-1 and ICAM1. Knockdown of MXRA7 alone could promoted the differentiation of THP-1 cells to macrophages. Furthermore, the MXRA7-knockdown THP-1 cells produced a more significant upregulation pattern with M1-type cytokines (TNF-α, IL-1β and IL-6) than with those M2-type molecules (TGF-β1 and IL-1RA) upon PMA stimulation, indicating that knockdown of MXRA7 facilitated THP-1 cells differentiation toward M1 macrophages. RNA sequencing analysis revealed the potential biological roles of MXRA7 in cell adhesion, macrophage and monocyte differentiation. Moreover, MXRA7 knockdown promoted the expression of NF-κB p52/p100, while PMA stimulation could increase the expression of NF-κB p52/p100 and activating MAPK signaling pathways in MXRA7 knockdown cells. In conclusion, MXRA7 affected the differentiation of THP-1 cells toward macrophages possibly through NF-κB signaling pathways.

Proteomics of urinary small extracellular vesicles in early diagnosis of kidney diseases in children-expectations and limitations

The primary function of the kidneys is to maintain systemic homeostasis (disruption of renal structure and function results in multilevel impairment of body function). Kidney diseases are characterized by a chronic, progressive course and may result in the development of chronic kidney disease (CKD). Evaluation of the composition of the proteome of urinary small extracellular vesicles (sEVs) as a so-called liquid biopsy is a promising new research direction. Knowing the composition of sEV could allow localization of cellular changes in specific sections of the nephron or the interstitial tissue before fixed changes, detectable only at an advanced stage of the disease, occur. Research is currently underway on the role of sEVs in the diagnosis and monitoring of many disease entities. Reports in the literature on the subject include: diabetic nephropathy, focal glomerulosclerosis in the course of glomerulopathies, renal fibrosis of various etiologies. Studies on pediatric patients are still few, involving piloting if small groups of patients without validation studies. Here, we review the literature addressing the use of sEV for diagnosis of the most common urinary disorders in children. We evaluate the clinical utility and define limitations of markers present in sEV as potential liquid biopsy.

Senescence-Associated Alterations in Matrisome of Mesenchymal Stem Cells

The process of aging is intimately linked to alterations at the tissue and cellular levels. Currently, the role of senescent cells in the tissue microenvironment is still being investigated. Despite common characteristics, different cell populations undergo distinctive morphofunctional changes during senescence. Mesenchymal stem cells (MSCs) play a pivotal role in maintaining tissue homeostasis. A multitude of studies have examined alterations in the cytokine profile that determine their regulatory function. The extracellular matrix (ECM) of MSCs is a less studied aspect of their biology. It has been shown to modulate the activity of neighboring cells. Therefore, investigating age-related changes in the MSC matrisome is crucial for understanding the mechanisms of tissue niche ageing. This study conducted a broad proteomic analysis of the matrisome of separated fractions of senescent MSCs, including the ECM, conditioned medium (CM), and cell lysate. This is the first time such an analysis has been conducted. It has been established that there is a shift in production towards regulatory molecules and a significant downregulation of the main structural and adhesion proteins of the ECM, particularly collagens, fibulins, and fibrilins. Additionally, a decrease in the levels of cathepsins, galectins, S100 proteins, and other proteins with cytoprotective, anti-inflammatory, and antifibrotic properties has been observed. However, the level of inflammatory proteins and regulators of profibrotic pathways increases. Additionally, there is an upregulation of proteins that can directly cause prosenescent effects on microenvironmental cells (SERPINE1, THBS1, and GDF15). These changes confirm that senescent MSCs can have a negative impact on other cells in the tissue niche, not only through cytokine signals but also through the remodeled ECM.

Construction and evaluation of a metabolic correlation diagnostic model for diabetes based on machine learning algorithms

BACKGROUND

Diabetes mellitus (DM) is a prevalent chronic disease marked by significant metabolic dysfunctions. Understanding its molecular mechanisms is vital for early diagnosis and treatment strategies.

METHODS

We used datasets GSE7014, GSE25724, and GSE156248 from the GEO database to build a diagnostic model for DM using Random Forest (RF) and LASSO regression models. GSE20966 served as a validation cohort. DM patients were classified into two subtypes for functional enrichment analysis. Expression levels of key diagnostic genes were validated using quantitative real-time PCR (qRT-PCR) on Peripheral Blood Mononuclear Cells (PBMCs) from DM patients and healthy controls, focusing on CXCL12 and PPP1R12B with GAPDH as the internal control.

RESULTS

After de-batching the datasets, we identified 131 differentially expressed genes (DEGs) between DM and control groups, with 70 up-regulated and 61 down-regulated. Enrichment analysis revealed significant down-regulation in the IL-12 signaling pathway, JAK signaling post-IL-12 stimulation, and the ferroptosis pathway in DM. Five genes (CXCL12, MXRA5, UCHL1, PPP1R12B, and C7) were identified as having diagnostic value. The diagnostic model showed high accuracy in both the training and validation cohorts. The gene set also enabled the subclassification of DM patients into groups with distinct functional traits. qRT-PCR results confirmed the bioinformatics findings, particularly the up-regulation of CXCL12 and PPP1R12B in DM patients.

CONCLUSION

Our study pinpointed seven energy metabolism-related genes differentially expressed in DM and controls, with five holding diagnostic value. Our model accurately diagnosed DM and facilitated patient subclassification, offering new insights into DM pathogenesis.

Mitochondrial reactive oxygen species impact human fibroblast responses to protracted γ-ray exposures

Purpose: Continuous exposure to ionizing radiation at a low dose rate poses significant health risks to humans on deep space missions, prompting the need for mechanistic studies to identify countermeasures against its deleterious effects. Mitochondria are a major subcellular locus of radiogenic injury, and may trigger secondary cellular responses through the production of reactive oxygen species (mtROS) with broader biological implications. Methods and Materials: To determine the contribution of mtROS to radiation-induced cellular responses, we investigated the impacts of protracted γ-ray exposures (IR; 1.1 Gy delivered at 0.16 mGy/min continuously over 5 days) on mitochondrial function, gene expression, and the protein secretome of human HCA2-hTERT fibroblasts in the presence and absence of a mitochondria-specific antioxidant mitoTEMPO (MT; 5 µM). Results: IR increased fibroblast mitochondrial oxygen consumption (JO2) and H2O2 release rates (JH2O2) under energized conditions, which corresponded to higher protein expression of NADPH Oxidase (NOX) 1, NOX4, and nuclear DNA-encoded subunits of respiratory chain Complexes I and III, but depleted mtDNA transcripts encoding subunits of the same complexes. This was associated with activation of gene programs related to DNA repair, oxidative stress, and protein ubiquination, all of which were attenuated by MT treatment along with radiation-induced increases in JO2 and JH2O2. IR also increased secreted levels of interleukin-8 and Type I collagens, while decreasing Type VI collagens and enzymes that coordinate assembly and remodeling of the extracellular matrix. MT treatment attenuated many of these effects while augmenting others, revealing complex effects of mtROS in fibroblast responses to IR. Conclusion: These results implicate mtROS production in fibroblast responses to protracted radiation exposure, and suggest potentially protective effects of mitochondrial-targeted antioxidants against radiogenic tissue injury in vivo.

Human dental pulp stem cells derived extracellular matrix promotes mineralization via Hippo and Wnt pathways

Extracellular matrix (ECM) is an intricate structure providing the microenvironment niche that influences stem cell differentiation. This study aimed to investigate the efficacy of decellularized ECM derived from human dental pulp stem cells (dECM_DPSCs) and gingival-derived mesenchymal stem cells (dECM_GSCs) as an inductive scaffold for osteogenic differentiation of GSCs. The proteomic analysis demonstrated that common and signature matrisome proteins from dECM_DPSCs and dECM_GSCs were related to osteogenesis/osteogenic differentiation. RNA sequencing data from GSCs reseeded on dECM_DPSCs revealed that dECM_DPSCs upregulated genes related to the Hippo and Wnt signaling pathways in GSCs. In the inhibitor experiments, results revealed that dECM_DPSCs superiorly promoted GSCs osteogenic differentiation, mainly mediated through Hippo and Wnt signaling. The present study emphasizes the promising translational application of dECM_DPSCs as a bio-scaffold rich in favorable regenerative microenvironment for tissue engineering.

Single-cell analysis reveals one cancer-associated fibroblasts subtype linked to metastasis in breast cancer: MXRA5 as a potential novel marker for prognosis

Cancer-associated fibroblasts (CAFs) are prevalent in the tumor microenvironment of breast cancer, comprising a group of cell subpopulations with spatial, phenotypic, and functional heterogeneity. Due to the lack of specific markers for CAF subpopulations, their specific mechanisms in breast cancer remain unclear. We identified eight distinct CAF phenotypes in breast cancer using multiple single-cell RNA sequencing datasets and determined distinct transcription factors (TFs) of CAFs through SCENIC analysis. Our study highlights one CAF subtype in breast cancer, FN1+CAF2, associated with metastasis and macrophage polarization. We observed elevated FN1 expression in the stromal tissue of breast cancer patients. Furthermore, FN1 knockdown in CAFs reduced the migration ability of breast cancer cells. We identified a regulatory gene, MXRA5, in CAF2, which may play crucial roles in breast cancer. Our results indicated upregulated MXRA5 expression in breast cancer tissues and CAFs from patients with lymph node metastasis in the following experiment. Overall, our study reveals that the FN1+CAF2 subtype is associated with metastasis and suggests that MXRA5 may be a novel marker mediating the effects of CAF2 on breast cancer metastasis. This study enriches our understanding of CAF heterogeneity and offers new insights for treating breast cancer metastasis.

Optimum Doses and Forms of Selenium Maintaining Reproductive Health via Regulating Homeostasis of Gut Microbiota and Testicular Redox, Inflammation, Cell Proliferation, and Apoptosis in Roosters

BACKGROUND

There is a U-shaped relationship between dietary selenium (Se) ingestion and optimal sperm quality.

OBJECTIVES

This study aimed to investigate the optimal dietary dose and forms of Se for sperm quality of breeder roosters and the relevant mechanisms.

METHODS

In experiment 1, 18-wk-old Jingbai laying breeder roosters were fed a Se-deficient base diet (BD, 0.06 mg Se/kg), or the BD + 0.1, 0.2, 0.3, 0.4, 0.5, or 1.0 mg Se/kg for 9 wk. In experiment 2, the roosters were fed the BD or the BD + sodium selenite (SeNa), seleno-yeast (SeY), or Se-nanoparticles (SeNPs) at 0.2 mg Se/kg for 9 wk.

RESULTS

In experiment 1, added dietary 0.2 and 0.3 mg Se/kg led to higher sperm motility and lower sperm mortality than the other groups at weeks 5, 7, and/or 9. Furthermore, added dietary 0.2-0.4 mg Se/kg produced better testicular histology and/or lower testicular 8-hydroxy-deoxyguanosine than the other groups. Moreover, integrated testicular transcriptomic and cecal microbiomic analysis revealed that inflammation, cell proliferation, and apoptosis-related genes and bacteria were dysregulated by Se deficiency or excess. In experiment 2, compared with SeNa, SeNPs slightly increased sperm motility throughout the experiment, whereas SeNPs slightly reduced sperm mortality compared with SeY at week 9. Both SeY and SeNPs decreased malondialdehyde in the serum than those of SeNa, and SeNPs led to higher glutathione peroxidase (GPX) and thioredoxin reductase activities and GPX1 and B-cell lymphoma 2 protein concentrations in the testis compared with SeY and SeNa.

CONCLUSIONS

The optimal dietary Se dose for reproductive health of breeder roosters is 0.25-0.35 mg Se/kg, and SeNPs displayed better effects on reproductive health than SeNa and SeY in laying breeder roosters. The optimal doses and forms of Se maintain reproductive health of roosters associated with regulation intestinal microbiota homeostasis and/or testicular redox balance, inflammation, cell proliferation, and apoptosis.

Advancing Biomarker Discovery and Therapeutic Targets in Duchenne Muscular Dystrophy: A Comprehensive Review

Mounting evidence underscores the intricate interplay between the immune system and skeletal muscles in Duchenne muscular dystrophy (DMD), as well as during regular muscle regeneration. While immune cell infiltration into skeletal muscles stands out as a prominent feature in the disease pathophysiology, a myriad of secondary defects involving metabolic and inflammatory pathways persist, with the key players yet to be fully elucidated. Steroids, currently the sole effective therapy for delaying onset and symptom control, come with adverse side effects, limiting their widespread use. Preliminary evidence spotlighting the distinctive features of T cell profiling in DMD prompts the immuno-characterization of circulating cells. A molecular analysis of their transcriptome and secretome holds the promise of identifying a subpopulation of cells suitable as disease biomarkers. Furthermore, it provides a gateway to unraveling new pathological pathways and pinpointing potential therapeutic targets. Simultaneously, the last decade has witnessed the emergence of novel approaches. The development and equilibrium of both innate and adaptive immune systems are intricately linked to the gut microbiota. Modulating microbiota-derived metabolites could potentially exacerbate muscle damage through immune system activation. Concurrently, genome sequencing has conferred clinical utility for rare disease diagnosis since innovative methodologies have been deployed to interpret the functional consequences of genomic variations. Despite numerous genes falling short as clinical targets for MD, the exploration of Tdark genes holds promise for unearthing novel and uncharted therapeutic insights. In the quest to expedite the translation of fundamental knowledge into clinical applications, the identification of novel biomarkers and disease targets is paramount. This initiative not only advances our understanding but also paves the way for the design of innovative therapeutic strategies, contributing to enhanced care for individuals grappling with these incapacitating diseases.

Balance between Pro- and Antifibrotic Proteins in Mesenchymal Stromal Cell Secretome Fractions Revealed by Proteome and Cell Subpopulation Analysis

Multipotent mesenchymal stromal cells (MSCs) regulate tissue repair through paracrine activity, with secreted proteins being significant contributors. Human tissue repair commonly results in fibrosis, where fibroblast differentiation into myofibroblasts is a major cellular mechanism. MSCs' paracrine activity can inhibit fibrosis development. We previously demonstrated that the separation of MSC secretome, represented by conditioned medium (CM), into subfractions enriched with extracellular vesicles (EV) or soluble factors (SF) boosts EV and SF antifibrotic effect. This effect is realized through the inhibition of fibroblast-to-myofibroblast differentiation in vitro. To unravel the mechanisms of MSC paracrine effects on fibroblast differentiation, we performed a comparative proteomic analysis of MSC secretome fractions. We found that CM was enriched in NF-κB activators and confirmed via qPCR that CM, but not EV or SF, upregulated NF-κB target genes (COX2, IL6, etc.) in human dermal fibroblasts. Furthermore, we revealed that EV and SF were enriched in TGF-β, Notch, IGF, and Wnt pathway regulators. According to scRNAseq, 11 out of 13 corresponding genes were upregulated in a minor MSC subpopulation disappearing in profibrotic conditions. Thus, protein enrichment of MSC secretome fractions and cellular subpopulation patterns shift the balance in fibroblast-to-myofibroblast differentiation, which should be considered in studies of MSC paracrine effects and the therapeutic use of MSC secretome.

Multiomic sequencing of paired primary and metastatic small bowel carcinoids

Background: Small bowel carcinoids are insidious tumors that are often metastatic when diagnosed. Limited mutation landscape studies of carcinoids indicate that these tumors have a relatively low mutational burden. The development of targeted therapies will depend upon the identification of mutations that drive the pathogenesis and metastasis of carcinoid tumors. Methods: Whole exome and RNA sequencing of 5 matched sets of normal tissue, primary small intestine carcinoid tumors, and liver metastases were investigated. Germline and somatic variants included: single nucleotide variants (SNVs), insertions/deletions (indels), structural variants, and copy number alterations (CNAs). The functional impact of mutations was predicted using Ensembl Variant Effect Predictor. Results: Large-scale CNAs were observed including the loss of chromosome 18 in all 5 metastases and 3/5 primary tumors. Certain somatic SNVs were metastasis-specific; including mutations in ATRX, CDKN1B, MXRA5 (leading to the activation of a cryptic splice site and loss of mRNA), SMARCA2, and the loss of UBE4B. Additional mutations in ATRX, and splice site loss of PYGL, leading to intron retention observed in primary and metastatic tumors. Conclusions: We observed novel mutations in primary/metastatic carcinoid tumor pairs, and some have been observed in other types of neuroendocrine tumors. We confirmed a previously observed loss of chromosome 18 and CDKN1B. Transcriptome sequencing added relevant information that would not have been appreciated with DNA sequencing alone. The detection of several splicing mutations on the DNA level and their consequences at the RNA level suggests that RNA splicing aberrations may be an important mechanism underlying carcinoid tumors.

Identifying immune cell infiltration and diagnostic biomarkers in heart failure and osteoarthritis by bioinformatics analysis

Heart failure (HF) and osteoarthritis (OA) are medical conditions that can significantly impact daily activities. Evidence has shown that HF and OA may share some pathogenic mechanisms. However, the underlying genomic mechanisms remain unclear. This study aimed to explore the underlying molecular mechanism and identify diagnostic biomarkers for HF and OA. With the cutoff criteria of fold change (FC) > 1.3 and P < .05, 920, 1500, 2195, and 2164 differentially expressed genes (DEGs) were identified in GSE57338, GSE116250, GSE114007, and GSE169077, respectively. After making the intersection of DEGs, we obtained 90 upregulated DEGs and 51 downregulated DEGs in HF datasets and 115 upregulated DEGs and 75 downregulated DEGs in OA datasets. Afterward, we conducted genome ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, protein-protein interaction (PPI) networks, and hub genes screening based on DEGs. Then, 4 common DEGs (fibroblast activation protein alpha [FAP], secreted frizzled-related protein 4 (SFRP4), Thy-1 cell surface antigen (THY1), matrix remodeling associated 5 [MXRA5]) between HF and OA were screened and validated in GSE5406 and GSE113825 datasets, based on which we established the support vector machine (SVM) models. The combined area under the receiver operating characteristic curve (AUC) of THY1, FAP, SFRP4, and MXRA5 in the HF training and test sets reached 0.949 and 0.928. While in the OA training set and test set, the combined AUC of THY1, FAP, SFRP4, and MXRA5 reached 1 and 1, respectively. The analysis of immune cells in HF revealed high levels of dendritic cell (DC), B cells, natural killer T cell (NKT), Type 1 regulatory T cell (Tr1), cytotoxic T cell (Tc), exhausted T cell (Tex), and mucosal-associated invariant T cell (MAIT), while displaying lower levels of monocytes, macrophages, NK, CD4 + T, gamma delta T (γδ T), T helper type 1 (Th1), T helper type 2 (Th2), and effector memory T cell (Tem). Moreover, the 4 common DEGs were positively correlated with DCs and B cells and negatively correlated with γδ T. In OA patients, the abundance of monocyte, macrophage, CD4 + naïve, and natural T regulatory cell (nTreg) was higher, while the infiltration of CD8 + T, γδ T, CD8 + naïve, and MAIT was lower. The expression of THY1 and FAP was significantly correlated with macrophage, CD8 + T, nTreg, and CD8 + naïve. SFRP4 was correlated with monocyte, CD8 + T, γδ T, CD4 + naïve, nTreg, CD8 + naïve and MAIT. MXRA5 was correlated with macrophage, CD8 + T, nTreg and CD8 + naïve. FAP, THY1, MXRA5, and SFRP4 may be diagnostic biomarkers for both HF and OA, and their correlation with immune cell infiltrations suggests shared immune pathogenesis.

Investigating the genetic basis of vertebrate dispersal combining RNA-seq, RAD-seq and quantitative genetics

Although animal dispersal is known to play key roles in ecological and evolutionary processes such as colonization, population extinction and local adaptation, little is known about its genetic basis, particularly in vertebrates. Untapping the genetic basis of dispersal should deepen our understanding of how dispersal behaviour evolves, the molecular mechanisms that regulate it and link it to other phenotypic aspects in order to form the so-called dispersal syndromes. Here, we comprehensively combined quantitative genetics, genome-wide sequencing and transcriptome sequencing to investigate the genetic basis of natal dispersal in a known ecological and evolutionary model of vertebrate dispersal: the common lizard, Zootoca vivipara. Our study supports the heritability of dispersal in semi-natural populations, with less variation attributable to maternal and natal environment effects. In addition, we found an association between natal dispersal and both variation in the carbonic anhydrase (CA10) gene, and in the expression of several genes (TGFB2, SLC6A4, NOS1) involved in central nervous system functioning. These findings suggest that neurotransmitters (serotonin and nitric oxide) are involved in the regulation of dispersal and shaping dispersal syndromes. Several genes from the circadian clock (CRY2, KCTD21) were also differentially expressed between disperser and resident lizards, supporting that the circadian rhythm, known to be involved in long-distance migration in other taxa, might affect dispersal as well. Since neuronal and circadian pathways are relatively well conserved across vertebrates, our results are likely to be generalisable, and we therefore encourage future studies to further investigate the role of these pathways in shaping dispersal in vertebrates.

Identification of matrix-remodeling associated 5 as a possible molecular oncotarget of pancreatic cancer

Pancreatic cancer has an extremely poor prognosis. Here we examined expression, potential functions and underlying mechanisms of MXRA5 (matrix remodeling associated 5) in pancreatic cancer. Bioinformatics studies revealed that MXRA5 transcripts are significantly elevated in pancreatic cancer tissues, correlating with the poor overall survival, high T-stage, N1 and pathologic stage of the patients. MXRA5 mRNA and protein expression is significantly elevated in microarray pancreatic cancer tissues and different pancreatic cancer cells. In primary and immortalized (BxPC-3 and PANC-1 lines) pancreatic cancer cells, shRNA-induced MXRA5 silencing or CRISPR/Cas9-mediated MXRA5 knockout suppressed cell survival, proliferation, migration, invasion, and epithelial-mesenchymal transition (EMT), while provoking cell apoptosis. Conversely, forced overexpression of MXRA5 further promoted pancreatic cancer cell progression and EMT. Bioinformatics studies and the protein chip analyses revealed that differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) in MXRA5-overexpressed primary pancreatic cancer cells were enriched in the PI3K-Akt-mTOR cascade. Indeed, Akt-mTOR activation in primary human pancreatic cancer cells was inhibited by MXRA5 shRNA or knockout, but was augmented following MXRA5 overexpression. In vivo, the growth of MXRA5 KO PANC-1 xenografts was largely inhibited in nude mice. Moreover, intratumoral injection of adeno-associated virus-packed MXRA5 shRNA potently inhibited primary pancreatic cancer cell growth in nude mice. Akt-mTOR activation was also largely inhibited in the MXRA5-depleted pancreatic cancer xenografts. Contrarily MXRA5 overexpression promoted primary pancreatic cancer cell growth in nude mice. Together, overexpressed MXRA5 is important for pancreatic cancer cell growth possibly through promoting EMT and Akt-mTOR activation. MXRA5 could be a potential therapeutic oncotarget for pancreatic cancer.

Proteomic profiling of extracellular vesicles in synovial fluid and plasma from Oligoarticular Juvenile Idiopathic Arthritis patients reveals novel immunopathogenic biomarkers

Introduction

New early low-invasive biomarkers are demanded for the management of Oligoarticular Juvenile Idiopathic Arthritis (OJIA), the most common chronic pediatric rheumatic disease in Western countries and a leading cause of disability. A deeper understanding of the molecular basis of OJIA pathophysiology is essential for identifying new biomarkers for earlier disease diagnosis and patient stratification and to guide targeted therapeutic intervention. Proteomic profiling of extracellular vesicles (EVs) released in biological fluids has recently emerged as a minimally invasive approach to elucidate adult arthritis pathogenic mechanisms and identify new biomarkers. However, EV-prot expression and potential as biomarkers in OJIA have not been explored. This study represents the first detailed longitudinal characterization of the EV-proteome in OJIA patients.

Methods

Fourty-five OJIA patients were recruited at disease onset and followed up for 24 months, and protein expression profiling was carried out by liquid chromatography-tandem mass spectrometry in EVs isolated from plasma (PL) and synovial fluid (SF) samples.

Results

We first compared the EV-proteome of SF vs paired PL and identified a panel of EV-prots whose expression was significantly deregulated in SF. Interaction network and GO enrichment analyses performed on deregulated EV-prots through STRING database and ShinyGO webserver revealed enrichment in processes related to cartilage/bone metabolism and inflammation, suggesting their role in OJIA pathogenesis and potential value as early molecular indicators of OJIA development. Comparative analysis of the EV-proteome in PL and SF from OJIA patients vs PL from age/gender-matched control children was then carried out. We detected altered expression of a panel of EV-prots able to differentiate new-onset OJIA patients from control children, potentially representing a disease-associated signature measurable at both the systemic and local levels with diagnostic potential. Deregulated EV-prots were significantly associated with biological processes related to innate immunity, antigen processing and presentation, and cytoskeleton organization. Finally, we ran WGCNA on the SF- and PL-derived EV-prot datasets and identified a few EV-prot modules associated with different clinical parameters stratifying OJIA patients in distinct subgroups.

Discussion

These data provide novel mechanistic insights into OJIA pathophysiology and an important contribution in the search of new candidate molecular biomarkers for the disease.

Transcriptomic differences between fibrotic and non-fibrotic testicular tissue reveal possible key players in Klinefelter syndrome-related testicular fibrosis

Klinefelter syndrome (KS; 47,XXY) affects 1-2 in 1000 males. Most men with KS suffer from an early germ cell loss and testicular fibrosis from puberty onwards. Mechanisms responsible for these processes remain unknown. Previous genomics studies on testis tissue from men with KS focused on germ cell loss, while a transcriptomic analysis focused on testicular fibrosis has not yet been performed. This study aimed to identify factors involved in the fibrotic remodelling of KS testes by analysing the transcriptome of fibrotic and non-fibrotic testicular tissue. RNA sequencing was performed to compare the genes expressed in testicular samples with (KS and testis atrophy) and without (Sertoli cell-only syndrome and fertile controls) fibrosis (n = 5, each). Additionally, differentially expressed genes (DEGs) between KS and testis atrophy samples were studied to reveal KS-specific fibrotic genes. DEGs were considered significant when p < 0.01 and log2FC > 2. Next, downstream analyses (GO and KEGG) were performed. Lastly, RNA in situ hybridization was performed to validate the results. The first analysis (fibrotic vs non-fibrotic) resulted in 734 significant DEGs (167 up- and 567 down-regulated). Genes involved in the extracellular structure organization (e.g. VCAM1) were found up-regulated. KEGG analysis showed an up-regulation of genes involved in the TGF-β pathway. The KS vs testis atrophy analysis resulted in 539 significant DEGs (59 up- and 480 down-regulated). Chronic inflammatory response genes were found up-regulated. The overlap of X-linked DEGs from the two analyses revealed three genes: matrix-remodelling associated 5 (MXRA5), doublecortin (DCX) and variable charge X-Linked 3B (VCX3B). RNA in situ hybridization showed an overexpression of VCAM1, MXRA5 and DCX within the fibrotic group compared with the non-fibrotic group. To summarize, this study revealed DEGs between fibrotic and non-fibrotic testis tissue, including VCAM1. In addition, X-linked fibrotic genes were revealed, e.g. MXRA5, DCX and VCX3B. Their potential role in KS-related testicular fibrosis needs further study.

Indirect Immobilised Jagged-1 Enhances Matrisome Proteins Associated with Osteogenic Differentiation of Human Dental Pulp Stem Cells: A Proteomic Study

The indirect immobilisation of Jagged-1 (Jagged-1) promoted osteogenic differentiation of human dental pulp cells (hDPs). Furthermore, the analysis of the Reactome pathway of RNA sequencing data indicates the upregulated genes involved with the extracellular matrix (ECM). Hence, our objective was to investigate the effects of Jagged-1 on proteomic profiles of human dental pulp stem cells (hDPSC). hDPSCs were cultured on the surface coated with human IgG Fc fragment (hFc) and the surface coated with rhJagged1/Fc recombinant protein-coated surface. Cells were differentiated to the osteogenic lineage using an osteogenic differentiation medium (OM) for 14 days, and cells cultured in a growth medium were used as a control. The protein component of the cultured cells was extracted into the cytosol, membrane, nucleus, and cytoskeletal compartment. Subsequently, the proteomic analysis was performed using liquid chromatography-tandem mass spectrometry (LC-MS). Metascape gene list analysis reported that Jagged-1 stimulated the expression of the membrane trafficking protein (DOP1B), which can indirectly improve osteogenic differentiation. hDPSCs cultured on Jagged-1 surface under OM condition expressed COL27A1, MXRA5, COL7A1, and MMP16, which played an important role in osteogenic differentiation. Furthermore, common matrisome proteins of all cellular components were related to osteogenesis/osteogenic differentiation. Additionally, the gene ontology categorised by the biological process of cytosol, membrane, and cytoskeleton compartments was associated with the biomineralisation process. The gene ontology of different culture conditions in each cellular component showed several unique gene ontologies. Remarkably, the Jagged-1_OM culture condition showed the biological process related to odontogenesis in the membrane compartment. In conclusion, the Jagged-1 induces osteogenic differentiation could, mainly through the regulation of protein in the membrane compartment.

Host gene expression is associated with viral shedding magnitude in blue-winged teals (Spatula discors) infected with low-path avian influenza virus

Intraspecific variation in host infectiousness affects disease transmission dynamics in human, domestic animal, and many wildlife host-pathogen systems including avian influenza virus (AIV); therefore, identifying host factors related to host infectiousness is important for understanding, controlling, and preventing future outbreaks. Toward this goal, we used RNA-seq data collected from low pathogenicity avian influenza virus (LPAIV)-infected blue-winged teal (Spatula discors) to determine the association between host gene expression and intraspecific variation in cloacal viral shedding magnitude, the transmissible fraction of virus. We found that host genes were differentially expressed between LPAIV-infected and uninfected birds early in the infection, host genes were differentially expressed between shed level groups at one-, three-, and five-days post-infection, host gene expression was associated with LPAIV infection patterns over time, and genes of the innate immune system had a positive linear relationship with cloacal viral shedding. This study provides important insights into host gene expression patterns associated with intraspecific LPAIV shedding variation and can serve as a foundation for future studies focused on the identification of host factors that drive or permit the emergence of high viral shedding individuals.

Inter-donor variability of extracellular matrix production in long-term cultures of human fibroblasts

Several tissue engineering approaches are based on the ability of mesenchymal cells to endogenously synthesize an extracellular matrix (ECM) in vitro, which can be seen as a form of biomaterial. Accordingly, the inter-donor variability of cell-assembled extracellular matrix (CAM) production is a key parameter to understand in order to progress towards clinical applications, especially for autologous strategies. In this study, CAMs were produced, under good manufacturing process conditions, from skin fibroblasts of 21 patients as part of a clinical trial to evaluate a tissue-engineered vascular graft. The inter-donor variability of CAM strength, thickness, hydroxyproline, and glycosaminoglycan was substantial (coefficient of variability of 33%, 19%, 24%, and 19%, respectively), but a significant correlation was observed between all four properties (Pearson r: 0.43 to 0.70; p-value ≤ 0.05). A CAM matrisome analysis, performed by mass spectrometry, revealed the presence of 70 ECM-related proteins. Our study shows that the relative abundance of 16 proteins (15 non-collagenous) correlated with CAM thickness. These proteins also correlated with CAM hydroxyproline content, as well as 21 other proteins that included fibrillar collagens and non-collagenous proteins. However, data demonstrated that only the relative abundance of type I collagen subunit alpha-1 was correlated to CAM strength. This study is the most extensive evaluation of CAM inter-donor variability to date and will help tissue engineers working with this type of biomaterial to design strategies that take into account this variability, especially for autologous tissue manufacturing.

Several tissue engineering approaches are based on the ability of mesenchymal cells to endogenously synthesize an extracellular matrix (ECM) in vitro, which can be seen as a form of biomaterial.

Gill and Liver Transcript Expression Changes Associated With Gill Damage in Atlantic Salmon (Salmo salar)

Gill damage represents a significant challenge in the teleost fish aquaculture industry globally, due to the gill's involvement in several vital functions and direct contact with the surrounding environment. To examine the local and systemic effects accompanying gill damage (which is likely to negatively affect gill function) of Atlantic salmon, we performed a field sampling to collect gill and liver tissue after several environmental insults (e.g., harmful algal blooms). Before sampling, gills were visually inspected and gill damage was scored; gill scores were assigned from pristine [gill score 0 (GS0)] to severely damaged gills (GS3). Using a 44K salmonid microarray platform, we aimed to compare the transcriptomes of pristine and moderately damaged (i.e., GS2) gill tissue. Rank Products analysis (5% percentage of false-positives) identified 254 and 34 upregulated and downregulated probes, respectively, in GS2 compared with GS0. Differentially expressed probes represented genes associated with functions including gill remodeling, wound healing, and stress and immune responses. We performed gill and liver qPCR for all four gill damage scores using microarray-identified and other damage-associated biomarker genes. Transcripts related to wound healing (e.g., neb and klhl41b) were significantly upregulated in GS2 compared with GS0 in the gills. Also, transcripts associated with immune and stress-relevant pathways were dysregulated (e.g., downregulation of snaclec 1-like and upregulation of igkv3) in GS2 compared with GS0 gills. The livers of salmon with moderate gill damage (i.e., GS2) showed significant upregulation of transcripts related to wound healing (i.e., chtop), apoptosis (e.g., bnip3l), blood coagulation (e.g., f2 and serpind1b), transcription regulation (i.e., pparg), and stress-responses (e.g., cyp3a27) compared with livers of GS0 fish. We performed principal component analysis (PCA) using transcript levels for gill and liver separately. The gill PCA showed that PC1 significantly separated GS2 from all other gill scores. The genes contributing most to this separation were pgam2, des, neb, tnnt2, and myom1. The liver PCA showed that PC1 significantly separated GS2 from GS0; levels of hsp70, cyp3a27, pparg, chtop, and serpind1b were the highest contributors to this separation. Also, hepatic acute phase biomarkers (e.g., serpind1b and f2) were positively correlated to each other and to gill damage. Gill damage-responsive biomarker genes and associated qPCR assays arising from this study will be valuable in future research aimed at developing therapeutic diets to improve farmed salmon welfare.

Single-Cell RNA-Seq Analysis Reveals Macrophage Involved in the Progression of Human Intervertebral Disc Degeneration

Intervertebral disc degeneration (IDD) has been considered as the primary pathological mechanism that underlies low back pain. Understanding the molecular mechanisms underlying human IDD is imperative for making strategies to treat IDD-related diseases. Herein, we report the molecular programs, lineage progression patterns, and paths of cellular communications during the progression of IDD using single-cell RNA sequencing (scRNA-seq) on nucleus pulposus (NP) cells from patients with different grades of IDD undergoing discectomy. New subtypes of cells and cell-type-specific gene signatures of the metabolic homeostatic NP cells (Met NPC), adhesive NP cells (Adh NPC), inflammatory response NP cells (IR NPC), endoplasmic reticulum stress NP cells (ERS NPC), fibrocartilaginous NP cells (Fc NPC), and CD70 and CD82+ progenitor NP cells (Pro NPC) were identified. In the late stage of IDD, the IR NPC and Fc NPC account for a large proportion of NPC. Importantly, immune cells including macrophages, T cells, myeloid progenitors, and neutrophils were also identified, and further analysis showed that significant intercellular interaction between macrophages and Pro NPC occurred via MIF (macrophage migration inhibitory factor) and NF-kB signaling pathways during the progression of IDD. In addition, dynamic polarization of macrophage M1 and M2 cell subtypes was found in the progression of IDD, and gene set functional enrichment analysis suggested a significant role of the macrophage polarization in regulating cell metabolism, especially the Pro NPC. Finally, we found that the NP cells in the late degenerative stage were mainly composed of the cell types related to inflammatory and endoplasmic reticulum (ER) response, and fibrocartilaginous activity. Our results provided new insights into the identification of NP cell populations at single-cell resolution and at the relatively whole-transcriptome scale, accompanied by cellular communications between immune cells and NP cells, and discriminative markers in relation to specific cell subsets. These new findings present clues for effective and functional manipulation of human IDD-related bioremediation and healthcare.

Three Months of Strength Training Changes the Gene Expression of Inflammation-Related Genes in PBMC of Older Women: A Randomized Controlled Trial

Here, we investigate changes in inflammation-related gene-expression in peripheral mononuclear blood cells (PBMC) by strength training. A total of 14 women aged ≥65 years were randomized into 3 months of either 3×/week intensive strength training (IST: 3×10 rep at 80% 1RM), strength endurance training (SET: 2×30 reps at 40% 1RM) or control (CON: 3×30 sec stretching). Differentially expressed genes (fold change ≤0.67 or ≥1.5) were identified by targeted RNA-sequencing of 407 inflammation-related genes. A total of 98 genes (n = 61 pro-inflammatory) were significantly affected. IST and SET altered 14 genes in a similar direction and 19 genes in the opposite direction. Compared to CON, IST changed the expression of 6 genes in the same direction, and 17 genes in the SET. Likewise, 18 and 13 genes were oppositely expressed for, respectively, IST and SET compared to CON. Changes in gene expression affected 33 canonical pathways related to chronic inflammation. None of the altered pathways overlapped between IST and SET. Liver X Receptor/Retinoid X Receptor Activation (LXR/RXR) and Triggering Receptor Expressed On Myeloid Cells 1 (TREM1) pathways were enriched oppositely in both training groups. We conclude that three months IST and SET can induce changes in CLIP-related gene expression in PBMC, but by affecting different genes and related pathways.

Cellular features of localized microenvironments in human meniscal degeneration: a single-cell transcriptomic study

Background

Musculoskeletal tissue degeneration impairs the life quality and function of many people. Meniscus degeneration is a major origin of knee osteoarthritis and a common threat to athletic ability, but its cellular mechanism remains elusive.

Methods

We built a cell atlas of 12 healthy or degenerated human meniscus samples from the inner and outer meniscal zones of 8 patients using scRNA-seq to investigate meniscal microenvironment homeostasis and its changes in the degeneration process and verified findings with immunofluorescent imaging.

Results

We identified and localized cell types in inner and outer meniscus and found new chondrocyte subtypes associated with degeneration. The observations suggested understandings on how cellular compositions, functions, and interactions participated in degeneration, and on the possible loop-like interactions among extracellular matrix disassembly, angiogenesis, and inflammation in driving the degeneration.

Conclusions

The study provided a rich resource reflecting variations in the meniscal microenvironment during degeneration and suggested new cell subtypes as potential therapeutic targets. The hypothesized mechanism could also be a general model for other joint degenerations.

Funding

The National Natural Science Foundation of China (81972123, 82172508, 62050178, 61721003), the National Key Research and Development Program of China (2021YFF1200901), Fundamental Research Funds for the Central Universities (2015SCU04A40); The Innovative Spark Project of Sichuan University (2018SCUH0034); Sichuan Science and Technology Program (2020YFH0075); Chengdu Science and Technology Bureau Project (2019-YF05-00090-SN); 1.3.5 Project for Disciplines of Excellence of West China Hospital Sichuan University (ZYJC21030, ZY2017301); 1.3.5 Project for Disciplines of Excellence - Clinical Research Incubation Project, West China Hospital, Sichuan University (2019HXFH039).

A single-cell transcriptome of mesenchymal stromal cells to fabricate bioactive hydroxyapatite materials for bone regeneration

The osteogenic microenvironment of bone-repairing materials plays a key role in accelerating bone regeneration but remains incompletely defined, which significantly limits the application of such bioactive materials. Here, the transcriptional landscapes of different osteogenic microenvironments, including three-dimensional (3D) hydroxyapatite (HA) scaffolds and osteogenic medium (OM), for mesenchymal stromal cells (MSCs) in vitro were mapped at single-cell resolution. Our findings suggested that an osteogenic process reminiscent of endochondral ossification occurred in HA scaffolds through sequential activation of osteogenic-related signaling pathways, along with inflammation and angiogenesis, but inhibition of adipogenesis and fibrosis. Moreover, we revealed the mechanism during OM-mediated osteogenesis involves the ZBTB16 and WNT signaling pathways. Heterogeneity of MSCs was also demonstrated. In vitro ossification of LRRC75A⁺ MSCs was shown to have better utilization of WNT-related ossification process, and PCDH10⁺ MSCs with superiority in hydroxyapatite-related osteogenic process. These findings provided further understanding of the cellular activity modulated by OM conditions and HA scaffolds, providing new insights for the improvement of osteogenic biomaterials. This atlas provides a blueprint for research on MSC heterogeneity and the osteogenic microenvironment of HA scaffolds and a database reference for the application of bioactive materials for bone regeneration.

Identification of N-glycoproteins of hip cartilage in patients with osteonecrosis of femoral head using quantitative glycoproteomics

N-glycosylation is a major post-translational modification of proteins and involved in many diseases, however, the state and role of N-glycosylation in cartilage degeneration of osteonecrosis of femoral head (ONFH) remain unclear. The aim of this study is to identify the glycoproteins of ONFH hip cartilage. Cartilage tissues were collected from nine patients with ONFH and nine individuals with traumatic femoral neck fracture. Cartilage glycoproteins were identified by glycoproteomics based on LC-MS/MS. The differentially N-glycoproteins including glycosites were identified in ONFH and controls. A total of 408 N-glycoproteins with 444 N-glycosites were identified in ONFH and control cartilage. Among them, 104 N-glycoproteins with 130 N-glycosites were significantly differential in ONFH and control cartilage, which including matrix-remodeling-associated protein 5, prolow-density lipoprotein receptor-related protein 1, clusterin and lysosome-associated membrane glycoprotein 2. Gene Ontology analysis revealed the significantly differential glycoproteins mainly belonged to protein metabolic process, single-multicellular organism process, proteolysis, biological adhesion and cell adhesion. KEGG pathway and protein-protein interaction analysis suggested that the significantly differential glycoproteins were associated with PI3K-Akt signalling pathway, ECM-receptor interaction, protein processing in the endoplasmic reticulum and N-glycan biosynthesis. This information provides substantial insight into the role of protein glycosylation in the development of cartilage degeneration of ONFH patients.

Equine Mesenchymal Stem/Stromal Cells Freeze-Dried Secretome (Lyosecretome) for the Treatment of Musculoskeletal Diseases: Production Process Validation and Batch Release Test for Clinical Use

In the last decades, it has been demonstrated that the regenerative therapeutic efficacy of mesenchymal stromal cells is primarily due to the secretion of soluble factors and extracellular vesicles, collectively known as secretome. In this context, our work described the preparation and characterization of a freeze-dried secretome (Lyosecretome) from adipose tissue-derived mesenchymal stromal cells for the therapy of equine musculoskeletal disorder. An intraarticular injectable pharmaceutical powder has been formulated, and the technological process has been validated in an authorized facility for veterinary clinical-use medicinal production. Critical parameters for quality control and batch release have been identified regarding (i) physicochemical properties; (ii) extracellular vesicle morphology, size distribution, and surface biomarker; (iii) protein and lipid content; (iv) requirements for injectable pharmaceutical dosage forms such as sterility, bacterial endotoxin, and Mycoplasma; and (v) in vitro potency tests, as anti-elastase activity and proliferative activity on musculoskeletal cell lines (tenocytes and chondrocytes) and mesenchymal stromal cells. Finally, proteins putatively responsible for the biological effects have been identified by Lyosecretome proteomic investigation: IL10RA, MXRA5, RARRES2, and ANXA1 modulate the inflammatory process RARRES2, NOD1, SERPINE1, and SERPINB9 with antibacterial activity. The work provides a proof-of-concept for the manufacturing of clinical-grade equine freeze-dried secretome, and prototypes are now available for safety and efficacy clinical trials in the treatment of equine musculoskeletal diseases

MXRA5 Is a Novel Immune-Related Biomarker That Predicts Poor Prognosis in Glioma

Background

Glioma is the most common primary intracranial tumor and is associated with poor prognosis. Identifying effective biomarkers for glioma is particularly important. MXRA5, a secreted glycoprotein, is involved in cell adhesion and extracellular matrix remodeling and has been reported to be expressed in many cancers. However, the role and mechanism of action of MXRA5 in gliomas remain unclear. This study was aimed at investigating the role of MXRA5 at the transcriptome level and its clinical prognostic value.

Methods

In this study, RNA microarray data of 301 glioma patients from the Chinese Glioma Genome Atlas (CGGA) were collected as a training cohort and RNA-seq data of 702 glioma samples from The Cancer Genome Atlas (TCGA) were used for validation. We analyzed the clinical and molecular characteristics as well as the prognostic value of MXRA5 in glioma. In addition, the expression level of MXRA was evaluated in 28 glioma tissue samples.

Results

We found that MXRA5 expression was significantly upregulated in high-grade gliomas and IDH wild-type gliomas compared to controls. Receiver operating characteristic (ROC) analysis showed that MXRA5 is a potential marker of the mesenchymal subtype of glioblastoma multiforme (GBM). We found that MXRA5 expression is highly correlated with immune checkpoint molecule expression levels and tumor-associated macrophage infiltration. High MXRA5 expression could be used as an independent indicator of poor prognosis in glioma patients.

Conclusion

Our study suggests that MXRA5 expression is associated with the clinicopathologic features and poor prognosis of gliomas. MXRA5 may play an important role in the immunosuppressive microenvironment of glioma. As a secreted glycoprotein, MXRA5 is a potential circulating biomarker for glioma, deserving further investigation.

Of mice and men - and guinea pigs?

This year marks the twentieth anniversary of the publication of the first draft of the human genome and its broad availability to the scientific community. In parallel, the annotation of the mouse genome led to the identification and analysis of countless genes by means of genetic manipulation. Today, when comparing both genomes, it might surprise that some genes are still seeking their respective homologs in either species. In this review, we aim at raising awareness for the remarkable differences between the researcher's favorite rodents, i.e., mice and rats, when it comes to the generation of rodent research models regarding genes with a particular delicate localization, namely the pseudoautosomal region on both sex chromosomes. Many of these genes are of utmost clinical relevance in humans and still miss a rodent disease model giving their absence in mice and rats or low sequence similarity compared to humans. The abundance of rodents within mammals prompted us to investigate different branches of rodents leading us to the re-discovery of the guinea pig as a mammalian research model for a distinct group of genes.

Intratumoral Fibrosis in Facilitating Renal Cancer Aggressiveness: Underlying Mechanisms and Promising Targets

Intratumoral fibrosis is a histologic manifestation of fibrotic tumor stroma. The interaction between cancer cells and fibrotic stroma is intricate and reciprocal, involving dysregulations from multiple biological processes. Different components of tumor stroma are implicated via distinct manners. In the kidney, intratumoral fibrosis is frequently observed in renal cell carcinoma (RCC). However, the underlying mechanisms remain largely unclear. In this review, we recapitulate evidence demonstrating how fibrotic stroma interacts with cancer cells and mechanisms shared between RCC tumorigenesis and renal fibrogenesis, providing promising targets for future studies.

TGF-β1 Signaling: Immune Dynamics of Chronic Kidney Diseases

Chronic kidney disease (CKD) is a major cause of morbidity and mortality worldwide, imposing a great burden on the healthcare system. Regrettably, effective CKD therapeutic strategies are yet available due to their elusive pathogenic mechanisms. CKD is featured by progressive inflammation and fibrosis associated with immune cell dysfunction, leading to the formation of an inflammatory microenvironment, which ultimately exacerbating renal fibrosis. Transforming growth factor β1 (TGF-β1) is an indispensable immunoregulator promoting CKD progression by controlling the activation, proliferation, and apoptosis of immunocytes via both canonical and non-canonical pathways. More importantly, recent studies have uncovered a new mechanism of TGF-β1 for de novo generation of myofibroblast via macrophage-myofibroblast transition (MMT). This review will update the versatile roles of TGF-β signaling in the dynamics of renal immunity, a better understanding may facilitate the discovery of novel therapeutic strategies against CKD.

Genome-wide meta-analysis identifies 127 open-angle glaucoma loci with consistent effect across ancestries

Primary open-angle glaucoma (POAG), is a heritable common cause of blindness world-wide. To identify risk loci, we conduct a large multi-ethnic meta-analysis of genome-wide association studies on a total of 34,179 cases and 349,321 controls, identifying 44 previously unreported risk loci and confirming 83 loci that were previously known. The majority of loci have broadly consistent effects across European, Asian and African ancestries. Cross-ancestry data improve fine-mapping of causal variants for several loci. Integration of multiple lines of genetic evidence support the functional relevance of the identified POAG risk loci and highlight potential contributions of several genes to POAG pathogenesis, including SVEP1, RERE, VCAM1, ZNF638, CLIC5, SLC2A12, YAP1, MXRA5, and SMAD6. Several drug compounds targeting POAG risk genes may be potential glaucoma therapeutic candidates.

Gene Expression Profiles Reveal Extracellular Matrix and Inflammatory Signaling in Radiation-Induced Premature Differentiation of Human Fibroblast in vitro

Purpose

Fibroblasts are considered to play a major role in the development of fibrotic reaction after radiotherapy and premature radiation-induced differentiation has been proposed as a cellular basis. The purpose was to relate gene expression profiles to radiation-induced phenotypic changes of human skin fibroblasts relevant for radiogenic fibrosis.

Materials and Methods

Exponentially growing or confluent human skin fibroblast strains were irradiated in vitro with 1–3 fractions of 4 Gy X-rays. The differentiated phenotype was detected by cytomorphological scoring and immunofluorescence microscopy. Microarray analysis was performed on Human Genome U133 plus2.0 microarrays (Affymetrix) with JMP Genomics software, and pathway analysis with Reactome R-package. The expression levels and kinetics of selected genes were validated with quantitative real-time PCR (qPCR) and Western blotting.

Results

Irradiation of exponentially growing fibroblast with 1 × 4 Gy resulted in phenotypic differentiation over a 5-day period. This was accompanied by downregulation of cell cycle-related genes and upregulation of collagen and other extracellular matrix (ECM)-related genes. Pathway analysis confirmed inactivation of proliferation and upregulation of ECM- and glycosaminoglycan (GAG)-related pathways. Furthermore, pathways related to inflammatory reactions were upregulated, and potential induction and signaling mechanisms were identified. Fractionated irradiation (3 × 4 Gy) of confluent cultures according to a previously published protocol for predicting the risk of fibrosis after radiotherapy showed similar downregulation but differences in upregulated genes and pathways.

Conclusion

Gene expression profiles after irradiation of exponentially growing cells were related to radiation-induced differentiation and inflammatory reactions, and potential signaling mechanisms. Upregulated pathways by different irradiation protocols may reflect different aspects of the fibrogenic process thus providing a model system for further hypothesis-based studies of radiation-induced fibrogenesis.

Proteome-wide and matrisome-specific alterations during human pancreas development and maturation

The extracellular matrix (ECM) is unique to each tissue and capable of guiding cell differentiation, migration, morphology, and function. The ECM proteome of different developmental stages has not been systematically studied in the human pancreas. In this study, we apply mass spectrometry-based quantitative proteomics strategies using N,N-dimethyl leucine isobaric tags to delineate proteome-wide and ECM-specific alterations in four age groups: fetal (18-20 weeks gestation), juvenile (5-16 years old), young adults (21-29 years old) and older adults (50-61 years old). We identify 3,523 proteins including 185 ECM proteins and quantify 117 of them. We detect previously unknown proteome and matrisome features during pancreas development and maturation. We also visualize specific ECM proteins of interest using immunofluorescent staining and investigate changes in ECM localization within islet or acinar compartments. This comprehensive proteomics analysis contributes to an improved understanding of the critical roles that ECM plays throughout human pancreas development and maturation.

TWEAK-Fn14 as a common pathway in the heart and the kidneys in cardiorenal syndrome

There is a complex relationship between cardiac and renal disease, often referred to as the cardiorenal syndrome. Heart failure adversely affects kidney function, and both acute and chronic kidney disease are associated with structural and functional changes to the myocardium. The pathological mechanisms and contributing interactions that surround this relationship remain poorly understood, limiting the opportunities for therapeutic intervention. The cytokine tumor necrosis factor-like weak inducer of apoptosis (TWEAK) and its receptor, fibroblast growth factor-inducible 14 (Fn14), are abundantly expressed in injured kidneys and heart. The TWEAK-Fn14 axis promotes responses that drive tissue injury such as inflammation, proliferation, fibrosis, and apoptosis, while restraining the expression of tissue protective factors such as the anti-aging factor Klotho and the master regulator of mitochondrial biogenesis peroxisome proliferator-activated receptor-γ coactivator-1α (PGC-1α). High levels of TWEAK induce cardiac remodeling, and promote inflammation, tubular and podocyte injury and death, fibroblast proliferation, and, ultimately, renal fibrosis. Accordingly, targeting the TWEAK-Fn14 axis is protective in experimental kidney and heart disease. TWEAK has also emerged as a biomarker of kidney damage and cardiovascular outcomes and has been successfully targeted in clinical trials. In this review, we update our current knowledge of the roles of the TWEAK-Fn14 axis in cardiovascular and kidney disease and its potential contribution to the cardiorenal syndrome. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Proteomic Architecture of Valvular Extracellular Matrix: FNDC1 and MXRA5 Are New Biomarkers of Aortic Stenosis

•

ECM proteins play an important role in maintaining the structural architecture and the mechanical behavior of the aortic valve.

•

Network analysis highlights a strong connection between metabolic markers and ECM proteins.

•

MXRA5 and FNDC1 were identified as new biomarkers of aortic stenosis in 2 independent cohorts

This study analyzed the expression of extracellular matrix (ECM) proteins during aortic valve calcification with mass spectrometry, and further validated in an independent human cohort using RNAseq data. The study reveals that valve calcification is associated with significant disruption in ECM and metabolic pathways, and highlights a strong connection between metabolic markers and ECM remodeling. It also identifies FNDC1 and MXRA5 as novel ECM biomarkers in calcified valves, electing them as potential targets in the development and progression of aortic stenosis.

Integrated analysis of mRNA and microRNA expression profiles reveals differential transcriptome signature in ischaemic and dilated cardiomyopathy induced heart failure

Cardiac remodelling is widely accepted as a common characteristic for many heart diseases, especially in heart failure (HF). Ischaemic cardiomyopathy (ICM) and dilated cardiomyopathy (DCM) are associated with cardiac remodelling. Both mRNA and microRNA are potential diagnostic markers and therapeutic targets of cardiac remodelling in HF. However, the mechanisms of microRNA-mRNA joint regulation in HF are still unclear. In this study, 3 gene expression profiles from patients with and without HF were analysed to harvest shared differentially expressed genes (microRNA and mRNA) with significant major biological function. Moreover, key genes highly related to ICM and DCM-induced HF were screened out through a Weighted Genes Co-Expression Network Analysis (WGCNA). Based on microRNA-mRNA analysis, several microRNAs and target genes were identified. Combined with pathway analysis, we found that miR-542-3p and its target gene CILP were likely involved in the regulation of TGF-β signalling pathway in ICM induced HF. Collectively, the microRNA-mRNA interaction network analysis revealed that miR-542-3p-CILP as mediator of TGF-β signalling pathway might be a new mechanism to mediate ICM induced HF. This study provides certain novel targets for diagnosis and therapeutic treatment of ICM- and DCM-induced HF.

Tunica-Specific Transcriptome of Abdominal Aortic Aneurysm and the Effect of Intraluminal Thrombus, Smoking, and Diameter Growth Rate

OBJECTIVE

There is no medical treatment to prevent abdominal aortic aneurysm (AAA) growth and rupture, both of which are linked to smoking. Our objective was to map the tunica-specific pathophysiology of AAA with consideration of the intraluminal thrombus, age, and sex, and to subsequently identify which mechanisms were linked to smoking and diameter growth rate. Approach and Results: Microarray analyses were performed on 246 samples from 76 AAA patients and 13 controls. In media and adventitia, there were 5889 and 2701 differentially expressed genes, respectively. Gene sets related to adaptive and innate immunity were upregulated in both tunicas. Media-specific gene sets included increased matrix disassembly and angiogenesis, as well as decreased muscle cell development, contraction, and differentiation. Genes implicated in previous genome-wide association studies were dysregulated in media. The intraluminal thrombus had a pro-proteolytic and proinflammatory effect on the underlying media. Active smoking resulted in increased inflammation, oxidative stress, and angiogenesis in all tissues and enriched lipid metabolism in adventitia. Processes enriched with active smoking in control aortas overlapped to a high extent with those differentially expressed between AAAs and controls. The AAA diameter growth rate (n=24) correlated with T- and B-cell expression in media, as well as lipid-related processes in the adventitia.

CONCLUSIONS

This tunica-specific analysis of gene expression in a large study enabled the detection of features not previously described in AAA disease. Smoking was associated with increased expression of aneurysm-related processes, of which adaptive immunity and lipid metabolism correlated with growth rate.

Deep Sequencing Analysis Identified a Specific Subset of Mutations Distinctive of Biphasic Malignant Pleural Mesothelioma

Malignant Pleural Mesothelioma (MPM) is a heterogeneous disease. Morphologically, three different phenotypes are distinguishable: epithelioid (e-), sarcomatoid (s-) and biphasic (biph-) MPM, the latest, being a mixture of e- and s-MPM cells. Being an intermediate entity, management of biph-MPM, remains debatable and controversial, with different guidelines recommending distinct approaches. Identification of biph-MPM associated genetic alterations, through deep sequencing analysis, may provide useful tools to understand these lesions. A retrospective cohort of 69 surgically resected MPMs, 39 biph-MPMs (56.5%) and 30 e-MPMs (43.5%) was selected. A separate set of 16 biph-MPM was used as validation set. Deep sequencing analysis on an MPM-specific custom panel (MPM_geneset) comprising 1041 amplicons spanning 34 genes was performed. A total of 588 variants and 5309 mutational events were detected. In total, 91.3% of MPMs showed at least one mutation and 76.8% showed co-occurrence of more than one alteration. Mutations in MXRA5 (p = 0.05) and NOD2 (p = 0.018) were significantly associated with biph-MPM both in the training and validation cohort and correlated with the extent of the sarcomatoid component. Mutations in NOD2 and XRCC6 correlated with patients’ survival. We demonstrated that biph-MPM are associated with a specific mutation set, and that genetic analysis at diagnosis may improve patients’ risk stratification.

Matrix Remodeling-Associated Protein 5 in Urinary Exosomes as a Potential Novel Marker of Obstructive Nephropathy in Children With Ureteropelvic Junction Obstruction

Recent investigations have described the use of urinary matrix remodeling-associated protein 5 (MXRA5) as a novel biomarker of kidney impairment in the setting of chronic kidney disease. In this study, we aimed to evaluate the possible clinical application of urinary MXRA5 as a useful non-invasive marker in the urine from the affected renal pelvis and bladder of children with ureteropelvic junction obstruction (UPJO). We conducted a prospective cohort study of patients aged <12 months with prenatally diagnosed unilateral UPJO who underwent dismembered pyeloplasty in 2018 or 2019, and a sex- and age-matched control group of healthy children. Blood urea nitrogen and creatinine levels were normal in all the patients. The whole urine and urinary exosomal concentrations of MXRA5 were measured by enzyme-linked immunosorbent assay. The correlations between bladder/renal pelvic MXRA5 levels and differential renal function (DRF) in the affected kidney were also determined. A total of 35 UPJO patients and 12 controls were enrolled in the study. There was no significant difference in whole-urine MXRA5 level between the controls and UPJO patients. However, the exosomal MXRA5 level was significantly lower in the controls than in patients with UPJO (p < 0.05). There were non-significant correlations between bladder and renal pelvis whole-urine MXRA5 levels and DRF (R 2 = 0.1115, p = 0.05 and R 2 = 0.3313, p = 0.0502, respectively). The strongest correlation was between exosomal MXRA5 level in the renal pelvis and DRF (R 2 = 0.8128, p < 0.0001). Urinary exosomal MXRA5 level was significantly higher in children with UPJO than controls. Higher urinary exosomal MXRA5 levels were significantly correlated with lower DRF in the affected kidney in children with UPJO.

Differential proteomic analysis of fetal and geriatric lumbar nucleus pulposus: immunoinflammation and age-related intervertebral disc degeneration

BACKGROUND

Intervertebral disc degeneration (IVDD) is a major cause of low back pain. Although the mechanism of degeneration remains unclear, aging has been recognized as a key risk factor for IVDD. Most studies seeking to identify IVDD-associated molecular alterations in the context of human age-related IVDD have focused only on a limited number of proteins. Differential proteomic analysis is an ideal method for comprehensively screening altered protein profiles and identifying the potential pathways related to pathological processes such as disc degeneration.

METHODS

In this study, tandem mass tag (TMT) labeling was combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS) for differential proteomic analysis of human fetal and geriatric lumbar disc nucleus pulposus (NP) tissue. Parallel reaction monitoring (PRM) and Western blotting (WB) techniques were used to identify target proteins. Bioinformatic analyses, including Gene Ontology (GO) annotation, domain annotation, pathway annotation, subcellular localization and functional enrichment analyses, were used to interpret the potential significance of the protein alterations in the mechanism of IVDD. Student's t-tests and two-tailed Fisher's exact tests were used for statistical analysis.

RESULTS

Six hundred forty five proteins were significantly upregulated and 748 proteins were downregulated in the geriatric group compared with the fetal group. Twelve proteins were verified to have significant differences in abundance between geriatric and fetal NP tissue; most of these have not been previously identified as being associated with human IVDD. The potential significance of the differentially expressed proteins in age-related IVDD was analyzed from multiple perspectives, especially with regard to the association of the immunoinflammatory response with IVDD.

CONCLUSIONS

Differential proteomic analysis was used as a comprehensive strategy for elucidating the protein alterations associated with age-related IVDD. The findings of this study will aid in the screening of new biomarkers and molecular targets for the diagnosis and therapy of IVDD. The results may also significantly enhance our understanding of the pathophysiological process and mechanism of age-related IVDD.

Identification and functional activity of matrix-remodeling associated 5 (MXRA5) in benign hyperplastic prostate

Objective: Benign prostatic hyperplasia (BPH) is a common condition in aging males. The current study aims to identify differentially expressed genes (DEGs) associated with BPH and to elucidate the role of matrix-remodeling associated 5 (MXRA5) protein and mitogen-activated protein kinase (MAPK) signaling pathways in BPH.

Results: A total of 198 DEGs and a number of related pathways were identified with MXRA5 being one of the most significantly altered DEGs. MXRA5 was upregulated in BPH samples and localized mostly in stroma. Knockdown of MXRA5 induced stromal cell cycle arrest instead of inhibiting apoptosis. Consistently, MXRA5 overexpression enhanced epithelial cell proliferation. In addition, phosphorylated ERK1/2 and p38, key members of the MAPK family, were strongly decreased with knockdown but increased with overexpression.

Conclusion: Our novel data demonstrates that upregulation of MXRA5 in the enlarged prostate could contribute to the development of BPH through increasing cell proliferation via the MAPK pathway. Thus, the MXRA5-MAPK system could be rediscovered as a new therapeutic target for treating BPH.

Methods: Microarray analysis and integrated bioinformatics were conducted. The expression and biologic functions of MXRA5 was investigated via RT-PCR, western-blot, immunofluorescence, flow cytometry and MTT assay. Finally, genes involved in regulation of the MAPK pathway were investigated.

The Impact of Transcription Factor Prospero Homeobox 1 on the Regulation of Thyroid Cancer Malignancy

Transcription factor Prospero homeobox 1 (PROX1) is continuously expressed in the lymphatic endothelial cells, playing an essential role in their differentiation. Many reports have shown that PROX1 is implicated in cancer development and acts as an oncoprotein or suppressor in a tissue-dependent manner. Additionally, the PROX1 expression in many types of tumors has prognostic significance and is associated with patient outcomes. In our previous experimental studies, we showed that PROX1 is present in the thyroid cancer (THC) cells of different origins and has a high impact on follicular thyroid cancer (FTC) phenotypes, regulating migration, invasion, focal adhesion, cytoskeleton reorganization, and angiogenesis. Herein, we discuss the PROX1 transcript and protein structures, the expression pattern of PROX1 in THC specimens, and its epigenetic regulation. Next, we emphasize the biological processes and genes regulated by PROX1 in CGTH-W-1 cells, derived from squamous cell carcinoma of the thyroid gland. Finally, we discuss the interaction of PROX1 with other lymphatic factors. In our review, we aimed to highlight the importance of vascular molecules in cancer development and provide an update on the functionality of PROX1 in THC biology regulation.

miR-30b facilitates preeclampsia through targeting MXRA5 to inhibit the viability, invasion and apoptosis of placental trophoblast cells

Preeclampsia (PE) may induce gestational failure threatening a significant number of pregnant women. Dysfunctional placental trophoblast cells have an important impact on PE progression. microRNAs (miRNAs) have been reported to participate in PE progression, whereas the mechanism that underlies miR-30b involved in PE progression and function of placental trophoblast cells remains poorly understood. Cell viability was investigated by cell counting kit-8 (CCK-8) assay. Cell apoptosis was detected by flow cytometry using Annexin V-FITC/propidium iodide (PI) staining. Cell invasion was analyzed by trans-well assay. The expression of miR-30b was measured by quantitative real-time polymerase chain reaction (qRT-PCR). The abundance of matrix-remodeling associated 5 (MXRA5) protein was detected by western blots (WB). The interaction between miR-30b and MXRA5 was investigated by bioinformatics analysis and luciferase activity assay. The effect of miR-30b and MXRA5 on mitogen-activated protein kinases (MAPK) pathway and invasion was evaluated by WB. Then we found miR-30b was highly expressed in PE and its overexpression inhibited cell viability and invasion while enhanced apoptosis in JEG-3 and HTR8/SVneo cells. Moreover, MXRA5 was targeted by miR-30b and MXRA5 restoration attenuated the effect of miR-30b on cell processes in HTR8/SVneo cells. Besides, both of miR-30b and MXRA5 were associated with MAPK pathway in HTR8/SVneo cells. Our data suggested miR-30b might contribute to PE through inhibiting cell viability, invasion while inducing apoptosis of placental trophoblast cells via MAPK pathway by targeting MXRA5. These indicated that miR-30b might be a novel biomarker for PE treatment.

Albumin downregulates Klotho in tubular cells

Background

Kidney tubular cells are the main sources of Klotho, a protein with phosphaturic action. Genetic Klotho deficiency causes premature cardiovascular aging in mice. Human chronic kidney disease (CKD) is characterized by acquired Klotho deficiency. Despite the lack of uremic toxin accumulation, Category G1 CKD [(normal glomerular filtration rate (GFR)] is already associated with decreased Klotho and with premature cardiovascular aging.

Methods

We have explored whether albuminuria, a criterion to diagnose CKD when GFR is normal, may directly decrease Klotho expression in human CKD, preclinical models and cultured tubular cells.

Results

In a CKD cohort, albuminuria correlated with serum phosphate after adjustment for GFR, age and sex. In this regard, urinary Klotho was decreased in patients with pathological albuminuria but preserved GFR. Proteinuria induced in rats by puromycin aminonucleoside and in mice by albumin overload was associated with interstitial inflammation and reduced total kidney Klotho messenger ribonucleic acid (mRNA) expression. Western blot disclosed reduced kidney Klotho protein in proteinuric rats and mice and immunohistochemistry localized the reduced kidney Klotho expression to tubular cells in proteinuric animals. In cultured murine and human tubular cells, albumin directly decreased Klotho mRNA and protein expression. This was inhibited by trichostatin A, an inhibitor of histone deacetylases, but unlike cytokine-induced Klotho downregulation, not by inhibitors of nuclear factor kappa-light-chain-enhancer of activated B cells.

Conclusions

In conclusion, albumin directly decreases Klotho expression in cultured tubular cells. This may explain, or at least contribute to, the decrease in Klotho and promote fibroblast growth factor 23 resistance in early CKD categories, as observed in preclinical and clinical proteinuric kidney disease.

Genome-wide ChIP-seq data with a transcriptome analysis reveals the groups of genes regulated by histone demethylase LSD1 inhibition in esophageal squamous cell carcinoma cells

Expression of genes is controlled by histone modification, histone acetylation and methylation, but abnormalities of these modifications have been observed in carcinogenesis and cancer development. The effect of the lysine-specific histone demethylase 1 (LSD1) inhibitor, a demethylating enzyme of histones, is thought to be caused by controlling the expression of genes. The aim of the present study is to elucidate the efficacies of the LSD1 inhibitor on the gene expression of esophageal cancer cell lines using chromatin immunoprecipitation (ChIP)-Seq. A comprehensive analysis of gene expression changes in esophageal squamous cell carcinoma (ESCC) cell lines induced by the LSD1 inhibitor NCL1 was clarified via analysis using microarray. In addition, ChIP-seq analysis was conducted using a SimpleChIP plus Enzymatic Chromatin IP kit. NCL1 strongly suppressed the proliferation of T.Tn and TE2 cells, which are ESCC cell lines, and further induced apoptosis. According to the combinatory analysis of ChIP-seq and microarray, 17 genes were upregulated, and 16 genes were downregulated in both cell lines. The comprehensive gene expression study performed in the present study is considered to be useful for analyzing the mechanism of the antitumor effect of the LSD1 inhibitor in patients with ESCC.

Matrix remodeling associated 7 promotes differentiation of bone marrow mesenchymal stem cells toward osteoblasts

The matrix remodeling associated 7 (MXRA7) gene had been ill-studied and its biology remained to be discovered. Inspired by our previous findings and public datasets concerning MXRA7, we hypothesized that the MXRA7 gene might be involved in bone marrow mesenchymal stem cells (BMSCs) functions related to bone formation, which was checked by utilizing in vivo or in vitro methodologies. Micro-computed tomography of MXRA7-deficient mice demonstrated retarded osteogenesis, which was reflected by shorter femurs, lower bone mass in both trabecular and cortical bones compared with wild-type (WT) mice. Histology confirmed the osteopenia-like feature including thinner growth plates in MXRA7-deficient femurs. Immunofluorescence revealed less osteoblasts in MXRA7-deficient femurs. Polymerase chain reaction or western blot analysis showed that when WT BMSCs were induced to differentiate toward osteoblasts or adipocytes in culture, MXRA7 messenger RNA or protein levels were significantly increased alongside osteoblasts induction, but decreased upon adipocytes induction. Cultured MXRA7-deficient BMSCs showed decreased osteogenesis upon osteogenic differentiation induction as reflected by decreased calcium deposition or lower expression of genes responsible for osteogenesis. When recombinant MXRA7 proteins were supplemented in a culture of MXRA7-deficient BMSCs, osteogenesis or gene expression was fully restored. Upon osteoblast induction, the level of active β-catenin or phospho-extracellular signal-regulated kinase in MXRA7-deficient BMSCs was decreased compared with that in WT BMSCs, and these impairments could be rescued by recombinant MXRA7 proteins. In adipogenesis induction settings, the potency of MXRA7-deficient BMSCs to differentiate into adipocytes was increased over the WT ones. In conclusion, this study demonstrated that MXRA7 influences bone formation via regulating the balance between osteogenesis and adipogenesis in BMSCs.

Integrative expression analysis identifies a novel interplay between CFTR and linc-SUMF1-2 that involves CF-associated gene dysregulation

Cystic fibrosis transmembrane regulator (CFTR) is a cyclic AMP-dependent Cl^- channel, and its dysfunction, due to CFTR gene mutations, causes the lethal inherited disorder cystic fibrosis (CF). To date, widespread dysregulation of certain coding genes in CF airway epithelial cells is well studied and considered as the driver of pulmonary abnormality. However, the involvement of non-coding genes, novel classes of functional RNAs with little or no protein-coding capacity, in the regulation of CF-associated gene dysregulation is poorly understood. Here, we utilized integrative analyses of human transcriptome array (HTA) and characterized 99 coding and 91 non-coding RNAs that are dysregulated in CFTR-defective CF bronchial epithelial cell line CFBE41o-. Among these genes, the expression level of linc-SUMF1-2, an intergenic non-coding RNA (lincRNA) whose function is unknown, was inversely correlated with that of WT-CFTR and consistently higher in primary human CF airway epithelial cells (DHBE-CF). Further integrative analyses under linc-SUMF1-knockdown condition determined MXRA5, SEMA5A, CXCL10, AK022877, CTGF, MYC, AREG and LAMB3 as both CFTR- and linc-SUMF1-2-dependent dysregulated gene sets in CF airway epithelial cells. Overall, our analyses reveal linc-SUMF1-2 as a dysregulated non-coding gene in CF as well as CFTR-linc-SUMF1-2 axis as a novel regulatory pathway involved in CF-associated gene dysregulation.